Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

General purpose algorithms for characterization of slow and fast phase nystagmus

In the overall aim for a better understanding of the vestibular and optokinetic systems and their roles in space motion sickness, the eye movement responses to various dynamic stimuli are measured. The vestibulo-ocular reflex (VOR) and the optokinetic response, as the eye movement responses are known, consist of slow phase and fast phase nystagmus. The specific objective is to develop software programs necessary to characterize the vestibulo-ocular and optokinetic responses by distinguishing between the two phases of nystagmus. The overall program is to handle large volumes of highly variable data with minimum operator interaction. The programs include digital filters, differentiation, identification of fast phases, and reconstruction of the slow phase with a least squares fit such that sinusoidal or psuedorandom data may be processed with accurate results. The resultant waveform, slow phase velocity eye movements, serves as input data to the spectral analysis programs previously developed for NASA to analyze nystagmus responses to pseudorandom angular velocity inputs.

Lessard, Charles S.↗
NASA NTRS1970

Dynamic characteristics of a variable-mass flexible missile

The general motion of a variable mass flexible missile with internal flow and aerodynamic forces is considered. The resulting formulation comprises six ordinary differential equations for rigid body motion and three partial differential equations for elastic motion. The simultaneous differential equations are nonlinear and possess time-dependent coefficients. The differential equations are solved by a semi-analytical method leading to a set of purely ordinary differential equations which are then solved numerically. A computer program was developed for the numerical solution and results are presented for a given set of initial conditions.

Meirovitch, L.↗
NASA NTRS1974

Finite element solution algorithm for incompressible fluid dynamics

A finite element solution algorithm is established for the two-dimensional Navier-Stokes equations governing the transient motion of a viscous incompressible fluid, i.e., hydrodynamics. Dependent variable transformation renders the differential equation description uniformly elliptic. The finite element algorithm is established using the Galerkin criterion on a local basis within the Method of Weighted Residuals. It is unconstrained with respect to system linearity, computational mesh uniformity or solution domain closure regularity. The finite element matrices are established using a linear 'natural coordinate function' description. Computational solutions using the COMOC computer program illustrate the various features of the algorithm including recirculating flows.

Baker, A. J.↗
DOE OSTI2015

Toward the credibility of Northeast United States summer precipitation projections in CMIP5 and NARCCAP simulations

Abstract Precipitation projections for the northeast United States and nearby Canada (Northeast) are examined for 15 Fifth Phase of the Coupled Model Intercomparison Project (CMIP5) models. A process‐based evaluation of atmospheric circulation features associated with wet Northeast summers is performed to examine whether credibility can be differentiated within the multimodel ensemble. Based on these evaluations, and an analysis of the interannual statistical properties of area‐averaged precipitation, model subsets were formed. Multimodel precipitation projections from each subset were compared to the multimodel projection from all of the models. Higher‐resolution North American Regional Climate Change Assessment Program (NARCCAP) regional climate models (RCMs) were subjected to a similar evaluation, grouping into subsets, and examination of future projections. CMIP5 models adequately simulate most large‐scale circulation features associated with wet Northeast summers, though all have errors in simulating observed sea level pressure and moisture divergence anomalies in the western tropical Atlantic/Gulf of Mexico. Relevant large‐scale processes simulated by the RCMs resemble those of their driving global climate models (GCMs), which are not always realistic. Future RCM studies could benefit from a process analysis of potential driving GCMs prior to dynamical downscaling. No CMIP5 or NARCCAP models were identified as clearly more credible, but six GCMs and four RCMs performed consistently better. Among the “Better” models, there is no consistency in the direction of future summer precipitation change. CMIP5 projections suggest that the Northeast precipitation response depends on the dynamics of the North Atlantic anticyclone and associated circulation and moisture convergence patterns, which vary among “Better” models. Even when model credibility cannot be clearly differentiated, examination of simulated processes provides important insights into their evolution under greenhouse warming.

Thibeault, Jeanne M.↗
NASA NTRS1993

CSI flight experiment projects of the Naval Research Laboratory

The Naval Research Laboratory (NRL) is involved in an active program of CSI flight experiments. The first CSI flight experiment of the Naval Research Laboratory, the Low Power Atmospheric Compensation Experiment (LACE) dynamics experiment, has successfully measured vibrations of an orbiting satellite with a ground-based laser radar. The observations, made on January 7, 8 and 10, 1991, represent the first ever measurements of this type. In the tests, a narrowband heterodyne CO2 laser radar, operating at a wavelength of 10.6 microns, detected vibration induced differential-Doppler signatures of the LACE satellite. Power spectral densities of forced oscillations and modal frequencies and damping rates of free-damped vibrations were obtained and compared with finite element structural models of the LACE system. Another manifested flight experiment is the Advanced Controls Technology Experiment (ACTEX) designed to demonstrate active and passive damping with piezo-electric (PZT) sensors and actuators. This experiment was developed under the management of the Air Force Phillips Laboratory with integration of the experiment at NRL. It is to ride as a secondary, or 'piggyback,' experiment on a future Navy satellite.

Fisher, Shalom↗
NASA NTRS1989

Real-time approximate optimal guidance laws for the advanced launch system

An approach to optimal ascent guidance for a launch vehicle is developed using an expansion technique. The problem is to maximize the payload put into orbit subject to the equations of motion of a rocket over a rotating spherical earth. It is assumed that the thrust and gravitational forces dominate over the aerodynamic forces. It is shown that these forces can be separated by a small parameter epsilon, where epsilon is the ratio of the atmospheric scale height to the radius of the earth. The Hamilton-Jacobi-Bellman or dynamic programming equation is expanded in a series where the zeroth-order term (epsilon = 0) can be obtained in closed form. The zeroth-order problem is that of putting maximum payload into orbit subject to the equations of motion of a rocket in a vacuum over a flat earth. The neglected inertial and aerodynamic terms are included in higher order terms of the expansion, which are determined from the solution of first-order linear partial differential equations requiring only quadrature integrations. These quadrature integrations can be performed rapidly, so that real-time approximate optimization can be used to construct the launch guidance law.

Speyer, Jason L.↗
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.↗
DOE OSTI2017

pyomo.dae: a modeling and automatic discretization framework for optimization with differential and algebraic equations

We describe pyomo.dae, an open source Python-based modeling framework that enables high-level abstract specification of optimization problems with differential and algebraic equations. The pyomo.dae framework is integrated with the Pyomo open source algebraic modeling language, and is available at http://www.pyomo.org. One key feature of pyomo.dae is that it does not restrict users to standard, predefined forms of differential equations, providing a high degree of modeling flexibility and the ability to express constraints that cannot be easily specified in other modeling frameworks. Other key features of pyomo.dae are the ability to specify optimization problems with high-order differential equations and partial differential equations, defined on restricted domain types, and the ability to automatically transform high-level abstract models into finite-dimensional algebraic problems that can be solved with off-the-shelf solvers. Moreover, pyomo.dae users can leverage existing capabilities of Pyomo to embed differential equation models within stochastic and integer programming models and mathematical programs with equilibrium constraint formulations. Collectively, these features enable the exploration of new modeling concepts, discretization schemes, and the benchmarking of state-of-the-art optimization solvers.

97 MATHEMATICS AND COMPUTING↗
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↗
DOE OSTI2017

Physics-based control-oriented modeling of the current density profile evolution in NSTX-Upgrade

Active control of the toroidal current density profile is among those plasma control milestones that the National Spherical Torus eXperiment-Upgrade (NSTX-U) program must achieve to realize its next-step operational goals. Motivated by the coupled, nonlinear, multivariable, distributed-parameter plasma dynamics, the first step towards control design is the development of a physics-based, control-oriented model for the current profile evolution in response to non-inductive current drives and heating systems. The evolution of the toroidal current density profile is closely related to the evolution of the poloidal magnetic flux profile, whose dynamics is modeled by a nonlinear partial differential equation (PDE) referred to as the magnetic-flux diffusion equation (MDE). The proposed control-oriented model predicts the spatial-temporal evolution of the current density profile by combining the nonlinear MDE with physics-based correlations obtained at NSTX-U for the electron density, electron temperature, and non-inductive current drives (neutral beams). The resulting first-principles-driven, control-oriented model is tailored for NSTX-U based on predictions by the time-dependent transport code TRANSP. Furthermore, main objectives and possible challenges associated with the use of the developed model for the design of both feedforward and feedback controllers are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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.↗
DOE OSTI2017

Managing time-substitutable electricity usage using dynamic controls

A predictive-control approach allows an electricity provider to monitor and proactively manage peak and off-peak residential intra-day electricity usage in an emerging smart energy grid using time-dependent dynamic pricing incentives. The daily load is modeled as time-shifted, but cost-differentiated and substitutable, copies of the continuously-consumed electricity resource, and a consumer-choice prediction model is constructed to forecast the corresponding intra-day shares of total daily load according to this model. This is embedded within an optimization framework for managing the daily electricity usage. A series of transformations are employed, including the reformulation-linearization technique (RLT) to obtain a Mixed-Integer Programming (MIP) model representation of the resulting nonlinear optimization problem. In addition, various regulatory and pricing constraints are incorporated in conjunction with the specified profit and capacity utilization objectives.

24 POWER TRANSMISSION AND DISTRIBUTION↗