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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

Development of a commercial enzymes system for lignocellulosic biomass saccharification

DSM Innovation Inc., in its four year effort was able to evaluate and develop its in-house DSM fungal cellulolytic enzymes system to reach enzyme efficiency mandates set by DoE Biomass program MYPP goals. DSM enzyme cocktail is uniquely active at high temperature and acidic pH, offering many benefits and product differentiation in 2G bioethanol production. Under this project, strain and process development, ratio optimization of enzymes, protein and genetic engineering has led to multitudes of improvement in productivity and efficiency making development of a commercial enzyme system for lignocellulosic biomass saccharification viable. DSM is continuing further improvement by additional biodiversity screening, protein engineering and overexpression of enzymes to continue to further lower the cost of enzymes for saccharification of biomass.

09 BIOMASS FUELS↗
NASA NTRS1975

Averaged initial Cartesian coordinates for long lifetime satellite studies

A set of initial Cartesian coordinates, which are free of ambiguities and resonance singularities, is developed to study satellite mission requirements and dispersions over long lifetimes. The method outlined herein possesses two distinct advantages over most other averaging procedures. First, the averaging is carried out numerically using Gaussian quadratures, thus avoiding tedious expansions and the resulting resonances for critical inclinations, etc. Secondly, by using the initial rectangular Cartesian coordinates, conventional, existing acceleration perturbation routines can be absorbed into the program without further modifications, thus making the method easily adaptable to the addition of new perturbation effects. The averaged nonlinear differential equations are integrated by means of a Runge Kutta method. A typical step size of several orbits permits rapid integration of long lifetime orbits in a short computing time.

Pines, S.↗
NASA NTRS1989

Experimental verification of dynamic simulation

The dynamics model here is a backhoe, which is a four degree of freedom manipulator from the dynamics standpoint. Two types of experiment are chosen that can also be simulated by a multibody dynamics simulation program. In the experiment, recorded were the configuration and force histories; that is, velocity and position, and force output and differential pressure change from the hydraulic cylinder, in the time domain. When the experimental force history is used as driving force in the simulation model, the forward dynamics simulation produces a corresponding configuration history. Then, the experimental configuration history is used in the inverse dynamics analysis to generate a corresponding force history. Therefore, two sets of configuration and force histories--one set from experiment, and the other from the simulation that is driven forward and backward with the experimental data--are compared in the time domain. More comparisons are made in regard to the effects of initial conditions, friction, and viscous damping.

Yae, K. Harold↗
NASA NTRS2012

Photogrammetry Measurements During a Tanking Test on the Space Shuttle External Tank, ET-137

On November 5, 2010, a significant foam liberation threat was observed as the Space Shuttle STS-133 launch effort was scrubbed because of a hydrogen leak at the ground umbilical carrier plate. Further investigation revealed the presence of multiple cracks at the tops of stringers in the intertank region of the Space Shuttle External Tank. As part of an instrumented tanking test conducted on December 17, 2010, a three dimensional digital image correlation photogrammetry system was used to measure radial deflections and overall deformations of a section of the intertank region. This paper will describe the experimental challenges that were overcome in order to implement the photogrammetry measurements for the tanking test in support of STS-133. The technique consisted of configuring and installing two pairs of custom stereo camera bars containing calibrated cameras on the 215-ft level of the fixed service structure of Launch Pad 39-A. The cameras were remotely operated from the Launch Control Center 3.5 miles away during the 8 hour duration test, which began before sunrise and lasted through sunset. The complete deformation time history was successfully computed from the acquired images and would prove to play a crucial role in the computer modeling validation efforts supporting the successful completion of the root cause analysis of the cracked stringer problem by the Space Shuttle Program. The resulting data generated included full field fringe plots, data extraction time history analysis, section line spatial analyses and differential stringer peak ]valley motion. Some of the sample results are included with discussion. The resulting data showed that new stringer crack formation did not occur for the panel examined, and that large amounts of displacement in the external tank occurred because of the loads derived from its filling. The measurements acquired were also used to validate computer modeling efforts completed by NASA Marshall Space Flight Center (MSFC).

Littell, Justin D.↗
NASA NTRS2020

Comparison of Raindrop Size Distribution between NASA’s S-Band Polarimetric Radar and Two-Dimensional Video Disdrometers

The ground-based-radar-derived raindrop size distribution (DSD) parameters—mass-weighted drop diameter Dmass and normalized intercept parameter NW—are the sole resource for direct validation of the National Aeronautics and Space Administration (NASA) Global Precipitation Measurement (GPM) mission Core Observatory satellite-based retrieved DSD. Both Dmass and NW are obtained from radar-measured reflectivity ZH and differential reflectivity ZDR through empirical relationships. This study uses existing relationships that were determined for the GPM ground validation (GV) program and directly compares the NASA S-band polarimetric radar (NPOL) observables of ZH and ZDR and derivedDmass and NW with those calculated by two-dimensional video disdrometer (2DVD). The joint NPOL and 2DVD datasets were acquired during three GPM GV field campaigns conducted in eastern Iowa, southern Appalachia, and western Washington State. The comparative study quantifies the level of agreement for ZH, ZDR, Dmass, and log(NW) at an optimum distance (15–40 km) from the radar as well as at distances greater than 60 km from radar and over mountainous terrain. Interestingly, roughly 10%–15% of the NPOL ZH–ZDR pairs were well outside the envelope of 2DVD-estimated ZH–ZDR pairs. The exclusion of these pairs improved the comparisons noticeably.

Ali Tokay↗
NASA NTRS2020

Comparison of Raindrop Size Distribution between NASA's S-Band Polarimetric Radar and Two-Dimensional Video Disdrometers

The ground-based-radar-derived raindrop size distribution (DSD) parameters—mass-weighted drop diameter D mass and normalized intercept parameter NW—are the sole resource for direct validation of the National Aeronautics and Space Administration (NASA) Global Precipitation Measurement (GPM ) mission Core Observatory satellite-based retrieved DSD. Both D mass and NW are obtained from radar-measured reflectivity Z Hand differential reflectivity ZDR through empirical relationships. This study uses existing relationships that were determined for the GPM ground validation (GV) program and directly compares the NASA S-band polarimetric radar (NPOL) observables of Z Hand Z DR and derived D mass and NW with those calculated by two-dimensional video disdrometer (2DVD). The joint NPOL and 2DVD datasets were acquired during three GPM GV field campaigns conducted in eastern Iowa, southern Appalachia, and western Washington State. The comparative study quantifies the level of agreement for ZH, ZDR, D mass, and log (NW) at an optimum distance (15–40 km) from the radar as well as at distances greater than 60 km from radar and over mountainous terrain. Interestingly, roughly 10%–15% of the NPOLZH–Z DR pairs were well outside the envelope of 2DVD-estimatedZH–Z DR pairs. The exclusion of these pairs improved the comparisons noticeably.

Ali Tokay↗
NASA NTRS2000

A Comparison of the Automated Meteorological Profiling System High Resolution Flight Element to the Kennedy Space Center 50 MHz Doppler Wind Profiler

Wind profile measurement and the simulation of aerodynamic loads on a launch vehicle play an important role in determining launch capability and post launch assessment of the vehicle's performance. To date, all United States range certified wind profile measurement systems have been based on balloon tracking. Since the 1960's, the standard used by the National Aeronautics and Space Administration and the Air Force at the Cape Canaveral Air Station (CCAS) for detailed wind profile measurements has been the radar tracked, aerodynamically stabilized Jimsphere balloon system. Currently, the Air Force is nearing certification and operational implementation of the Automated Meteorological Profiling System (AMPS) at CCAS and Vandenburg Air Force Base (VAFB). AMPS uses the Global Positioning System for tracking the Jimsphere balloon. It is anticipated that the AMPS/Jimsphere, named the High Resolution Flight Element (HRFE), will have equivalent, or better resolution than the radar tracked Jimsphere, especially when the balloon is far downrange, at a low elevation angle. By the 1980's, the development of Doppler Wind Profilers (DWP) had become sufficiently advanced to justify an experimental measurement program at Kennedy Space Center (KSC). In 1989 a 50 MHz DWP was installed at KSC. In principal, the 50 MHz DWP has the capability to track the evolution of wind profile dynamics within 5 minutes of a launch. Because of fundamental differences in the measurement technique, there is a significant time and space differential between 50 MHz DWP and HRFE wind profiles. This paper describes a study to quantify these differences from a sample of 50 MHz DWP/HRFE pairs obtained during the AMPS certification test program.

Roberts, Barry C.↗
DOE OSTI2020

Coupled Time-Lapse Full-Waveform Inversion for Subsurface Flow Problems Using Intrusive Automatic Differentiation

We describe a novel framework for estimating subsurface properties, such as rock permeability and porosity, from time-lapse observed seismic data by coupling full-waveform inversion (FWI), subsurface flow processes, and rock physics models. For the inverse modeling, we handle the back propagation of gradients by an intrusive automatic differentiation strategy that offers three levels of user control: (1) At the wave physics level, we adopted the discrete adjoint method in order to use our existing high-performance FWI code; (2) at the rock physics level, we used built-in automatic differentiation operators from the TensorFlow backend; (3) at the flow physics level, we implemented customized partial differential equation (PDE) operators for the multiphase flow equations. The three-level coupled inversion strategy strikes a good balance between computational efficiency and programming efforts, and when the gradients are chained together, it constitutes a coupled inverse system. Our numerical experiments demonstrate that the three-level coupled inverse problem is superior in terms of accuracy to a traditional decoupled inversion strategy. Additionally, our method is able to simultaneously invert for parameters in empirical relationships such as the rock physics models. Our proposed inverted model can be used for reservoir performance prediction and reservoir management/optimization purposes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

CHARGED-CURRENT DOUBLE-DIFFERENTIAL CROSS SECTION MEASUREMENT OF ELECTRON NEUTRINOS IN LOW RECOIL REGION ON HYDROCARBON AT ∼7 GEV

Measuring neutrino mixing angles, especially the CP violation phase and mass hierarchy, is one of the priorities of the neutrino community and long baseline neutrino oscillation programs, such as T2K and NOvA, provided the most precise measurements to date. One of the dominating uncertainties of the measurements is the neutrino interaction model due to the difficulty of modeling the nuclear effects from first principles. As a result, neutrino scattering experiments such as MINERvA were conducted to provide data as input to reduce this uncertainty. This thesis presents a charged-current double-differential cross section measurement of electron neutrinos using the MINERvA detector in terms of available energy (a proxy of energy transfer), three-momentum transfer, and/or transverse momentum of the final state electron. This measurement focuses on the transition region between quasi-elastic and resonance processes, where the multi-nucleon knockout process dominates. The result is compared to a modified GENIE model and similar measurements using muon neutrino samples. In addition, this measurement reproduces the photon-like excess seen by previous MINERvA electron neutrino measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Exponential Decay of Sensitivity in Graph-Structured Nonlinear Programs

We study solution sensitivity for nonlinear programs (NLPs) whose structures are induced by graphs. These NLPs arise in many applications such as dynamic optimization, stochastic optimization, optimization with partial differential equations, and network optimization. We illustrate that for a given pair of nodes, the sensitivity of the primal-dual solution at one node against a data perturbation at the other node decays exponentially with respect to the distance between these two nodes on the graph. In other words, the solution sensitivity decays as one moves away from the perturbation point. This result, which we call exponential decay of sensitivity, holds under the strong second-order sufficiency condition and the linear independence constraint qualification. We also present conditions under which the decay rate remains uniformly bounded; this allows us to characterize the sensitivity behavior of NLPs defined over subgraphs of infinite graphs. The theoretical developments are illustrated with numerical examples.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1980

Steady, Oscillatory, and Unsteady Subsonic and Supersonic Aerodynamics, production version (SOUSSA-P 1.1). Volume 1: Theoretical manual

Recent developments of the Green's function method and the computer program SOUSSA (Steady, Oscillatory, and Unsteady Subsonic and Supersonic Aerodynamics) are reviewed and summarized. Applying the Green's function method to the fully unsteady (transient) potential equation yields an integro-differential-delay equation. With spatial discretization by the finite-element method, this equation is approximated by a set of differential-delay equations in time. Time solution by Laplace transform yields a matrix relating the velocity potential to the normal wash. Premultiplying and postmultiplying by the matrices relating generalized forces to the potential and the normal wash to the generalized coordinates one obtains the matrix of the generalized aerodynamic forces. The frequency and mode-shape dependence of this matrix makes the program SOUSSA useful for multiple frequency and repeated mode-shape evaluations.

Morino, L.↗
NASA NTRS1994

(abstract) Optimal Low Thrust Trajectories Using Differential Inclusion Concepts

Low thrust propulsion systems typically have their greatest benefit for high energy missions or missions with large post-launch maneuver requirements. Missions which have been examined include main belt asteroid rendezvous, comet rendezvous, outer planet and Mercury orbiters, Pluto flyby, and solar probe missions. Low thrust mission design software used to determine these trajectories is based on two distinct formulations of the optimal control problem: the indirect and direct methods. The traditional approach (indirect) is to use the calculus of variations to obtain first order necessary conditions on the states and costates. In contrast, direct methods are conceptually different in that no explicit integration takes place. A direct method based on differential inclusion concepts has been developed and successfully used to compute low thrust trajectories. This new approach removes explicit control dependence from the problem thereby reducing the dimension of the parameter space for the nonlinear programming problem. Also when compared to other direct methods, fewer nonlinear constraints are required to represent the dynamics of the problem.

low↗
NASA NTRS1976

Skylab results - Review and outlook

In the history of space flight, Project Skylab marks the transition from short manned flights to long-time manned orbital operations, from single-purpose spacecraft to multipurpose space stations, and from a time of space exploration to a time of space utilization. In the present paper, some of the Skylab scientific programs in the areas of solar physics, stellar astronomy, space physics, earth observations, life sciences, and materials sciences are reviewed. Visible, UV, and X-ray observations showed the sun's surface to be far more active and more differentiated than was earlier assumed. The unexpected apparition of the comet Kohoutek provided valuable information on the composition of cometary nuclei and tails. Pictures taken with earth-looking cameras have provided insight into specific data analysis and assessment problems. Results of the materials science program offer much promise for future work.

Stuhlinger, E.↗
DOE OSTI2016

Nucleon resonance structure studies via exclusive KY electroproduction

Studying the structure of excited nucleon states employing the electroproduction of exclusive reactions is an important avenue for exploring the nature of the non-perturbative strong interaction. The electrocouplings of $N^*$ states in the mass range below 1.8~GeV have been determined from analyses of CLAS $\pi N$, $\eta N$, and $\pi \pi N$ data. This work has made it clear that consistent results from independent analyses of several exclusive channels with different couplings and non-resonant backgrounds but the same $N^*$ electro excitation amplitudes, is essential to have confidence in the extracted results. In terms of hadronic coupling, many high-lying $N^*$ states preferentially decay through the $\pi \pi N$ channel instead of $\pi N$. Data from the $KY$ channels will therefore be critical to provide an independent analysis to compare the extracted electrocouplings for the high-lying $N^*$ states against those determined from the $\pi N$ and $\pi \pi N$ channels. Lastly, a program to study excited $N^*$ state structure in both non-strange and strange exclusive electroproduction channels using CLAS12 will measure differential cross sections and polarization observables to be used as input to extract the $\gamma_vNN^*$ electrocoupling amplitudes for the most prominent $N^*$ states in the range of invariant energy $W$ up 3~GeV in the virtually unexplored domain of momentum transfers $Q^2$ up to 12~GeV$^2$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

PythonFOAM: In-situ data analyses with OpenFOAM and Python

Here, we outline the development of a general-purpose Python-based data analysis tool for OpenFOAM. Our implementation relies on the construction of OpenFOAM applications that have bindings to data analysis libraries in Python. Double precision data in OpenFOAM is cast to a NumPy array using the NumPy C-API and Python modules may then be used for arbitrary data analysis and manipulation on flow-field information. We highlight how the proposed wrapper may be used for an in-situ online singular value decomposition (SVD) implemented in Python and accessed from the OpenFOAM solver PimpleFOAM. Here, 'in-situ' refers to a programming paradigm that allows for a concurrent computation of the data analysis on the same computational resources utilized for the partial differential equation solver. In addition, to demonstrate parallel deployments, we deploy a distributed SVD, which collects snapshot data across the ranks of a distributed simulation to compute the global left singular vectors. Crucially, both OpenFOAM and Python share the same message passing interface (MPI) communicator for this deployment which allows Python objects and functions to exchange NumPy arrays across ranks. Subsequently, we provide scaling assessments of this distributed SVD on multiple nodes of Intel Broadwell and KNL architectures for canonical test cases such as the large eddy simulations of a backward facing step and a channel flow at friction Reynolds number of 395. Finally, we demonstrate the deployment of a deep neural network for compressing the flow-field information using an autoencoder to demonstrate an ability to use state-of-the-art machine learning tools in the Python ecosystem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Energy Dependent Fission Product Yields (2nd year Milestone and Deliverables)

Introduction to the project: One of the Laboratory missions is to provide consistent, high-precision fission product; yield data critical for testing fission models and maintaining the safety and security of the nation’s nuclear weapons stockpile. However, high-quality, energy-differential fission product yield data is missing for certain actinides and neutron energies important to constrain the new U.S. Nuclear Data Program evaluation effort; It has been shown that the reactor antineutrino anomaly may be at least partially caused by roughly 20 fission products. The fission product yield data is missing or incomplete for many of these isotopes, thus it is necessary to accurately determine these values to better constrain the anomaly; Aitor Bracho is measuring very short-lived (seconds to minutes) fission product yields of 235 U and 239 Pu using monoenergetic neutron beams at E n = 60 and 560 keV; Aitor Bracho is using a direct approach utilizing a state-of-the-art rabbit transfer system, superior HPGe detector, and digital acquisition systems for fission decay measurements. The goals of this project: lop experimental capabilities and data analysis techniques to carry out the gamma-ray spectra analysis necessary for fission product yield calculation; Provide high-precision and energy-dependent fission product data supporting fission theory, neutrino physics, and applied physics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1973

Monitoring vegetation cover on mine dumps with ERTS-1 imagery: Some initial results

ERTS-1 imagery is being used in an attempt to differentiate between mine dumps having varying degrees of vegetative cover. At this stage it is clear that the various mine dumps can be located and identified. Differences in vegetative cover can be seen and measured. Patterns of vegetative growth, some characteristic to particular dumps, can also be seen. It is therefore tentatively concluded that mine dumps can be differentiated with respect to their vegetative cover on the imagery received to date. Subsequent imagery showing seasonal variations should facilitate this program.

Gilbertson, B. P.↗