Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical optimal control”

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 667 records · Page 37

Optimal control of spin-stabilized spacecraft with telescoping appendages

The control of a spin-stabilized spacecraft consisting of a rigid central hub and one or two movable offset telescoping booms (with end masses) is considered. The equations of rotational motion are linearized about either of two desired final states. A control law for the boom end mass position is sought such that a quadratic cost functional involving the weighted components of angular velocity plus the control is minimized when the final time is unspecified and involves the solution of the matrix Riccati algebraic equation. For three axis control more than one offset boom (orthogonal to each other) is required. For two-axis control with a single boom offset from a symmetrical hub, an analytic solution is obtained; when this system is used for nutation decay the time constant is one order of magnitude smaller than previously achieved using nonoptimal control logic. For the general case results are obtained numerically.

Bainum, P. M.↗

Evolving Design Criteria for Very Large Aperture Space-Based Telescopes and Their Influence on the Need for Integrated Tools in the Optimization Process

NASA's Advanced Mirror Technology Development (AMTD) program has been developing the means to design and build the future generations of space based telescopes. With the nearing completion of the James Webb Space Telescope (JWST), the astrophysics community is already starting to define the requirements for follow on observatories. The restrictions of available launch vehicles and the possibilities of planned future vehicles have fueled the competition between monolithic primaries (with better optical quality) and segmented primaries (with larger apertures, but with diffraction, costs and figure control issues). Regardless of the current shroud sizes and lift capacities, these competing architectures share the need for rapid design tools. As part of the AMTD program a number of tools have been developed and tested to speed up the design process. Starting with the Arnold Mirror Modeler (which creates Finite Element Models (FEM) for structural analysis) and now also feeds these models into thermal stability analyses. They share common file formats and interchangeable results. During the development of the program, numerous trade studies were created for 4 meter and 8 meter monolithic primaries, complete with support systems. Evaluation of these results has led to a better understanding of how the specification drives the results. This paper will show some of the early trade studies for typical specification requirements such as lowest mirror bending frequency and suspension system lowest frequency. The results use representative allowable stress values for each mirror substrate material and construction method and generic material properties. These studies lead to some interesting relationships between feasible designs and the realities of actually trying to build these mirrors. Much of the traditional specifications were developed for much smaller systems, where the mass and volume of the primary where a small portion of the overall satellite. JWST shows us that as the aperture grows, the primary takes up the majority of the mass and volume and the established rules need to be adjusted. For example, a small change in lowest frequency requirement can change the cost by millions of dollars.

Space-based telescopes↗

A man in the loop trajectory optimization program (MILTOP)

An interactive trajectory optimization program is developed for use in initial fixing of launch configurations. The program is called MILTOP for Man-In-the-Loop-Trajectory Optimization-Program. The program is designed to facilitate quick look studies using man-machine decision combinations to reduce the time required to solve a given problem. MILTOP integrates the equations of motion of a point-mass in 3-Dimensions with drag as the only aerodynamic force present. Any point in time at which an integration step terminates, may be used as a decision-break-point, with complete user control over all variables and routines at this point. Automatic phases are provided for different modes of control: vertical rise, pitch-over, gravity turn, chi-freeze and control turn. Stage parameters are initialized from a separate routine so the user may fly as many stages as his problem demands. The MILTOP system uses both interactively on storage scope consoles, or in batch mode with numerical output on the live printer.

Reinfields, J.↗

Evolving Design Criteria for Very Large Aperture Space Based Telescopes and Their Influence on the Need for Integrated Tools in the Optimization Process

NASA's Advanced Mirror Technology Development (AMTD) program has been developing the means to design and build the future generations of space based telescopes. With the nearing completion of the James Webb Space Telescope (JWST), the astrophysics community is already starting to define the requirements for follow-on observatories. The restrictions of available launch vehicles and the possibilities of planned future vehicles have fueled the competition between monolithic primaries (with better optical quality) and segmented primaries (with larger apertures, but with diffraction, costs and figure control issues). Regardless of the current shroud sizes and lift capacities, these competing architectures share the need for rapid design tools. As part of the AMTD program a number of tools have been developed and tested to speed up the design process. Starting with the Arnold Mirror Modeler (which creates Finite Element Models (FEM) for structural analysis) and now also feeds these models into thermal stability analyses. They share common file formats and interchangeable results. During the development of the program, numerous trade studies were created for 4-meter and 8-meter monolithic primaries, complete with support systems. Evaluation of these results has led to a better understanding of how the specification drives the results. This paper will show some of the early trade studies for typical specification requirements such as lowest mirror bending frequency and suspension system lowest frequency. The results use representative allowable stress values for each mirror substrate material and construction method and generic material properties. These studies lead to some interesting relationships between feasible designs and the realities of actually trying to build these mirrors. Much of the traditional specifications were developed for much smaller systems, where the mass and volume of the primary where a small portion of the overall satellite. JWST shows us that as the aperture grows, the primary takes up the majority of the mass and volume and the established rules need to be adjusted. For example, a small change in lowest frequency requirement can change the cost by millions of dollars. The paper uses numerous trade studies created during the software development phase of the Arnold Mirror Modeler to illustrate the influences of system specifications on the design space. The future telescopes will require better performance, stability and documented feasibility to meet the hurdles of today's budget and schedules realities. AMTD is developing the tools, but the basic system planning mentality also has to adopt to the requirements of these very large and complex physical structures.

Arnold, William R., Sr.↗

A Study into Validating a Coupled Method of Characteristics and Direct Simulation Monte Carlo Method Against Empirical Data

The following will outline the methodology and results of validating a coupled Method of Characteristics (MOC) and Direct Simulation Monte Carlo (DSMC) method. This research focused specifically on modeling plume impingement, induced by Reaction Control System (RCS) thrusters that flew on the National Aeronautics and Space Administration’s (NASA’s) space shuttle Discovery. For each simulation, the continuum portion of the RCS thruster was simulated using MOC for solving hyperbolic Partial Differential Equations (PDEs) and computed with the NASA code, Reacting and Multi-phase Program (RAMP). The solution was then implemented as a starting condition into the NASA DSMC code, Direct Simulation and Monte Carlo Analysis Code (DAC). Typically, DSMC models rely on code-to-code validation for fidelity. The significance of this research is in its ability to validate its models against empirical data. Prior to computing solutions for these simulations, the mesh size and structure were optimized and many variants of DSMC input parameters were iterated on in order to acquire a reliable, mesh-independent, fully optimized numerical solution. This research will discuss the mathematical formulation of MOC for nozzle flow and DSMC for rarefied gases. Additionally, it will provide an explanation of how to implement these mathematical concepts into the two solvers: RAMP and DAC. Ultimately, this research will demonstrate that the overall process illustrated produces results in good agreement with empirical data. As a consequence, the methodology presented is granted an increased level of confidence and will greatly contribute to the aerospace industry and its effort in understanding and predicting rarefied flow fields.

Direct Simulation Monte Carlo Analysis Code↗

A Study into Validating A Coupled Method of Characteristics And Direct Simulation Monte Carlo Method Against Empirical Data

The following will outline the methodology and results of validating a coupled Method of Characteristics (MOC) and Direct Simulation Monte Carlo (DSMC) method. This research focused specifically on modeling plume impingement, induced by Reaction Control System (RCS) thrusters that flew on the National Aeronautics and Space Administration’s (NASA’s) space shuttle Discovery. For each simulation, the continuum portion of the RCS thruster was simulated using MOC for solving hyperbolic Partial Differential Equations (PDEs) and computed with the NASA code, Reacting and Multi-phase Program (RAMP). The solution was then implemented as a starting condition into the NASA DSMC code, Direct Simulation and Monte Carlo Analysis Code (DAC). Typically, DSMC models rely on code-to-code validation for fidelity. The significance of this research is in its ability to validate its models against empirical data. Prior to computing solutions for these simulations, the mesh size and structure were optimized and many variants of DSMC input parameters were iterated on in order to acquire a reliable, mesh-independent, fully optimized numerical solution. This research will discuss the mathematical formulation of MOC for nozzle flow and DSMC for rarefied gases. Additionally, it will provide an explanation of how to implement these mathematical concepts into the two solvers: RAMP and DAC. Ultimately, this research will demonstrate that the overall process illustrated produces results in good agreement with empirical data. As a consequence, the methodology presented is granted an increased level of confidence and will greatly contribute to the aerospace industry and its effort in understanding and predicting rarefied flow fields.

DAC↗

A Simple Two Aircraft Conflict Resolution Algorithm

Conflict detection and resolution methods are crucial for distributed air-ground traffic management in which the crew in, the cockpit, dispatchers in operation control centers sad and traffic controllers in the ground-based air traffic management facilities share information and participate in the traffic flow and traffic control functions. This paper describes a conflict detection, and a conflict resolution method. The conflict detection method predicts the minimum separation and the time-to-go to the closest point of approach by assuming that both the aircraft will continue to fly at their current speeds along their current headings. The conflict resolution method described here is motivated by the proportional navigation algorithm, which is often used for missile guidance during the terminal phase. It generates speed and heading commands to rotate the line-of-sight either clockwise or counter-clockwise for conflict resolution. Once the aircraft achieve a positive range-rate and no further conflict is predicted, the algorithm generates heading commands to turn back the aircraft to their nominal trajectories. The speed commands are set to the optimal pre-resolution speeds. Six numerical examples are presented to demonstrate the conflict detection, and the conflict resolution methods.

Chatterji, Gano B.↗

Optimal guidance and propulsion control for transatmospheric vehicles

Problems associated with on-board trajectory optimization and with the synthesis of guidance laws are addressed for ascent to LEO of an air-breathing, single-stage-to-orbit vehicle. A multimode propulsion system is assumed which incorporates turbojet, ramjet, scramjet, and rocket engines. An energy-state approximation is applied to a four-state dynamic model for flight of a point mass over a spherical nonrotating earth. An algorithm for generating fuel-optimal climb profiles is derived via singular perturbation theory. This algorithm results from application of the minimum principle to a low-order dynamic model that includes general functional dependence on angle of attack and a normal component of thrust. Switching conditions are derived which, under appropriate assumptions, govern optimal transition from one propulsion mode to another. The use of bank angle to modulate the magnitude of the vertical component of lift is shown to improve the index performance. Numerical results illustrate the nature of the resulting fuel-optimal climb paths.

Corban, J. E.↗

Integrated Aerodynamic and Control System Design of Oblique Wing Aircraft

An efficient high speed aircraft design must achieve a high lift to drag ratio at transonic and supersonic speeds. In 1952 Dr. R. T. Jones proved that for any flight Mach number minimum drag at a fixed lift is achieved by an elliptic wing planform with an appropriate oblique sweep angle. Since then, wind tunnel tests and numerical flow models have confirmed that the compressibility drag of oblique wing aircraft is lower than similar symmetrical sweep designs. At oblique sweep angles above thirty degrees the highly asymmetric planform gives rise to aerodynamic and inertia couplings which affect stability and degrade the aircraft's handling qualities. In the case of the NASA-Rockwell Oblique Wing Research Aircraft, attempts to improve the handling qualities by implementing a stability augmentation system have produced unsatisfactory results because of an inherent lack of controllability in the proposed design. The present work focuses on improving the handling qualities of oblique wing aircraft by including aerodynamic configuration parameters as variables in the control system synthesis to provide additional degrees of freedom with which to further decouple the aircraft's response. Handling qualities are measured using a quadratic cost function identical to that considered in optimal control problems, but the controller architecture is not restricted to full state feedback. An optimization procedure is used to simultaneously solve for the aircraft configuration and control gains which maximize a handling qualities measure, while meeting imposed constraints on trim. In some designs wing flexibility is also modeled and reduced order controllers are implemented. Oblique wing aircraft synthesized by this integrated design method show significant improvement in handling qualities when compared to the originally proposed closed loop aircraft. The integrated design synthesis method is then extended to show how handling qualities may be traded for other types of mission performance (drag, weight, etc.). Examples are presented which show how performance can be maximized while maintaining a desired level of handling quality.

Morris, Stephen James↗

Automatically Finding the Control Variables for Complex System Behavior

Testing large-scale systems is expensive in terms of both time and money. Running simulations early in the process is a proven method of finding the design faults likely to lead to critical system failures, but determining the exact cause of those errors is still time-consuming and requires access to a limited number of domain experts. It is desirable to find an automated method that explores the large number of combinations and is able to isolate likely fault points. Treatment learning is a subset of minimal contrast-set learning that, rather than classifying data into distinct categories, focuses on finding the unique factors that lead to a particular classification. That is, they find the smallest change to the data that causes the largest change in the class distribution. These treatments, when imposed, are able to identify the factors most likely to cause a mission-critical failure. The goal of this research is to comparatively assess treatment learning against state-of-the-art numerical optimization techniques. To achieve this, this paper benchmarks the TAR3 and TAR4.1 treatment learners against optimization techniques across three complex systems, including two projects from the Robust Software Engineering (RSE) group within the National Aeronautics and Space Administration (NASA) Ames Research Center. The results clearly show that treatment learning is both faster and more accurate than traditional optimization methods.

Gay, Gregory↗

System identification methods for aircraft flight control development and validation

System-identification methods compose a mathematical model, or series of models, from measurements of inputs and outputs of dynamic systems. The extracted models allow the characterization of the response of the overall aircraft or component subsystem behavior, such as actuators and on-board signal processing algorithms. This paper discusses the use of frequency-domain system-identification methods for the development and integration of aircraft flight-control systems. The extraction and analysis of models of varying complexity from nonparametric frequency-responses to transfer-functions and high-order state-space representations is illustrated using the Comprehensive Identification from FrEquency Responses (CIFER) system-identification facility. Results are presented for test data of numerous flight and simulation programs at the Ames Research Center including rotorcraft, fixed-wing aircraft, advanced short takeoff and vertical landing (ASTOVL), vertical/short takeoff and landing (V/STOL), tiltrotor aircraft, and rotor experiments in the wind tunnel. Excellent system characterization and dynamic response prediction is achieved for this wide class of systems. Examples illustrate the role of system-identification technology in providing an integrated flow of dynamic response data around the entire life-cycle of aircraft development from initial specifications, through simulation and bench testing, and into flight-test optimization.

Tischler, Mark B.↗

Trajectory Optimization: OTIS 4

The latest release of the Optimal Trajectories by Implicit Simulation (OTIS4) allows users to simulate and optimize aerospace vehicle trajectories. With OTIS4, one can seamlessly generate optimal trajectories and parametric vehicle designs simultaneously. New features also allow OTIS4 to solve non-aerospace continuous time optimal control problems. The inputs and outputs of OTIS4 have been updated extensively from previous versions. Inputs now make use of objectoriented constructs, including one called a metastring. Metastrings use a greatly improved calculator and common nomenclature to reduce the user s workload. They allow for more flexibility in specifying vehicle physical models, boundary conditions, and path constraints. The OTIS4 calculator supports common mathematical functions, Boolean operations, and conditional statements. This allows users to define their own variables for use as outputs, constraints, or objective functions. The user-defined outputs can directly interface with other programs, such as spreadsheets, plotting packages, and visualization programs. Internally, OTIS4 has more explicit and implicit integration procedures, including high-order collocation methods, the pseudo-spectral method, and several variations of multiple shooting. Users may switch easily between the various methods. Several unique numerical techniques such as automated variable scaling and implicit integration grid refinement, support the integration methods. OTIS4 is also significantly more user friendly than previous versions. The installation process is nearly identical on various platforms, including Microsoft Windows, Apple OS X, and Linux operating systems. Cross-platform scripts also help make the execution of OTIS and post-processing of data easier. OTIS4 is supplied free by NASA and is subject to ITAR (International Traffic in Arms Regulations) restrictions. Users must have a Fortran compiler, and a Python interpreter is highly recommended.

Riehl, John P.↗

Hyper-X Stage Separation Trajectory Validation Studies

An independent twelve degree-of-freedom simulation of the X-43A separation trajectory was created with the Program to Optimize Simulated trajectories (POST II). This simulation modeled the multi-body dynamics of the X-43A and its booster and included the effect of two pyrotechnically actuated pistons used to push the vehicles apart as well as aerodynamic interaction forces and moments between the two vehicles. The simulation was developed to validate trajectory studies conducted with a 14 degree-of-freedom simulation created early in the program using the Automatic Dynamic Analysis of Mechanics Systems (ADAMS) simulation software. The POST simulation was less detailed than the official ADAMS-based simulation used by the Project, but was simpler, more concise and ran faster, while providing similar results. The increase in speed provided by the POST simulation provided the Project with an alternate analysis tool. This tool was ideal for performing separation control logic trade studies that required the running of numerous Monte Carlo trajectories.

Tartabini, Paul V.↗

Predictive and Neural Predictive Control of Uncertain Systems

Accomplishments and future work are:(1) Stability analysis: the work completed includes characterization of stability of receding horizon-based MPC in the setting of LQ paradigm. The current work-in-progress includes analyzing local as well as global stability of the closed-loop system under various nonlinearities; for example, actuator nonlinearities; sensor nonlinearities, and other plant nonlinearities. Actuator nonlinearities include three major types of nonlineaxities: saturation, dead-zone, and (0, 00) sector. (2) Robustness analysis: It is shown that receding horizon parameters such as input and output horizon lengths have direct effect on the robustness of the system. (3) Code development: A matlab code has been developed which can simulate various MPC formulations. The current effort is to generalize the code to include ability to handle all plant types and all MPC types. (4) Improved predictor: It is shown that MPC design using better predictors that can minimize prediction errors. It is shown analytically and numerically that Smith predictor can provide closed-loop stability under GPC operation for plants with dead times where standard optimal predictor fails. (5) Neural network predictors: When neural network is used as predictor it can be shown that neural network predicts the plant output within some finite error bound under certain conditions. Our preliminary study shows that with proper choice of update laws and network architectures such bound can be obtained. However, much work needs to be done to obtain a similar result in general case.

Kelkar, Atul G.↗

Hypersonic Vehicle Trajectory Optimization and Control

Two classes of neural networks have been developed for the study of hypersonic vehicle trajectory optimization and control. The first one is called an 'adaptive critic'. The uniqueness and main features of this approach are that: (1) they need no external training; (2) they allow variability of initial conditions; and (3) they can serve as feedback control. This is used to solve a 'free final time' two-point boundary value problem that maximizes the mass at the rocket burn-out while satisfying the pre-specified burn-out conditions in velocity, flightpath angle, and altitude. The second neural network is a recurrent network. An interesting feature of this network formulation is that when its inputs are the coefficients of the dynamics and control matrices, the network outputs are the Kalman sequences (with a quadratic cost function); the same network is also used for identifying the coefficients of the dynamics and control matrices. Consequently, we can use it to control a system whose parameters are uncertain. Numerical results are presented which illustrate the potential of these methods.

Balakrishnan, S. N.↗

Construction of a Piezoresistive Neural Sensor Array

The construction of a piezoresistive - piezoelectric sensor (or actuator) array is proposed using 'neural' connectivity for signal recognition and possible actuation functions. A closer integration of the sensor and decision functions is necessary in order to achieve intrinsic identification within the sensor. A neural sensor is the next logical step in development of truly 'intelligent' arrays. This proposal will integrate 1-3 polymer piezoresistors and MLC electroceramic devices for applications involving acoustic identification. The 'intelligent' piezoresistor -piezoelectric system incorporates printed resistors, composite resistors, and a feedback for the resetting of resistances. A model of a design is proposed in order to simulate electromechanical resistor interactions. The goal of optimizing a sensor geometry for improving device reliability, training, & signal identification capabilities is the goal of this work. At present, studies predict performance of a 'smart' device with a significant control of 'effective' compliance over a narrow pressure range due to a piezoresistor percolation threshold. An interesting possibility may be to use an array of control elements to shift the threshold function in order to change the level of resistance in a neural sensor array for identification, or, actuation applications. The proposed design employs elements of: (1) conductor loaded polymers for a 'fast' RC time constant response; and (2) multilayer ceramics for actuation or sensing and shifting of resistance in the polymer. Other material possibilities also exist using magnetoresistive layered systems for shifting the resistance. It is proposed to use a neural net configuration to test and to help study the possible changes required in the materials design of these devices. Numerical design models utilize electromechanical elements, in conjunction with structural elements in order to simulate piezoresistively controlled actuators and changes in resistance of sensors. The construction of these devices may show significant improvement in ability to interrogate signals and in the control of effective compliance. This work focuses on the development a variety of series/parallel interconnected piezoresistive control elements for the neural sensing function.

Carlson, W. B.↗

Shuttle Entry Imaging Using Infrared Thermography

During the Columbia Accident Investigation, imaging teams supporting debris shedding analysis were hampered by poor entry image quality and the general lack of information on optical signatures associated with a nominal Shuttle entry. After the accident, recommendations were made to NASA management to develop and maintain a state-of-the-art imagery database for Shuttle engineering performance assessments and to improve entry imaging capability to support anomaly and contingency analysis during a mission. As a result, the Space Shuttle Program sponsored an observation campaign to qualitatively characterize a nominal Shuttle entry over the widest possible Mach number range. The initial objectives focused on an assessment of capability to identify/resolve debris liberated from the Shuttle during entry, characterization of potential anomalous events associated with RCS jet firings and unusual phenomenon associated with the plasma trail. The aeroheating technical community viewed the Space Shuttle Program sponsored activity as an opportunity to influence the observation objectives and incrementally demonstrate key elements of a quantitative spatially resolved temperature measurement capability over a series of flights. One long-term desire of the Shuttle engineering community is to calibrate boundary layer transition prediction methodologies that are presently part of the Shuttle damage assessment process using flight data provided by a controlled Shuttle flight experiment. Quantitative global imaging may offer a complementary method of data collection to more traditional methods such as surface thermocouples. This paper reviews the process used by the engineering community to influence data collection methods and analysis of global infrared images of the Shuttle obtained during hypersonic entry. Emphasis is placed upon airborne imaging assets sponsored by the Shuttle program during Return to Flight. Visual and IR entry imagery were obtained with available airborne imaging platforms used within DoD along with agency assets developed and optimized for use during Shuttle ascent to demonstrate capability (i.e., tracking, acquisition of multispectral data, spatial resolution) and identify system limitations (i.e., radiance modeling, saturation) using state-of-the-art imaging instrumentation and communication systems. Global infrared intensity data have been transformed to temperature by comparison to Shuttle flight thermocouple data. Reasonable agreement is found between the flight thermography images and numerical prediction. A discussion of lessons learned and potential application to a potential Shuttle boundary layer transition flight test is presented.

Horvath, Thomas↗

The Future of in-Situ Sequencing-Based Microbial Monitoring: Development of a Shelf-Stable Method for Artemis and Beyond

Microbial monitoring onboard the International Space Station (ISS) is essential for assessing the efficiency of the Environmental Control and Life Support Systems (ECLSS) and providing insight into potential risk to both crew and spacecraft. Historically, this monitoring required the need to culture organisms onboard, return these cultures to Earth, and then complete the identifications, a process that would take months. Over the past decade, and through numerous payloads, advances in molecular biology have enabled in-flight microbial identifications using nanopore sequencing. The swab-to-sequencer method resulting from these efforts was transitioned from research to operations for microbial monitoring under the Crew Health Care Systems (CHeCS) BioMole. Collectively, these accomplishments have propelled the swab-to-sequencer method to be selected as the Microbial Surface Monitor (MSM) for Gateway, as well as a payload on Artemis IV. However, the lack of cold stowage availability for Artemis requires modifications to the entire method due to the thermal instability of the reagents required for sample preparation. To achieve this, new development, optimization, and validations were undertaken. Key considerations included enzyme concentration, buffer compatibility, and equal or enhanced sensitivity and specificity. At each step, thorough side-by-side comparisons with the current ISS method were performed. The development of a robust shelf-stable method will ensure continued sequencing-based microbial monitoring for Artemis and beyond, providing data in near real-time, enhancing risk response time, and yielding clear insight into the microbiome of spacecraft.

Christian G Mena↗