Search NASA⌕ Search

SEARCH · Search NASA

Results for “differentiable 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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

GenASiS Mathematics: Object-oriented manifolds, operations, and solvers for large-scale physics simulations (version 2)

We report GenASiS Mathematics provides modern Fortran classes furnishing extensible object-oriented functionality for the solution of fields governed by selected partial differential equations. The initial release included extensible object-oriented implementations of simple meshes and the evolution of generic conserved currents thereon. This revision - Version 2 of Mathematics - includes significant reorganization and streamlining of these classes, higher-order reconstruction by a different method, a Poisson solver, coarsening to avoid Courant time step limitations near coordinate singularities, and the offloading of computational kernels to GPUs.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1991

Generating Three-Dimensional Grids About Anything

Three-Dimensional Grids About Anything by Poisson's Equation (3DGRAPE) computer program designed to make computational grids in or about almost any shape. Generated by solution of Poisson's differential equations in three dimensions. Program automatically finds its own values for inhomogeneous terms giving near-orthogonality and controlled grid-cell height at boundaries. Grids generated applied to both viscous and inviscid aerodynamic problems, and to problems in other areas of fluid dynamics. Written in 100 percent FORTRAN 77.

Sorenson, Reese L.↗
DOE OSTI2017

Stan : A Probabilistic Programming Language

Stan is a probabilistic programming language for specifying statistical models. A Stan program imperatively defines a log probability function over parameters conditioned on specified data and constants. As of version 2.14.0, Stan provides full Bayesian inference for continuous-variable models through Markov chain Monte Carlo methods such as the No-U-Turn sampler, an adaptive form of Hamiltonian Monte Carlo sampling. Penalized maximum likelihood estimates are calculated using optimization methods such as the limited memory Broyden-Fletcher-Goldfarb-Shanno algorithm. Stan is also a platform for computing log densities and their gradients and Hessians, which can be used in alternative algorithms such as variational Bayes, expectation propagation, and marginal inference using approximate integration. To this end, Stan is set up so that the densities, gradients, and Hessians, along with intermediate quantities of the algorithm such as acceptance probabilities, are easily accessible. Stan can also be called from the command line using the cmdstan package, through R using the rstan package, and through Python using the pystan package. All three interfaces support sampling and optimization-based inference with diagnostics and posterior analysis. rstan and pystan also provide access to log probabilities, gradients, Hessians, parameter transforms, and specialized plotting.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1978

Lidar remote sensing of tropospheric pollutants and trace gases - Programs of NASA Langley Research Center

NASA Langley Research Center is engaged in a number of in-house and contracted programs to develop, evaluate, and apply lidar techniques to remote measurements of pollutant gases, trace molecules, and aerosols. The differential absorption lidar (DIAL) research programs include the development and evaluation of a UV DIAL system for remote measurements of SO2 and O3 concentrations and aerosol dispersion in urban and power-plant stack plumes, a near-IR DIAL system for vertical measurements of water vapor concentrations, and a high-power tunable IR DIAL system which will measure the concentration of gas species that have absorption features in the 1.4- to 4.4-micron wavelength range. The lidar programs which deal with aerosol measurements include the evaluation of a mobile ruby lidar system for use in plume dispersion model verification and the development of a high-spectral-resolution lidar system for measurement of aerosol extinction and backscattering coefficients. System characteristics and measurement sensitivities are discussed, and comparisons are made between simulations and experimental results.

Browell, E. V.↗
DOE OSTI2021

A solution framework for linear PDE-constrained mixed-integer problems

Abstract We present a general numerical solution method for control problems with state variables defined by a linear PDE over a finite set of binary or continuous control variables. We show empirically that a naive approach that applies a numerical discretization scheme to the PDEs to derive constraints for a mixed-integer linear program (MILP) leads to systems that are too large to be solved with state-of-the-art solvers for MILPs, especially if we desire an accurate approximation of the state variables. Our framework comprises two techniques to mitigate the rise of computation times with increasing discretization level: First, the linear system is solved for a basis of the control space in a preprocessing step. Second, certain constraints are just imposed on demand via the IBM ILOG CPLEX feature of a lazy constraint callback. These techniques are compared with an approach where the relations obtained by the discretization of the continuous constraints are directly included in the MILP. We demonstrate our approach on two examples: modeling of the spread of wildfire and the mitigation of water contamination. In both examples the computational results demonstrate that the solution time is significantly reduced by our methods. In particular, the dependence of the computation time on the size of the spatial discretization of the PDE is significantly reduced.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

NO 2 oxidation reactivity and burning mode of diesel particulates

The NO 2 oxidation kinetics and burning mode for diesel particulate from light-duty and medium-duty engines fueled with either ultra low sulfur diesel or soy methyl ester biodiesel blends have been investigated and are shown to be significantly different from oxidation by O 2 . Oxidation kinetics were measured using a flow-through packed bed microreactor for temperature programmed reactions and isothermal differential pulsed oxidation reactions. The burning mode was evaluated using the same reactor system for flowing BET specific surface area measurements and HR-TEM with fringe analysis to evaluate the nanostructure of the nascent and partially oxidized particulates. The low activation energy measured, specific surface area progression with extent of oxidation, HR-TEM images and difference plots of fringe length and tortuosity paint a consistent picture of higher reactivity for NO 2 , which reacts indiscriminately immediately upon contact with the surface, leading to the Zone I or shrinking core type oxidation. In comparison, O 2 oxidation is shown to have relatively lower reactivity, preferentially attacking highly curved lamella, which are more reactive due to bond strain, and short lamella, which have a higher proportion of more reactive edge sites. Furthermore, this preferential oxidation leads to Zone II type oxidation, where solid phase diffusion of oxygen via pores contributes significantly to slowing the overall oxidation rate, by comparison.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2011

Peridynamics with LAMMPS : a user guide.

Peridynamics is a nonlocal extension of classical continuum mechanics. The discrete peridynamic model has the same computational structure as a molecular dynamics model. This document provides a brief overview of the peridynamic model of a continuum, then discusses how the peridynamic model is discretized within LAMMPS. An example problem is also included.

99 GENERAL AND MISCELLANEOUS↗
NASA NTRS1975

Flight instrumentation specification for parameter identification: Program user's guide

A set of four digital computer programs is presented which can be used to investigate the effects of instrumentation errors on the accuracy of aircraft and helicopter stability-and-control derivatives identified from flight test data. The programs assume that the differential equations of motion are linear and consist of small perturbations about a quasi-steady flight condition. It is also assumed that a Newton-Raphson optimization technique is used for identifying the estimates of the parameters. Flow charts and printouts are included.

Mohr, R. L.↗
NASA NTRS1975

An adaptive maneuvering logic computer program for the simulation of one-to-one air-to-air combat. Volume 2: Program description

A detailed description is presented of the computer programs in order to provide an understanding of the mathematical and geometrical relationships as implemented in the programs. The individual sbbroutines and their underlying mathematical relationships are described, and the required input data and the output provided by the program are explained. The relationship of the adaptive maneuvering logic program with the program to drive the differential maneuvering simulator is discussed.

Burgin, G. H.↗
NASA NTRS1978

Estimate of Orbiter static aeroelasticity properties via FLEXSTAB

The theory underlying the use of rigid-model experimental data in FLEXSTAB aeroelasticity calculations is summarized. A description is given of the aerodynamic and structural models used in the calculation of Space Shuttle Orbiter aeroelasticity properties during the re-entry phase. Experience gained with the use of pressure, nonlinear force and moment, and area-ratio correction methods in the Orbiter application is described and typical program output data are included. The inability of the program to correctly estimate differential control surface pressures and to accept external corrections for these terms is shown to be a major limitation in the Orbiter application.

Cavin, R. K., III↗
NASA NTRS1995

Numerical thermal analyses of heat exchangers for the stirling engine application

The Regenerator, Cooler and Heater for the NASA Space Power Research Engine (SPRE) have been analyzed in detail for laminar, incompressible and oscillatory flow conditions. Each component has been analyzed independently and in detail with the regenerator being modeled as two-parallel-plates channel with a solid wall. The ends of the channel are exposed to two reservoir maintained at different temperature thus providing an axial temperature gradient along the channel. The cooler and heater components have been modeled as circular pipes with isothermal walls. Two different types of thermal boundary conditions have been investigated for the cooler and heater, namely, symmetric and asymmetric temperature inflow. In symmetric temperature inflow the flow enters the channel with the same temperature in throughout the velocity cycle whereas, in asymmetric temperature inflow the flow enters with a different temperature in each half cycle. The study was conducted over a wide range of Maximum Reynolds number (RE(max) varying from 75 to 60000, Valensi number (Va) from 2.5 to 800, and relative amplitude of fluid displacement (A(sub r) from 0.357 to 1.34. A two dimensional Finite volume method based on the SIMPLE algorithm was used to solve the governing partial differential equations. Post processing programs were developed to effectively describe the heat transfer mechanism under oscillatory flows. The computer code was validated by comparing with existing analytical solutions for oscillating flows. The thermal field have been studied with the help of temperature contour and three dimensional plots. The instantaneous friction factor, wall heat flux and heat transfer coefficient have been examined. It has been concluded that in general, the frictional factor and heat transfer coefficient are higher under oscillatory flow conditions when the Valensi number is high. Also, the thermal efficiency decreases for lower A(r) values. Further, the usual steady state definition for the heat transfer coefficient does not seem to be valid.

Kannapareddy, Mohan Raj↗
NASA NTRS1961

Fourier Series Operating Package

This report presents a computer program for multiplying, adding, differentiating, integrating, "barring" and scalarly multiplying "literal" Fourier series as such, and for extracting the coefficients of specified terms.

Charnow, Milton L.↗
NASA NTRS1999

Demonstration of Automatically-Generated Adjoint Code for Use in Aerodynamic Shape Optimization

Gradient-based optimization requires accurate derivatives of the objective function and constraints. These gradients may have previously been obtained by manual differentiation of analysis codes, symbolic manipulators, finite-difference approximations, or existing automatic differentiation (AD) tools such as ADIFOR (Automatic Differentiation in FORTRAN). Each of these methods has certain deficiencies, particularly when applied to complex, coupled analyses with many design variables. Recently, a new AD tool called ADJIFOR (Automatic Adjoint Generation in FORTRAN), based upon ADIFOR, was developed and demonstrated. Whereas ADIFOR implements forward-mode (direct) differentiation throughout an analysis program to obtain exact derivatives via the chain rule of calculus, ADJIFOR implements the reverse-mode counterpart of the chain rule to obtain exact adjoint form derivatives from FORTRAN code. Automatically-generated adjoint versions of the widely-used CFL3D computational fluid dynamics (CFD) code and an algebraic wing grid generation code were obtained with just a few hours processing time using the ADJIFOR tool. The codes were verified for accuracy and were shown to compute the exact gradient of the wing lift-to-drag ratio, with respect to any number of shape parameters, in about the time required for 7 to 20 function evaluations. The codes have now been executed on various computers with typical memory and disk space for problems with up to 129 x 65 x 33 grid points, and for hundreds to thousands of independent variables. These adjoint codes are now used in a gradient-based aerodynamic shape optimization problem for a swept, tapered wing. For each design iteration, the optimization package constructs an approximate, linear optimization problem, based upon the current objective function, constraints, and gradient values. The optimizer subroutines are called within a design loop employing the approximate linear problem until an optimum shape is found, the design loop limit is reached, or no further design improvement is possible due to active design variable bounds and/or constraints. The resulting shape parameters are then used by the grid generation code to define a new wing surface and computational grid. The lift-to-drag ratio and its gradient are computed for the new design by the automatically-generated adjoint codes. Several optimization iterations may be required to find an optimum wing shape. Results from two sample cases will be discussed. The reader should note that this work primarily represents a demonstration of use of automatically- generated adjoint code within an aerodynamic shape optimization. As such, little significance is placed upon the actual optimization results, relative to the method for obtaining the results.

Green, Lawrence↗