Search NASA⌕ Search

SEARCH · Search NASA

Results for “Differential 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Measurement and Characterization of Helicopter Noise in Steady-State and Maneuvering Flight

A special acoustic flight test program was performed on the Bell 206B helicopter outfitted with an in-flight microphone boom/array attached to the helicopter while simultaneous acoustic measurements were made using a linear ground array of microphones arranged to be perpendicular to the flight path. Air and ground noise measurements were made in steady-state longitudinal and steady turning flight, and during selected dynamic maneuvers. Special instrumentation, including direct measurement of the helicopter s longitudinal tip-path-plane (TPP) angle, Differential Global Positioning System (DGPS) and Inertial Navigation Unit (INU) measurements, and a pursuit guidance display were used to measure important noise controlling parameters and to make the task of flying precise operating conditions and flight track easier for the pilot. Special care was also made to test only in very low winds. The resulting acoustic data is of relatively high quality and shows the value of carefully monitoring and controlling the helicopter s performance state. This paper has shown experimentally, that microphones close to the helicopter can be used to estimate the specific noise sources that radiate to the far field, if the microphones are positioned correctly relative to the noise source. Directivity patterns for steady, turning flight were also developed, for the first time, and connected to the turning performance of the helicopter. Some of the acoustic benefits of combining normally separated flight segments (i.e. an accelerated segment and a descending segment) were also demonstrated.

Schmitz, Fredric H.↗
NASA NTRS1978

Numerical method for solution of systems of non-stationary spatially one-dimensional nonlinear differential equations

A computational scheme and a standard program is proposed for solving systems of nonstationary spatially one-dimensional nonlinear differential equations using Newton's method. The proposed scheme is universal in its applicability and its reduces to a minimum the work of programming. The program is written in the FORTRAN language and can be used without change on electronic computers of type YeS and BESM-6. The standard program described permits the identification of nonstationary (or stationary) solutions to systems of spatially one-dimensional nonlinear (or linear) partial differential equations. The proposed method may be used to solve a series of geophysical problems which take chemical reactions, diffusion, and heat conductivity into account, to evaluate nonstationary thermal fields in two-dimensional structures when in one of the geometrical directions it can take a small number of discrete levels, and to solve problems in nonstationary gas dynamics.

Morozov, S. K.↗
NASA NTRS1990

Reduced-order filtering for flexible space structures

There is a need for feedback control of the large flexible space structures which are going to be increasingly important in the future of the space program. These structures are very lightly damped, and vibrations may persist for a long time when the system is disturbed unless an active feedback control strategy is used to damp out the vibrations. The system is best described by a partial differential equation description, but the more common approach is to use a large set of second order differential equations, where a large number of modes must be retained if the mathematics is to provide an adequate description of the dynamical process. Sensors, such as accelerometers and rate gyros, may provide data to the feedback controller so that it may respond appropriately to control the system. The data from the sensors is not perfect, but is subject to noise, called measurement noise, and the dynamical process itself is subject to disturbances referred to as process noise. Filtering the sensor signals to remove the measurement noise, and using the resulting state estimates to control the system are investigated.

Sims, Craig S.↗
NASA NTRSDate not supplied

A Verification Framework for Runtime Assurance of Autonomous UAS

Runtime Assurance (RTA) is a design-time architecture for safety-critical systems where an internal monitor acts upon detecting a violation of a property. The simplex architecture is an instance of RTA, where the action taken is to hand control of the overall system to a trusted controller when an untrusted one violates a safety property. Simplex RTA is emerging as a method for allowing AI/ML and other unverified software to be integrated into safety-critical applications like aircraft. To this end, the American Society for Testing and Materials (ASTM) and NASA have each published guidelines on the use of RTA in such systems. In the simplex RTA framework, a system has an advanced controller (AC) and a reversionary controller (RC). The system is allowed to operate with the AC until a runtime monitor detects that some property has been violated and then the RC takes over. Assuming that the sample rate of the monitor will detect improper functioning with enough time for the RC to correct the impending problem, and that the RC is trusted, the system will operate as intended. This use of the simplex RTA framework can allow for the integration of untrusted, but possibly more performant, controllers in a safe way. This paper presents a formalization of a simplex RTA framework in the Prototype Verification System (PVS) theorem prover using an embedding of differential dynamic logic (DDL) called Plaidypvs. A novel feature of this framework is that it can be instantiated at different levels of abstraction. This feature allows for the formal verification of a system with an untrusted black box component, such as an AI/ML controller. This paper does not address the many difficulties in deploying RTA in an industrial-level system. Instead, the focus is on the formal verification of the simplex RTA framework in the language of hybrid programs. Hybrid programs are programs that include both discrete and continuous dynamics and can be used to model complex cyber-physical systems. Plaidypvs is a tool that enables formalization of hybrid programs in the PVS theorem prover. Plaidypvs enables the verification of the general simplex RTA framework and then, by specializing some components of the hybrid program, verifying instances of the framework while treating the untrusted component as a black box. A selection of Unmanned Aircraft Systems (UAS) operations are shown as instances of the general RTA framework in PVS. This offers the benefit of design time verification of relevant safety properties to the system, and it also gives requirements on the sample rate of sensors that determine the time interval in which the ‘switch’ property of the RTA framework is checked.

PVS↗
NASA NTRS1989

Prototype secondary mirror assembly design for the space infrared telescope facility

The authors describe the design of a liquid helium temperature prototype secondary mirror assembly (PSMA) under developement for the NASA Space Infrared Telescope Facility (SIRTF) program. The SIRTF secondary mirror assembly must operate below 4 K and provide highly precise two-axis dynamic tilting (chopping) in addition to the conventional functions needed by the SIRTF observatory. The PSMA design uses a fused quartz mirror kinematically attached at its center to an aluminum cruciform. The mirror/cruciform assembly is driven in tilt about its combined center of mass using a unique flexure pivot and a four-actuator control system with feedback provided by pairs of differential position sensors. The voice coil actuators are mounted on a second flexture-pivoted mass, enhancing servo system stability and isolating the telescope from vibration-induced disturbances. The mirror/cruciform and the reaction mass are attached to opposite sides of an aluminum mounting plate whose dimensional characteristics are nominally identical to those of the aluminum flexure pivot material. The mounting plate is connected to the outer housing by a six-degree-of-freedom focus and centering mechanism using pivoted actuation levers driven by lead screw/harmonic drive/stepper motor assemblies.

Stier, M.↗
NASA NTRS1989

A vectorized algorithm for 3D dynamics of a tethered satellite

Equations of motion characterizing the three dimensional motion of a tethered satellite during the retrieval phase are studied. The mathematical model involves an arbitrary number of point masses connected by weightless cords. Motion occurs in a gravity gradient field. The formulation presented accounts for general functions describing support point motion, rate of tether retrieval, and arbitrary forces applied to the point masses. The matrix oriented program language MATLAB is used to produce an efficient vectorized formulation for computing natural frequencies and mode shapes for small oscillations about the static equilibrium configuration; and for integrating the nonlinear differential equations governing large amplitude motions. An example of time response pertaining to the skip rope effect is investigated.

Wilson, Howard B.↗
NASA NTRS1990

A computational algorithm for spacecraft control and momentum management

Developments in the area of nonlinear control theory have shown how coordinate changes in the state and input spaces of a dynamical system can be used to transform certain nonlinear differential equations into equivalent linear equations. These techniques are applied to the control of a spacecraft equipped with momentum exchange devices. An optimal control problem is formulated that incorporates a nonlinear spacecraft model. An algorithm is developed for solving the optimization problem using feedback linearization to transform to an equivalent problem involving a linear dynamical constraint and a functional approximation technique to solve for the linear dynamics in terms of the control. The original problem is transformed into an unconstrained nonlinear quadratic program that yields an approximate solution to the original problem. Two examples are presented to illustrate the results.

Dzielski, John↗
NASA NTRS2004

Further Improvement in 3DGRAPE

3DGRAPE/AL:V2 denotes version 2 of the Three-Dimensional Grids About Anything by Poisson's Equation with Upgrades from Ames and Langley computer program. The preceding version, 3DGRAPE/AL, was described in Improved 3DGRAPE (ARC-14069) NASA Tech Briefs, Vol. 21, No. 5 (May 1997), page 66. These programs are so named because they generate volume grids by iteratively solving Poisson's Equation in three dimensions. The grids generated by the various versions of 3DGRAPE have been used in computational fluid dynamics (CFD). The main novel feature of 3DGRAPE/AL:V2 is the incorporation of an optional scheme in which anisotropic Lagrange-based trans-finite interpolation (ALBTFI) is coupled with exponential decay functions to compute and blend interior source terms. In the input to 3DGRAPE/AL:V2 the user can specify whether or not to invoke ALBTFI in combination with exponential-decay controls, angles, and cell size for controlling the character of grid lines. Of the known programs that solve elliptic partial differential equations for generating grids, 3DGRAPE/AL:V2 is the only code that offers a combination of speed and versatility with most options for controlling the densities and other characteristics of grids for CFD.

Alter, Stephen↗
NASA NTRS1969

Dynamic characteristics of a variable-mass flexible missile: Dynamics of a two-stage variable-mass flexible rocket

The dynamic characteristics of two-stage slender elastic body were investigated. The first stage, containing a solid-fuel rocket, possesses variable mass while the second stage, envisioned as a flexible case, contains packaged instruments of constant mass. The mathematical formulation was in terms of vector equations of motion transformed by a variational principle into sets of scalar differential equations in terms of generalized coordinates. Solutions to the complete equations were obtained numerically by means of finite difference techniques. The problem has been programmed in the FORTRAN 4 language and solved on an IBM 360/50 computer. Results for limited cases are presented showing the nature of the solutions.

Meirovitch, L.↗
NASA NTRSDate not supplied

Progress of 2.05 μm Fiber Laser Development for CO2 DIAL Measurements of Martian CO2 and Pressure

The Decadal Survey by the US National Academy of Sciences in the United States and NASA’s Science Mission Directorate (SMD) Science Plan both require global CO2 observations for the Martian atmospheric pressure, dynamics and CO2 variations. However, there are significant observational gaps in polar regions and during nighttime for global air pressure and CO2 on Mars. Therefore, recently we proposed a new concept of Martian differential absorption lidar (DIAL) operating in the 2.05 μm CO2 absorption band for global, including poles, atmospheric CO2 and pressure observations day and night. Based on the concept, we are awarded to develop 2.05 μm fiber lasers by the NASA Planetary Instrument Concepts for the Advancement of Solar System Observations (PICASSO) Program. The laser is designed to be an all-fiber master oscillator and power amplifier (MOPA) system with laser out put of ~3 mJ at a repetition frequency of 2 kHz. The primary master oscillator (PMO) is locked to the center (2.0504280 μm) of the selected CO2 absorption line. The frequency of a second MO (SMO) is locked to that of PMO. This SMO frequency is adjustable and switched between the online (2.05044156 μm) and offline (2.05050812 μm) wavelengths. The online wavelength is optimally selected so that the CO2 absorption optical depth (AOD) is ~1.1 at 3 km and the measurement in the low Martian atmosphere has largest signal-to-noise ratio. The online wavelength can also be adjusted to a line slope location where AOD is larger to observe atmospheric pressure at higher altitudes. We will present more details about this project and instrument development.

CO2↗
NASA NTRS1995

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.↗
NASA NTRSDate not supplied

Development of 2.05 µm Fiber Lasers for CO2 DIAL Lidar Measuring Martian CO2 and Pressure

The Decadal Survey by the National Academy of Sciences in the United States and NASA’s Science Mission Directorate (SMD) Science Plan both require global CO 2 observations for the Martian atmospheric pressure, dynamics and chemistry. However, there are significant observational gaps in polar regions and during nighttime for global air pressure and CO 2 on Mars. Therefore, recently we proposed a new concept of Martian differential absorption lidar (DIAL) operating in the 2.05 µm CO 2 absorption band for global, including poles, atmospheric CO 2 and pressure observations day and night [1]. Based on the concept, we are awarded to develop 2.05 µm fiber lasers by the NASA Planetary Instrument Concepts for the Advancement of Solar System Observations (PICASSO) Program. The laser is designed to be an all-fiber master oscillator and power amplifier (MOPA) system with laser output of ~3 mJ at a repetition frequency of 2 kHz. The primary master oscillator (PMO) is locked to the center (2.0504280 µm) of the selected CO 2 absorption line. The frequency of a second MO (SMO) is locked to that of PMO. This SMO frequency is adjustable and switched between the online (2.05044156 µm) and offline (2.05050812 µm) wavelengths. The online wavelength is optimally selected so that the CO 2 absorption optical depth (AOD) is ~1.1 at 3 km and the measurement in the low Martian atmosphere has largest signal-to-noise ratio. The online wavelength can also be adjusted to a line slope location where AOD is larger to observe atmospheric pressure at higher altitudes. We will present more detail about this project and instrument development at the conference.

Zhaoyan Liu↗
NASA NTRS1984

Practical aspects of modeling aircraft dynamics from flight data

The purpose of parameter estimation, a subset of system identification, is to estimate the coefficients (such as stability and control derivatives) of the aircraft differential equations of motion from sampled measured dynamic responses. In the past, the primary reason for estimating stability and control derivatives from flight tests was to make comparisons with wind tunnel estimates. As aircraft became more complex, and as flight envelopes were expanded to include flight regimes that were not well understood, new requirements for the derivative estimates evolved. For many years, the flight determined derivatives were used in simulations to aid in flight planning and in pilot training. The simulations were particularly important in research flight test programs in which an envelope expansion into new flight regimes was required. Parameter estimation techniques for estimating stability and control derivatives from flight data became more sophisticated to support the flight test programs. As knowledge of these new flight regimes increased, more complex aircraft were flown. Much of this increased complexity was in sophisticated flight control systems. The design and refinement of the control system required higher fidelity simulations than were previously required.

Iliff, K. W.↗
NASA NTRS1994

Parallel Calculation of Sensitivity Derivatives for Aircraft Design using Automatic Differentiation

Sensitivity derivative (SD) calculation via automatic differentiation (AD) typical of that required for the aerodynamic design of a transport-type aircraft is considered. Two ways of computing SD via code generated by the ADIFOR automatic differentiation tool are compared for efficiency and applicability to problems involving large numbers of design variables. A vector implementation on a Cray Y-MP computer is compared with a coarse-grained parallel implementation on an IBM SP1 computer, employing a Fortran M wrapper. The SD are computed for a swept transport wing in turbulent, transonic flow; the number of geometric design variables varies from 1 to 60 with coupling between a wing grid generation program and a state-of-the-art, 3-D computational fluid dynamics program, both augmented for derivative computation via AD. For a small number of design variables, the Cray Y-MP implementation is much faster. As the number of design variables grows, however, the IBM SP1 becomes an attractive alternative in terms of compute speed, job turnaround time, and total memory available for solutions with large numbers of design variables. The coarse-grained parallel implementation also can be moved easily to a network of workstations.

Bischof, c. H.↗
NASA NTRS1993

Approximate optimal guidance for the advanced launch system

A real-time guidance scheme for the problem of maximizing the payload into orbit subject to the equations of motion for a rocket over a spherical, non-rotating earth is presented. An approximate optimal launch guidance law is developed based upon an asymptotic expansion of the Hamilton - Jacobi - Bellman or dynamic programming equation. The expansion is performed in terms of a small parameter, which is used to separate the dynamics of the problem into primary and perturbation dynamics. For the zeroth-order problem the small parameter is set to zero and a closed-form solution to the zeroth-order expansion term of Hamilton - Jacobi - Bellman equation is obtained. Higher-order terms of the expansion include the effects of the neglected perturbation dynamics. These higher-order terms are determined from the solution of first-order linear partial differential equations requiring only the evaluation of quadratures. This technique is preferred as a real-time, on-line guidance scheme to alternative numerical iterative optimization schemes because of the unreliable convergence properties of these iterative guidance schemes and because the quadratures needed for the approximate optimal guidance law can be performed rapidly and by parallel processing. Even if the approximate solution is not nearly optimal, when using this technique the zeroth-order solution always provides a path which satisfies the terminal constraints. Results for two-degree-of-freedom simulations are presented for the simplified problem of flight in the equatorial plane and compared to the guidance scheme generated by the shooting method which is an iterative second-order technique.

Feeley, T. S.↗
NASA NTRS2002

Microgravity: Molecular Dynamics Simulations at the NCCS Probe the Behavior of Liquids in Low Gravity

The life of the very small, whether in something as complicated as a human cell or as simple as a drop of water, is of fundamental scientific interest: By knowing how a tiny amount of material reacts to changes in its environment, scientists maybe able to answer questions about how a bulk of material would react to comparable changes. NASA is in the forefront of computational research into a broad range of basic scientific questions about fluid dynamics and the nature of liquid boundary instability. For example, one important issue for the space program is how drops of water and other materials will behave in the low-gravity environment of space and how the low gravity will affect the transport and containment of these materials. Accurate prediction of this behavior is among the aims of a set of molecular dynamics experiments carried out on the NCCSs Cray supercomputers. In conventional computational studies of materials, matter is treated as continuous - a macroscopic whole without regard to its molecular parts - and the behavior patterns of the matter in various physical environments are studied using well-established differential equations and mathematical parameters based on physical properties such as compressibility density, heat capacity, and vapor pressure of the bulk material.

Source record↗
NASA NTRS1993

Inversion of gravity and bathymetry in oceanic regions for long-wavelength variations in upper mantle temperature and composition

Long-wavelength variations in geoid height, bathymetry, and SS-S travel times are all relatable to lateral variations in the characteristic temperature and bulk composition of the upper mantle. The temperature and composition are in turn relatable to mantle convection and the degree of melt extraction from the upper mantle residuum. Thus the combined inversion of the geoid or gravity field, residual bathymetry, and seismic velocity information offers the promise of resolving fundamental aspects of the pattern of mantle dynamics. The use of differential body wave travel times as a measure of seismic velocity information, in particular, permits resolution of lateral variations at scales not resolvable by conventional global or regional-scale seismic tomography with long-period surface waves. These intermediate scale lengths, well resolved in global gravity field models, are crucial for understanding the details of any chemical or physical layering in the mantle and of the characteristics of so-called 'small-scale' convection beneath oceanic lithosphere. In 1991 a three-year project to the NASA Geophysics Program was proposed to carry out a systematic inversion of long-wavelength geoid anomalies, residual bathymetric anomalies, and differential SS-S travel time delays for the lateral variation in characteristic temperature and bulk composition of the oceanic upper mantle. The project was funded as a three-year award, beginning on 1 Jan. 1992.

Solomon, Sean C.↗
NASA NTRS1991

On the use of maximum likelihood estimation for the assembly of Space Station Freedom

Distributed parameter models of the Solar Array Flight Experiment, the Mini-MAST truss, and Space Station Freedom assembly are discussed. The distributed parameter approach takes advantage of (1) the relatively small number of model parameters associated with partial differential equation models of structural dynamics, (2) maximum-likelihood estimation using both prelaunch and on-orbit test data, (3) the inclusion of control system dynamics in the same equations, and (4) the incremental growth of the structural configurations. Maximum-likelihood parameter estimates for distributed parameter models were based on static compliance test results and frequency response measurements. Because the Space Station Freedom does not yet exist, the NASA Mini-MAST truss was used to test the procedure of modeling and parameter estimation. The resulting distributed parameter model of the Mini-MAST truss successfully demonstrated the approach taken. The computer program PDEMOD enables any configuration that can be represented by a network of flexible beam elements and rigid bodies to be remodeled.

Taylor, Lawrence W., Jr.↗