SEARCH · Search NASA

Results for “differentiable simulation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Examining ICARUS Cosmic Muon Signal Shapes

In the search for new physics, such as sterile neutrinos, we must compare our experimental data to the case where this new physics does not exist, which is provided by simulations. However, our detectors and simulations are not perfect, so we need to be able to differentiate imperfections in our simulations, detector effects, and unknown unknowns from new physics. Thus, we need to quantify ICARUS detector systematic uncertainties; in other words, we need to know how much difference between experiment and simulation we can expect due to only detector systematics, so when we see differences greater than this, we can be confident they are due to new physics. To calculate this uncertainty, we study the signals produced in ICARUS by cosmic muons, since these muons are well understood. Ideally, we want to fit these waveforms as Gaussians and compare the fits from experimental data to fits from simulations to calculate the uncertainties, but first, we need to know how accurately these curve s can be described by Gaussians. We examined the peak and the full width at half maximum (FHWM) of signals from simulations, and we produced plots of the distribution of peaks and FHWMs for these signals. We further studied how the peaks and FHWMs varied depending on where the signal came from in the detector. Ultimately, by comparing the actual distribution of peaks and FHWMs to the distribution predicted by the Gaussian fits, we hope to determine how well these signals can be described as Gaussians.

Patino, Nicolas↗

DOE OSTI2025

End-to-end differentiable digital twin for the IOTA/FAST facility

As the design complexity of modern accelerators grows, there is more interest in using controllable-fidelity simulations that have fast execution time or can yield additional insights about accelerator state. One notable example of additional information are gradients of physical observables with respect to design parameters produced by differentiable simulations. The IOTA/FAST facility has recently begun a program to implement and experimentally validate an end-to-end digital twin to serve as a virtual accelerator test stand, allowing for rapid prototyping of new software and experiments with minimal beam time costs. In this contribution we will discuss our plans and progress. Specifically, we will cover the selection and benchmarking of both physics and ML codes for linac and ring simulation, the development of generic interfaces between surrogate and physics-based sections, and presenting the control interface as either a deterministic event loop or a fully asynchronous EPICS soft input/output controller. We will also discuss challenges in model calibration and uncertainty quantification, as well as future plans to implement larger proton accelerators like PIPII and Booster.

Kuklev, N. [Fermilab]↗

NASA NTRSDate not supplied

Training the Powered-Lift Evaluation Pilot

This poster describes a project to prepare pilots for a study assessing novel aircraft automation concepts for electric Vertical Takeoff and Landing (eVTOL) aircraft using NASA’s Vertical Motion Simulator (VMS). By exploring the operational and learning challenges related to transitioning between forward flight and vertical landing, we seek to establish baselines of pilot workload and aircraft handling qualities across varying atmospheric conditions and automation states. The simulated eVTOL design differentiates flight control allocations as a function of airspeed across four speed ranges as the vehicle transitions between fully thrust-borne lift and wing-borne lift. As speed increases, side stick controls command: translational ground speeds, vertical and lateral acceleration, vertical rate, vertical flight path angle, and bank angle. This novel approach to flight control allocation creates a significant learning challenge for pilots. Since initial eVTOL aircraft may have limitations on hover capabilities, automation and flight guidance cues also vary with airspeed to provide efficient landing profiles while still providing cues suitable for cruise flight. The NASA team prepared the study pilots to follow these flight guidance cues along curved Required Navigation Performance (RNP) approaches and along 6o and 12o glide paths to energy-efficient assistive-hover landing and goarounds. The pre-VMS preparation sought to prepare pilots from diverse levels of experience and background. To do this, NASA researchers designed and developed a fixed-based, large field-ofview simulator with terrain, structures, and air traffic. With one day of combined classroom learning and skill development in the fixedbase simulator, pilots were largely able to fly the simulated eVTOL in the VMS with sufficient mastery to provide handling quality assessments using the Cooper-Harper Handling Qualities Rating and workload assessments through the Bedford Workload Scale.

AAM↗

DOE OSTI2024

Online randomized interpolative decomposition with a posteriori error estimator for temporal PDE data reduction

Traditional low-rank approximation is a powerful tool for compressing large data matrices that arise in simulations of partial differential equations (PDEs), but suffers from high computational cost and requires several passes over the PDE data. The compressed data may also lack interpretability thus making it difficult to identify feature patterns from the original data. Here, to address these issues, we present an online randomized algorithm to compute the interpolative decomposition (ID) of large-scale data matrices in situ. Compared to previous randomized IDs that used the QR decomposition to determine the column basis, we adopt a streaming ridge leverage score-based column subset selection algorithm that dynamically selects proper basis columns from the data and thus avoids an extra pass over the data to compute the coefficient matrix of the ID. In particular, we adopt a single-pass error estimator based on the non-adaptive Hutch++ algorithm to provide real-time error approximation for determining the best coefficients. As a result, our approach only needs a single pass over the original data and thus is suitable for large and high-dimensional matrices stored outside of core memory or generated in PDE simulations. A strategy to improve the accuracy of the reconstructed data gradient, when desired, within the ID framework is also presented. We provide numerical experiments on turbulent channel flow and ignition simulations, and on the NSTX Gas Puff Image dataset, comparing our algorithm with the offline ID algorithm to demonstrate its utility in real-world applications.

Column subset selection↗

NASA NTRS2024

Observing System Simulation Experiments Exploring Potential Spaceborne Deployment Options for a Differential Absorption Radar Measuring Marine Surface Pressures

A new technology for remote measurements of marine surface pressure has been proposed, employing a V-band differential absorption radar and a radiometric temperature sounder to calculate the total column atmospheric mass. Observing System Simulation Experiments (OSSEs) are performed to evaluate the potential impact of Spaceborne Marine Surface Pressure (SMSP) on Numerical Weather Prediction. These experiments build on prior efforts (Privé, McLinden, et al., 2023, https://doi.org/10.16993/tellusa.3254), but with an updated version of the OSSE framework and with more sophisticated simulation of the SMSP observations and a longer experiment period. Several different instrument configurations are compared, including both scanning and non-scanning orbits. SMSP impacts are calculated for analysis quality and forecast skill, and a forecast sensitivity observation impact tool is employed to place SMSP observations in context with the global observing network. The effects of rain contamination on observation quality are explored. Different magnitudes of simulated SMSP observation error are tested in the context of data assimilation to show the range of potential behaviors. Overall, SMSP observations are found to be most beneficial in the southern hemisphere extratropics, with statistically significant forecast improvements for the first 72 hr of the forecast. A constellation of four non-scanning SMSP satellites is found to outperform a single scanning instrument with a 250 km wide swath.

N C Prive↗

DOE OSTI2024

A discontinuous piecewise polynomial generalized moving least squares scheme for robust finite element analysis on arbitrary grids

A variational approach is developed with a meshless discretization to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. Here, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality, which we call discontinuous piecewise polynomial generalized moving least squares (DPP-GMLS). At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework, providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5 x improvement in accuracy and allows for taking an 8 x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

QuadConv: Quadrature-based convolutions with applications to non-uniform PDE data compression

We present a new convolution layer for deep learning architectures which we call QuadConv — an approximation to continuous convolution via quadrature. Our operator is developed explicitly for use on non-uniform, mesh-based data, and accomplishes this by learning a continuous kernel that can be sampled at arbitrary locations. Moreover, the construction of our operator admits an efficient implementation which we detail and construct. As an experimental validation of our operator, we consider the task of compressing partial differential equation (PDE) simulation data from fixed meshes. Here, we show that QuadConv can match the performance of standard discrete convolutions on uniform grid data by comparing a QuadConv autoencoder (QCAE) to a standard convolutional autoencoder (CAE). Further, we show that the QCAE can maintain this accuracy even on non-uniform data. In both cases, QuadConv also outperforms alternative unstructured convolution methods such as graph convolution.

Compression↗

DOE OSTI2022

A Differentiable Hydrology Approach for Modeling With Time-Varying Parameters

In the analysis of observational data with mathematical models, it is often desirable to let one or more model parameters vary with time to account for changing environmental conditions or to offer more flexibility. A standard approach for dynamic models is to allow for the number of parameters to grow as a function of time, leading to a high-dimensional inference problem for long data records. We propose using Hamiltonian Monte Carlo, a gradient-based Markov chain Monte Carlo (MCMC) method, to perform Bayesian inference for time-varying parameters (TVP) in hydrology models. As derivatives of model error functions with regard to parameters are not available in closed form, we implemented the GR4J rainfall-runoff models in Theano/PyMC3, and Jax/NumPyro, allowing for the application of automatic differentiation techniques. A simulation experiment assessing the viability of these methods for recovering underlying temporal variation in parameters indicates that a discrete Gaussian random walk-based prior appears to be best suited across a range of scenarios while the continuous-time Gaussian process compares relatively poorly. In an analysis of 20 years of daily real-world streamflow records from the Model Parameter Estimation Experiment (MOPEX) hydrology data set with increasing data sparsity and at multiple temporal resolutions, we find that TVP inference with gradient-based MCMC is a flexible and reliable approach for analyzing hydrology models with dynamic parameter sets.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Efficient six-dimensional phase space reconstructions from experimental measurements using generative machine learning

Next-generation accelerator concepts, which hinge on the precise shaping of beam distributions, demand equally precise diagnostic methods capable of reconstructing beam distributions within six-dimensional position-momentum spaces. However, the characterization of intricate features within six-dimensional beam distributions using current diagnostic techniques necessitates a substantial number of measurements, using many hours of valuable beam time. Novel phase space reconstruction techniques are needed to reduce the number of measurements required to reconstruct detailed, high-dimensional beam features in order to resolve complex beam phenomena and as a feedback in precision beam shaping applications. In this study, we present a novel approach to reconstructing detailed six-dimensional phase space distributions from experimental measurements using generative machine learning and differentiable beam dynamics simulations. We demonstrate that this approach can be used to resolve six-dimensional phase space distributions from scratch, using basic beam manipulations and as few as 20 two-dimensional measurements of the beam profile. We also demonstrate an application of the reconstruction method in an experimental setting at the Argonne Wakefield Accelerator, where it is able to reconstruct the beam distribution and accurately predict previously unseen measurements 75× faster than previous methods.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Simulating Open Quantum Systems Using Hamiltonian Simulations

We present a novel method to simulate the Lindblad equation, drawing on the relationship between Lindblad dynamics, stochastic differential equations, and Hamiltonian simulations. We derive a sequence of unitary dynamics in an enlarged Hilbert space that can approximate the Lindblad dynamics up to an arbitrarily high order. This unitary representation can then be simulated using a quantum circuit that involves only Hamiltonian simulation and tracing out the ancilla qubits. There is no need for additional postselection in measurement outcomes, ensuring a success probability of one at each stage. Our method can be directly generalized to the time-dependent setting. We provide numerical examples that simulate both time-independent and time-dependent Lindbladian dynamics with accuracy up to the third order. Published by the American Physical Society 2024

Ding, Zhiyan (ORCID:000000018863403X)↗

DOE OSTI2023

Solving multiphysics-based inverse problems with learned surrogates and constraints

Abstract Solving multiphysics-based inverse problems for geological carbon storage monitoring can be challenging when multimodal time-lapse data are expensive to collect and costly to simulate numerically. We overcome these challenges by combining computationally cheap learned surrogates with learned constraints. Not only does this combination lead to vastly improved inversions for the important fluid-flow property, permeability, it also provides a natural platform for inverting multimodal data including well measurements and active-source time-lapse seismic data. By adding a learned constraint, we arrive at a computationally feasible inversion approach that remains accurate. This is accomplished by including a trained deep neural network, known as a normalizing flow, which forces the model iterates to remain in-distribution, thereby safeguarding the accuracy of trained Fourier neural operators that act as surrogates for the computationally expensive multiphase flow simulations involving partial differential equation solves. By means of carefully selected experiments, centered around the problem of geological carbon storage, we demonstrate the efficacy of the proposed constrained optimization method on two different data modalities, namely time-lapse well and time-lapse seismic data. While permeability inversions from both these two modalities have their pluses and minuses, their joint inversion benefits from either, yielding valuable superior permeability inversions and CO 2 plume predictions near, and far away, from the monitoring wells.

Yin, Ziyi (ORCID:0000000250248771)↗

DOE OSTI2020

Regression Based Approach for Robust Finite Element Analysis on Arbitrary Grids. LDRD Final Report

This report summarizes the work performed under a one-year LDRD project aiming to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. In this project, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality. At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5x improvement in accuracy and allows for taking an 8x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method. The report concludes with a brief summary of ongoing projects and collaborations that utilize or extend the products of this work.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1984

Temperature and humidity profiles in the atmosphere from spaceborne lasers: A feasibility study

Computer simulations of the differential absorption lidar technique in a space craft for the purpose of temperature and humidity profiling indicate: (1) Current technology applied to O2 and H2O lines in the .7 to .8 micrometers wavelength band gives sufficiently high signal-to-noise ratios (up to 50 for a single pulse pair) if backscattering by aerosol particles is high, i.e. profiling accurate to 2 K for temperature and 10% for humidity should be feasible within the turbid lower troposphere in 1 km layers and with an averaging over approximately 100 pulses. (2) The impact of short term fluctuations in aerosol particle concentration is too big for a one laser system. Only a two laser system firing at a time lag of about 1 millisecond can surmount these difficulties. (3) The finite width of the laser line and the quasi-random shift of this line introduce tolerable, partly systematic errors.

Grassl, H.↗

NASA NTRS1983

Differential absorption lidar technique for measurement of the atmospheric pressure profile

A new two-wavelength lidar technique for remotely measuring the pressure profile using the trough absorption region between two strong lines in the oxygen A band is described. The theory of integrated vertical path, differential ranging, and horizontal-path pressure measurements is given, with methods to desensitize and correct for temperature effects. The properties of absorption troughs are described and shown to reduce errors due to laser frequency jitter by up to two orders of magnitude. A general analysis, including laser bandwidth effects, demonstrates that pressure measurements with an integrated-vertical-path technique are typically fifty times more accurate than with a differential ranging technique. Simulations show 0.1-0.3 percent accuracy for ground and Shuttle-based pressure-profile and surface-pressure experiments.

Korb, C. L.↗

NASA NTRS1987

The loss and depth of CO2 ice in comet nuclei

An analytical model was developed to simulate the material differentiation of a cometary nucleus composed of water ice, putative unclathrated CO2 ice and silicate dust in specified proportions. Selective sublimation of any free CO2 ice present in a comet would produce a surface layer of water ice and dust overlying the original CO2 rich material. This surface layer reduces the temperature of buried CO2 ice and restricts the outflow of gaseous CO2. On each orbit, water sublimation at smaller heliocentric distances temporarily reduces the thickness of the water ice and dust layer and liberates dust. The model includes the effects of nucleus rotation, arbitrary orientation of the rotation axis, latitude, heat conduction into the deep interior of the nucleus and restriction of CO2 gas outflow by the water ice and dust layer. The effects of the permeability of the surface water ice layer, the nucleus rotation rate, and the latitude were investigated. Comparison of these and similar results with observations could yield information regarding the permeability and chemical composition of cometary material and suggest sampling strategies to minimize fractionation effects.

Fanale, F. P.↗

NASA NTRS1988

Partitioning of regular computation on multiprocessor systems

Problem partitioning of regular computation over two dimensional meshes on multiprocessor systems is examined. The regular computation model considered involves repetitive evaluation of values at each mesh point with local communication. The computational workload and the communication pattern are the same at each mesh point. The regular computation model arises in numerical solutions of partial differential equations and simulations of cellular automata. Given a communication pattern, a systematic way to generate a family of partitions is presented. The influence of various partitioning schemes on performance is compared on the basis of computation to communication ratio.

Lee, Fung Fung↗

NASA NTRS1990

Partitioning of regular computation on multiprocessor systems

Lee, Fung F.↗

NASA NTRS1990

The electron beam instability and turbulence theories

Extensions and practical applications of recent observations of electron beam-plasma interactions are investigated for the range of turbulence theories, extending from quasi-linear to strong turbulence theory, which have been developed on the basis of the Langmuir-wave excitation model. Electron foreshock observations have indicated that linear instability theory must encompass the excitation of waves whose frequencies are substantially different from those of the plasma frequency; the point of departure for such extensions should be a quantitative test of existing theories, and particle simulations conducive to such testing are presented. A step-by-step addition of physical considerations is used in such simulation studies to differentiate among nonlinear turbulence effects.

Dum, C. T.↗