SEARCH · Engineering Papers

Results for “Differentiable physics”

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 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2022

The role of geographical spreaders in infectious pattern formation and front propagation speeds

The pattern formation and spatial spread of infectious populations are investigated using a kernel-based Susceptible–Infectious–Recovered (SIR) model applicable across a wide range of basic reproduction numbers R O . The goal is to examine the role of geographical spreaders on transient spatial pattern formation of infectious populations and the associated maximum invasive front speeds c max . In the simulations conducted here, geographical spreaders are defined as a portion of the infected population Φ experiencing high mobility between identical communities. The spatial organization of the infected population and c max are determined when the infections are randomly initiated in space within multiple communities. For small but finite , scaling analysis and numerical simulations in 1-dimension suggest that when the spreading kernel is Gaussian-shaped, where is the inverse of the infectious duration. This finding for agrees with a diffusion-based representation of mobility in 1-D. Numerical simulations in 2-D across wide-ranging suggest that , the variance of the spatial kernel describing mobility of long-distance geographical spreaders across communities, determines the spatial organization of infections across communities. When (long-distance mobility, where is the minimum spatial extent defining adjacent communities), the infectious population will experience a transient but spatially coherent pattern with a wavelength that can be derived from the spreading kernel properties. Moreover, the 2-D simulations for the bounded kernel suggest that attainment of is also dictated by but the magnitude is not sensitive to unlike diffusion-based models.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

Stiff neural ordinary differential equations

Neural Ordinary Differential Equations (ODEs) are a promising approach to learn dynamical models from time-series data in science and engineering applications. This work aims at learning neural ODEs for stiff systems, which are usually raised from chemical kinetic modeling in chemical and biological systems. We first show the challenges of learning neural ODEs in the classical stiff ODE systems of Robertson’s problem and propose techniques to mitigate the challenges associated with scale separations in stiff systems. We then present successful demonstrations in stiff systems of Robertson’s problem and an air pollution problem. The demonstrations show that the usage of deep networks with rectified activations, proper scaling of the network outputs as well as loss functions, and stabilized gradient calculations are the key techniques enabling the learning of stiff neural ODEs. The success of learning stiff neural ODEs opens up possibilities of using neural ODEs in applications with widely varying time-scales, such as chemical dynamics in energy conversion, environmental engineering, and life sciences.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Energy transport between heat baths with oscillating temperatures

Not Available

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Physics Guided Simulation of Electrostatic Discharge: Technical Report

Triboelectrically-charged objects may create threshold sparks, electrostatic discharge (ESD) events, to equilibrate charge between themselves and other relatively charged objects. ESD events exhibit many complex physical phenomena. They are a nexus of several fields of physics with disparate characteristic scales: plasma physics, chemical kinetics, hydrodynamics, circuit models, etc. These scales can span many orders of magnitude from the varied collisions thermalizing information within a plasma on the $\mathcal{O}(fs/ps)$ to the physical size of the plasma channel on the $\mathcal{O}(100µm)$, to the speed of a nonlinear hydrodynamic wave propagating at $\mathcal{O}(µm, ns)$. These threshold ESD events may occur in situations of programmatic importance, delivering energy and power profiles to a “victim load” generating deleterious consequences. To predict and mitigate these consequences we must answer questions about the spark’s energy budget: how much energy goes into producing the spark channel; how much gets radiated away; how much energy is advected away into the hydrodynamics; and how much energy is delivered to a victim load. An ESD simulation toolset has been created and evolved in order to answer these questions. An appropriate, physics guided implementation for simulation can be done by gaining insight into its constituent physics and leveraging that intuition to choose a suitable numerical operator. We examine in detail the chemical kinetics, circuit discharge, and hydrodynamics to deter mine dominant regimes, values, timescales, and interactions to uncover the underlying physical dynamics. We also examine and propose model reduction schemes for high-dimensional chemical kinetics. We use past and current work with experimentally validated and theoretically-verified hydrodynamics to calculate applicability limits of the non-ionizing strong shock limit. We quantify the energy budget from a hydrodynamic perspective and demonstrate that a significant fraction of the stored energy is “earmarked” for hydrodynamic advection as an energy terminus. Lastly, we combine the constituent physics of an ESD event (chemical kinetics, circuit model, and hydrodynamics) into a cohesive, actionable toolset and obtain promising results from an isothermal test case. We then propose a viable, modular evolution of the ESD toolset based upon the performed examination of the physics uncovering dominant physical scales and the stiffness of the compositional differential system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1971

Similarity and generalized finite-difference solutions of parabolic partial differential equations.

Techniques are presented for obtaining generalized finite-difference solutions to partial differential equations of the parabolic type. It is shown that the advantages of similarity in the solution of similar problems are generally not lost if the solution to the original partial differential equations is effected in the physical plane by finite-difference methods. The analysis results in a considerable saving in computational effort in the solution of both similar and nonsimilar problems. Several examples, including both the heat-conduction equation and the boundary-layer equations, are given. The analysis also provides a practical means of estimating the accuracy of finite-difference solutions to parabolic equations.

Clausing, A. M.↗

NASA NTRS1975

Mixed boundary value problems in mechanics

Certain boundary value problems were studied over a domain D which may contain the point at infinity and may be multiply connected. Contours forming the boundary are assumed to consist of piecewise smooth arcs. Mixed boundary value problems are those with points of flux singularity on the boundary; these are points on the surface, either side of which at least one of the differential operator has different behavior. The physical system was considered to be described by two quantities, the potential and the flux type quantities. Some of the examples that were illustrated included problems in potential theory and elasticity.

Erdogan, F.↗

NASA NTRS1977

Cryogenic temperature control by means of energy storage materials

An investigation was conducted to study the concept of thermal control by means of physical or chemical reaction heats for applications involving the storage of cryogens during long-term space voyages. The investigation included some preliminary experimental tests of energy storage material (ESM) effectiveness. The materials considered can store and liberate large amounts of thermal energy by means of mechanisms such as sensible heat, heat of fusion, and physical or chemical reaction heat. A differential thermal analysis was utilized in the laboratory tests. Attention is given to the evaluation of cryogenic ESM thermal control concepts, the experimental determination of phase change materials characteristics, and adsorption ESMs. It is found that an ESM shield surrounded by multiple layer insulation provides the best protection for a cryogen store.

Grodzka, P. G.↗

NASA NTRS1979

A fast semi-implicit algorithm for problems of mixed type

Certain physical processes are modeled by partial differential equations which are parabolic over part of the domain and elliptic over the remainder. A family of semi-implicit algorithms which are well suited to initial-boundary value problems of this mixed type is discussed. One important feature of these algorithms is the use of an approximate inverse for the solution of the implicit linear system. A strong error analysis results in an estimate of the total error as a function of approximate inverse error e and time step h.

Frederickson, P. O.↗

NASA NTRS1960

The Interactions Between Nitrogen and Oxygen Molecules

Lippincott's delta-function model for atomic interactions is analyzed, both physically and mathematically, and extended, by differentiation between K- and L-shell electrons and the introduction of a variable parameter in the expression for the delta-function strength, to cover homonuclear molecules more complex than hydrogen. In addition, modifications are made which allow treatments of diatomic, heteronuclear molecules. This theory, in conjunction with a reasonably extensive study of resonance, dispersion, and configuration interaction phenomena, as well as the use of simple quantum mechanical arguments, is then applied to the N2-N2, N2-O2, and O2-O2 interactions.

Meador, Willard E., Jr.↗

DOE OSTI2022

Comparative critical mass calculations for NNL and ENDF/B-VIII.0 - Zirconium hydride thermal neutron scattering laws

Zirconium hydride (ZrH{sub x}) is a moderator material for TRIGA reactors and historical space reactor systems, such as SNAP-10A. Thermal neutron scattering laws (TSL) for two phases of this material, δ and ε, have been previously evaluated by Naval Nuclear Laboratory (NNL) and submitted to the National Nuclear Data Center (NNDC) for inclusion in the US national ENDF/B-VIII.1 nuclear data library. In contrast to the current ENDF/B-VIII.0 TSL evaluations, which consider only a single phase, the new evaluations are derived from separate ab initio calculations for both phases and include coherent elastic effects of the zirconium sublattice. To estimate the impact of these changes to the TSL evaluation of this material, comparative critical mass calculations were performed with MC21 for homogenous mixtures of high-enriched uranium (HEU) and ZrH{sub x} in bare and water reflected sphere configurations. These calculations yield an impact on the estimated critical mass as a function of {sup 235}U loading density with maximum differences as large as 1% - 5% for over-moderated thermal spectrum systems. Consequently, the NNL TSL evaluations are anticipated to have a small impact on criticality calculations of thermal reactor systems regardless of the material phase. Nevertheless, characteristic differences exist in the predicted thermal spectra as function of energy for the two sets of TSL evaluations, which are attributed to difference in the underlying phonon density of states of hydrogen bound in ZrH{sub x}. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Development of a phonon-based sampling method for thermal neutron scattering data

Simulations of reactor systems require access to accurate nuclear data. For many systems, thermal neutron scattering data can have large effects on the eigenvalue and neutron flux distributions. Inelastic thermal neutron scattering can excite or de-excite vibrational, rotational, and translational modes in a material, so thermal scattering evaluations are often obtained by summing over the number of phonons created/destroyed by a scattering event. In recent years, the thermal scattering cross sections and angular distributions have greatly improved in accuracy, but the format in which this data is delivered to simulation codes has remained virtually unchanged. Thermal scattering data is typically either compiled into large tables and sorted by incoming neutron energy, outgoing neutron energy, scattering angle, and material temperature, or represented as cumulative distribution functions of momentum exchange or energy exchange. Either method can be quite memory intensive when fine bins are used. In an effort to decrease the amount of space that processed thermal scattering data requires, an alternate format is proposed. The phonon-based sampling method introduced here can sample the number of phonons excited for each collision, the change in neutron energy, and the scattering angle while avoiding pre-computed angular bins and limiting the amount of data that is dependent on incoming energy. Through this method, the generation and storage of large interpolation tables is avoided, which could have benefits in both memory storage and accuracy. While the initial implementation of this method is slower than current alternatives, it is significantly more resistant to grid coarseness errors and has good potential for improvement. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1990

Galaxy interactions and strength of nuclear activity

Analysis of data in the literature for differential velocities and projected separations of nearby Seyfert galaxies with possible companions shows a clear difference in projected separations between type 1's and type 2's. This kinematic difference between the two activity classes reinforces other independent evidence that their different nuclear characteristics are related to a non-nuclear physical distinction between the two classes. The differential velocities and projected separations of the galaxy pairs in this sample yield mean galaxy masses, sizes, and mass to light ratios which are consistent with those found by the statistical methods of Karachentsev. Although the galaxy sample discussed here is too small and too poorly defined to provide robust support for these conclusions, the results strongly suggest that nuclear activity in Seyfert galaxies is associated with gravitational perturbations from companion galaxies, and that there are physical distinctions between the host companions of Seyfert 1 and Seyfert 2 nuclei which may depend both on the environment and the structure of the host galaxy itself.

Simkin, S. M.↗

NASA NTRS1995

Investigation of Benchmark $k$ eff Sensitivity and Uncertainty for 239 Pu fission in Specific Energy Ranges

Nuclear data at intermediate energies (from 1 to 100s of keV) are evaluated based on scarce differential data and theory unable to capture physics’ expected structure. There is also a lack of integral data. This is a known deficiency and is challenging to address. Calculated effective multiplication factor, k eff , values for intermediate energy experiments are ~25× further from experiment than for fast energies and are often well outside the experimental uncertainties. The goal of the PARADIGM (PARallel Approach of Differential and InteGral Measurements) project is to significantly re duce the uncertainties of intermediate energy nuclear data for 239 Pu. To this end, PARADIGM simultaneously optimizes experiments at both the Los Alamos Neutron Science Center (LANSCE) and National Criticality Experiments Research Center (NCERC). The combined set of data will inform new intermediate-energy nuclear data. By execution of differential and integral experiments, establishment of new theory, and undertaking nuclear data evaluation in parallel, the timeline to deliver improved nuclear data to users will be reduced significantly that is to three years. For the PARADIGM project, it was decided to optimize an integral experiment for two neutron energy ranges, within the full intermediate energy range. The low energy range goes from 1 to 30 keV, while the higher energy range goes from 30 to 600 keV. This work focuses on nuclear data sensitivities and uncertainties for 239 Pu fission for existing experiments in the International Criticality Safety Benchmark Evaluation Project (ICSBEP). When designing new experiments, it is important to understand what benchmarks currently exist. For a more traditional experiment design (in which a specific application model(s) exists), comparisons would be made between the application model(s) and existing benchmarks. For PARADIGM, there is no specific application model, but instead the specific nuclear data reaction and energy ranges of interest can be explored for existing benchmarks.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRSDate not supplied

Uncertainty Quantification using Deep Ensembles for Decision Making in Cyber-Physical-Human Systems

In this paper and its companion, Differential Equation Approximation Using Gradient-Boosted Quantile Regression, Robison et al., we examine an approach to quantifying model uncertainty with the aim of increasing the trustworthiness of computational models in human-machine interactions. In Differential Equation Approximation Using Gradient-Boosted Quantile Regression, we focus on gradient-boosted decision trees, while in this one, we give more details about deep ensembles. Uncertainty quantification is crucial for building trustworthy autonomous decision-making agents in human-machine teams. There are two types of uncertainties: aleatoric and epistemic. The former is related to the inherent stochasticity (noise) of the process, whereas the latter is associated with the lack of knowledge or representation capability of models, such as neural networks. By lack of knowledge, we mean the model’s inability to accurately predict outputs for all possible inputs. The aleatory uncertainty can be estimated fairly easily with, for example, filters, whereas epistemic uncertainty is challenging to compute. This paper uses deep ensembles to quantify both aleatory and epistemic uncertainty. It can act as an uncertainty-aware surrogate transition model for decision-making frameworks. "Uncertainty-aware" means that the surrogate transition model should make predictions along with confidence in those predictions. In the context of decision-making, the transition models are ordinary differential equations (ODEs). Since ODEs can be simulated to make one-step or multi-step predictions, a good surrogate model for them should perform reasonably well in both modes. In a multi-step approach, the trajectory sampling method TS∞ was used to propagate uncertainty over multiple steps. The cartpole dynamical system was selected to demonstrate the ability of deep ensembles as good surrogate transition models for decision-making frameworks. The deep ensembles modeled the dynamics of cartpole ODEs and made uncertainty-aware predictions in single-step and multi-step transition modes.

CPH systems↗

DOE OSTI2020

A physics-informed operator regression framework for extracting data-driven continuum models

The application of deep learning toward discovery of data-driven models requires careful application of inductive biases to obtain a description of physics which is both accurate and robust. We present here a framework for discovering continuum models from high fidelity molecular simulation data. Our approach applies a neural network parameterization of governing physics in modal space, allowing a characterization of differential operators while providing structure which may be used to impose biases related to symmetry, isotropy, and conservation form. Here, we demonstrate the effectiveness of our framework for a variety of physics, including local and nonlocal diffusion processes and single and multiphase flows. For the flow physics we demonstrate this approach leads to a learned operator that generalizes to system characteristics not included in the training sets, such as variable particle sizes, densities, and concentration.

42 ENGINEERING↗

DOE OSTI2021

Method and apparatus for shape-based energy analysis of solids

A computer-readable medium stores instructions including storing a model of a physical structure and defining a mesh for the model. The mesh includes finite elements, each defined by a respective set of edges. The instructions include, for each finite element, identifying a governing differential equation and a set of complementary functions that exactly satisfy the governing differential equation. The instructions include determining an applied physical stimulus for the physical structure. The instructions include generating an energy optimization model that minimizes a difference between internal energy and external energy of the finite elements in response to the applied physical stimulus. The model includes a matrix of respective scalar multipliers for the complementary functions for each finite element. The instructions include transforming the matrix and calculating a physical parameter of interest. The instructions include, in response to the physical parameter not satisfying a design parameter, updating the model of the physical structure.

Spears, Robert E.↗

DOE OSTI2021

Structure–Capacitance Relationships of Graphene/Ionic Liquid Electrolyte Double Layers

The differential capacitance profile of electrochemical interfaces reflects the physical properties of the double layer. For carbon electrodes and ionic liquid based electrolytes, these capacitance profiles are not fully understood. In this work, we utilize constant voltage molecular dynamics simulations to compute differential capacitance profiles of ionic liquids [BMIm + ][BF 4 - ] and [BMIm + ][TFSI - ] mixed with acetonitrile and 1,2-dichloroethane, at model graphene electrodes. We find that both pure and 10% mole fraction ionic liquid electrolytes exhibit camel shaped capacitance profiles with two peaks on either side of a minimum centered at the potential of zero charge. This profile shape results from electric field induced rearrangement of ion structure within the inner layer closest to the electrode interface. At low potential, the ionic liquid inner layer is concentrated with non-polar tri-fluoromethyl and butyl functional groups of the anions and cations, corresponding to the minimum of the capacitance profiles. With increasing voltage, electrostatic interactions of polar/charged functional groups with the electrode surface compete with these non-polar interactions, leading to ion rearrangement that increases the inner layer charge density and results in higher capacitance. After the ion restructuring is complete, the response saturates and capacitance diminishes. The presence of organic solvent significantly changes the composition of the inner layer. For example, strong non-polar interactions between dichloroethane molecules and the graphene surface substantially block ion/electrode contact at moderate potentials. Overall, our simulations highlight the dynamic nature of the inner region of organic electrolyte double layers, and the sensitive dependence on electrolyte composition and applied voltage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗