Optimal control of a flexible launch vehicle.
Linear optimal control technique for flexible- booster control system design, showing drift minimum model with matrix transformations for closed-loop dynamics
OPTIMAL CONTROL↗
SEARCH · Search NASA
Results for “dynamic optimization”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Figure control concepts for segmented reflector telescopes
Control analysis activities related to the development of figure control technologies for large space telescopes with precision segmented actively controlled primary reflectors are reviewed. The Precision Segmented Reflector (PSR) configuration is described along with the development and use of the PSR models. Geometric and dynamic models, characterization of figure estimation errors and optimal sensor placement, and the development of quasi-static and dynamic control concepts are outlined. The structure of a quasi-static controller is presented, that utilizes edge sensor measurements to estimate displacement errors at the actuator level and generates figure correction commands. The second approach considers decentralization of the reflector system at the panel level, while the third considers decentralization at the actuator level.
Boussalis, D.↗
Practical and Optimal Sequential Bayesian Experimental Design for Complex Systems Incorporating Human Experimenter Preferences (Final Scientific/Technical Report)
Experiments are indispensable for developing models of complex systems. Carefully designed experiments can provide substantial savings for these expensive data-acquisition opportunities. However, designs based on heuristics are often suboptimal for systems with multiphysics, nonlinear dynamics, and uncertain and noisy environments. Optimal experimental design, while leveraging predictive models, seeks to systematically quantify and maximize the value of experiments. In this project, we focused on the design of multiple experiments, where current approaches are largely suboptimal: batch-design does not adapt to new data acquired during the experiment campaign (no feedback), and greedy/myopic design ignores future dynamics and consequences (no lookahead). We developed the mathematical framework and computational methods for sequential optimal experimental design (sOED) for complex systems. We enabled tractable model-based sOED in a rigorous manner through novel algorithms based on reinforcement learning, and investigated the effects of human experimenters on the design process. Our methods are fully Bayesian, able to quantify and update uncertainty in a principled manner. The traits aimed by our approach—mathematical rigor and optimality, human effects and uncertainty quantification, computational practicality—are crucial for elevating the standards of artificial intelligence (AI) to support decision-making in scientific domains, and contribute toward trust and realistic adoption of AI in experimental design practice.
97 MATHEMATICS AND COMPUTING↗
Feedforward equilibrium trajectory optimization with GSPulse
One of the common tasks required for designing new plasma scenarios or evaluating capabilities of a tokamak is to design the desired equilibria using a Grad-Shafranov (GS) equilibrium solver. However, most standard equilibrium solvers are time-independent and do not include dynamic effects such as plasma current flux consumption, induced vessel currents, or voltage constraints. Another class of tools, plasma equilibrium evolution simulators, do include time-dependent effects. These are generally structured to solve the forward problem of evolving the plasma equilibrium given feedback-controlled voltages. In this work, we introduce GSPulse, a novel algorithm for equilibrium trajectory optimization, that is more akin to a pulse planner than a pulse simulator. GSPulse includes time-dependent effects and solves the inverse problem: given a user-specified set of target equilibrium shapes, as well as limits on the coil currents and voltages, the optimizer returns trajectories of the voltages, currents, and achievable equilibria. This task is useful for scoping performance of a tokamak and exploring the space of achievable pulses. The computed equilibria satisfy both Grad-Shafranov force balance and axisymmetric circuit dynamics. The optimization is performed by restructuring the free-boundary equilibrium evolution equations into a form where it is computationally efficient to optimize the entire dynamic sequence. GSPulse can solve for hundreds of equilibria simultaneously within a few minutes. GSPulse has been validated against NSTX-U and MAST-U experiments and against SPARC feedback control simulations, and is being used to perform scenario design for SPARC. The computed trajectories can be used as feedforward inputs that are connected to the feedback controller to inform and improve feedback performance. The code for GSPulse is available open-source at github.com/jwai-cfs/GSPulse_public.
equilibrium↗
Dynamic analysis of fully constrained Cable-Driven Parallel Robots for automated prefabricated component installation
This paper presents a dynamic analysis and validation framework to assess a fully constrained six-anchor Cable-Driven Parallel Robot (CDPR) for automated installation of prefabricated facade components. Compared with conventional eight-anchor systems, the six-anchor configuration simplifies setup and reduces cost, but it also reduces control authority, shrinks the wrench-feasible workspace, and tightens orientation limits. Consequently, it is unclear a priori whether dynamically feasible trajectories exist to move the end effector from pickup to the facade. A constrained trajectory optimization is formulated to enforce the system dynamics, cable-tension bounds, and pose/velocity limits, and the framework is evaluated in simulation at three levels: (i) an idealized reference model, (ii) a lab-scale prototype model incorporating measured anchor misalignments and identified damping, and (iii) a full-scale three-story building model with load decomposition for structural feasibility checks. Across these scenarios, the analysis shows that optimal, constraint-satisfying trajectories exist that move the end effector from pickup to installation while maintaining a near-plumb, level orientation at the final pose. Collectively, this multi-scale dynamic analysis and validation framework supports the deployment readiness of the six-anchor CDPR and provides a prototype-based sensitivity case study of how measured anchor placement deviations affect feasibility.
CDPR↗
Gate assisted turn-off thyristor with cathode shunts and dynamic gate
A 1,000-V, 200-A gate-assisted turn-off thyristor (GATT) is described, whose design features include an interdigitated shunted cathode, a dynamic gate, a means for optimizing the carrier lifetime level, and a bypass diode. The device physics of gate-assisted turn-off are reviewed. Based on this, improvements in the design are described. It is shown that a prime failure mode can be eliminated and that the gate-assist signal voltage can be substantially decreased by employing a shunted cathode emitter. The test data show excellent turn-on characteristics due to the dynamic gate and the long perimeter of the edge of the main cathode. Turn-off times as short as 3 microsec are obtained. The combination of controlling the carrier lifetime with a precisely controlled and easily variable irradiation dose of high energy electrons with gate assist current provides for simple, precision tailoring of the device characteristics to the intended application.
Schlegel, E. S.↗
Optimization Study of the Ames 0.5 Two-Stage Light Gas Gun
There is a need for more faithful simulation of space debris impacts on various space vehicles. Space debris impact velocities can range up to 14 km/sec and conventional two-stage light gas guns with moderately heavy saboted projectiles are limited to launch velocities of 7-8 km/sec. Any increases obtained in the launch velocities will result in more faithful simulations of debris impacts. It would also be valuable to reduce the maximum gun and projectile base pressures and the gun barrel erosion rate. In this paper, the results of a computational fluid dynamics (CFD) study designed to optimize the performance of the NASA Ames 0.5' gun by systematically varying seven gun operating parameters are reported. Particularly beneficial effects were predicted to occur if (1) the piston mass was decreased together with the powder mass and the hydrogen fill pressure and (2) the pump tube length was decreased. The optimum set of changes in gun operating conditions were predicted to produce an increase in muzzle velocity of 0.7-1.0 km/sec, simultaneously with a substantial decrease in gun erosion. Preliminary experimental data have validated the code predictions. Velocities of up to 8.2 km/sec with a 0.475 cm diameter saboted aluminum sphere have been obtained, along with large reductions in gun erosion rates.
Bogdanoff, D. W.↗
The Analysis and Design of Low Boom Configurations Using CFD and Numerical Optimization Techniques
The use of computational fluid dynamics (CFD) for the analysis of sonic booms generated by aircraft has been shown to increase the accuracy and reliability of predictions. CFD takes into account important three-dimensional and nonlinear effects that are generally neglected by modified linear theory (MLT) methods. Up to the present time, CFD methods have been primarily used for analysis or prediction. Some investigators have used CFD to impact the design of low boom configurations using trial and error methods. One investigator developed a hybrid design method using a combination of Modified Linear Theory (e.g. F-functions) and CFD to provide equivalent area due to lift driven by a numerical optimizer to redesign or modify an existing configuration to achieve a shaped sonic boom signature. A three-dimensional design methodology has not yet been developed that completely uses nonlinear methods or CFD. Constrained numerical optimization techniques have existed for some time. Many of these methods use gradients to search for the minimum of a specified objective function subject to a variety of design variable bounds, linear and nonlinear constraints. Gradient based design optimization methods require the determination of the objective function gradients with respect to each of the design variables. These optimization methods are efficient and work well if the gradients can be obtained analytically. If analytical gradients are not available, the objective gradients or derivatives with respect to the design variables must be obtained numerically. To obtain numerical gradients, say, for 10 design variables, might require anywhere from 10 to 20 objective function evaluations. Typically, 5-10 global iterations of the optimizer are required to minimize the objective function. In terms of using CFD as a design optimization tool, the numerical evaluation of gradients can require anywhere from 100 to 200 CFD computations per design for only 10 design variables. If one CFD computation requires an hour of computational time on a Cray computer, one can see that the use of constrained numerical optimization quickly becomes impractical.Hence, in order to practically couple a numerical design optimization technique with a CFD method, the CFD method must be extremely efficient with running times on the order of only minutes. The CFD Euler code developed under NASA sponsorship and referred to as MIM3D-SB for the most part fulfills these efficiency requirements. Analysis of wing- body configurations can be computed in a matter of a few minutes. The present study will concentrate on the feasibility of the use of this CFD code in conjunction with a numerical design optimization technique for the sonic boom reduction of candidate HSCT configurations. A preliminary supersonic aircraft design system has been established that utilizes the numerical design optimization code NPSOL developed at Stanford University coupled with the supersonic NUM3D-SB CFD code. Many questions still need to be answered in regard to using CFD and numerical optimizers as design tools. There are difficulties related to both the CFD codes and the numerical optimizers. Numerical optimizers can converge to a local minima rather than a global minima. This behavior is largely a function of the initial guess in the design space. The optimizer also is searching for a minimum of the function in terms of its derivative without any regard to the actual function value. Numerically (i.e. CFD) determined gradients can also generate spurious numerical local minima. In addition, for the sonic boom problem, grid fineness will also determine the accuracy of the final design solution. Design optimization methods work well on problems defined by continuous objective functions. The sonic boom signature design problem is not necessarily defined by a continuous objective function. The signature can have a variety of shapes; i.e. from N-wave to multiple shocks. The far-field or ground signature may not transition continuously from one shape to another and hence, may exhibit discontinuous behavior. This is also a source of difficulty in using design optimization methods.In the following sections, several low boom and one reference aircraft configuration will be analyzed to predict their sonic boom signature characteristics. Modifications to some of these designs will also be presented to demonstrate the feasibility of using CFD as a design tool and to demonstrate the feasibility of designing shaped sonic boom signatures. Design modifications to some configurations will be presented to demonstrate the feasibility of achieving shaped signatures with reduced levels and not necessarily to represent realistic or aerodynamically efficient design modifications. Fuselage volume or camber are used as design variables in order to have a minimal effect on the primary wing aerodynamics. The paper will also seek to demonstrate whether a hybrid or ramped signature is feasible to achieve. For the low-boom configurations, the CFD predicted signatures will be compared qualitatively to their MLT design signatures.
Siclari, Michael J.↗
The poststall nonlinear dynamics and control of an F-18: A preliminary investigation
The successful high angle of attack (HAOA) operation of fighter aircraft will necessarily require the introduction of a new onboard control methodology that address the nonlinearity of the system when flown at the stall/poststall limits of the craft's flight envelope. As a precursor to this task, a researcher endeavored to familarize himself with the dynamics of one specific aircraft, the F-18, when it is flown at HAOA. This was accomplished by conducting a number of real time flight sorties using the NASA-Langley Research Center's F-18 simulator, which was operated with a pilot in the loop. In addition to developing a first hand familarity with the aircraft's dynamic characteristic at HAOA, work was also performed to identify the input/output operational footprint of the F-18's control surfaces. This investigator proposes to employ the nonlinear models of the plant identified this summer in a subsequent research effort that will make it possible to fly the F-18 effectively at poststall angles of attack. The controller design used there will rely on a new technique proposed by this investigator that provides for the automatic generation of online optimal control solutions for nonlinear dynamic systems.
Patten, William N.↗
Frequency-Domain Deconvolution for Flight Dynamics Applications
A deconvolution method is presented for estimating input data from measured output data and a model of the dynamic process involved. The method uses an optimal Wiener filter for separating the measured data into signal and noise components, and a high-accuracy Fourier transform for inverting the model dynamics in the frequency domain. The method is an extension of optimal Fourier smoothing, and uses a technique to enhance the contrast between the signal and noise spectra in designing the Wiener filter. The deconvolution method was applied to simulation and flight test data for the purposes of removing unwanted distortions introduced by signal-conditioning filters and sensor dynamics, and for reconstructing turbulence inputs from measured sensor data. Results indicated hat the method performs well given good signal-to-noise levels and accurate models of the dynamic process.
Grauer, Jared A.↗
Optimal compensator structure for linear time-invariant plant with inaccessible states
The problem is considered of designing an optimal linear time-invariant dynamic compensator for the regulation of an n-th order linear time-invariant plant with m independent outputs. The initial plant state is characterized by its first and second moments, and the cost is usual quadratic infinite-time penalty on the state and control, averaged over the initial plant and compensator states. The compensator is based on a minimal-order Luenberger observer and consequently has fixed dimension n-m. Necessary and sufficient conditions are derived for optimality of the compensator gains. The optimal compensator is shown to be unique if the plant has a particular canonical form, and, in general, for any arbitrary plant, the class of all optimal compensators is precisely determined.
Blanvillain, P. J. P.↗
A reduced-dynamic technique for precise orbit determination
Observations of the Global Positioning System (GPS) will enable a reduced-dynamic technique for achieving subdecimeter orbit determination of earth-orbiting satellites. With this technique, information on the transition between satellite states at different observing times is furnished by both a formal dynamic model and observed satellite positional change (which is inferred kinematically from continuous GPS carrier-phase data). The relative weighting of dynamic and kinematic information can be freely varied. Covariance studies show that in situations where observing geometry is poor and the dynamic model is good, the model dominates determination of the state transition; where the dynamic model is poor and the geometry strong, carrier phase governs the determination of the transition. When neither kinematic nor dynamic information is clearly superior, the reduced-dynamic combination of the two can substantially improve the orbit-determination solution. Guidelines are given here for selecting a near-optimal weighting for the reduced-dynamic solution, and sensitivity of solution accuracy to this weighting is examined.
Wu, S. C.↗
Universal energy-speed-accuracy trade-offs in driven nonequilibrium systems
The connection between measure theoretic optimal transport and dissipative nonequilibrium dynamics provides a language for quantifying nonequilibrium control costs, leading to a collection of thermodynamic speed limits, which rely on the assumption that the target probability distribution is perfectly realized. This is almost never the case in experiments or numerical simulations, so here we address the situation in which the external controller is imperfect. We obtain a lower bound for the dissipated work in generic nonequilibrium control problems that (1) is asymptotically tight and (2) matches the thermodynamic speed limit in the case of optimal driving. Along with analytically solvable examples, we refine this imperfect driving notion to systems in which the controlled degrees of freedom are slow relative to the nonequilibrium relaxation rate, and identify independent energy contributions from fast and slow degrees of freedom. Furthermore, we develop a strategy for optimizing minimally dissipative protocols based on optimal transport flow matching, a generative machine learning technique. Furthermore, this latter approach ensures the scalability of both the theoretical and computational framework we put forth. Crucially, we demonstrate that we can compute the terms in our bound numerically using efficient algorithms from the computational optimal transport literature and that the protocols we learn saturate the bound.
59 BASIC BIOLOGICAL SCIENCES↗
Trajectory optimization and guidance law development for national aerospace plane applications
The problem of onboard trajectory optimization for an airbreathing, single-stage-to-orbit vehicle is examined. A simple model representative of the aerospace plane concept, including a dual-mode propulsion system composed of scramjet and rocket engines, is presented. Consideration is restricted to hypersonic flight within the atmosphere. An energy state approximation is used in a four-state model for flight of a point mass in a vertical plane. Trajectory constraints, including those of dynamic pressure and aerodynamic heating, are initially ignored. Singular perturbation methods are applied in solving the optimal control problem of minimum fuel climb. The resulting reduced solution for the energy state dynamics provides an optimal altitude profile dependent on energy level and control for rocket thrust. A boundary-layer analysis produces an approximate lift control solution in feedback form and accounts for altitude and flight path angle dynamics. The reduced solution optimal climb path is presented for the unconstrained case and the case for which a maximum dynamic pressure constraint is enforced.
Calise, A. J.↗
Performance of Optimized Actuator and Sensor Arrays in an Active Noise Control System
Experiments have been conducted in NASA Langley's Acoustics and Dynamics Laboratory to determine the effectiveness of optimized actuator/sensor architectures and controller algorithms for active control of harmonic interior noise. Tests were conducted in a large scale fuselage model - a composite cylinder which simulates a commuter class aircraft fuselage with three sections of trim panel and a floor. Using an optimization technique based on the component transfer functions, combinations of 4 out of 8 piezoceramic actuators and 8 out of 462 microphone locations were evaluated against predicted performance. A combinatorial optimization technique called tabu search was employed to select the optimum transducer arrays. Three test frequencies represent the cases of a strong acoustic and strong structural response, a weak acoustic and strong structural response and a strong acoustic and weak structural response. Noise reduction was obtained using a Time Averaged/Gradient Descent (TAGD) controller. Results indicate that the optimization technique successfully predicted best and worst case performance. An enhancement of the TAGD control algorithm was also evaluated. The principal components of the actuator/sensor transfer functions were used in the PC-TAGD controller. The principal components are shown to be independent of each other while providing control as effective as the standard TAGD.
Palumbo, D. L.↗
Performance of Optimized Actuator and Sensor Arrays in an Active Noise Control System
Experiments have been conducted in NASA Langley s Acoustics and Dynamics Laboratory to determine the effectiveness of optimized actuator/sensor architectures and controller algorithms for active control of harmonic interior noise. Tests were conducted in a large scale fuselage model - a composite cylinder which simulates a commuter class aircraft fuselage with three sections of trim panel and a floor. Using an optimization technique based on the component transfer functions, combinations of 4 out of 8 piezoceramic actuators and 8 out of 462 microphone locations were evaluated against predicted performance. A combinatorial optimization technique call tabu search was employed to select the optimum transducer arrays. Three test frequencies represent the cases of a strong acoustic and strong structural response, a weak acoustic and strong structural response and a strong acoustic and weak structural response. Noise reduction was obtained using a Time Averaged/Gradient Descent (TAGD) controller. Results indicate that the optimization technique successfully predicted best and worst case performance. An enhancement of the TAGD control algorithm was also evaluated. The principal components of the actuator/sensor transfer functions were used in the PC-TAGD controller. The principal components are shown to be independent of each other while providing control as effective as the standard TAGD.
Palumbo, D. L.↗
Reward based optimization of resonance-enhanced piezoresponse spectroscopy
Dynamic spectroscopies in scanning probe microscopy (SPM) are critical for probing material properties, such as force interactions, mechanical properties, polarization switching, electrochemical reactions, and ionic dynamics. However, the practical implementation of these measurements is constrained by the need to balance imaging time and data quality. Signal to noise requirements favor long acquisition times and high frequencies to improve signal fidelity. However, these are limited on the low end by contact resonant frequency and photodiode sensitivity and on the high end by the time needed to acquire high-resolution spectra or the propensity for sample degradation under high field excitation over long times. The interdependence of key parameters such as instrument settings, acquisition times, and sampling rates makes manual tuning labor-intensive and highly dependent on user expertise, often yielding operator-dependent results. These limitations are prominent in techniques like dual amplitude resonance tracking in piezoresponse force microscopy that utilize multiple concurrent feedback loops for topography and resonance frequency tracking. Here, a reward-driven workflow is proposed that automates the tuning process, adapting experimental conditions in real time to optimize data quality. Furthermore, this approach significantly reduces the complexity and time required for manual adjustments and can be extended to other SPM spectroscopic methods, enhancing overall efficiency and reproducibility.
47 OTHER INSTRUMENTATION↗
Optimal interpolation and the Kalman filter
The estimation theory of stochastic-dynamic systems is described and used in a numerical study of optimal interpolation. The general form of data assimilation methods is reviewed. The Kalman-Bucy, KB filter, and optimal interpolation (OI) filters are examined for effectiveness in performance as gain matrices using a one-dimensional form of the shallow-water equations. Control runs in the numerical analyses were performed for a ten-day forecast in concert with the OI method. The effects of optimality, initialization, and assimilation were studied. It was found that correct initialization is necessary in order to localize errors, especially near boundary points. Also, the use of small forecast error growth rates over data-sparse areas was determined to offset inaccurate modeling of correlation functions near boundaries.
Cohn, S.↗