Search NASA⌕ Search

SEARCH · Search NASA

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 109 records · Page 6

Numerical simulations of cometary dust

Most observations of comets are done photometrically or spectrophotometrically. The interpretation of the aperture-averaged flux is relatively simple for an isotropic, radially expanding coma of infinite extent - the canonical model. However, the interpretation of the observations is not so clear when the motion of the dust is affected by radiation pressure, or when the emission is time-varying and anisotropic. For example, in a sample of CCD images of 10 comets, Jewitt and Meech (1987, Ap.J. 317, 992) found that the photometric profiles of only three comets were consistent, within the observational errors, with the profiles predicted form the canonical model. Photometric observations with large apertures, however, seem to suggest that the canonical model may be quite adequate (c.f. Osip, Schleicher, and Millis, 1992, Icarus 98, 115). The dust itself is characterized by a size distribution, with size dependencies on the expansion velocity, the scattered and thermal radiation, the response to radiation pressure, and probably the density. How good then are the approximations normally used in determining the production rates of the dust when these effects are present? As part of a program to better understand the dynamics of cometary dust and gas, a computer program has been developed which numerically simulates the emission of both dust and gas from a tilted rotating nucleus.

Lien, David J.↗

Dynamics Explorer 1: Energetic Ion Composition Spectrometer (EICS)

The Energetic Ion Composition Spectrometer (EICS) experiment was selected as part of the Dynamics Explorer (DE) Program. One of the primary goals of the DE program was to investigate in detail the plasma physics processes responsible for energizing thermal (approximately 1 eV) ionospheric ions and transporting them to the earth's plasma sheet and distant polar cap. The results of the EICS data analysis (including support of other investigators) and of the archiving efforts supported by this contract are summarized in this document. Also reported are some aspects of our operational support activities.

Shelley, E. G.↗

Alternate space station freedom configuration considerations to accommodate solar dynamic power

The results of a technical audit of the Space Station Freedom Program conducted by the Program Director was announced in early 1989 and included a proposal to use solar dynamic power generation systems to provide primary electrical energy for orbital flight operations rather than photovoltaic solar array systems. To generate the current program baseline power of 75 kW, two or more solar concentrators approximately 50 feet in diameter would be required to replace four pairs of solar arrays whose rectangular blanket size is approximately 200 feet by 30 feet. The photovoltaic power system concept uses solar arrays to generate electricity that is stored in nickel-hydrogen batteries. The proposed concept uses the solar concentrator dishes to reflect and focus the Sun's energy to heat helium-xenon gas to drive electricity generating turbines. The purpose here is to consider the station configuration issues for incorporation of solar dynamic power system components. Key flight dynamic configuration geometry issues are addressed and an assembly sequence scenario is developed.

Deryder, L. J.↗

Deep Learning Explicit Differentiable Predictive Control Laws for Buildings

We present a differentiable predictive control (DPC) methodology for learning constrained control laws for unknown nonlinear systems. DPC poses an approximate solution to multiparametric programming problems emerging from explicit nonlinear model predictive control (MPC). Contrary to approximate MPC, DPC does not require supervision by an expert controller. Instead, a system dynamics model is learned from a small dataset of recorded observations of the perturbed system's dynamics and the control law is optimized offline by interaction with the learned system model. The DPC method is based on two sequential steps, i) system identification using a constrained neural state-space model, and ii) optimization of an explicit control law parametrized by another neural network in closed-loop simulation with the identified neural state-space model. The combination of a differentiable closed-loop system and penalty methods for constraint handling of system outputs and inputs allows us to optimize the control law's parameters directly by backpropagating economic MPC loss through the learned system model. By incorporating domain knowledge and leveraging established techniques from optimal control, our method leverages deep neural networks as nonlinear function approximators for system identification and control while avoiding concomitant costs of intractably large datasets, and computationally expensive over-parametrized models. The scalability, data efficiency, and constrained optimal control capability of the proposed DPC method are demonstrated in simulation using a multi-zone building emulator.

Drgona, Jan↗

Optimization of thermal protection systems for the space shuttle vehicle. Volume 1: Final report

A study performed to continue development of computational techniques for the Space Shuttle Thermal Protection System is reported. The resulting computer code was used to perform some additional optimization studies on several TPS configurations. The program was developed in Fortran 4 for the CDC 6400, and it was converted to Fortran 5 to be used for the Univac 1108. The computational methodology is developed in modular fashion to facilitate changes and updating of the techniques and to allow overlaying the computer code to fit into approximately 131,000 octal words of core storage. The program logic involves subroutines which handle input and output of information between computer and user, thermodynamic stress, dynamic, and weight/estimate analyses of a variety of panel configurations. These include metallic, ablative, RSI (with and without an underlying phase change material), and a thermodynamic analysis only of carbon-carbon systems applied to the leading edge and flat cover panels. Two different thermodynamic analyses are used. The first is a two-dimensional, explicit precedure with variable time steps which is used to describe the behavior of metallic and carbon-carbon leading edges. The second is a one-dimensional implicity technique used to predict temperature in the charring ablator and the noncharring RSI. The latter analysis is performed simply by suppressing the chemical reactions and pyrolysis of the TPS material.

Source record↗

Empirical relation between induced velocity, thrust, and rate of descent of a helicopter rotor as determined by wind-tunnel tests on four model rotors

The empirical relation between the induced velocity, thrust, and rate of vertical descent of a helicopter rotor was calculated from wind tunnel force tests on four model rotors by the application of blade-element theory to the measured values of the thrust, torque, blade angle, and equivalent free-stream rate of descent. The model tests covered the useful range of C(sub t)/sigma(sub e) (where C(sub t) is the thrust coefficient and sigma(sub e) is the effective solidity) and the range of vertical descent from hovering to descent velocities slightly greater than those for autorotation. The three bladed models, each of which had an effective solidity of 0.05 and NACA 0015 blade airfoil sections, were as follows: (1) constant-chord, untwisted blades of 3-ft radius; (2) untwisted blades of 3-ft radius having a 3/1 taper; (3) constant-chord blades of 3-ft radius having a linear twist of 12 degrees (washout) from axis of rotation to tip; and (4) constant-chord, untwisted blades of 2-ft radius. Because of the incorporation of a correction for blade dynamic twist and the use of a method of measuring the approximate equivalent free-stream velocity, it is believed that the data obtained from this program are more applicable to free-flight calculations than the data from previous model tests.

Castles, Walter, Jr.↗

Empirical Relation Between Induced Velocity, Thrust, and Rate of Descent of a Helicopter Rotor as Determined by Wind-tunnel Tests on Four Model Rotors

The empirical relation between the induced velocity, thrust, and rate of vertical descent of a helicopter rotor was calculated from wind tunnel force tests on four model rotors by the application of blade-element theory to the measured values of the thrust, torque, blade angle, and equivalent free-stream rate of descent. The model tests covered the useful range of C(sub t)/sigma(sub e) (where C(sub t) is the thrust coefficient and sigma(sub e) is the effective solidity) and the range of vertical descent from hovering to descent velocities slightly greater than those for autorotation. The three bladed models, each of which had an effective solidity of 0.05 and NACA 0015 blade airfoil sections, were as follows: (1) constant-chord, untwisted blades of 3-ft radius; (2) untwisted blades of 3-ft radius having a 3/1 taper; (3) constant-chord blades of 3-ft radius having a linear twist of 12 degrees (washout) from axis of rotation to tip; and (4) constant-chord, untwisted blades of 2-ft radius. Because of the incorporation of a correction for blade dynamic twist and the use of a method of measuring the approximate equivalent free-stream velocity, it is believed that the data obtained from this program are more applicable to free-flight calculations than the data from previous model tests.

Castles, Walter, Jr↗

Optimum Design of High-Speed Prop-Rotors

An integrated multidisciplinary optimization procedure is developed for application to rotary wing aircraft design. The necessary disciplines such as dynamics, aerodynamics, aeroelasticity, and structures are coupled within a closed-loop optimization process. The procedure developed is applied to address two different problems. The first problem considers the optimization of a helicopter rotor blade and the second problem addresses the optimum design of a high-speed tilting proprotor. In the helicopter blade problem, the objective is to reduce the critical vibratory shear forces and moments at the blade root, without degrading rotor aerodynamic performance and aeroelastic stability. In the case of the high-speed proprotor, the goal is to maximize the propulsive efficiency in high-speed cruise without deteriorating the aeroelastic stability in cruise and the aerodynamic performance in hover. The problems studied involve multiple design objectives; therefore, the optimization problems are formulated using multiobjective design procedures. A comprehensive helicopter analysis code is used for the rotary wing aerodynamic, dynamic and aeroelastic stability analyses and an algorithm developed specifically for these purposes is used for the structural analysis. A nonlinear programming technique coupled with an approximate analysis procedure is used to perform the optimization. The optimum blade designs obtained in each case are compared to corresponding reference designs.

Chattopadhyay, Aditi↗

Multi-Stage Stochastic Programming to Joint Economic Dispatch for Energy and Reserve With Uncertain Renewable Energy

To address the uncertain renewable energy in the day-ahead optimal dispatch of energy and reserve, a multi-stage stochastic programming model is established in this paper to minimize the expected total costs. The uncertainties over the multiple stages are characterized by a scenario tree and the optimal dispatch scheme is cast as a decision tree which guarantees the flexibility to decide the reasonable outputs of generation and the adequate reserves accounting for different realizations of renewable energy. Most importantly, to deal with the “Curse of Dimensionality” of stochastic programming, stochastic dual dynamic programming (SDDP) is employed, which decomposes the original problem into several sub-problems according to the stages. Specifically, the SDDP algorithm performs forward pass and backward pass repeatedly until the convergence criterion is satisfied. At each iteration, the original problem is approximated by creating a linear piecewise function. Besides, an improved convergence criterion is adopted to narrow the optimization gaps. The results on the IEEE 118-bus system and real-life provincial power grid show the effectiveness of the proposed model and method.

economic dispatch↗

GPR_calculator: An on-the-fly surrogate model to accelerate massive nudged elastic band calculations

We present GPR_calculator, a package based on Python and C++ programming languages to build an on-the-fly surrogate model using Gaussian Process Regression (GPR) to approximate computationally expensive electronic structure calculations. The key idea is to dynamically train a GPR model during the simulation that can accurately predict energies and forces with uncertainty quantification. When the uncertainty is high, the costly electronic structure calculation is performed to obtain the ground truth data, which is then used to update the GPR model. To illustrate the effectiveness of GPR_calculator, we demonstrate its application in Nudged Elastic Band (NEB) simulations of surface diffusion and reactions, achieving 3-10 times acceleration compared to pure ab initio calculations. The source code is available at https://github.com/MaterSim/GPR_calculator.

Gaussian process regression↗

A Deep Generative Model for Non-Intrusive Identification of EV Charging Profiles

The proliferation of electric vehicles (EVs) brings environmental benefits and technical challenges to power grids. An identification algorithm which can accurately extract individual EV charging profiles out of widely available smart meter measurements has attracted great interests. This paper proposes a non-intrusive identification framework for EV charging profile extraction, which is driven by deep generative models (DGM). First, the proposed DGM is designed as a representation layer embedded into the Markov process and used to model the joint probability distribution of available time-series data. A novel contribution is to approximate posterior distributions by neural networks whose parameters are obtained by variational inference and supervised learning. Second, the EV charging status is inferred from the DGM via dynamic programming. Lastly, the desired EV charging profile can be reconstructed by the rated power of EV models and inferred status. Compared with the benchmark Hidden Markov Models, the proposed framework can better handle noise in data with less computational complexity and better overall accuracy performances with smaller recall. The proposed framework is validated by numerical experiments on the Pecan Street dataset.

33 ADVANCED PROPULSION SYSTEMS↗

Graphical models for optimal power flow

Optimal power flow (OPF) is the central optimization problem in electric power grids. Although solved routinely in the course of power grid operations, it is known to be strongly NP-hard in general, and weakly NP-hard over tree networks. In this paper, we formulate the optimal power flow problem over tree networks as an inference problem over a tree-structured graphical model where the nodal variables are low-dimensional vectors. We adapt the standard dynamic programming algorithm for inference over a tree-structured graphical model to the OPF problem. Combining this with an interval discretization of the nodal variables, we develop an approximation algorithm for the OPF problem. Further, we use techniques from constraint programming (CP) to perform interval computations and adaptive bound propagation to obtain practically efficient algorithms. Compared to previous algorithms that solve OPF with optimality guarantees using convex relaxations, our approach is able to work for arbitrary tree-structured distribution networks and handle mixed-integer optimization problems. Further, it can be implemented in a distributed message-passing fashion that is scalable and is suitable for “smart grid” applications like control of distributed energy resources. In conclusion, numerical evaluations on several benchmark networks show that practical OPF problems can be solved effectively using this approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A design approach for improving the performance of single-grid planar retarding potential analyzers

Planar retarding potential analyzers (RPAs) have a long flight history and have been included on numerous spaceflight missions including Dynamics Explorer, the Defense Meteorological Satellite Program, and the Communications/Navigation Outage Forecast System. RPAs allow for simultaneous measurement of plasma composition, density, temperature, and the component of the velocity vector normal to the aperture plane. Internal conductive grids are used to approximate ideal potential planes within the instrument, but these grids introduce perturbations to the potential map inside the RPA and cause errors in the measurement of the parameters listed above. A numerical technique is presented herein for minimizing these grid errors for a specific mission by varying the depth and spacing of the grid wires. The example mission selected concentrates on plasma dynamics near the sunset terminator in the equatorial region. The international reference ionosphere model is used to discern the average conditions expected for this mission, and a numerical model of the grid-particle interaction is used to choose a grid design that will best fulfill the mission goals.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Verification of Java Programs using Symbolic Execution and Invariant Generation

Software verification is recognized as an important and difficult problem. We present a norel framework, based on symbolic execution, for the automated verification of software. The framework uses annotations in the form of method specifications an3 loop invariants. We present a novel iterative technique that uses invariant strengthening and approximation for discovering these loop invariants automatically. The technique handles different types of data (e.g. boolean and numeric constraints, dynamically allocated structures and arrays) and it allows for checking universally quantified formulas. Our framework is built on top of the Java PathFinder model checking toolset and it was used for the verification of several non-trivial Java programs.

Pasareanu, Corina↗

Synchronous acceleration with tapered dielectric-lined waveguides

Here, we present a general concept to accelerate non-relativistic charged particles. Our concept employs an adiabatically-tapered dielectric-lined waveguide which supports accelerating phase velocities for synchronous acceleration. We propose an ansatz for the transient field equations, show it satisfies Maxwell's equations under an adiabatic approximation and find excellent agreement with a finite-difference time-domain computer simulation. The fields were implemented into the particle-tracking program {\sc astra} and we present beam dynamics results for an accelerating field with a 1-mm-wavelength and peak electric field of 100~MV/m. The numerical simulations indicate that a $\sim 200$-keV electron beam can be accelerated to an energy of $\sim10$~MeV over $\sim 10$~cm. The novel scheme is also found to form electron beams with parameters of interest to a wide range of applications including, e.g., future advanced accelerators, and ultra-fast electron diffraction.

43 PARTICLE ACCELERATORS↗

Geologic Mapping of Athabasca Valles

Two factors drive us to map the Athabasca Valles area in unusual detail: (1) the extremely well-preserved and exposed surface morphologies and (2) the extensive high resolution imaging. In particular, the near-complete CTX coverage of Athabasca Valles proper and the extensive coverage of its surroundings have been invaluable. The mapping has been done exclusively in ArcGIS, using individual CTX, THEMIS VIS, and MOC frames overlying the THEMIS IR daytime basemap. MOLA shot points and gridded DTMs are also included. It was found that CTX images processed through ISIS are almost always within 300 m of the MOLA derived locations, and usually within tens of meters, with no adjustments to camera pointing. THEMIS VIS images appear to be systematically shifted to the southwest of their correct positions and MOC images are often kilometers off. The good SNR and minimal artifacts make the CTX images vastly more useful than the THEMIS VIS or MOC images. The bulk of the mapping was done at 1:50,000 scale on CTX images. In more complex areas, mapping at 1:24,000 proved necessary. The CTX images were usually simultaneously viewed on a second monitor using the ISIS3 qview program to display the full dynamic range of the CTX data. Where CTX data was not available, mapping was often done at 1:100,000 and most contacts are mapped as approximate.

Keszthelyi, L. P.↗

Dynamic loads analysis system (DYLOFLEX) summary. Volume 1: Engineering formulation

The DYLOFLEX computer program system expands the aeroelastic cycle from that in the FLEXSTAB computer program system to include dynamic loads analyses involving active controls. Two aerodynamic options exist within DYLOFLEX. The analyst can formulate the problem with unsteady aerodynamics calculated using the doublet lattice method or with quasi-steady aerodynamics formulated from either FLEXSTAB or doublet lattice steady state aerodynamics with unsteady effects approximated by indicial lift growth functions. The equations of motion are formulated assuming straight and level flight and small motions. Loads are calculated using the force summation technique. DYLOFLEX consists of nine standalone programs which can be linked with each other by magnetic files used to transmit the required data between programs.

Miller, R. D.↗