Search NASA⌕ Search

SEARCH · Search NASA

Results for “solution 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 451 records · Page 25

A hybrid approach to near-optimal launch vehicle guidance

This paper evaluates a proposed hybrid analytical/numerical approach to launch-vehicle guidance for ascent to orbit injection. The feedback-guidance approach is based on a piecewise nearly analytic zero-order solution evaluated using a collocation method. The zero-order solution is then improved through a regular perturbation analysis, wherein the neglected dynamics are corrected in the first-order term. For real-time implementation, the guidance approach requires solving a set of small dimension nonlinear algebraic equations and performing quadrature. Assessment of performance and reliability are carried out through closed-loop simulation for a vertically launched 2-stage heavy-lift capacity vehicle to a low earth orbit. The solutions are compared with optimal solutions generated from a multiple shooting code. In the example the guidance approach delivers over 99.9 percent of optimal performance and terminal constraint accuracy.

Leung, Martin S. K.↗

A Physics-Based Work-Energy Formulation for Real-Time Trajectory Guidance of A Lunar Lander

Throughout the years, many researchers have calculated and optimized trajectory solutions for lunar landing systems by employing sophisticated mathematical methods, that include: Hamilton’s Principle of Variation, Pontryagin’s maximum principle, and well known convex-optimization techniques among others. Many of these approaches typically require expensive computational resources to achieve convergence in the solution. In an effort to reduce complexity and the computational load required to obtain real-time guidance commands, a simple physics-based work-energy approach has been formulated. This approach is based on the dissipation of the mechanical energy of the vehicle to its final desired energy state required to achieve a safe landing. The rocket engine(s) employed during landing (among other maneuvers) dissipates mechanical energy by both doing work against the velocity vector of the vehicle (thus defining the trajectory path), and by jettisoning mass. Therefore, by solving the energy dissipation problem at every step of the maneuver, a much simpler formulation that naturally and quickly attains convergence is obtained. This formulation is not limited to approach, landing, and divert maneuvers, but in principle it can be employed during de-orbiting, braking burn, ascent, as well as orbit insertion.

Guidance↗

A Physics-Based Work-Energy Formulation for Real-Time Trajectory Guidance of a Lunar Lander

Throughout the years, many researchers have calculated and optimized trajectory solutions for lunar landing systems by employing sophisticated mathematical methods, that include: Hamilton’s Principle of Variation, Pontryagin’s maximum principle, and well known convex-optimization techniques among others. Many of these approaches typically require expensive computational resources to achieve convergence in the solution. In an effort to reduce complexity and the computational load required to obtain real-time guidance commands, a simple physics-based work-energy approach has been formulated. This approach is based on the dissipation of the mechanical energy of the vehicle to its final desired energy state required to achieve a safe landing. The rocket engine(s) employed during landing (among other maneuvers) dissipates mechanical energy by both doing work against the velocity vector of the vehicle (thus defining the trajectory path), and by jettisoning mass. Therefore, by solving the energy dissipation problem at every step of the maneuver, a much simpler formulation that naturally and quickly attains convergence is obtained. This formulation is not limited to approach, landing, and divert maneuvers, but in principle it can be employed during de-orbiting, braking burn, ascent, as well as orbit insertion.

Guidance↗

Affine Generalized Inverse for Optimal Control Allocation

This research is a follow on to the "Optimal Control Prediction Method for Control Allocation" paper in which the Prediction Method iterative algorithm was introduced. Previously, the Prediction Method was shown to provide optimal control allocation solutions over the entire Attainable Moment Set for the Moore-Penrose and the generalized (weighted) inverse. As an extension to the Prediction Method, this paper introduces a family of Moore Penrose Affine Generalized Inverses, applicable for all moments, which compute control allocation solutions using a constant matrix and fixed null-space vector. The Moore-Penrose Affine Generalized Inverse is proven to yield equivalent solutions to those of the Prediction Method and therefore is guaranteed to yield Moore-Penrose optimal control allocation solutions. While the Prediction Method is applicable for any moment along an a priori specified moment direction, the Affine Generalized Inverse is shown to yield optimal control allocation solutions in a neighborhood of the given moment which is not restricted to a specified moment direction. Furthermore, the Affine Generalized Inverse is shown to provide the time derivative of optimal control allocation solutions and to facilitate maintaining solutions within control effector rate limitations. The Moore-Penrose Affine Generalized Inverse is broadened to encompass any arbitrary (weighted) Affine Generalized Inverse. Finally, a method of creating a moment lookup table is outlined to utilize the Affine Generalized Inverse as an offline control allocation solution for all moments in the Attainable Moment Set.

Acheson, Michael J.↗

Determination of optimal control using imbedding of the terminal conditions

The utility of a method to solve optimal control problems is discussed with emphasis placed on the solutions of singular optimal control problems. The method involves imbedding the optimization problem of interest in a family of optimization problems which are parameterized by their terminal conditions. The solution to the optimization problem is obtained by collapsing the terminal conditions of the family of problems onto those of the original problem, while appropriately modifying the control function to account for the changing terminal conditions. The utility of this method is illustrated by obtaining the solution to an aerial attack problem where the attacking vehicle must maneuver subject to control limitations so as to place the enemy in a certain relative position for launching a missile in the shortest possible time.

Montgomery, R. C.↗

Numerical analysis of discrete fractional integrodifferential structural dampers

This paper develops solution algorithms enabling the handling of the dynamic response of nonlinear structures contained discretely attached dampers modeled by fractional integrodifferential operators of the Grunwald-Liouville-Riemann type. The development consists of two levels of formulation, namely: (1) numerical approximations of fractional operators and, (2) the establishment of global level implicit schemes enabling the solution to nonlinear structural formulations. To generalize the overall results, error estimates are derived for the fractional operator approximation algorithm. These enable an ongoing optimization of solution efficiency for a given error tolerance. To benchmark the scheme, the results of several numerical experiments are presented. These illustrate the numerical characteristics of the overall formulation.

Padovan, Joseph↗

A complete analytical solution for the inverse instantaneous kinematics of a spherical-revolute-spherical (7R) redundant manipulator

Using a method based upon resolving joint velocities using reciprocal screw quantities, compact analytical expressions are generated for the inverse solution of the joint rates of a seven revolute (spherical-revolute-spherical) manipulator. The method uses a sequential decomposition of screw coordinates to identify reciprocal screw quantities used in the resolution of a particular joint rate solution, and also to identify a Jacobian null-space basis used for the direct solution of optimal joint rates. The results of the screw decomposition are used to study special configurations of the manipulator, generating expressions for the inverse velocity solution for all non-singular configurations of the manipulator, and identifying singular configurations and their characteristics. Two functions are therefore served: a new general method for the solution of the inverse velocity problem is presented; and complete analytical expressions are derived for the resolution of the joint rates of a seven degree of freedom manipulator useful for telerobotic and industrial robotic application.

Podhorodeski, R. P.↗

Modeling of Multicomponent Trace Contaminant Adsorption on Activated Carbon

Modeling of multicomponent trace contaminant adsorption on activated carbon for spacecraft life support applications presents challenges due to the very low contaminant concentrations in cabin air coupled with very strong adsorption of heavier components such as cyclic polydimethylsiloxanes. Both gas-phase and solid-phase micropore mass transfer resistances can be important, and competitive interactions between adsorbed components can result in behavior such as roll-up where lighter components displaced by heavier components have a higher concentration at the outlet than at the inlet. This paper describes the further development of a dynamic trace contaminant adsorption model based on Ideal Adsorbed Solution Theory (IAST). Alternative 1-dimensional and 2-dimensional versions of the model are described with different mass transfer resistance and driving force assumptions. Optimization of solution methods for speed and stability is also described. Mass balance equations are scaled to provide sensitivity over the wide range of component concentrations in cabin air. Model parameters are estimated from isotherm and breakthrough test data on single adsorbate and binary adsorbate systems. Multicomponent results are compared against predictions of a heritage model used extensively by NASA for trace contaminant control system (TCCS) carbon bed sizing. Finally, efforts to implement the model in a second software platform for wider access are described.

Kevin E Lange↗

Machine Learning for Dynamic Test Sensor Placement

There are multiple different algorithms to perform modal test sensor placement optimization: effective independence, residual kinetic energy, iterative Guyan reduction, genetic algorithms, or a brute-force methodology. However, any of these methods may be computationally expensive, especially for structural models with a large number of degrees of freedom. Given the high-cost and the need to optimize the solution, modal sensor placement is a great application for machine learning (ML) algorithms. In this paper, we will apply ML algorithms to determine the optimal sensor locations for simple and complex structures. We will also discuss the benefits and drawbacks of using machine learning over other sensor placement algorithms.

Kelsey Buckles↗

Efficient single-level solution of hierarchical problems in structural optimization

A demonstration is made of the use of single-level hierarchical solution techniques for engineering design, which, unlike the formulation stage-applied multilevel techniques, are applied at the solution stage. A single-level decomposition technique is developed which reduces computational costs and memory requirements without incurring the disadvantages of multilevel optimization. Savings of 75 percent are obtained in a case involving 500 design variables and 2408 constraints.

Thareja, Rajiv R.↗

Optimal boost trajectories for the shuttle vehicle.

Solutions of the optimal boost problem for the MSC April 1970 base line shuttle vehicle show significant payload increases may be achieved by flying optimal pitch profiles. Results consistent with aerodynamic and acceleration constraints indicate payload gains of approximately 2.0 per cent of the total injected weight may be achieved. This corresponds to a real payload increase of about 20 per cent. Trajectory parameters for optimal launches under various conditions and a standard gravity turn launch are presented along with a brief discussion of the problem modeling and assumptions. The data are used to support a number of observations and conclusions and to provide some insight into the methods by which payload gains are achieved.

Sullivan, H. C.↗

Dual approximations in optimal control

A dual approximation for the solution to an optimal control problem is analyzed. The differential equation is handled with a Lagrange multiplier while other constraints are treated explicitly. An algorithm for solving the dual problem is presented.

Hager, W. W.↗

Combined optimal control and estimation.

Combined optimization problem, equivalent to dual control problem, considering determination of optimal control policies for plant under random disturbances, using iterative equations

CONTROL SYSTEM↗

A linear stochastic model of the human operator

A linear stochastic model of the human operator is developed and applied to the problem of piloted control of an aircraft. The pilot and aircraft are modeled as linear time-invariant systems containing both process and measurement noise. The loop closure by the pilot is determined by formulating the problem as an optimal stochastic control problem. The solution to the optimal control problem yields not only the pilot's optimal control output which he uses to control the vehicle, but also the optimal combination of his observations of the vehicle states upon which the pilot bases his control. A method is presented so that, using experimental pilot vehicle data, the cost functional which is minimized in the optimal control problem will be numerically equal to the pilot rating that the pilot would associate with the given vehicle and task.

Durrett, J. C.↗

On optimizing the treatment of exchange perturbations

A method using the zeroth plus first order wave functions, obtained by optimizing the basic equation used in exchange perturbation treatments, is utilized in an attempt to determine the exact energy and wave function in the exchange process. Attempts to determine the first order perturbation solution by optimizing the sum of the first and second order energies were unsuccessful.

Hirschfelder, J. O.↗

Experiments on high speed ejectors

Experimental studies were conducted to investigate the flow and the performance of thrust augmenting ejectors for flight Mach numbers in the range of 0.5 to 0.8, primary air stagnation pressures up to 107 psig (738 kPa), and primary air stagnation temperatures up to 1250 F (677 C). The experiment verified the existence of the second solution ejector flow, where the flow after complete mixing is supersonic. Thrust augmentation in excess of 1.2 was demonstrated for both hot and cold primary jets. The experimental ejector performed better than the corresponding theoretical optimal first solution ejector, where the mixed flow is subsonic. Further studies are required to realize the full potential of the second solution ejector. The research program was started by the Flight Dynamics Research Corporation (FDRC) to investigate the characteristic of a high speed ejector which augments thrust of a jet at high flight speeds.

Wu, J. J.↗

Automated Generation of Finite-Element Meshes for Aircraft Conceptual Design

This paper presents a novel approach for automated generation of fully connected finite-element meshes for all internal structural components and skins of a given wing-body geometry model, controlled by a few conceptual-level structural layout parameters. Internal structural components include spars, ribs, frames, and bulkheads. Structural layout parameters include spar/rib locations in wing chordwise/spanwise direction and frame/bulkhead locations in longitudinal direction. A simple shell thickness optimization problem with two load conditions is used to verify versatility and robustness of the automated meshing process. The automation process is implemented in ModelCenter starting from an OpenVSP geometry and ending with a NASTRAN 200 solution. One subsonic configuration and one supersonic configuration are used for numerical verification. Two different structural layouts are constructed for each configuration and five finite-element meshes of different sizes are generated for each layout. The paper includes various comparisons of solutions of 20 thickness optimization problems, as well as discussions on how the optimal solutions are affected by the stress constraint bound and the initial guess of design variables.

Li, Wu↗