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 451 records · Page 25

Differentiable Multiphysics Codes: A Breakthrough Technology for Simulation and Computing

This document summarizes the findings of a strategic planning exercise commissioned by the Weapons Simulation and Computing, Computational Physics (WSC/CP) program at the Lawrence Livermore National Laboratory (LLNL) in FY24. During the year, the committee met with multiple stakeholder communities to gather input, opinions, suggestions and concerns which have been incorporated throughout this document. The key findings from this exercise are summarized: • The development of multiphysics modelling and simulation (mod/sim) codes and software technologies, their deployment on exascale compute platforms, and their broad adoption across the NNSA is a major success of the Advanced Simulation and Computing (ASC) program and the Exascale Computing Project (ECP). Sustained investment in these core technologies is essential. • Today’s state of the art involves running ensembles of O(100K) simulations to perform uncertainty quantification (UQ) and design studies using multiple statistical methods such as Bayesian optimization to understand sensitivities of our models and explore parameterized design spaces. Even with exascale computing, we are practically limited to O(10) parameters in these studies since the number of simulations required to sample the space scales exponentially with the number of design parameters. • The data from these simulation ensembles is increasingly being used to train machine learned (ML) surrogates (or reduced order models, ROMs) which can then be used for optimization or real time design exploration. However, the trained surrogates are still limited in the number of parameters they can represent due to the sampling limitations previously noted. • Augmenting our suite of integrated multiphysics simulation codes, both current and emerging, with the ability to compute gradients (solution derivatives) of arbitrary simulation outputs with respect to (some or all) simulation inputs would be a breakthrough technology, opening the door to a new era of efficient and automated inverse design based on verified and validated mod/sim capabilities. • This capability, which we refer to as differentiable multiphysics codes (DMCs), would revolutionize both UQ and optimization studies by breaking the curse of dimensionality that presently limits our “gradient-free” ensemble based computing approach. A similar breakthrough occurred in the AI/ML community once the ability to compute gradients of arbitrary loss functions using back-propagation became commonplace. Gradient information from the multiphysics codes can also be used to dramatically improve the efficiency and scale of training of ML/ROM surrogates for rapid assessments. • Achieving this in our suite of codes will be a grand challenge, similar to the amount of effort that was required to transition from CPU to GPU computing. It will require buy-in from the entire WSC/CP program and beyond, including all integrated codes, physics and engineering models, third-party library dependencies and performance portability abstractions. It will also require investment in research and development of numerical methods for computing adjoints of coupled physics across multiple adaptively refined moving meshes and of stochastic (Monte Carlo) and mesh free (SPH) methods. • New software and numerical techniques, largely pioneered by the AI/ML community, make this feasible. Chief among these is automatic differentiation (AD), the ability to employ AD at point-wise locations in a physics calculation (instead of traditional black-box approaches) and the ability to perform “back-propagation in time” (or reverse mode AD) for non-linear partial differential equations (PDEs). Fundamentally, the conclusion of this strategic planning exercise is that the time is right to undertake a large scale effort in WSC, centered on the existing integrated codes, to continue the natural evolution of mod/sim in the age of AI/ML. Instead of attempting to replace mod/sim with purely data driven AI/ML models, we believe the key to success is to integrate AI/ML by building on top of the decades of hard-won knowledge and the verified/validated multiphysics modelling capability that is the hallmark of the ASC program.

97 MATHEMATICS AND COMPUTING↗

Psyche Early Project Verification & Validation Planning Development

The Psyche mission to the asteroid (16) Psyche was selected as the fourteenth mission in the Discovery program in January 2017. The Psyche mission will determine if (16) Psyche is the core of a larger differentiated body. As part of the development of this mission a Verification and Validation (V&V) engineer was assigned early in the project’s design Phase B. This paper will discuss some of the strategies that the Psyche team is using to take full advantage of the early planning for V&V on flight projects and specifically how the Psyche mission is approaching these tasks. This paper will discuss the effects of having a V&V mindset on: 1) The Psyche requirements development process, and how focusing not only on how the team will verify these requirements but also on developing the tools necessary to track and monitor that verification feeds back into the requirement development process. 2) The Psyche testbed development, and how using a verification mindset is useful for identifying holes in the testbed development process, including the required testbed speed and how to think about testbed certification early in the process while encouraging trades and developing relationships with the testbed team. 3) Developing new V&V Tools for Psyche, and how developing tools early in the project development process means that they can influence the development of other requirement and scheduling tools. 4) Verification Activity Planning, which is typically done to a preliminary level during Phase B of the project, and is important for understanding the major testing that is needed to ensure that the system as built represents the design. On Psyche we are not only focused on bringing this planning to a preliminary level, we are also developing a V&V focused schedule to help us de-conflict V&V activities that may require similar resources early on in the program’s development. 5) Validation Planning, for which we are working with the testbed and model development teams to ensure that their models can be effectively validated and that the plans to do so are in place. In this paper we will describe how the Psyche mission is approaching each of these V&V areas and identify lessons that can be taken by other space missions trying to decide how much effort should be put into V&V early in the project lifecycle.

Solish, Benjamin↗

Computer program documentation for the dynamic analysis of a noncontacting mechanical face seal

A computer program is presented which achieves a numerical solution for the equations of motion of a noncontacting mechanical face seal. The flexibly-mounted primary seal ring motion is expressed by a set of second order differential equations for three degrees of freedom. These equations are reduced to a set of first order equations and the GEAR software package is used to solve the set of first order equations. Program input includes seal design parameters and seal operating conditions. Output from the program includes velocities and displacements of the seal ring about the axis of an inertial reference system. One example problem is described.

Auer, B. M.↗

Modeling the missile-launch tube problem in DYSCO

DYSCO is a versatile, general purpose dynamic analysis program which assembles equations and solves dynamics problems. The executive manages a library of technology modules which contain routines that compute the matrix coefficients of the second order ordinary differential equations of the components. The executive performs the coupling of the equations of the components and manages the solution of the coupled equations. Any new component representation may be added to the library if, given the state vector, a FORTRAN program can be written to compute M, C, K, and F. The problem described demonstrates the generality of this statement.

Berman, Alex↗

Computer simulation program is adaptable to industrial processes

The Reaction kinetics ablation program /REKAP/, developed to simulate ablation of various materials, provides mathematical formulations for computer programs which can simulate certain industrial processes. The programs are based on the use of nonsymmetrical difference equations that are employed to solve complex partial differential equation systems.

Schultz, F. E.↗

Effect of gravity on methane-air combustion

Analytical and numerical techniques dealing with the theoretical description of the influence of zero and reduced gravitational acceleration on diffusion flames, with a view to improving understanding of fires in space vehicles, were developed in support of experimental work performed in this area. This was done in order to confirm qualitative understanding of the process, to determine the quantitative accuracy of numerical predictions, and to establish a mathematical model of the process for subsequent use as a predictive and exploratory tool. The following results were accomplished: (1) derivation of differential equations and boundary conditions describing the system, (2) details of the computations, using a FORTRAN computer program, for calculating the flow and heat and mass transfer in two dimensions (both steady and unsteady). It was shown that the experimental behavior can be reproduced with fair accuracy, provided that the time step is sufficiently short.

Elghobashi, S.↗

Conductive spacecraft materials development program

The objectives of this program are to provide design criteria, techniques, materials, and test methods to ensure control of absolute and differential charging of spacecraft surfaces. The control of absolute and differential charging of spacecraft cannot be effected without the development of new and improved or modified materials or techniques that will provide electrical continuity over the surface of the spacecraft. The materials' photoemission, secondary emission, thermooptical, physical, and electrical properties in the space vacuum environment both in the presence and absence of electrical stress and ultraviolet, electron, and particulate radiation, are important to the achievement of charge control. The materials must be stable or have predictable response to exposure to the space environment for long periods of time. The materials of interest include conductive polymers, paints, transparent films and coatings as well as fabric coating interweaves.

Lehn, W. L.↗

Solid polymer electrolyte water electrolysis preprototype subsystem

Hardware and controls developed for an electrolysis demonstration unit for use with the life sciences payload program and in NASA's regenerative life support evaluation program are described. Components discussed include: the electrolysis module; power conditioner; phase separator-pump and hydrogen differential regulator; pressure regulation of O2, He, and N2; air-cooled heat exchanger; water accumulator; fluid flow sight gage assembly; catalytic O2/H2 sensor; gas flow sensors; low voltage power supply; 100 Amp DC contactor assembly; and the water purifier design.

Source record↗

Research through simulation

The design of the computer operating system at Langley Research Center allows for concurrent support of time-critical simulations and background analytical computing on the same machine. Signal path interconnections between computing hardware and flight simulation hardware is provided to allow up to six simulation programs to be in operation at one time. Capabilities and research applications are discussed for the: (1) differential maneuvering simulator; (2) visual motion simulator; (3) terminal configured vehicle simulator; (4) general aviation aircraft simulator; (5) general purpose fixed based simulator; (6) transport simulator; (7) digital fly by wire simulator; (8) general purpose fighter simulator; and (9) the roll-up cockpit. The visual landing display system and graphics display system are described and their simulator support applications are listed.

Copeland, J. L.↗

High voltage space plasma interactions

Two primary problems resulted from plasma interactions; one of concern to operations in geosynchronous orbit (GEO), the other in low orbits (LEO). The two problems are not the same. Spacecraft charging has become widely recognized as a problem, particularly for communications satellites operating in GEO. The very thin thermal plasmas at GEO are insufficient to bleed off voltage buildups due to higher energy charged particle radiation collected on outer surfaces. Resulting differential charging/discharging causes electrical transients, spurious command signals and possible direct overload damage. An extensive NASA/Air Force program has been underway for several years to address this problem. At lower altitudes, the denser plasmas of the plasmasphere/ionosphere provide sufficient thermal current to limit such charging to a few volts or less. Unfortunately, these thermal plasma currents which solve the GEO spacecraft charging problem can become large enough to cause just the opposite problem in LEO.

Mccoy, J. E.↗

Analytical approach to peel stresses in bonded composite stiffened panels

A closed-form solution was obtained for the stresses and displacements of two bonded beams. A system of two fourth-order and two second-order differential equations with the associated boundary equations was determined using a variational work approach. A FORTRAN computer program was devised to solve for the eigenvalues and eigenvectors of this system and to calculate the coefficients from the boundary conditions. The results were then compared with NASTRAN finite-element solutions and shown to agree closely.

Barkey, Derek A.↗

User's guide to PMESH: A grid-generation program for single-rotation and counterrotation advanced turboprops

A detailed operating manual is presented for a grid generating program that produces 3-D meshes for advanced turboprops. The code uses both algebraic and elliptic partial differential equation methods to generate single rotation and counterrotation, H or C type meshes for the z - r planes and H type for the z - theta planes. The code allows easy specification of geometrical constraints (such as blade angle, location of bounding surfaces, etc.), mesh control parameters (point distribution near blades and nacelle, number of grid points desired, etc.), and it has good runtime diagnostics. An overview is provided of the mesh generation procedure, sample input dataset with detailed explanation of all input, and example meshes.

Warsi, Saif A.↗

Microphysical and radiative characteristics of convective clouds during COHMEX

This paper has two purposes: first, to document the structure and evolution of two strong thunderstorms in Alabama using radar multiparameter data; and second, to relate the inferred microphysics to the resulting upwelling microwave radiance observed concurrently by high altitude aircraft. These measurements were collected during the COHMEX field program in the summer of 1986. The radar analysis includes a description of the parameters reflectivity Z, differential reflectivity Z(DR), linear depolarization ratio LDR, and hail signal HS for two thunderstorm cases on July 11, 1986. The simultaneous aircraft data include passive microwave brightness temperature, T(B), measurements at four frequencies ranging from 18 to 183 GHz, as well as visible and infrared data. The remote radar observations reveal the existence of large ice particles within the storms, which is likely to have caused the observed low microwave brightness temperatures. By relating the evolution of the radar measureables to the microwave T(B)s it has been found that knowledge of the storm microphysics and its evolution is important to adequately understand the microwave T(B)s.

Fulton, Richard↗

Continuous Optimization on Constraint Manifolds

This paper demonstrates continuous optimization on the differentiable manifold formed by continuous constraint functions. The first order tensor geodesic differential equation is solved on the manifold in both numerical and closed analytic form for simple nonlinear programs. Advantages and disadvantages with respect to conventional optimization techniques are discussed.

Dean, Edwin B.↗

Three-Dimensional Finite Element Ablative Thermal Response and Thermostructural Design of Thermal Protection Systems

A finite element ablation and thermal response program is presented for simulation of three-dimensional transient thermostructural analysis. The three-dimensional governing differential equations and finite element formulation are summarized. A novel probabilistic design methodology for thermal protection systems is presented. The design methodology is an eight step process beginning with a parameter sensitivity study and is followed by a deterministic analysis whereby an optimum design can determined. The design process concludes with a Monte Carlo simulation where the probabilities of exceeding design specifications are estimated. The design methodology is demonstrated by applying the methodology to the carbon phenolic compression pads of the Crew Exploration Vehicle. The maximum allowed values of bondline temperature and tensile stress are used as the design specifications in this study.

Dec, John A.↗

First Temperature Profile of a Stellar Flare Using Differential Chromatic Refraction

We present the first derivation of a stellar flare temperature profile from single-band photometry. Stellar flare DWF 030225.574−545707.45129 was detected in 2015 by the Dark Energy Camera as part of the Deeper, Wider, Faster program. The brightness (Δm g = −6.12) of this flare, combined with the high air mass (1.45 ≲ X ≲ 1.75) and blue filter (DES g, 398–548 nm) in which it was observed, provided ideal conditions to measure the zenithward apparent motion of the source due to differential chromatic refraction (DCR) and, from that, infer the effective temperature of the event. We model the flare’s spectral energy distribution as a blackbody to produce the constraints on flare temperature and geometric properties derived from single-band photometry. We additionally demonstrate how simplistic assumptions on the flaring spectrum, as well as on the evolution of flare geometry, can result in solutions that overestimate the effective temperature. Exploiting DCR enables studying chromatic phenomena with ground-based astrophysical surveys, and stellar flares on M dwarfs are a particularly enticing target for such studies due to their ubiquity across the sky and the heightened color contrast between their red quiescent photospheres and the blue flare emission. Our novel method will enable similar temperature constraints for a large sample of objects in upcoming photometric surveys like the Vera C. Rubin Legacy Survey of Space and Time.

Clarke, Riley W. [University of Delaware, Newark, ↗

Program for standard statistical distributions

Development of procedure to describe frequency distributions involved in statistical theory is discussed. Representation of frequency distributions by first order differential equation is presented. Classification of various types of distributions based on Pearson parameters is analyzed.

Falls, L. W.↗

DOE/NNSA Nuclear Criticality Safety Program – Brief Overview [Slides]

This presentation discusses nuclear data measurements and evaluation work for the NCSP. The objective of the NCSP is to provide measured and evaluated thermal, resonance, unresolved resonance, and fast region cross section data to address the priority NCSP nuclear data needs. It also discusses its vision which is to address multiple Nuclear Data 5- and 10-year goals and attributes identified in the NCSP. The intended final product is for rigorous ENDF/B evaluations to be produced from cross section measurements and analyses. It also discusses NCSP integral experiments and differential experiments in addition to NCSP support of SCALE in the future.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗