Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximate dynamic programming problem”

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 91 records · Page 5

Techniques for developing approximate optimal advanced launch system guidance

An extension to the authors' previous technique used to develop a real-time guidance scheme for the Advanced Launch System is presented. The approach is to construct an optimal guidance law based upon an asymptotic expansion associated with small physical parameters, epsilon. The trajectory of a rocket modeled as a point mass is considered with the flight restricted to an equatorial plane while reaching an orbital altitude at orbital injection speeds. The dynamics of this problem can be separated into primary effects due to thrust and gravitational forces, and perturbation effects which include the aerodynamic forces and the remaining inertial forces. An analytic solution to the reduced-order problem represented by the primary dynamics is possible. The Hamilton-Jacobi-Bellman or dynamic programming equation is expanded in an asymptotic series where the zeroth-order term (epsilon = 0) can be obtained in closed form.

Feeley, Timothy S.↗

A computer program for the geometrically nonlinear static and dynamic analysis of arbitrarily loaded shells of revolution, theory and users manual

A digital computer program known as SATANS (static and transient analysis, nonlinear, shells) for the geometrically nonlinear static and dynamic response of arbitrarily loaded shells of revolution is presented. Instructions for the preparation of the input data cards and other information necessary for the operation of the program are described in detail and two sample problems are included. The governing partial differential equations are based upon Sanders' nonlinear thin shell theory for the conditions of small strains and moderately small rotations. The governing equations are reduced to uncoupled sets of four linear, second order, partial differential equations in the meridional and time coordinates by expanding the dependent variables in a Fourier sine or cosine series in the circumferential coordinate and treating the nonlinear modal coupling terms as pseudo loads. The derivatives with respect to the meridional coordinate are approximated by central finite differences, and the displacement accelerations are approximated by the implicit Houbolt backward difference scheme with a constant time interval. The boundaries of the shell may be closed, free, fixed, or elastically restrained. The program is coded in the FORTRAN 4 language and is dimensioned to allow a maximum of 10 arbitrary Fourier harmonics and a maximum product of the total number of meridional stations and the total number of Fourier harmonics of 200. The program requires 155,000 bytes of core storage.

Ball, R. E.↗

Rotordynamic Analysis of the SSME Turbopumps Using Reduced Models

Alternative methods for the rotor-dynamic and sensitivity analysis of large rotor systems are examined. The methods are assessed for their ability to utilize accurate models of reduced size along with effective procedures for describing the dynamic behavior of the systems. Frequency response-based techniques are developed for determining the steady state response to imbalance of the space shuttle main engine (SSME) turbopumps and the related eigenvalue problem. The rotor and housing are represented by reduced receptances associated with their coupling points. The housing may be described by all of its normal modes within a frequency range of interest. The effects of truncated higher and lower modes are accounted for in an approximate manner. A procedure is described for determining the sensitivity of the coupling forces to changes in the coupling elements and rotor speed of the turbopump systems. In addition, an eigenvalue sensitivity analysis technique is adopted for application to the systems. Computer programs were developed for the numerical implementation of the impedance and eigenvalue sensitivity formulated in this study.

Noah, S. T.↗

Down-to-earth supercomputing at NASA Ames

The requirements for computers for aerospace vehicle design problems, in particular fluid dynamics, are examined. The levels of approximation to Navier-Stokes equations are: (1) linearized inviscid form, (2) nonlinear inviscid equation, (3) Reynolds-averaged Navier-Stokes equation, (4) large eddy simulation with models for subgrid scale turbulence, and (5) full Navier-Stokes equations to calculate laminar/turbulent transition and turbulence dissipation. The Cray 2 with 256 megaword memory and 4.1 nsec clock and the Cray X-MP/48 with a 128 megaword memory and a 9.5 nsec clock of the Ames facility are described, and an example revealing their applicability is presented. The use of zonal methods to simulate aerodynamic flow fields numerically is discussed. Consideration is given to the development of a program for simulation of aerodynamic flow fields, the compilers, symbolic debuggers, the scheduling of multiple CPUs, and the analysis and storage of the data.

Stevens, Kenneth G., Jr.↗

Tractable minor-free generalization of planar zero-field Ising models

In this work, we present a new family of zero-field Ising models over N binary variables/spins obtained by consecutive 'gluing' of planar and O(1)-sized components and subsets of at most three vertices into a tree. The polynomial time algorithm of the dynamic programming type for solving exact inference (computing partition function) and exact sampling (generating i.i.d. samples) consists of sequential application of an efficient (for planar) or brute-force (for O(1)-sized) inference and sampling to the components as a black box. To illustrate the utility of the new family of tractable graphical models, we first build a polynomial algorithm for inference and sampling of zero-field Ising models over K 33 -minor-free topologies and over K 5 -minor-free topologies—both of which are extensions of the planar zero-field Ising models—which are neither genus- nor treewidth-bounded. Second, we empirically demonstrate an improvement in the approximation quality of the NP-hard problem of inference over the square-grid Ising model in a node-dependent nonzero 'magnetic' field.

97 MATHEMATICS AND COMPUTING↗

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↗

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.↗

SIM PlanetQuest: Science with the Space Interferometry Mission

SIM - the Space Interferometry Mission - will perform precision optical astrometry on objects as faint as R magnitude 20. It will be the first space-based astrometric interferometer, operating in the optical band with a 10-m baseline. The Project is managed by the Jet Propulsion Laboratory, California Institute of Technology, in close collaboration with two industry partners, Lockheed Martin Missiles and Space, and TRW Inc., Space and Electronics Group. Launch of SIM is currently planned for 2009. In its wide-angle astrometric mode, SIM will yield 4 microarcsecond absolute position and parallax measurements. Astrometric planet searches will be done in a narrow-angle mode, with an accuracy of 4 microarcseconds or better in a single measurement. As a pointed rather than a survey instrument, SIM will maintain.its astrometric accuracy down to the faintest, magnitudes, opening up the opportunity for astrometry of active galactic nuclei to better than 10 pas. SIM will define a new astrometric reference frame, using a grid of approximately 1500 stars with positions accurate to 4 microarcseconds. The SIM Science Team comprises the Principal Investigators of ten Key Projects, and five Mission Scientists contributing their expertise to specific areas of the mission. Their science programs cover a wide range of topics in Galactic and extragalactic astronomy. They include: searches for low-mass planets - including analogs to our own solar system - tlie formation and dynamics of our Galaxy, calibration of the cosmic distance scale, and fundamental stellar astrophysics. All of the science observing on SIM is competitively awarded; the Science Team programs total about 40% of the total available, and the remainder will be assigned via future NASA competitions. This report is a compilation of science summaries by members of the Science Team, and it illustrates the wealth of scientific problems that microarcsecond-precision astrometry can contribute to. More information on SIM, including copies of this report, may be obtained from the project web site, at http://sim. jpl.nasa.gov.

Unwin, Stephen↗

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.↗

En route Spacing Tool: Efficient Conflict-free Spacing to Flow-Restricted Airspace

This paper describes the Air Traffic Management (ATM) problem within the U.S. of flow-restricted en route airspace, an assessment of its impact on airspace users, and a set of near-term tools and procedures to resolve the problem. The FAA is committed, over the next few years, to deploy the first generation of modem ATM decision support tool (DST) technology under the Free-Flight Phase-1 (FFp1) program. The associated en route tools include the User Request Evaluation Tool (URET) and the Traffic Management Advisor (TMA). URET is an initial conflict probe (ICP) capability that assists controllers with the detection and resolution of conflicts in en route airspace. TMA orchestrates arrivals transitioning into high-density terminal airspace by providing controllers with scheduled times of arrival (STA) and delay feedback advisories to assist with STA conformance. However, these FFPl capabilities do not mitigate the en route Miles-In-Trail (MIT) restrictions that are dynamically applied to mitigate airspace congestion. National statistics indicate that en route facilities (Centers) apply Miles-In-Trail (MIT) restrictions for approximately 5000 hours per month. Based on results from this study, an estimated 45,000 flights are impacted by these restrictions each month. Current-day practices for implementing these restrictions result in additional controller workload and an economic impact of which the fuel penalty alone may approach several hundred dollars per flight. To mitigate much of the impact of these restrictions on users and controller workload, a DST and procedures are presented. The DST is based on a simple derivative of FFP1 technology that is designed to introduce a set of simple tools for flow-rate (spacing) conformance and integrate them with conflict-probe capabilities. The tool and associated algorithms are described based on a concept prototype implemented within the CTAS baseline in 1995. A traffic scenario is used to illustrate the controller's use of the tool, and potential display options are presented for future controller evaluation.

Green, S.↗

Overvoltage Prevention and Curtailment Reduction Using Adaptive Droop-Based Supplementary Control in Smart Inverters

Recent developments in the renewable energy sector have seen an unprecedented growth in residential photovoltaic (PV) installations. However, high PV penetration levels often lead to overvoltage problems in low-voltage (LV) distribution feeders. Smart inverter control such as active power curtailment (APC)-based overvoltage control can be implemented to overcome these challenges. The APC technique utilizes a constant droop-based approach which curtails power rigidly, which can lead to significant energy curtailment in the LV distribution feeders. In this paper, different variations of the APC technique with linear, quadratic, and exponential droops have been analyzed from the point-of-view of energy curtailment for a LV distribution network in North America. Further, a combinatorial approach using various droop-based APC methods in conjunction with adaptive dynamic programming (ADP) as a supplementary control scheme has also been proposed. The proposed approach minimizes energy curtailment in the LV distribution network by adjusting the droop gains. Simulation results depict that ADP in conjunction with exponential droop reduces the energy curtailment to approximately 50% compared to using the standard linear droop.

42 ENGINEERING↗

Comparison of some optimal control methods for the design of turbine blades

This paper attempts a comparative study of some numerical methods for the optimal control design of turbine blades whose vibration characteristics are approximated by Timoshenko beam idealizations with shear and incorporating simple boundary conditions. The blade was synthesized using the following methods: (1) conjugate gradient minimization of the system Hamiltonian in function space incorporating penalty function transformations, (2) projection operator methods in a function space which includes the frequencies of vibration and the control function, (3) epsilon-technique penalty function transformation resulting in a highly nonlinear programming problem, (4) finite difference discretization of the state equations again resulting in a nonlinear program, (5) second variation methods with complex state differential equations to include damping effects resulting in systems of inhomogeneous matrix Riccatti equations some of which are stiff, (6) quasi-linear methods based on iterative linearization of the state and adjoint equation. The paper includes a discussion of some substantial computational difficulties encountered in the implementation of these techniques together with a resume of work presently in progress using a differential dynamic programming approach.

Desilva, B. M. E.↗

Progress on the Use of Combined Analog and Photon Counting Detection for Raman Lidar

The Atmospheric Radiation Measurement (ARM) program Raman Lidar (CARL) was upgraded in 2004 with a new data system that provides simultaneous measurements of both the photomultiplier analog output voltage and photon counts. The so-called merge value added procedure (VAP) was developed to combine the analog and count-rate signals into a single signal with improved dynamic range. Earlier versions of this VAP tended to cause unacceptably large biases in the water vapor mixing ratio during the daytime as a result of improper matching between the analog and count-rate signals in the presence of elevated solar background levels. We recently identified several problems and tested a modified version of the merge VAP by comparing profiles of water vapor mixing ratio derived from CARL with simultaneous sonde data over a six month period. We show that the modified merge VAP significantly reduces the daytime bias, and results in mean differences that are within approximately 1% for both nighttime and daytime measurements.

Newsom, Rob↗

Dynamic Ride-Matching for Large-Scale Transportation Systems

Efficient dynamic ride-matching (DRM) in large-scale transportation systems is a key driver in transport simulations to yield answers to challenging problems. Although the DRM problem is simple to solve, it quickly becomes a computationally challenging problem in large-scale transportation system simulations. Therefore, this study thoroughly examines the DRM problem dynamics and proposes an optimization-based solution framework to solve the problem efficiently. To benefit from parallel computing and reduce computational times, the problem’s network is divided into clusters utilizing a commonly used unsupervised machine learning algorithm along with a linear programming model. Then, these sub-problems are solved using another linear program to finalize the ride-matching. At the clustering level, the framework allows users adjusting cluster sizes to balance the trade-off between the computational time savings and the solution quality deviation. A case study in the Chicago Metropolitan Area, U.S., illustrates that the framework can reduce the average computational time by 58% at the cost of increasing the average pick up time by 26% compared with a system optimum, that is, non-clustered, approach. Another case study in a relatively small city, Bloomington, Illinois, U.S., shows that the framework provides quite similar results to the system-optimum approach in approximately 62% less computational time.

33 ADVANCED PROPULSION SYSTEMS↗

Planar dynamics of a uniform beam with rigid bodies affixed to the ends

The planar dynamics of a uniform elastic beam subject to a variety of geometric and natural boundary conditions and external excitations were analyzed. The beams are inextensible and capable of small transverse bending deformations only. Classical beam vibration eigenvalue problems for a cantilever with tip mass, a cantilever with tip body and an unconstrained beam with rigid bodies at each are examined. The characteristic equations, eigenfunctions and orthogonality relations for each are derived. The forced vibration of a cantilever with tip body subject to base acceleration is analyzed. The exact solution of the governing nonhomogeneous partial differential equation with time dependent boundary conditions is presented and compared with a Rayleigh-Ritz approximate solution. The arbitrary planar motion of an elastic beam with rigid bodies at the ends is addressed. Equations of motion are derived for two modal expansions of the beam deflection. The motion equations are cast in a first order form suitable for numerical integration. Selected FORTRAN programs are provided.

Storch, J.↗

Minimal norm constrained interpolation

In computational fluid dynamics and in CAD/CAM, a physical boundary is usually known only discreetly and most often must be approximated. An acceptable approximation preserves the salient features of the data such as convexity and concavity. In this dissertation, a smooth interpolant which is locally concave where the data are concave and is locally convex where the data are convex is described. The interpolant is found by posing and solving a minimization problem whose solution is a piecewise cubic polynomial. The problem is solved indirectly by using the Peano Kernal theorem to recast it into an equivalent minimization problem having the second derivative of the interpolant as the solution. This approach leads to the solution of a nonlinear system of equations. It is shown that Newton's method is an exceptionally attractive and efficient method for solving the nonlinear system of equations. Examples of shape-preserving interpolants, as well as convergence results obtained by using Newton's method are also shown. A FORTRAN program to compute these interpolants is listed. The problem of computing the interpolant of minimal norm from a convex cone in a normal dual space is also discussed. An extension of de Boor's work on minimal norm unconstrained interpolation is presented.

Irvine, L. D.↗

Which Way is Up? Lessons Learned from Space Shuttle Sensorimotor Research

The Space Shuttle Program provided the opportunity to examine sensorimotor adaptation to space flight in unprecedented numbers of astronauts, including many over multiple missions. Space motion sickness (SMS) severity was highly variable across crewmembers. SMS generally lasted 2-3 days in-flight with approximately 1/3 of crewmembers experiencing moderate to severe symptoms, and decreased incidence in repeat flyers. While SMS has proven difficult to predict from susceptibility to terrestrial analogs, symptoms were alleviated by medications, restriction of early activities, maintaining familiar orientation with respect to the visual environment and maintaining contact cues. Adaptive changes were also reflected by the oculomotor and perceptual disturbances experienced early inflight and by the perceptual and motor coordination problems experienced during re-entry and landing. According to crew self-reports, systematic head movements performed during reentry, as long as paced within one's threshold for motion tolerance, facilitated the early readaptation process. The Shuttle provided early postflight crew access to document the initial performance decrements and time course of recovery. These early postflight measurements were critical to inform the program of risks associated with extending the duration of Shuttle missions. Neurological postflight deficits were documented using a standardized subjective rating by flight surgeons. Computerized dynamic posturography was also implemented as a quantitative means of assessing sensorimotor function to support crew return-to-duty assessments. Towards the end of the Shuttle Program, more emphasis has been placed on mapping physiological changes to functional performance. Future commercial flights will benefit from pre-mission training including exposures to launch and entry G transitions and sensorimotor adaptability assessments. While SMS medication usage will continue to be refined, non-pharmacological countermeasures (e.g., sensory aids) will have both space and Earth-based applications. Early postflight field tests are recommended to provide the evidence base for best practices for future commercial flight programs. Learning Objective: Overview of the Space Shuttle Program regarding adaptive changes in sensorimotor function, including what was learned from research, what was implemented for medical operations, and what is recommended for commercial flights.

Wood, S. J.↗