Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Approximate dynamic programming”

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 235 records · Page 13

Advanced launch system trajectory optimization using suboptimal control

The maximum-final mass trajectory of a proposed configuration of the Advanced Launch System is presented. A model for the two-stage rocket is given; the optimal control problem is formulated as a parameter optimization problem; and the optimal trajectory is computed using a nonlinear programming code called VF02AD. Numerical results are presented for the controls (angle of attack and velocity roll angle) and the states. After the initial rotation, the angle of attack goes to a positive value to keep the trajectory as high as possible, returns to near zero to pass through the transonic regime and satisfy the dynamic pressure constraint, returns to a positive value to keep the trajectory high and to take advantage of minimum drag at positive angle of attack due to aerodynamic shading of the booster, and then rolls off to negative values to satisfy the constraints. Because the engines cannot be throttled, the maximum dynamic pressure occurs at a single point; there is no maximum dynamic pressure subarc. To test approximations for obtaining analytical solutions for guidance, two additional optimal trajectories are computed: one using untrimmed aerodynamics and one using no atmospheric effects except for the dynamic pressure constraint. It is concluded that untrimmed aerodynamics has a negligible effect on the optimal trajectory and that approximate optimal controls should be able to be obtained by treating atmospheric effects as perturbations.

Shaver, Douglas A.↗

NASA's Earth Science Use of Commercially Availiable Remote Sensing Datasets: Cover Image

The cover image incorporates high resolution stereo pairs acquired from the DigitalGlobe(R) QuickBird sensor. It shows a digital elevation model of Meteor Crater, Arizona at approximately 1.3 meter point-spacing. Image analysts used the Leica Photogrammetry Suite to produce the DEM. The outside portion was computed from two QuickBird panchromatic scenes acquired October 2006, while an Optech laser scan dataset was used for the crater s interior elevations. The crater s terrain model and image drape were created in a NASA Constellation Program project focused on simulating lunar surface environments for prototyping and testing lunar surface mission analysis and planning tools. This work exemplifies NASA s Scientific Data Purchase legacy and commercial high resolution imagery applications, as scientists use commercial high resolution data to examine lunar analog Earth landscapes for advanced planning and trade studies for future lunar surface activities. Other applications include landscape dynamics related to volcanism, hydrologic events, climate change, and ice movement.

Underwood, Lauren W.↗

The multiple-function multi-input/multi-output digital controller system for the AFW wind tunnel model

A real-time multiple-function digital controller system was developed for the Active Flexible Wing (AFW) Program. The digital controller system (DCS) allowed simultaneous execution of two control laws: flutter suppression and either roll trim or a rolling maneuver load control. The DCS operated within, but independently of, a slower host operating system environment, at regulated speeds up to 200 Hz. It also coordinated the acquisition, storage, and transfer of data for near real-time controller performance evaluation and both open- and closed-loop plant estimation. It synchronized the operation of four different processing units, allowing flexibility in the number, form, functionality, and order of control laws, and variability in selection of sensors and actuators employed. Most importantly, the DCS allowed for the successful demonstration of active flutter suppression to conditions approximately 26 percent (in dynamic pressure) above the open-loop boundary in cases when the model was fixed in roll and up to 23 percent when it was free to roll. Aggressive roll maneuvers with load control were achieved above the flutter boundary.

Hoadley, Sherwood T.↗

An MHD Study of the Interaction Between the Solar Wind and the Interstellar Medium

The overall objective of this research program is to obtain a better understanding of the interaction between the solar wind and the interstellar medium through the use of numerical solutions of the time-dependent magnetohydrodynamic (MHD) equations. The simulated results have been compared with observations where possible and with the results from previous analytic and numerical studies. The primary accomplishment of this project has been the development of codes for 2-D models in both spherical and cylindrical coordinates and the application of the codes to the solar wind/interstellar medium interaction. Computations have been carried out for both a relatively simple gas-dynamic interaction and a flow-aligned interstellar magnetic field. The results have been shown to compare favorably with models that use more approximations and to modify and extend the previous results as would be expected. The simulations have also been used along with a data analysis study to provide a quantitative estimate of the distance to the termination and bow shocks. Some of the specific topics that have been studied are: (1) gas dynamic models of the solar wind/interstellar medium interaction, (2) termination shock response to large-scale solar wind fluctuations, and (3) distances to the termination shock and heliopause. The main results from each of these studies are summarized. The results were published in three papers which are included as attachments.

Steinolfson, R. S.↗

Ground-to-Flight Handling Qualities Comparisons for a High Performance Airplane

A flight test program was conducted in conjunction with a ground-based piloted simulation study to enable a comparison of handling qualities ratings for a variety of maneuvers between flight and simulation of a modern high performance airplane. Specific objectives included an evaluation of pilot-induced oscillation (PIO) tendencies and a determination of maneuver types which result in either good or poor ground-to-flight pilot handling qualities ratings. A General Dynamics F-16XL aircraft was used for the flight evaluations, and the NASA Langley Differential Maneuvering Simulator was employed for the ground based evaluations. Two NASA research pilots evaluated both the airplane and simulator characteristics using tasks developed in the simulator. Simulator and flight tests were all conducted within approximately a one month time frame. Maneuvers included numerous fine tracking evaluations at various angles of attack, load factors and speed ranges, gross acquisitions involving longitudinal and lateral maneuvering, roll angle captures, and an ILS task with a sidestep to landing. Overall results showed generally good correlation between ground and flight for PIO tendencies and general handling qualities comments. Differences in pilot technique used in simulator evaluations and effects of airplane accelerations and motions are illustrated.

Brandon, Jay M.↗

Adaptation of primate vestibuloocular reflex to altered peripheral vestibular inputs. I. Frequency-specific recovery of horizontal VOR after inactivation of the lateral semicircular canals

1. The adaptive plasticity of the vestibuloocular reflex (VOR) following a selective lesion of the peripheral vestibular organs was investigated in rhesus monkeys whose lateral semicircular canals were inactivated by plugging of the canal lumen in both ears. Gain and phase of horizontal, vertical, and torsional slow-phase eye velocity were determined from three-dimensional eye movement recordings obtained acutely after the plugging operation, as well as in regular intervals up to 10 mo later. 2. Acutely after plugging, horizontal VOR was minimal during yaw rotation with gains of < 0.1 at all frequencies. Horizontal VOR gain gradually increased over time, reaching gains of 0.4-0.5 for yaw oscillations at 1.1 Hz approximately 5 mo after lateral canal inactivation. This response recovery was strongly frequency dependent: horizontal VOR gains were largest at the highest frequency tested and progressively decreased for lower frequencies. Below approximately 0.1 Hz, no consistent horizontal VOR could be elicited even 10 mo after plugging. 3. The frequency-dependent changes in gain paralleled changes in horizontal VOR phase. Below approximately 0.1-0.05 Hz large phase leads were present, similarly as in semicircular canal primary afferents. Smaller phase leads were also present at higher frequencies, particularly at 1.1 Hz (the highest frequency tested). 4. Consistent with the afferent-like dynamics of the adapted horizontal VOR, per- and postrotatory horizontal responses to constant-velocity yaw rotations were short lasting. Time constants of the slow-phase eye velocity envelope of the horizontal postrotatory nystagmus were approximately 2 s. Nonetheless, a consistent horizontal optokinetic afternystagmus was evoked in plugged animals. 5. A torsional component that was absent in intact animals was consistently present during yaw rotation acutely after lateral canal inactivation and remained approximately constant thereafter. The frequency response characteristics of this torsional component resembled those of the adapted horizontal slow-phase responses: gain decreased and large phase leads were introduced at frequencies below approximately 0.05-0.1 Hz. Torsional responses elicited by roll oscillations in supine position, on the other hand, were indistinguishable in their dynamics from intact animals. No consistent vertical nystagmus was elicited during yaw rotation. 6. Our results show that there is a slow, frequency-specific recovery of horizontal VOR after selective inactivation of the lateral semicircular canals. Both the spatial organization and the dynamic properties of the adapted VOR responses are distinctly different from responses in intact animals, suggesting complex changes in the underlying vestibuloocular circuitry.

NASA Program Space Physiology and Countermeasures↗

Assessing Huygens Probe Entry, Descent, and Landing at Titan Simulation using Dragonfly Atmosphere Model

Dragonfly is a New Frontiers Program mission that will deliver a rotorcraft to Saturn's moon, Titan. This mission follows Huygens as the previous mission that successfully landed a vehicle on Titan. A flight mechanics simulation of Dragonfly's Entry, Descent, and Landing sequence has been developed using the Program to Optimize Simulated Trajectories II. The simulation incorporates several subsystem models, including aerodynamics, gravity, and mass properties, to fully capture the multi-body six degree of freedom dynamics. Among all the subsystem models that inform the Entry, Descent, and Landing dynamics, the atmosphere model of Titan is a critical component. The atmosphere model characterizes the density, temperature, pressure, and winds that the entry vehicle experiences during the descent. This impacts several aspects of the descent such as the peak heating, aerodynamics, parachute release conditions, the dynamics of the vehicle and parachutes, and the landing ellipse. In the course of developing Dragonfly, an updated model of the Titan atmosphere has been created corresponding to Dragonfly's arrival in the mid-2030s, approximately one Titan year after the Huygens mission successfully landed a probe on Titan. This work leverages previous work done to investigate Huygens EDL sequence to assess the atmosphere model developed for Dragonfly. This is done by utilizing the updated Titan atmosphere model, the Dragonfly atmosphere model, within the Huygens POST2-based flight simulation with the goal of characterizing the differences between the atmospheric models and assessing how the Dragonfly atmosphere model impacts Huygens entry dynamics.

Huygens↗

Assessing Huygens Probe Entry, Descent, and Landing at Titan Simulation using Dragonfly Atmosphere Model

Dragonfly is a New Frontiers Program mission that will deliver a rotorcraft to Saturn's moon, Titan. This mission follows Huygens as the previous mission that successfully landed a vehicle on Titan. A flight mechanics simulation of Dragonfly's Entry, Descent, and Landing sequence has been developed using the Program to Optimize Simulated Trajectories II. The simulation incorporates several subsystem models, including aerodynamics, gravity, and mass properties, to fully capture the multi-body six degree of freedom dynamics. Among all the subsystem models that inform the Entry, Descent, and Landing dynamics, the atmosphere model of Titan is a critical component. The atmosphere model characterizes the density, temperature, pressure, and winds that the entry vehicle experiences during the descent. This impacts several aspects of the descent such as the peak heating, aerodynamics, parachute release conditions, the dynamics of the vehicle and parachutes, and the landing ellipse. In the course of developing Dragonfly, an updated model of the Titan atmosphere has been created corresponding to Dragonfly's arrival in the mid-2030s, approximately one Titan year after the Huygens mission successfully landed a probe on Titan. This work leverages previous work done to investigate Huygens EDL sequence to assess the atmosphere model developed for Dragonfly. This is done by utilizing the updated Titan atmosphere model, the Dragonfly atmosphere model, within the Huygens POST2-based flight simulation with the goal of characterizing the differences between the atmospheric models and assessing how the Dragonfly atmosphere model impacts Huygens entry dynamics.

Entry Descent Landing↗

The multiple-function multi-input/multi-output digital controller system for the AFW wind-tunnel model

A real time multiple-function digital controller system was developed for the Active Flexible Wing (AFW) Program. The digital controller system (DCS) allowed simultaneous execution of two control laws: flutter suppression and either roll trim or a rolling maneuver load control. The DCS operated within, but independently of, a slower host operating system environment, at regulated speeds up to 200 Hz. It also coordinated the acquisition, storage, and transfer of data for near real time controller performance evaluation and both open- and closed-loop plant estimation. It synchronized the operation of four different processing units, allowing flexibility in the number, form, functionality, and order of control laws, and variability in the selection of the sensors and actuators employed. Most importantly, the DCS allowed for the successful demonstration of active flutter suppression to conditions approximately 26 percent (in dynamic pressure) above the open-loop boundary in cases when the model was fixed in roll and up to 23 percent when it was free to roll. Aggressive roll maneuvers with load control were achieved above the flutter boundary. The purpose here is to present the development, validation, and wind tunnel testing of this multiple-function digital controller system.

Hoadley, Sherwood T.↗

An MHD study of the interaction between the solar wind and the interstellar medium

The overall objective of this research program is to obtain a better understanding of the interaction between the solar wind and the interstellar medium through the use of numerical solutions of the time-dependent magnetohydrodynamic (MHD) equations. The simulated results will be compared with observations where possible and with the results from previous analytic and numerical studies. The primary progress during the first two years has been to develop codes for 2-D models in both spherical and cylindrical coordinates and to apply them to the solar wind-interstellar medium interaction. Computations have been carried out for both a relatively simple gas-dynamic interaction and a flow-aligned interstellar magnetic field. The results have been shown to compare favorably with models that use more approximations and to modify and extend the previous results as would be expected. Work has also been initiated on the development of a 3-D MHD code in spherical coordinates.

Steinolfson, R. S.↗

Transonic aeroelastic analysis of launch vehicle configurations

A numerical study of the aeroelastic stability of typical launch vehicle configurations in transonic flight is performed. Recent computational fluid dynamics techniques are used to simulate the transonic aerodynamic flow fields, as opposed to relying on experimental data for the unsteady aerodynamic pressures. The flow solver is coupled to an appropriate structural representation of the vehicle. The aerodynamic formulation is based on the thin layer approximation to the Reynolds-Averaged Navier-Stokes equations, where the account for turbulent mixing is done by the two-layer Baldwin and Lomax algebraic eddy viscosity model. The structural-dynamic equations are developed considering free-free flexural vibration of an elongated beam with variable properties and are cast in modal form. Aeroelastic analyses are performed by integrating simultaneously in the two sets of equations. By tracing the growth or decay of a perturbed oscillation, the aeroelastic stability of a given constant configuration can be ascertained. The method is described in detail, and results that indicate its application are presented. Applications include some validation cases for the algorithm developed, as well as the study of configurations known to have presented flutter programs in the past.

Filgueirasdeazevedo, Joao Luiz↗

Full Multigrid Flow Solver

FMG3D (full multigrid 3 dimensions) is a pilot computer program that solves equations of fluid flow using a finite difference representation on a structured grid. Infrastructure exists for three dimensions but the current implementation treats only two dimensions. Written in Fortran 90, FMG3D takes advantage of the recursive subroutine feature, dynamic memory allocation, and structured-programming constructs of that language. FMG3D supports multi-block grids with three types of block-to-block interfaces: periodic, C-zero, and C-infinity. For all three types, grid points must match at interfaces. For periodic and C-infinity types, derivatives of grid metrics must be continuous at interfaces. The available equation sets are as follows: scalar elliptic equations, scalar convection equations, and the pressure-Poisson formulation of the Navier-Stokes equations for an incompressible fluid. All the equation sets are implemented with nonzero forcing functions to enable the use of user-specified solutions to assist in verification and validation. The equations are solved with a full multigrid scheme using a full approximation scheme to converge the solution on each succeeding grid level. Restriction to the next coarser mesh uses direct injection for variables and full weighting for residual quantities; prolongation of the coarse grid correction from the coarse mesh to the fine mesh uses bilinear interpolation; and prolongation of the coarse grid solution uses bicubic interpolation.

Mineck, Raymond E.↗

Kalman Filtering with Inequality Constraints for Turbofan Engine Health Estimation

Kalman filters are often used to estimate the state variables of a dynamic system. However, in the application of Kalman filters some known signal information is often either ignored or dealt with heuristically. For instance, state variable constraints (which may be based on physical considerations) are often neglected because they do not fit easily into the structure of the Kalman filter. This paper develops two analytic methods of incorporating state variable inequality constraints in the Kalman filter. The first method is a general technique of using hard constraints to enforce inequalities on the state variable estimates. The resultant filter is a combination of a standard Kalman filter and a quadratic programming problem. The second method uses soft constraints to estimate state variables that are known to vary slowly with time. (Soft constraints are constraints that are required to be approximately satisfied rather than exactly satisfied.) The incorporation of state variable constraints increases the computational effort of the filter but significantly improves its estimation accuracy. The improvement is proven theoretically and shown via simulation results. The use of the algorithm is demonstrated on a linearized simulation of a turbofan engine to estimate health parameters. The turbofan engine model contains 16 state variables, 12 measurements, and 8 component health parameters. It is shown that the new algorithms provide improved performance in this example over unconstrained Kalman filtering.

Simon, Dan↗

Characterization of the Scale Model Acoustic Test Overpressure Environment using Computational Fluid Dynamics

The Scale Model Acoustic Test (SMAT) is a 5% scale test of the Space Launch System (SLS), which is currently being designed at Marshall Space Flight Center (MSFC). The purpose of this test is to characterize and understand a variety of acoustic phenomena that occur during the early portions of lift off, one being the overpressure environment that develops shortly after booster ignition. The pressure waves that propagate from the mobile launcher (ML) exhaust hole are defined as the ignition overpressure (IOP), while the portion of the pressure waves that exit the duct or trench are the duct overpressure (DOP). Distinguishing the IOP and DOP in scale model test data has been difficult in past experiences and in early SMAT results, due to the effects of scaling the geometry. The speed of sound of the air and combustion gas constituents is not scaled, and therefore the SMAT pressure waves propagate at approximately the same speed as occurs in full scale. However, the SMAT geometry is twenty times smaller, allowing the pressure waves to move down the exhaust hole, through the trench and duct, and impact the vehicle model much faster than occurs at full scale. The DOP waves impact portions of the vehicle at the same time as the IOP waves, making it difficult to distinguish the different waves and fully understand the data. To better understand the SMAT data, a computational fluid dynamics (CFD) analysis was performed with a fictitious geometry that isolates the IOP and DOP. The upper and lower portions of the domain were segregated to accomplish the isolation in such a way that the flow physics were not significantly altered. The Loci/CHEM CFD software program was used to perform this analysis.

Nielsen, Tanner↗

Volume 35, AMT-1 Cruise Report and Preliminary Results

This report documents the scientific activities on board the Royal Research Ship (RRS) 'James Clark Ross' during the irst Atlantic Meridional Transect (AMT-1), 21 September to 24 October 1995. The ship sailed from Grimsby (England) for Montevideo (Uruguay) and then continued on to Stanley (Falkland Islands). The primary objective of the AMT program is to investigate basic biological processes in the open Atlantic Ocean over very broad spatial scales. For AMT-1, the meridional range covered was approximately 50 deg N to 50 deg S or nearly 8,000 nmi. The measurements to be taken during the AMT cruises are fundamental for the calibration, validation, and continuing understanding of remotely sensed observations of biological oceanography. They are also important for understanding plankton community structure over latitudinal scales and the role of the world ocean in global carbon cycles. During AMT-1 a variety of instruments were used to map the physical, chemical, and biological structure of the upper 200 m of the water column. Ocean color measurements were made using state-of-the-art sensors, whose calibration was traceable to the highest international standards. New advances in fluorometry were used to measure photosynthetic activity, which was then used to further interpret primary productivity. A unique set of samples and data were collected for the planktonic assemblages that vary throughout the range of the transect. These data will yield new interpretations on community composition and their role in carbon cycling. While the various provinces of the Atlantic Ocean were being crossed, the partial pressure of CO2 was related to biological productivity. This comparison revealed the areas of drawdown of atmospheric CO2 and how these areas relate to the surrounding biological productivity. These data, plus the measurements of light attenuation and phytoplankton optical properties, will be used as a primary input for basin-scale biological productivity models to help develop ecosystem dynamics models which will be important for improving the forecasting abilities of modelers. The AMT program is also attempting to meet the needs of international agencies in their implementation of Sensor Intercomparison and Merger for Biological and Interdisciplinary Ocean Studies (SIMBIOS), a program to develop a methodology and operational capability to combine data products from the various ocean color satellite missions.

Hooker, Stanford B.↗

Propellant Slosh Analysis for the Solar Dynamics Observatory

The Solar Dynamics Observatory (SDO) mission, part of the Living With a Star program, is a geosynchronous satellite with tight pointing requirements. Due to a large amount of liquid propellant, a detailed slosh analysis is required to ensure the tight pointing budget can be satisfied. Much of the high fidelity slosh analysis and simulation has been performed via computational fluid dynamics. Even though this method of simulation is very accurate, it requires significant computational effort and specialized knowledge, limiting the ability of the SDO project to access fluid dynamics simulations at will. Furthermore, it is very difficult to incorporate most of these models into simulations of the overall spacecraft and its environment. Ultimately, the effects of the propellant slosh on the attitude stability and pointing performance of the entire spacecraft are of great interest to attitude control engineers. Equivalent mechanical models, such as models that approximate the fluid slosh effects by analogy to the movements of a point-mass pendulum, are important tools in simulating propellant slosh dynamics as part of the entire attitude determination and control system. This paper describes some of the current methods used to analyze and model slosh. It focuses on equivalent mechanical models and their incorporation into control-based analysis tools such as Simulink. The SDO mission is used as the case study for this work.

Mason, Paul A. C.↗

Unsteady transonic flow calculations for realistic aircraft configurations

A transonic unsteady aerodynamic and aeroelasticity code has been developed for application to realistic aircraft configurations. The new code is called CAP-TSD which is an acronym for Computational Aeroelasticity Program - Transonic Small Disturbance. The CAP-TSD code uses a time-accurate approximate factorization (AF) algorithm for solution of the unsteady transonic small-disturbance equation. The AF algorithm is very efficient for solution of steady and unsteady transonic flow problems. It can provide accurate solutions in only several hundred time steps yielding a significant computational cost savings when compared to alternative methods. The new code can treat complete aircraft geometries with multiple lifting surfaces and bodies including canard, wing, tail, control surfaces, launchers, pylons, fuselage, stores, and nacelles. Applications are presented for a series of five configurations of increasing complexity to demonstrate the wide range of geometrical applicability of CAP-TSD. These results are in good agreement with available experimental steady and unsteady pressure data. Calculations for the General Dynamics one-ninth scale F-16C aircraft model are presented to demonstrate application to a realistic configuration. Unsteady results for the entire F-16C aircraft undergoing a rigid pitching motion illustrated the capability required to perform transonic unsteady aerodynamic and aeroelastic analyses for such configurations.

Batina, John T.↗

Unsteady transonic flow calculations for realistic aircraft configurations

A transonic unsteady aerodynamic and aeroelasticity code has been developed for application to realistic aircraft configurations. The new code is called CAP-TSD which is an acronym for Computational Aeroelasticity Program - Transonic Small Disturbance. The CAP-TSD code uses a time-accurate approximate factorization (AF) algorithm for solution of the unsteady transonic small-disturbance equation. The AF algorithm is very efficient for solution of steady and unsteady transonic flow problems. It can provide accurate solutions in only several hundred time steps yielding a significant computational cost savings when compared to alternative methods. The new code can treat complete aircraft geometries with multiple lifting surfaces and bodies including canard, wing, tail, control surfaces, launchers, pylons, fuselage, stores, and nacelles. Applications are presented for a series of five configurations of increasing complexity to demonstrate the wide range of geometrical applicability of CAP-TSD. These results are in good agreement with available experimental steady and unsteady pressure data. Calculations for the General Dynamics one-ninth scale F-16C aircraft model are presented to demonstrate application to a realistic configuration. Unsteady results for the entire F-16C aircraft undergoing a rigid pitching motion illustrated the capability required to perform transonic unsteady aerodynamic and aeroelastic analyses for such configurations.

Batina, John T.↗