Automatic control
Numerical solution of optimal control and programming problems
FEEDBACK CONTROL SYSTEM↗
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.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
On the use of co-ordinate stretching in the numeral computation of high frequency scattering
The scattering of the sound of a jet engine by an airplane fuselage is modeled by solving the axially symmetric Helmholtz equation exterior to a long thin ellipsoid. The integral equation method based on the single layer potential formulation is used. A family of coordinate systems on the body is introduced and an algorithm is presented to determine the optimal coordinate system. Numerical results verify that the optimal choice enables the solution to be computed with a grid that is coarse relative to the wavelength.
Bayliss, A.↗
Multiabsorber Transition-Edge Sensors for X-Ray Astronomy
We are developing arrays of position-sensitive microcalorimeters for future x-ray astronomy applications. These position-sensitive devices commonly referred to as hydras consist of multiple x-ray absorbers, each with a different thermal coupling to a single-transition-edge sensor microcalorimeter. Their development is motivated by a desire to achieve very large pixel arrays with some modest compromise in performance. We report on the design, optimization, and first results from devices with small pitch pixels (<75 μm) being developed for a high-angular and energy resolution imaging spectrometer for Lynx. The Lynx x-ray space telescope is a flagship mission concept under study for the National Academy of Science 2020 decadal survey. Broadband full-width-half-maximum (FWHM) resolution measurements on a 9-pixel hydra have demonstrated ΔE(FWHM) = 2.23 ± 0.14 eV at Al-Kα, ΔE(FWHM) = 2.44 ± 0.29 eV at Mn-Kα, and ΔE(FWHM) = 3.39 ± 0.23 eV at Cu-Kα. Position discrimination is demonstrated to energies below <1 keV and the device performance is well-described by a finite-element model. Results from a prototype 20-pixel hydra with absorbers on a 50-μm pitch have shown ΔE(FWHM) = 3.38 ± 0.20 eV at Cr-Kα1. We are now optimizing designs specifically for Lynx and extending the number of absorbers up to 25/hydra. Numerical simulation suggests optimized designs could achieve ∼3 eV while being compatible with the bandwidth requirements of the state-of-the art multiplexed readout schemes, thus making a 100,000 pixel microcalorimeter instrument a realistic goal.
Smith, Stephen J.↗
Chebyshev Polynomials Are Not Always Optimal
The authors are concerned with the problem of finding among all polynomials of degree at most n and normalized to be 1 at c the one with minimal uniform norm on Epsilon. Here, Epsilon is a given ellipse with both foci on the real axis and c is a given real point not contained in Epsilon. Problems of this type arise in certain iterative matrix computations, and, in this context, it is generally believed and widely referenced that suitably normalized Chebyshev polynomials are optimal for such constrained approximation problems. In this note, the authors show that this is not true in general. Moreover, the authors derive sufficient conditions which guarantee that Chebyshev polynomials are optimal. Also, some numerical examples are presented.
CHEBYSHEV APPROXIMATION↗
Numerical Methodology For Full-Body TPS Sizing And Optimization For Access-To-Space Vehicles
This paper presents details of the methodology and numerical procedures developed at NASA Ames for full-body TEXAS sizing and optimization for Access to Space vehicle concepts. The core of the procedures is a robust implicit solver for one dimensional transient heat conduction in reusable multilayer TEXAS stackups. The solver includes an arbitrary number of material layers, contact resistances between materials, temperature and pressure dependent material and surface properties, numerous boundary-condition options, and self-adaptive time stepping. The solver is coupled with the Access-to-Space material database of 23 candidate TPS and structural materials and a thermal-environment database obtained from trajectory-based fullbody Navier-Stokes computations of the external flowfield. The thermal environment and material response are coupled through the use of T-type heat transfer coefficients. TIPS sizing and weight optimization are performed at every surface point on the vehicle based on sizing constraints which include material temperature limits, maximum backwall temperature, and cumulative interior heat flux. Typical results are presented for a lifting body concept with 10000 surface points, which required 35 minutes to compute on an SGI Indigo 2.
Milos, Frank S.↗
Development of an adaptive hp-version finite element method for computational optimal control
In this research effort, the usefulness of hp-version finite elements and adaptive solution-refinement techniques in generating numerical solutions to optimal control problems has been investigated. Under NAG-939, a general FORTRAN code was developed which approximated solutions to optimal control problems with control constraints and state constraints. Within that methodology, to get high-order accuracy in solutions, the finite element mesh would have to be refined repeatedly through bisection of the entire mesh in a given phase. In the current research effort, the order of the shape functions in each element has been made a variable, giving more flexibility in error reduction and smoothing. Similarly, individual elements can each be subdivided into many pieces, depending on the local error indicator, while other parts of the mesh remain coarsely discretized. The problem remains to reduce and smooth the error while still keeping computational effort reasonable enough to calculate time histories in a short enough time for on-board applications.
Hodges, Dewey H.↗
Parametric study of optimized liquid-hydrogen thermal protection systems for nuclear interplanetary spacecraft. Volume 3 - Numerical data Final report
Data tables for optimal thermal protection of liquid hydrogen storage in nuclear interplanetary spacecraft
Barry, D. G.↗
Structure-preserving neural networks for the regularized entropy-based closure of a linear, kinetic, radiative transport equation
The main challenge of large-scale numerical simulation of radiation transport is the high memory and computation time requirements of discretization methods for kinetic equations. In this work, we derive and investigate a neural network-based approximation to the entropy-based closure method to accurately compute the solution of the multi-dimensional moment system with a low memory footprint and competitive computational time. We extend methods developed for the standard entropy-based closure to the regularized entropy-based closures. The main idea is to interpret structure-preserving neural network approximations of the regularized entropy-based closure as a two-stage approximation to the original entropy-based closure. We conduct a numerical analysis of this approximation and investigate optimal parameter choices. Our numerical experiments demonstrate that the method has a much lower memory footprint than traditional methods with competitive computation times and simulation accuracy. The code and all trained networks are provided on GitHub.
entropy closure↗
Development of an hp-version finite element method for computational optimal control
The purpose of this research effort was to begin the study of the application of hp-version finite elements to the numerical solution of optimal control problems. Under NAG-939, the hybrid MACSYMA/FORTRAN code GENCODE was developed which utilized h-version finite elements to successfully approximate solutions to a wide class of optimal control problems. In that code the means for improvement of the solution was the refinement of the time-discretization mesh. With the extension to hp-version finite elements, the degrees of freedom include both nodal values and extra interior values associated with the unknown states, co-states, and controls, the number of which depends on the order of the shape functions in each element. One possible drawback is the increased computational effort within each element required in implementing hp-version finite elements. We are trying to determine whether this computational effort is sufficiently offset by the reduction in the number of time elements used and improved Newton-Raphson convergence so as to be useful in solving optimal control problems in real time. Because certain of the element interior unknowns can be eliminated at the element level by solving a small set of nonlinear algebraic equations in which the nodal values are taken as given, the scheme may turn out to be especially powerful in a parallel computing environment. A different processor could be assigned to each element. The number of processors, strictly speaking, is not required to be any larger than the number of sub-regions which are free of discontinuities of any kind.
Hodges, Dewey H.↗
Computer simulation for time optimal or energy optimal attitude control of spin-stabilized spacecraft.
An original technique for determining the optimal magnetic torque strategy for control of the attitude of spin stabilized spacecraft is presented. By employing Lagrange multipliers and the Calculus of Variations, optimal control equations are derived which define minimum time and minimum energy attitude maneuvers. Computer program algorithms to numerically solve these optimal control equations are also described. The performance of this technique is compared with a commonly employed planning method.
Woolley, R. D.↗
Development of an hp-version finite element method for computational optimal control
The purpose of this research effort is to develop a means to use, and to ultimately implement, hp-version finite elements in the numerical solution of optimal control problems. The hybrid MACSYMA/FORTRAN code GENCODE was developed which utilized h-version finite elements to successfully approximate solutions to a wide class of optimal control problems. In that code the means for improvement of the solution was the refinement of the time-discretization mesh. With the extension to hp-version finite elements, the degrees of freedom include both nodal values and extra interior values associated with the unknown states, co-states, and controls, the number of which depends on the order of the shape functions in each element.
Hodges, Dewey H.↗
Numerical simulation of the helium gas spin-up channel performance of the relativity gyroscope
The dependence of the spin-up system efficiency on each geometrical parameter of the spin-up channel and the exhaust passage of the Gravity Probe-B (GPB) is individually investigated. The spin-up model is coded into a computer program which simulates the spin-up process. Numerical results reveal optimal combinations of the geometrical parameters for the ultimate spin-up performance. Comparisons are also made between the numerical results and experimental data. The experimental leakage rate can only be reached when the gap between the channel lip and the rotor surface increases beyond physical limit. The computed rotating frequency is roughly twice as high as the measured ones although the spin-up torques fairly match.
Karr, Gerald R.↗
Minimum fuel spacecraft reorientation
Fuel optimal solutions for the reorientation of an inertially symmetric rigid spacecraft with independent 3-axes controls are investigated. All possible optimal control strategies are identified. These include bangbang solutions, finite order singular arcs and infinite order singular arcs. Higher order necessary conditions for optimality of finite order singular arcs are presented. Numerical examples of fuel optimal solutions with fixed maneuver time are presented involving all of the theoretically possible control logics.
Seywald, H.↗
Design optimization of lightweight automotive seatback through additive manufacturing compression overmolding of metal polymer composites
With the growing demand for enhanced automotive fuel efficiency and environmental sustainability, there is a need for lightweighting automotive components through innovative design and manufacturing processes. Here, this study leverages a combination of numerical iterative design optimization and hybrid additive manufacturing–compression molding (AM-CM) technique for metal polymer composites to lightweight an automotive seatback. The AM-CM process enables robust mechanical interlocking between metals and composites, boasting high stiffness and strength with low overall density. Replacing metallic components with such metal polymer composites allows for comparable mechanical performance while significantly reducing the overall weight. First, the automotive seatback design space is reduced to critical load carrying regions using topology optimization and high stress concentration areas are identified using finite element analysis. Next, a lightweight metal polymer subcomponent is designed for a high stress concentration region. The full seatback frame with spatially heterogeneous material-specific design is then iteratively optimized to enable enhanced stiffness with minimal weight. Overall, the automotive seatback frame designed with location-specific metal, polymer, and metal polymer composite materials weighs 20% less than the metal-only design while exhibiting similar stiffness.
36 MATERIALS SCIENCE↗
Comparison of polynomial approximations and artificial neural nets for response surfaces in engineering optimization
Engineering optimization problems involve minimizing some function subject to constraints. In areas such as aircraft optimization, the constraint equations may be from numerous disciplines such as transfer of information between these disciplines and the optimization algorithm. They are also suited to problems which may require numerous re-optimizations such as in multi-objective function optimization or to problems where the design space contains numerous local minima, thus requiring repeated optimizations from different initial designs. Their use has been limited, however, by the fact that development of response surfaces randomly selected or preselected points in the design space. Thus, they have been thought to be inefficient compared to algorithms to the optimum solution. A development has taken place in the last several years which may effect the desirability of using response surfaces. It may be possible that artificial neural nets are more efficient in developing response surfaces than polynomial approximations which have been used in the past. This development is the concern of the work.
Carpenter, William C.↗
Global optimality of extremals: An example
The question of the existence and location of Darboux points is crucial for minimally sufficient conditions for global optimality and for computation of optimal trajectories. A numerical investigation is presented of the Darboux points and their relationship with conjugate points for a problem of minimum fuel, constant velocity, and horizontal aircraft turns to capture a line. This simple second order optimal control problem shows that ignoring the possible existence of Darboux points may play havoc with the computation of optimal trajectories.
Kreindler, E.↗
Computational experiments in the optimal slewing of flexible structures
Numerical experiments on the problem of moving a flexible beam are discussed. An optimal control problem is formulated and transcribed into a form which can be solved using semi-infinite optimization techniques. All experiments were carried out on a SUN 3 microcomputer.
Baker, T. E.↗
Optimization and analysis of large chemical kinetic mechanisms using the solution mapping method - Combustion of methane
A method of systematic optimization, solution mapping, as applied to a large-scale dynamic model is presented. The basis of the technique is parameterization of model responses in terms of model parameters by simple algebraic expressions. These expressions are obtained by computer experiments arranged in a factorial design. The developed parameterized responses are then used in a joint multiparameter multidata-set optimization. A brief review of the mathematical background of the technique is given. The concept of active parameters is discussed. The technique is applied to determine an optimum set of parameters for a methane combustion mechanism. Five independent responses - comprising ignition delay times, pre-ignition methyl radical concentration profiles, and laminar premixed flame velocities - were optimized with respect to thirteen reaction rate parameters. The numerical predictions of the optimized model are compared to those computed with several recent literature mechanisms. The utility of the solution mapping technique in situations where the optimum is not unique is also demonstrated.
Frenklach, Michael↗