Search NASA⌕ Search

SEARCH · Search NASA

Results for “(non)linear 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.

At least 163 records · Page 9

Squash-Box Feasibility Driven Differential Dynamic Programming

Recently, Differential Dynamic Programming (DDP) and other similar algorithms have become the solvers of choice when performing non-linear Model Predictive Control (nMPC) with modern robotic devices. The reason is that they have a lower computational cost per iteration when compared with off-the-shelf Non-Linear Programming (NLP) solvers, which enables its online operation. However, they cannot handle constraints, and are known to have poor convergence capabilities. In this paper, we propose a method to solve the optimal control problem with control bounds through a squashing function (i.e., a sigmoid, which is bounded by construction). It has been shown that a naive use of squashing functions damage the convergence rate. To tackle this, we first propose to add a quadratic barrier that avoids the difficulty of the plateau produced by the sigmoid. Second, we add an outer loop that adapts both the sigmoid and the barrier; it makes the optimal control problem with the squashing function converge to the original control-bounded problem. To validate our method, we present simulation results for different types of platforms including a multi-rotor, a biped, a quadruped and a humanoid robot.

Navarro, Angel Santamaria↗

Towards fault-tolerant optimal control

The paper considers the design of fault-tolerant controllers that may endow systems with dynamic reliability. Results for jump linear quadratic Gaussian control problems are extended to include random jump costs, trajectory discontinuities, and a simple case of non-Markovian mode transitions.

Chizeck, H. J.↗

Mission Operations Planning with Preferences: An Empirical Study

This paper presents an empirical study of some nonexhaustive approaches to optimizing preferences within the context of constraint-based, mixed-initiative planning for mission operations. This work is motivated by the experience of deploying and operating the MAPGEN (Mixed-initiative Activity Plan GENerator) system for the Mars Exploration Rover Mission. Responsiveness to the user is one of the important requirements for MAPGEN, hence, the additional computation time needed to optimize preferences must be kept within reasonabble bounds. This was the primary motivation for studying non-exhaustive optimization approaches. The specific goals of rhe empirical study are to assess the impact on solution quality of two greedy heuristics used in MAPGEN and to assess the improvement gained by applying a linear programming optimization technique to the final solution.

Bresina, John L.↗

Optimal Aerosol Parameterization for Remote Sensing Retrievals

We have developed a new algorithm for the retrieval of aerosol and gases from SAGE It1 solar transmission measurements. This algorithm improves upon the NASA operational algorithm in several key aspects, including solving the problem non-linearly and incorporating a new methodology for separating the contribution of aerosols and gases. In order to extract aerosol information we have built a huge database of aerosol models for both stratospheric and tropospheric aerosols, and polar stratospheric cloud particles. This set of models allows us to calculate a vast range of possible extinction spectra for aerosols. and from these, derive a set of eigenvectors which then provide the basis set used in our inversion algorithm. Our aerosol algorithm and retrievals are described in several articles (listed in References Section) published under this grant. In particular they allow us to analyze the spectral properties of aerosols and PSCs and ultimately derive their microphysical properties. We have found some considerable differences between our spectra and the ones derived from the SAGE III operational algorithm. These are interesting as they provide an independent check on the validity of published aerosol data and, in particular, on their associated uncertainties. In order to understand these differences, we are assembling independent aerosol data from other sources with which to make comparisons. We have carried out extensive comparisons of our ozone retrievals with both SAGE III and independent lidar, ozonesonde, and satellite measurements (Polyakov et al., 2004). These show very good agreement throughout the stratosphere and help to quantify differences which can be attributed to natural variation in ozone versus that produced by algorithmic differences. In the mid - upper stratosphere, agreement with independent data was generally within 5 - 20%. but in the lower stratosphere the differences were considerably larger. We believe that a large proportion of this discrepancy in the lower stratosphere is attributable to natural variation, and is also seen in comparisons between lidar and ozonesonde measurements. NO2 profiles obtained with our algorithm were compared to those obtained through the SAGE III operational algorithm and exhibited differences of 20 - 40%. Our retrieved profiles agree with the HALOE NO2 measurements significantly better than those of the operational retrieval. In other work (described below), we are extending our aerosol retrievals into the infrared regime and plan to perform retrievals from combined uv-visible-infrared spectra. This work will allow us to use the spectra to derive the size and composition of aerosols, and we plan to employ our algorithms in the analysis of PSC spectra. We are presently also developing a limb-scattering algorithm to retrieve aerosol data from limb measurements of solar scattered radiation.

Newchurch, Michael J.↗

Three-dimensional Finite Element Formulation and Scalable Domain Decomposition for High Fidelity Rotor Dynamic Analysis

This paper has two objectives. The first objective is to formulate a 3-dimensional Finite Element Model for the dynamic analysis of helicopter rotor blades. The second objective is to implement and analyze a dual-primal iterative substructuring based Krylov solver, that is parallel and scalable, for the solution of the 3-D FEM analysis. The numerical and parallel scalability of the solver is studied using two prototype problems - one for ideal hover (symmetric) and one for a transient forward flight (non-symmetric) - both carried out on up to 48 processors. In both hover and forward flight conditions, a perfect linear speed-up is observed, for a given problem size, up to the point of substructure optimality. Substructure optimality and the linear parallel speed-up range are both shown to depend on the problem size as well as on the selection of the coarse problem. With a larger problem size, linear speed-up is restored up to the new substructure optimality. The solver also scales with problem size - even though this conclusion is premature given the small prototype grids considered in this study.

Datta, Anubhav↗

The linear and non-linear aerodynamics of three-surface aircraft concepts

It is noted that most modern fighter aircraft rely on vortex interaction to provide lift enhancement at maneuvering angles of attack. It is shown that the close-coupled horizontal canard in a three-surface configuration provides a control surface which in addition to its other control functions, can be used to optimize this vortex interaction. Attention is given to a study intended to provide a detailed understanding of the aerodynamics of this type of configuration. The discussion examines the results of this investigation and hypotheses are presented to explain the linear and nonlinear aerodynamic phenomena observed.

Agnew, J. W.↗

Generating Far-Infrared Radiation By Two-Wave Mixing

Far-infrared radiation 1 to 6 GHz generated by two-wave mixing in asymmetrically grown GaAs/AlxGa1-xAs multiple-quantum-well devices. Two near-infrared semiconductor diode lasers phase-locked. Outputs amplified, then combined in semiconductor nonlinear multiple-quantum-well planar waveguide. Necessary to optimize design of device with respect to three factors: high degree of confinement of electromagnetic field in nonlinear medium to maximize power density, phase matching to extend length of zone of interaction between laser beams in non-linear medium, and nonlinear susceptibility. Devices used as tunable local oscillators in heterodyne-detection radiometers.

Borenstain, Shmuel↗

Controls-structures-interaction dynamics during RCS control of the Orbiter/SRMS/SSF configuration

During the assembly flights of the Space Station Freedom (SSF), the Orbiter will either dock with the SSF and retract to the final berthed position, or will grapple the SSF using the Shuttle Remote Manipulator System (SRMS) and maneuver the SRMS coupled vehicles to their final berthed position. The SRMS method is expected to take approximately one to one and a half hours to complete and require periodic attitude corrections by either the Orbiter or the SSF reaction control system (RCS) or continuous control by a control moment gyro (CMG) system with RCS desaturation as required. Free drift of the attached vehicles is not currently thought to be acceptable because the desired system attitude will quickly deteriorate due to unbalanced gravity gradient and aerodynamic torques resulting in power generation problems, thermodynamic control problems, and communications problems. This paper deals with the simulation and control of the SRMS during trunnion/latch interaction dynamics and during RCS maneuvers. The SRMS servo drive joints have highly non-linear elastic characteristics which tend to degrade sensitive control strategies. In addition the system natural frequencies are extremely low and depend on the drive joint deflections and SRMS geometric position. The lowest mean period of oscillation for the Orbiter/SRMS/SSF(MB6) system in brakes hold mode positioned near the final berthed position is approximately 120 seconds. A detailed finite element model of the SRMS has been developed and used in a newly developed SRMS systems dynamics simulation to investigate the non-linear transient response dynamics of the Orbiter/SRMS/SSF systems. The present SRMS control strategy of brakes only recommended by the Charles Draper Labs is contrasted with a robust controller developed by the authors. The robust controller uses an optimal inear quadratic regulator (LQR) to optimally place the closed-loop poles of a multivariable continuous-time system within the common region of an open sector with the sector angle plus or minus 45 degrees from the negative real axis, and the left-hand side of a parallel to the imaginary axis in the complex s-plane. This guarantees that the critical damping ratio for the desired control modes is equal to or in excess of 0.707. The matrix sign function is used for solving the Riccati equations which appear in the controller design procedure. Fast and stable algorithms have recently been developed for the computation of the matrix sign function. Simulation results are given which demonstrate the potential CSI involvement for the current SRMS control system and the proposed control system.

Schliesing, J. A.↗

On the application of deterministic optimization methods to stochastic control problems.

A technique is presented by which one can apply the Minimum Principle of Pontryagin to stochastic optimal control problems formulated around linear systems with Gaussian noises and general cost criteria. Using this technique, the stochastic nature of the problem is suppressed but for two expectation operations, the optimization being essentially deterministic. The technique is applied to systems with quadratic and non-quadratic costs to illustrate its use.

Kramer, L. C.↗

Neural networks: What non-linearity to choose

Neural networks are now one of the most successful learning formalisms. Neurons transform inputs (x(sub 1),...,x(sub n)) into an output f(w(sub 1)x(sub 1) + ... + w(sub n)x(sub n)), where f is a non-linear function and w, are adjustable weights. What f to choose? Usually the logistic function is chosen, but sometimes the use of different functions improves the practical efficiency of the network. The problem of choosing f as a mathematical optimization problem is formulated and solved under different optimality criteria. As a result, a list of functions f that are optimal under these criteria are determined. This list includes both the functions that were empirically proved to be the best for some problems, and some new functions that may be worth trying.

Kreinovich, Vladik YA.↗

Pterodactyl: Guidance and Control of a Symmetric Deployable Entry Vehicle using an Aerodynamic Control System

The NASA-funded Pterodactyl project seeks to advance the state-of-the-art for varying entry vehicle types by developing unconventional guidance and control technologies for Deployable Entry Vehicles (DEVs) that can be applied to different entry vehicle configurations. Prior work by the authors [1–5] involved developing both traditional and novel integrated guidance and control solutions for a Pterodactyl Baseline Vehicle (PBV), a variant of an asymmetric DEV called the Lifting Nano ADEPT (LNA) [6]. In the prior studies, two different guidance schemes were designed and implemented for the PBV: (i) traditional bank angle guidance developed using the Fully Numerical Predictor-Corrector Entry Guidance (FNPEG) and (ii) novel angle of attack and sideslip (α - β) guidance developed using FNPEG with Uncoupled Range Control [4]. Using Linear Quadratic Regulator (LQR) optimal control methods with state-feedback integral control designs, these guidance trajectories were designed to be tracked using (i) a conventional propulsive entry vehicle control hardware architecture - reaction control systems (RCS) and (ii) novel non-propulsive entry vehicle control systems - aerodynamic flap control system (FCS) and moving mass control system (MMCS) [1]. The novel FCS and MMCS architectures were designed to track α - β guidance commands while the RCS was designed to track bank angle commands. It was discovered that the asymmetric DEV, the PBV, experienced a non-zero induced roll moment due to sideslip that the FCS and MMCS architectures had limited capability to trim out. These two architectures were designed to provide independent angle of attack and sideslip commands with limited consideration for roll moment generation to trim. As a result, for the PBV, the FCS and MMCS configurations as designed, were limited in providing the control authority needed to track an α - β guidance trajectory [1]. These results form the motivation for the work presented in this paper - utilizing an aerodynamic control system to track α - β guidance commands for a symmetric DEV with the expectation that a symmetric entry vehicle will have zero or significantly reduced roll moment due to sideslip that the FCS can handle when tracking an α - β guidance trajectory. To demonstrate the feasibility of a novel guidance and control architecture on a DEV, we utilize a symmetric DEV, the PBV-II, for (i) the novel α - β guidance development using FNPEG with Uncoupled Range Control and (ii) LQR control design using eight aerodynamic control surfaces. This paper demonstrates that the novel uncoupled α - β guidance tracking can be achieved using aerodynamic control surfaces on a symmetric deployable entry vehicle configuration.

Wendy A Okolo↗

Pterodactyl: Control Architectures Development for Integrated Control Design of a Mechanically Deployed Entry Vehicle

The need to return high mass payloads is driving the development of a new class of vehicles, Deployable Entry Vehicles (DEV) for which feasible and optimized control architectures have not been developed. The Pterodactyl project, seeks to advance the current state-of-the-art for entry vehicles by developing a design, test, and build capability for DEVs that can be applied to various entry vehicle configurations. This paper details the efforts on the NASA-funded Pterodactyl project to investigate multiple control techniques for the Lifting Nano-ADEPT (LNA) DEV. We design and implement multiple control architectures on the LNA and evaluate their performance in achieving varying guidance commands during entry.First we present an overview of DEVs and the Lifting Nano-ADEPT (LNA), along with the physical LNA configuration that influences the different control designs. Existing state-of-the-art for entry vehicle control is primarily propulsive as reaction control systems (RCS) are widely employed. In this work, we analyze the feasibility of using both propulsive control systems such as RCS to generate moments, and non-propulsive control systems such as aerodynamic control surfaces and internal moving mass actuations to shift the LNA center of gravity and generate moments. For these diverse control systems, we design different multi-input multi-output (MIMO) state-feedback integral controllers based on linear quadratic regulator (LQR) optimal control methods. The control variables calculated by the controllers vary, depending on the control system being utilized and the outputs to track for the controller are either the (i) bank angle or the (ii) angle of attack and sideslip angle as determined by the desired guidance trajectory. The LQR control design technique allows the relative allocation of the control variables through the choice of the weighting matrices in the cost index. Thus, it is easy to (i) specify which and how much of a control variable to use, and (ii) utilize one control design for different control architectures by simply modifying the choice of the weighting matrices.By providing a comparative analysis of multiple control systems, configurations, and performance, this paper and the Pterodactyl project as a whole will help entry vehicle system designers and control systems engineers determine suitable control architectures for integration with DEVs and other entry vehicle types.

Okolo, Wendy A.↗

An integrated approach to the synthesis of geometrically non-linear structures

An integrated approach to the minimum weight design of geometrically nonlinear three-dimensional truss structures with geometric imperfections, subject to inequality constraints on static displacements, stresses, local buckling and cross sectional areas, is investigated. The integrated structural synthesis problem involves design and response quantities as independent variables and equilibrium equations, describing the finite element model, as equality constraints. The nonlinear structural analysis and the optimization are thus merged together into a single process. A computer program developed to compute the constraint values and analytical gradients is coupled with a generalized reduced gradient algorithm to solve the integrated problem. Numerical results for a geometrically nonlinear shallow dome example problem are presented for various types of imperfections. Furthermore, it is found that the algorithm is capable of detecting and guarding against system as well as element elastic instability using equilibrium information only, that is, without imposing system and local buckling inequality constraints.

Smaoui, H.↗

Motion Cueing Algorithm Modification for Improved Turbulence Simulation

Atmospheric turbulence cueing produced by flight simulator motion systems has been less than satisfactory because the turbulence profiles have been attenuated by the motion cueing algorithms. Cardullo and Ellor initially addressed this problem by directly porting the turbulence model output to the motion system. Reid and Robinson addressed the problem by employing a parallel aircraft model, which is only stimulated by the turbulence inputs and adding a filter specially designed to pass the higher turbulence frequencies. There have been advances in motion cueing algorithm development at the Man-Machine Systems Laboratory, at SUNY Binghamton. In particular, the system used to generate turbulence cues has been studied. The Reid approach, implemented by Telban and Cardullo, was employed to augment the optimal motion cueing algorithm installed at the NASA LaRC Simulation Laboratory, driving the Visual Motion Simulator. In this implementation, the output of the primary flight channel was added to the output of the turbulence channel and then sent through a non-linear cueing filter. The cueing filter is an adaptive filter; therefore, it is not desirable for the output of the turbulence channel to be augmented by this type of filter. The likelihood of the signal becoming divergent was also an issue in this design. After testing on-site it became apparent that the architecture of the turbulence algorithm was generating unacceptable cues. As mentioned above, this cueing algorithm comprised a filter that was designed to operate at low bandwidth. Therefore, the turbulence was also filtered, augmenting the cues generated by the model. If any filtering is to be done to the turbulence, it will utilize a filter with a much higher bandwidth, above the frequencies produced by the aircraft response to turbulence. The authors have developed an implementation wherein only the signal from the primary flight channel passes through the nonlinear cueing filter. This paper discusses three new algorithms. Testing shows that the new methods provide the pilot with a more realistic sensation of turbulence; the cues are not attenuated by algorithm. Results of offline testing show the credibility of the models. Offline test verification was based primarily on the evaluation of the power spectral density of the outputs and the time response.

Ercole, Anthony V.↗

A Neural Network Aero Design System for Advanced Turbo-Engines

An inverse design method calculates the blade shape that produces a prescribed input pressure distribution. By controlling this input pressure distribution the aerodynamic design objectives can easily be met. Because of the intrinsic relationship between pressure distribution and airfoil physical properties, a Neural Network can be trained to choose the optimal pressure distribution that would meet a set of physical requirements. Neural network systems have been attempted in the context of direct design methods. From properties ascribed to a set of blades the neural network is trained to infer the properties of an 'interpolated' blade shape. The problem is that, especially in transonic regimes where we deal with intrinsically non linear and ill posed problems, small perturbations of the blade shape can produce very large variations of the flow parameters. It is very unlikely that, under these circumstances, a neural network will be able to find the proper solution. The unique situation in the present method is that the neural network can be trained to extract the required input pressure distribution from a database of pressure distributions while the inverse method will still compute the exact blade shape that corresponds to this 'interpolated' input pressure distribution. In other words, the interpolation process is transferred to a smoother problem, namely, finding what pressure distribution would produce the required flow conditions and, once this is done, the inverse method will compute the exact solution for this problem. The use of neural network is, in this context, highly related to the use of proper optimization techniques. The optimization is used essentially as an automation procedure to force the input pressure distributions to achieve the required aero and structural design parameters. A multilayered feed forward network with back-propagation is used to train the system for pattern association and classification.

Sanz, Jose M.↗

Convergence Estimates for Multidisciplinary Analysis and Optimization

A quantitative analysis of coupling between systems of equations is introduced. This analysis is then applied to problems in multidisciplinary analysis, sensitivity, and optimization. For the sensitivity and optimization problems both multidisciplinary and single discipline feasibility schemes are considered. In all these cases a "convergence factor" is estimated in terms of the Jacobians and Hessians of the system, thus it can also be approximated by existing disciplinary analysis and optimization codes. The convergence factor is identified with the measure for the "coupling" between the disciplines in the system. Applications to algorithm development are discussed. Demonstration of the convergence estimates and numerical results are given for a system composed of two non-linear algebraic equations, and for a system composed of two PDEs modeling aeroelasticity.

Arian, Eyal↗

Comparing Consider-Covariance Analysis with Sigma-Point Consider Filter and Linear-Theory Consider Filter Formulations

Recent literature in applied estimation theory reflects growing interest in the sigma-point (also called unscented ) formulation for optimal sequential state estimation, often describing performance comparisons with extended Kalman filters as applied to specific dynamical problems [c.f. 1, 2, 3]. Favorable attributes of sigma-point filters are described as including a lower expected error for nonlinear even non-differentiable dynamical systems, and a straightforward formulation not requiring derivation or implementation of any partial derivative Jacobian matrices. These attributes are particularly attractive, e.g. in terms of enabling simplified code architecture and streamlined testing, in the formulation of estimators for nonlinear spaceflight mechanics systems, such as filter software onboard deep-space robotic spacecraft. As presented in [4], the Sigma-Point Consider Filter (SPCF) algorithm extends the sigma-point filter algorithm to the problem of consider covariance analysis. Considering parameters in a dynamical system, while estimating its state, provides an upper bound on the estimated state covariance, which is viewed as a conservative approach to designing estimators for problems of general guidance, navigation and control. This is because, whether a parameter in the system model is observable or not, error in the knowledge of the value of a non-estimated parameter will increase the actual uncertainty of the estimated state of the system beyond the level formally indicated by the covariance of an estimator that neglects errors or uncertainty in that parameter. The equations for SPCF covariance evolution are obtained in a fashion similar to the derivation approach taken with standard (i.e. linearized or extended) consider parameterized Kalman filters (c.f. [5]). While in [4] the SPCF and linear-theory consider filter (LTCF) were applied to an illustrative linear dynamics/linear measurement problem, in the present work examines the SPCF as applied to nonlinear sequential consider covariance analysis, i.e. in the presence of nonlinear dynamics and nonlinear measurements. A simple SPCF for orbit determination, exemplifying an algorithm hosted in the guidance, navigation and control (GN&C) computer processor of a hypothetical robotic spacecraft, was implemented, and compared with an identically-parameterized (standard) extended, consider-parameterized Kalman filter. The onboard filtering scenario examined is a hypothetical spacecraft orbit about a small natural body with imperfectly-known mass. The formulations, relative complexities, and performances of the filters are compared and discussed.

Lisano, Michael E.↗

GOES-17 Magnetometer On-Orbit Calibration Design and Results

The Geostationary Operational Environmental Satellite-R (GOES-R) is a series of four NOAA weather satellites. Each spacecraft flies a Magnetometer instrument which includes two boom mounted fluxgate sensors located 6.3 and 8.5 meters from the spacecraft. Sensor readings are corrupted by both internal (bias and scale factor) and external (spacecraft emitted fields and misalignments) error sources. In order to estimate the bias, an on-orbit calibration is performed, slewing the spacecraft about multiple axes and then performing a fit of the data to determine zero offset and misalignment terms. The first satellite in the series, GOES-16 was launched in November 2016. Its magnetometer calibration maneuver was found to be inadequate providing insufficient rotation about orthogonal axes. The uncertainty in the estimated biases was found to be in the 3-5 nanotesla range, exceeding the magnetometer accuracy requirement of 1.0 nanotesla. Based on lessons learned, the maneuver was redesigned to provide multiple 360 degree rotations about two orthogonal axes, providing optimal visibility into bias and misalignment. Varying ambient fields over the duration of the calibration maneuvers, corrupts the ability to accurately determine parameters of interest. To address this issue, rather than assuming a constant field for the duration of the maneuver, a non-linear least squares filter was employed that uses a series of spline functions to better fit the ambient field. The number of splines can be varied to minimize the least squares filter residuals. Simulations were performed using 360 days of GOES-16 data to bound the uncertainty in the biases for given residuals of the fit to the ambient field. Following the launch of the second satellite in the series, GOES-17 in March 2018, the new calibration maneuver and enhanced bias estimation algorithm were exercised twice. Based on the residuals of the fit, the uncertainty of the estimated zero offset was reduced by an order of magnitude from 3-5 nanotesla to approximately 0.4 nanotesla per axis. This paper provides details of the maneuver design and the innovative algorithm used to estimate the GOES-17 biases and misalignments.

GOES-17↗