SEARCH · Search NASA

Results for “Numerical 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 541 records · Page 30

All sources NASA NTRS DOE OSTI

NASA NTRS1973

Optimal detumbling of a large manned spacecraft using an internal moving mass

The use of a movable mass control system to stabilize an arbitrarily tumbling asymmetric vehicle about the maximum inertia axis is considered. A first-order gradient optimization technique is used to minimize angular velocity components along the intermediate and minimum inertia axes. This method permits a wide range of initial guesses for mass position history. Motion of the control mass is along a linear track fixed in the vehicle. The control variable is taken as mass acceleration with respect to body coordinates. Motion is limited to defined quantities and a penalty function is used to insure a given range of positions. Numerical solutions of the optimization equations verify that minimum time detumbling is achieved with the largest permissible movable mass, length of linear track, and positions of the mass on the two coordinates perpendicular to the linear motion.

Kunciw, B. G.↗

NASA NTRS1984

A rotor optimization using regression analysis

The design and development of helicopter rotors is subject to the many design variables and their interactions that effect rotor operation. Until recently, selection of rotor design variables to achieve specified rotor operational qualities has been a costly, time consuming, repetitive task. For the past several years, Kaman Aerospace Corporation has successfully applied multiple linear regression analysis, coupled with optimization and sensitivity procedures, in the analytical design of rotor systems. It is concluded that approximating equations can be developed rapidly for a multiplicity of objective and constraint functions and optimizations can be performed in a rapid and cost effective manner; the number and/or range of design variables can be increased by expanding the data base and developing approximating functions to reflect the expanded design space; the order of the approximating equations can be expanded easily to improve correlation between analyzer results and the approximating equations; gradients of the approximating equations can be calculated easily and these gradients are smooth functions reducing the risk of numerical problems in the optimization; the use of approximating functions allows the problem to be started easily and rapidly from various initial designs to enhance the probability of finding a global optimum; and the approximating equations are independent of the analysis or optimization codes used.

Giansante, N.↗

NASA NTRS1991

An optimal control strategy for collision avoidance of mobile robots in non-stationary environments

An optimal control formulation of the problem of collision avoidance of mobile robots in environments containing moving obstacles is presented. Collision avoidance is guaranteed if the minimum distance between the robot and the objects is nonzero. A nominal trajectory is assumed to be known from off-line planning. The main idea is to change the velocity along the nominal trajectory so that collisions are avoided. Furthermore, time consistency with the nominal plan is desirable. A numerical solution of the optimization problem is obtained. Simulation results verify the value of the proposed strategy.

Kyriakopoulos, K. J.↗

NASA NTRS1972

Aeroelastic optimization of a panel in high Mach number supersonic flow.

Solution for a least-weight skin thickness distribution for a panel with a flutter parameter constraint. This panel weighs less than any similar constant thickness panel, but has the same critical supersonic panel flutter parameter. The panel rests on simple supports and is of sandwich construction. The span to chord ratio is large enough that the inertial, elastic, and aerodynamic behavior is one-dimensional. The Mach number is great enough that the aerodynamic forces acting on the upper panel surface may be accurately described by quasi-steady, linearized, supersonic aerodynamic theory. The final optimum design is obtained from theoretical and numerical methods adapted from optimal control theory. The results of this investigation show that the optimal panel thickness distribution is symmetric about the panel chord midpoint. Compared to a reference panel with constant thickness, optimum panels are found to be nearly 12% lighter.

Weisshaar, T. A.↗

NASA NTRS1990

Direct handling of equality constraints in multilevel optimization

In recent years there have been several hierarchic multilevel optimization algorithms proposed and implemented in design studies. Equality constraints are often imposed between levels in these multilevel optimizations to maintain system and subsystem variable continuity. Equality constraints of this nature will be referred to as coupling equality constraints. In many implementation studies these coupling equality constraints have been handled indirectly. This indirect handling has been accomplished using the coupling equality constraints' explicit functional relations to eliminate design variables (generally at the subsystem level), with the resulting optimization taking place in a reduced design space. In one multilevel optimization study where the coupling equality constraints were handled directly, the researchers encountered numerical difficulties which prevented their multilevel optimization from reaching the same minimum found in conventional single level solutions. The researchers did not explain the exact nature of the numerical difficulties other than to associate them with the direct handling of the coupling equality constraints. The coupling equality constraints are handled directly, by employing the Generalized Reduced Gradient (GRG) method as the optimizer within a multilevel linear decomposition scheme based on the Sobieski hierarchic algorithm. Two engineering design examples are solved using this approach. The results show that the direct handling of coupling equality constraints in a multilevel optimization does not introduce any problems when the GRG method is employed as the internal optimizer. The optimums achieved are comparable to those achieved in single level solutions and in multilevel studies where the equality constraints have been handled indirectly.

Renaud, John E.↗

NASA NTRS1985

An Integrated Control and Minimum Mass Structural Optimization Algorithm for Large Space Structures

A new approach is discussed for solving dual structural control optimization problems for high-order flexible space structures, where reduced-order structural models are employed and minimum mass designs are sought. For a given initial structural design, a quadratic control cost is minimized subject to a constant-mass constraint. The sensitivity of the optimal control cost with respect to the structural design variables is then determined and used to obtain successive structural redesigns, using a constrained gradient optimization algorithm. This process is repeated until the constrained control cost sensitivity becomes negligible. The minimum mass design is obtained by solving a sequence of neighboring optimal constant mass designs, where the sequence of optimal performance indices has a minimum at the optimal minimum mass design. A numerical example is presented which demonstrates that this new approach effectively addresses the problem of dual optimization for potentially very high-order structures.

Messac, A.↗

NASA NTRS1989

Recent experience in simultaneous control-structure optimization

To show the feasibility of simultaneous optimization as design procedure, low order problems were used in conjunction with simple control formulations. The numerical results indicate that simultaneous optimization is not only feasible, but also advantageous. Such advantages come at the expense of introducing complexities beyond those encountered in structure optimization alone, or control optimization alone. Examples include: larger design parameter space, optimization may combine continuous and combinatoric variables, and the combined objective function may be nonconvex. Future extensions to include large order problems, more complex objective functions and constraints, and more sophisticated control formulations will require further research to ensure that the additional complexities do not outweigh the advantages of simultaneous optimization. Some areas requiring more efficient tools than currently available include: multiobjective criteria and nonconvex optimization. Efficient techniques to deal with optimization over combinatoric and continuous variables, and with truncation issues for structure and control parameters of both the model space as well as the design space need to be developed.

Salama, M.↗

NASA NTRS1979

Leading-edge slat optimization for maximum airfoil lift

A numerical procedure for determining the position (horizontal location, vertical location, and deflection) of a leading edge slat that maximizes the lift of multielement airfoils is presented. The structure of the flow field is calculated by iteratively coupling potential flow and boundary layer analysis. This aerodynamic calculation is combined with a constrained function minimization analysis to determine the position of a leading edge slat so that the suction peak on the nose of the main airfoil is minized. The slat position is constrained by the numerical procedure to ensure an attached boundary layer on the upper surface of the slat and to ensure negligible interaction between the slat wake and the boundary layer on the upper surface of the main airfoil. The highest angle attack at which this optimized slat position can maintain attached flow on the main airfoil defines the optimum slat position for maximum lift. The design method is demonstrated for an airfoil equipped with a leading-edge slat and a trailing edge, single-slotted flap. The theoretical results are compared with experimental data, obtained in the Ames 40 by 80 Foot Wind Tunnel, to verify experimentally the predicted slat position for maximum lift. The experimentally optimized slat position is in good agreement with the theoretical prediction, indicating that the theoretical procedure is a feasible design method.

Olson, L. E.↗

NASA NTRS1992

Optimal motion planning for collision avoidance of mobile robots in non-stationary environments

An optimal control formulation of the problem of collision avoidance of mobile robots moving in general terrains containing moving obstacles is presented. A dynamic model of the mobile robot and the dynamic constraints are derived. Collision avoidance is guaranteed if the minimum distance between the robot and the object is nonzero. A nominal trajectory is assumed to be known from off-line planning. The main idea is to change the velocity along the nominal trajectory so that collisions are avoided. Time consistency with the nominal plan is desirable. A numerical solution of the optimization problem is obtained. A perturbation control type of approach is used to update the optimal plan. Simulation results verify the value of the proposed strategy.

Kyriakopoulos, K. J.↗

NASA NTRS1991

Doppler lidar for measurement of atmospheric wind fields

Measurements of wind fields in the earth's troposphere with daily global coverage is widely considered as a significant advance for forecasting and transport studies. For optimal use by NWP (Numerical Weather Prediction) models the horizontal and vertical resolutions should be approximately 100 km and 1 km, respectively. For boundary layer studies vertical resolution of a few hundred meters seems essential. Earth-orbiting Doppler lidar has a unique capability to measure global winds in the troposphere with the high vertical resolution required. The lidar approach depends on transmission of pulses with high spectral purity and backscattering from the atmospheric aerosol particles or layered clouds to provide a return signal. Recent field measurement campaigns using NASA research aircraft have resulted in collection of aerosol and cloud data which can be used to optimize the Doppler lidar instrument design and measurement strategy.

Menzies, Robert T.↗

NASA NTRS1985

Formulation and application of Russell's method

It is shown that the numerical technique of Russell's momentum approach can be derived by using Hamilton's principle and Vance's numerical scheme. It results in a set of first order differnce equations for solving the angular velocities. The numerical examples show that the method is reliable. The algorithm is modified next to perform the analysis of N-body systems with closed loop topology. To increase the formulation flexibility, the equations of motion are represented by using Cartesian coordinates and Lagrange multipliers. The algorithm consists of two parts, Vance's scheme and an unconstrained minimization. The Vance's scheme is used to find the angular velocities, and the unconstrained minimization is applied to provide the correct angular displacements. The proposed scheme is further extended to find the design sensitivity of an N-body system with closed loop configuration, and to carry out the design optimization as well. The numerical example of a small-scaled mechanical system is presented to verify the proposed formulation.

Hou, J. W.↗

NASA NTRS2002

Trajectories for High Specific Impulse High Specific Power Deep Space Exploration

Flight times and deliverable masses for electric and fusion propulsion systems are difficult to approximate. Numerical integration is required for these continuous thrust systems. Many scientists are not equipped with the tools and expertise to conduct interplanetary and interstellar trajectory analysis for their concepts. Several charts plotting the results of well-known trajectory simulation codes were developed and are contained in this paper. These charts illustrate the dependence of time of flight and payload ratio on jet power, initial mass, specific impulse and specific power. These charts are intended to be a tool by which people in the propulsion community can explore the possibilities of their propulsion system concepts. Trajectories were simulated using the tools VARITOP and IPOST. VARITOP is a well known trajectory optimization code that involves numerical integration based on calculus of variations. IPOST has several methods of trajectory simulation; the one used in this paper is Cowell's method for full integration of the equations of motion. An analytical method derived in the companion paper was also evaluated. The accuracy of this method is discussed in the paper.

Polsgrove, Tara↗

NASA NTRS1984

Design of an optimal output feedback control system with modal insensitivity

This paper deals with the design of an output feedback controller which results in selected modal insensitivity, and at the same time optimizes a quadratic performance index representative of desired system performance for nominal plant parameter values. The approach taken here is to characterize the class of attainable eigenvectors for a given set of eigenvalues (distinct or non-distinct) which lie in a subspace called the 'Modal Insensitivity Subspace'. A constraint is established on the feedback matrix which results in modal insensitivity. Necessary conditions for optimality subject to the constraint on the feedback matrix are given. This forms the basis for a numerical algorithm to compute the optimal feedback gain which analyzed for convergence. To illustrate the procedure, a design is carried out using the lateral dynamics of an L-1011 aircraft.

Raman, K. V.↗

NASA NTRS2005

Minimum impulse transfers to rotate the line of apsides

While an optimal scenario for the general two-impulse transfer between coplanar orbits is not known, there are optimal scenarios for various special cases. We consider in-plane rotations of the line of apsides. Numerical comparisons with a trajectory optimization program support the claim that the optimal deltaV required by two impulses is about half that required by a single impulse, regardless of semi-major axes. We observe that this estimate becomes more conservative with larger angles of rotation and eccentricities, and thus also present a more accurate two-impulse rotation deltaV estimator.

orbit rotation↗

DOE OSTI2025

A Scalable Interior‐Point Gauss–Newton Method for PDE‐Constrained Optimization With Bound Constraints

Here, we present a scalable approach to solve a class of partial differential equation (PDE)‐constrained optimization problems with bound constraints. This approach utilizes a robust full‐space interior‐point (IP)‐Gauss–Newton optimization method. To cope with the poorly‐conditioned IP‐Gauss–Newton saddle‐point linear systems that need to be solved approximately, once per optimization step, we propose two spectrally related preconditioners. These preconditioners leverage the limited informativeness of data in regularized PDE‐constrained optimization problems. A block Gauss–Seidel preconditioner is proposed for the GMRES‐based solution of the IP‐Gauss–Newton linear systems. It is shown, for a large‐class of PDE‐ and bound‐constrained optimization problems, that the spectrum of the block Gauss–Seidel preconditioned IP‐Gauss–Newton matrix is asymptotically independent of discretization and is not impacted by the ill‐conditioning that notoriously plagues interior‐point methods. We exploit symmetry of the IP‐Gauss–Newton linear systems and propose a regularization and log‐barrier Hessian preconditioner for the preconditioned conjugate gradient (PCG)‐based solution of the equivalent IP‐Gauss–Newton–Schur complement linear systems. The eigenvalues of the block Gauss–Seidel preconditioned IP‐Gauss–Newton matrix, that are not equal to one, are identical to the eigenvalues of the regularization and log‐barrier Hessian preconditioned Schur complement matrix. The scalability of the approach is demonstrated on two example problems. The numerical solution of these optimization problems is shown to require a discretization independent number of IP‐Gauss–Newton linear solves. Furthermore, the linear systems are solved in a discretization and IP ill‐conditioning independent number of preconditioned Krylov subspace iterations. The parallel scalability of the preconditioner, achieved via algebraic multigrid component solvers when applicable, and the aforementioned algorithmic scalability permits a parallel scalable means to compute solutions of a large class of PDE‐ and bound‐constrained problems.

PDE-constrained optimization↗

DOE OSTI2025

Learning with Adaptive Conservativeness for Distributionally Robust Optimization: Incentive Design for Voltage Regulation

Information asymmetry between the Distribution System Operator (DSO) and Distributed Energy Resource Aggregators (DERAs) obstructs designing effective incentives for voltage regulation. To capture this effect, we employ a Stackelberg game-theoretic framework, where the DSO seeks to overcome the information asymmetry and refine its incentive strategies by learning from DERA behavior over multiple iterations. We introduce a model-based online learning algorithm for the DSO, aimed at inferring the relationship between incentives and DERA responses. Given the uncertain nature of these responses, we also propose a distributionally robust incentive design model to control the probability of voltage regulation failure and then reformulate it into a convex problem. This model allows the DSO to periodically revise distribution assumptions on uncertain parameters in the decision model of the DERA. Finally, we present a gradient-based method that permits the DSO to adaptively modify its conservativeness level, measured by the size of a Wasserstein metric-based ambiguity set, according to historical voltage regulation performance. The effectiveness of our proposed method is demonstrated through numerical experiments.

adaptation models↗

NASA NTRS1992

The role of under-determined approximations in engineering and science application

There is currently a great deal of interest in using response surfaces in the optimization of aircraft performance. The objective function and/or constraint equations involved in these optimization problems may come from numerous disciplines such as structures, aerodynamics, environmental engineering, etc. In each of these disciplines, the mathematical complexity of the governing equations usually dictates that numerical results be obtained from large computer programs such as a finite element method program. Thus, when performing optimization studies, response surfaces are a convenient way of transferring information from the various disciplines to the optimization algorithm as opposed to bringing all the sundry computer programs together in a massive computer code. Response surfaces offer another advantage in the optimization of aircraft structures. A characteristic of these types of optimization problems is that evaluation of the objective function and response equations (referred to as a functional evaluation) can be very expensive in a computational sense. Because of the computational expense in obtaining functional evaluations, the present study was undertaken to investigate under-determinined approximations. An under-determined approximation is one in which there are fewer training pairs (pieces of information about a function) than there are undetermined parameters (coefficients or weights) associated with the approximation. Both polynomial approximations and neural net approximations were examined. Three main example problems were investigated: (1) a function of one design variable was considered; (2) a function of two design variables was considered; and (3) a 35 bar truss with 4 design variables was considered.

Carpenter, William C.↗

NASA NTRS1979

Study and evaluation of impulse mass spectrometers for ion analysis in the D and E regions of the ionosphere

Theoretical and numerical analyses were made of planar, cylindrical and spherical electrode time-of-flight mass spectrometers in order to optimize their operating conditions. A numerical analysis of potential barrier gating in time-of-flight spectrometers was also made. The results were used in the design of several small mass spectrometers. These were constructed and tested in a laboratory space simulator. Detailed experimental studies of a miniature cylindrical electrode time of flight mass spectrometer and of a miniature hemispherical electrode time of flight mass spectrometer were made. The extremely high sensitivity of these instruments and their ability to operate at D region pressures with an open source make them ideal instruments for D region ion composition measurements.

Kendall, B. R.↗