Engineering Papers⌕ Search

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.

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

At least 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Multiparticle integral and differential correlation functions

This paper formalizes the use of integral and differential cumulants for measurements of multiparticle event-by-event transverse momentum fluctuations, rapidity fluctuations, as well as net-charge fluctuations. This enables the introduction of multiparticle balance functions, defined based on differential correlation functions (factorial cumulants), that suppress two- and three-prong resonance decays effects and enable measurements of underlying long-range correlations obeying quantum number conservation constraints. These multiparticle balance functions satisfy simple sum rules determined by quantum number conservation. It is additionally shown that these multiparticle balance functions arise as an intrinsic component of high-order net-charge cumulants. This implies that the magnitude of these cumulants, measured in a specific experimental acceptance, is strictly constrained by charge conservation and primarily determined by the rapidity and momentum width of these balance functions. The paper also presents techniques to reduce the computation time of differential correlation functions up to order n = 10 based on the methods of moments. Published by the American Physical Society 2024

Physics↗
DOE OSTI2022

Divergence of Velocity Fields in Electrochemical Systems

The passage of current through electrochemical systems results in the development of concentration gradients in the electrolytic phase that can be modeled using concentrated solution theory. Application of this theory requires knowledge of three concentration-dependent transport coefficients, which are often taken to be conductivity, diffusion coefficient, and the cation transference number with respect to the solvent velocity. The governing diffusion equation for molar concentration contains two additional terms - the thermodynamic factor which is related to activity of the electrolytic species and the solvent velocity. The main advance in this paper is the derivation of an expression for the divergence of the solvent velocity. Solving this equation requires knowledge of the partial molar volume of the electrolyte. Analogous expressions are derived for the mass average and molar average velocities. These velocities occur naturally in the diffusion equation if concentration is expressed as weight fraction and mole fraction of the electrolytic species, respectively.

diffusion↗
DOE OSTI2025

DONUT: Physics-aware Machine Learning for Real-time X-ray Nanodiffraction Analysis

SF-25-088 Coherent X-ray scattering techniques are critical for investigating the fundamental structural properties of materials at the nanoscale. While advancements have made these experiments more accessible, real-time analysis remains a significant bottleneck, often hindered by artifacts and computational demands. In scanning X-ray nanodiffraction microscopy, which is widely used to spatially resolve structural heterogeneities, this challenge is compounded by the convolution of the divergent beam with the sample’s local structure. To address this, we introduce DONUT (Diffraction with Optics for Nanobeam by Unsupervised Training), a physics-aware neural network designed for the rapid and automated analysis of nanobeam diffraction data. By incorporating a differentiable geometric diffraction model directly into its architecture, DONUT learns to predict crystal lattice strain and orientation in real-time. Crucially, this is achieved without reliance on labeled datasets or pre-training, overcoming a fundamental limitation for supervised machine learning in X-ray science. We demonstrate experimentally that DONUT accurately extracts all features within the data over 200 times more efficiently than conventional fitting methods.

Zhou, Tao [Argonne National Laboratory (ANL), Argo↗
DOE OSTI2025

Electrically‐Driven Metal‐Insulator Transitions Emerging from Localizing Current Density and Temperature

Negative differential resistance (NDR) is a key electronic response enabling two‐terminal artificial neurons that can be achieved through different physical phenomena, including phase‐homogeneous current density and temperature (electro‐thermal) localizations and spatially‐localized metal‐insulator phase transitions (MITs). These two effects have been observed to occur sequentially in select electrically‐biased transition metal oxides. However, it is unknown why and under what conditions localizing behaviors precede MITs, particularly as a function of device length scale. To this end, the interplay between phase‐homogeneous electro‐thermal localizations and MITs is investigated in a 3D multiphysics simulation of a lateral thin film device, using the material properties of the prototype MIT material VO 2 . These findings demonstrate that the MIT is nucleated through dynamically localizing current density and temperature. A critical device width (≈0.7 µm in this study) is identified, below which both the electrically‐induced electro‐thermal and phase inhomogeneities cease to appear. It is demonstrated that the formation of spatial inhomogeneities directly relates to device dimensions, and demonstrate the decoupling of NDR from the MIT through device scaling relationships. These results provide insight into the material phenomena underlying the material's electrical responses, clarifying conditions under which spatial inhomogeneities form in electrically‐biased MIT materials.

artificial neuron↗
DOE OSTI2024

Cosmic neutrino decoupling and its observable imprints: insights from entropic-dual transport

Abstract Very different processes characterize the decoupling of neutrinos to form the cosmic neutrino background (CνB) and the much later decoupling of photons from thermal equilibrium to form the cosmic microwave background (CMB). The CνB emerges from the fuzzy, energy-dependent neutrinosphere and encodes the physics operating in the early universe in the temperature rangeT∼ 10 MeV toT∼ 10 keV. This is the epoch where beyond Standard Model (BSM) physics, especially in the neutrino sector, may be influential in setting the light element abundances, the necessarily distorted fossil neutrino energy spectra, and other light particle energy density contributions. Here we use techniques honed in extensive CMB studies to analyze the CνB as calculated in detailed neutrino energy transport and nuclear reaction simulations of the protracted weak decoupling and primordial nucleosynthesis epochs. Our moment method, relative entropy, and differential visibility approach can leverage future high precision CMB and light element primordial abundance measurements to provide new insights into the CνB and any BSM physics it encodes. We demonstrate that the evolution of the energy spectrum of the CνB throughout the weak decoupling epoch is accurately captured in the Standard Model by only three parameters per species, a non-trivial conclusion given the deviation from thermal equilibrium and the impact of the decrease of electron-positron pairs. Furthermore, we can interpret each of the three parameters as physical characteristics of a non-equilibrium system. Though the treatment presented here makes some simplifying assumptions including ignoring neutrino flavor oscillations, the success of our compact description within the Standard Model motivates its use also in BSM scenarios. We further demonstrate how observations of primordial light element abundances can be used to place constraints on the CνB energy spectrum, deriving response functions that can be applied for general deviations from a thermal spectrum. Combined with the description of those deviations that we develop here, our methods provide a convenient and powerful framework to constrain the impact of BSM physics on the CνB.

Astronomy & Astrophysics↗
DOE OSTI2024

Stellarator equilibrium axis-expansion to all orders in distance from the axis for arbitrary plasma beta

A systematic theory of the asymptotic expansion of the magnetohydrostatics (MHS) equilibrium in the distance from the magnetic axis is developed to include arbitrary smooth currents near the magnetic axis. Compared with the vacuum and the force-free system, an additional magnetic differential equation must be solved to obtain the pressure-driven currents. It is shown that there exist variables in which the rest of the MHS system closely mimics the vacuum system. Thus, a unified treatment of MHS fields is possible. The mathematical structure of the near-axis expansions to arbitrary order is examined carefully to show that the double-periodicity of physical quantities in a toroidal domain can be satisfied order by order. The essential role played by the leading-order Birkhoff–Gustavson normal form in solving the magnetic differential equations is highlighted. Several explicit examples of vacuum, force-free and MHS equilibrium in different geometries are presented.

Physics↗
DOE OSTI2021

An explicitly solvable energy-conserving algorithm for pitch-angle scattering in magnetized plasmas

In this work, we develop an Explicitly Solvable Energy-Conserving (ESEC) algorithm for the Stochastic Differential Equation (SDE) describing the pitch-angle scattering process in magnetized plasmas. The Cayley transform is used to calculate both the deterministic gyromotion and stochastic scattering, affording the algorithm to be explicitly solvable and exactly energy conserving. An unusual property of the SDE for pitch-angle scattering is that its coefficients diverge at the zero velocity and do not satisfy the global Lipschitz condition. Consequently, when standard numerical methods, such as the Euler-Maruyama (EM), are applied, numerical convergence is difficult to establish. For the proposed ESEC algorithm, its energy-preserving property enables us to overcome this obstacle. We rigorously prove that the ESEC algorithm is order 1/2 strongly convergent. This result is confirmed by detailed numerical studies. For the case of pitch-angle scattering in a magnetized plasma with a constant magnetic field, the numerical solution is benchmarked against the analytical solution, and excellent agreements are found.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Physics informed deep neural network embedded in a chemical transport model for the Amazon rainforest

Secondary organic aerosols (SOA) are fine particles in the atmosphere, which interact with clouds, radiation and affect the Earth’s energy budget. SOA formation involves chemistry in gas phase, aqueous aerosols, and clouds. Simulating these chemical processes involve solving a stiff set of differential equations, which are computationally expensive steps for three-dimensional chemical transport models. Deep neural networks (DNNs) are universal function approximators that could be used to represent the complex nonlinear changes in aerosol physical and chemical processes; however, key challenges such as generalizability to extended time periods, preservation of mass balance, simulating sparse model outputs, and maintaining physical constraints have limited their use in atmospheric chemistry. Here, we develop an approach of using a physics-informed DNN that overcomes previous such challenges and demonstrates its applicability for the chemical formation processes of isoprene epoxydiol SOA (IEPOX-SOA) over the Amazon rainforest. The DNN is trained with data generated by simulating IEPOX-SOA over the entire atmospheric column, using the Weather Research and Forecasting Model coupled with Chemistry (WRF-Chem). The trained DNN is then embedded within WRF-Chem to replace the computationally expensive default solver of IEPOX-SOA formation. The trained DNN predictions generalizes well with the default model simulation of the IEPOX-SOA mass concentrations and its size distribution (20 size bins) over several days of simulations in both dry and wet seasons. The embedded DNN reduces the computational expense of WRF-Chem by a factor of 2. Our approach shows promise in terms of application to other computationally expensive chemistry solvers in climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Structure and dynamics of magneto-inertial, differentially rotating laboratory plasmas

We present a detailed characterization of the structure and evolution of differentially rotating plasmas driven on the MAGPIE pulsed-power generator (1.4 MA peak current, 240 ns rise time). The experiments were designed to simulate physics relevant to accretion discs and jets on laboratory scales. A cylindrical aluminium wire array Z pinch enclosed by return posts with an overall azimuthal off-set angle was driven to produce ablation plasma flows that propagate inwards in a slightly off-radial trajectory, injecting mass, angular momentum and confining ram pressure to a rotating plasma column on the axis. However, the plasma is free to expand axially, forming a collimated, differentially rotating axial jet that propagates at ${\approx }100\,{\rm km}\,{\rm s}^{-1}$ . The density profile of the jet corresponds to a dense shell surrounding a low-density core, which is consistent with the centrifugal barrier effect being sustained along the jet's propagation. We show analytically that, as the rotating plasma accretes mass, conservation of mass and momentum implies plasma radial growth scaling as $r \propto t^{1/3}$ . As the characteristic moment of inertia increases, the rotation velocity is predicted to decrease and settle on a characteristic value ${\approx }20\,{\rm km}\,{\rm s}^{-1}$ . We find that both predictions are in agreement with Thomson scattering and optical self-emission imaging measurements.

Valenzuela-Villaseca, V. (ORCID:0000000167721441)↗
DOE OSTI2023

Redox and mineral controls on Fe and Ti isotopic fractionations during calc-alkaline magmatic differentiation

Titanium and Fe isotopic compositions of lavas from a calc-alkaline differentiation suite and corresponding mineral separates from the Rindjani Volcano, Indonesia show that Fe and Ti isotopic fractionations between minerals and melts are lower than those recorded in other suites at all stages of differentiation. Here, the limited isotopic fractionation for Ti is likely due to low-Ti magnetite and clinopyroxene being the dominant carriers of Ti in Rindjani lavas, as these minerals are thought to have limited equilibrium Ti isotopic fractionation relative to silicate magmas. Other magmatic differentiation suites controlled by removal of Ti-rich magnetite and characterized by a lesser role of clinopyroxene have larger Ti isotopic fractionations. This effect is an indirect consequence of the elevated Fe 3+ /Fe 2+ ratio of calc-alkaline magmas such as Rindjani, which promotes Fe 3+ incorporation into magnetite at the expense of Fe 2+ -Ti 4+ pairs, such that increased oxygen fugacity will subdue Ti isotopic fractionation in global magmatic series. Similarly, we find negligible Fe isotopic fractionation in Rindjani bulk rocks and mineral separates, unlike previous studies. This is also likely due to the oxidized nature of the Rindjani differentiation suite, which leads to similar Fe 3+ /Fe 2+ ratios in melt and minerals and decreases overall mineral-melt Fe fractionation factors. Paired Ti and Fe isotopic analyses may therefore represent a powerful tool to assess oxygen fugacity during differentiation, independent from Fe 3+ determinations of erupted samples.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

How to identify different new neutrino oscillation physics scenarios at DUNE

Next generation neutrino oscillation experiments are expected to measure the remaining oscillation parameters with very good precision. They will have unprecedented capabilities to search for new physics that modify oscillations. DUNE, with its broad band beam, good particle identification, and relatively high energies will provide an excellent environment to search for new physics. If deviations from the standard three-flavor oscillation picture are seen however, it is crucial to know which new physics scenario is found so that it can be verified elsewhere and theoretically understood. We investigate several benchmark new physics scenarios by looking at existing long-baseline accelerator neutrino data from NOvA and T2K and determine at what sensitivity DUNE can differentiate among them. We consider sterile neutrinos and both vector and scalar non-standard neutrino interactions, all with new complex phases, the latter of which could conceivably provide absolute neutrino mass scale information. We find that, in many interesting cases, DUNE will have good model discrimination. We also perform a new fit to NOvA and T2K data with scalar NSI.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Electron emission from bromouracil and uracil induced by protons and radiosensitization

Absolute double differential cross sections (DDCS) of electrons emitted from uracil and 5-bromouracil (BrU) in collisions with protons of energy 200 keV have been measured for various forward and backward emission angles over wide range of electron energies. The measured DDCS are compared with the continuum distorted wave-eikonal initial state (CDW-EIS) calculations. The optimized structure of the BrU was estimated along with the population analysis of all the occupied orbitals using a self-consistent field density. A comparison between the measured DDCS data for the two molecules show that the cross section of low energy electrons emitted from BrU is substantially larger than that for uracil. The BrU-to-uracil DDCS ratios obtained from the present measurements indicate an enhancement of the electron emission by a factor which is as large as 2.0 to 2.5. These electrons being the major agent for damaging the DNA/RNA of the malignant tissues, the present results are expected to provide an important input for the radiosensitization effect in hadron therapy. It is noteworthy to mention that the CDW-EIS calculations for Coulomb ionization cannot predict such enhancement. A large angular asymmetry is observed for uracil with a broad structure, which is absent in case of BrU.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Reliable extrapolation of deep neural operators informed by physics or sparse observations

Deep neural operators can learn nonlinear mappings between infinite-dimensional function spaces via deep neural networks. As promising surrogate solvers of partial differential equations (PDEs) for real-time prediction, deep neural operators such as deep operator networks (DeepONets) provide a new simulation paradigm in science and engineering. Pure data-driven neural operators and deep learning models, in general, are usually limited to interpolation scenarios, where new predictions utilize inputs within the support of the training set. However, in the inference stage of real-world applications, the input may lie outside the support, i.e., extrapolation is required, which may result to large errors and unavoidable failure of deep learning models. Here, we address this challenge of extrapolation for deep neural operators. First, we systematically investigate the extrapolation behavior of DeepONets by quantifying the extrapolation complexity, via the 2-Wasserstein distance between two function spaces and propose a new strategy of bias–variance trade-off for extrapolation with respect to model capacity. Subsequently, we develop a complete workflow, including extrapolation determination, and we propose five reliable learning methods that guarantee a safe prediction under extrapolation by requiring additional information—the governing PDEs of the system or sparse new observations. The proposed methods are based on either fine-tuning a pre-trained DeepONet or multifidelity learning. We demonstrate the effectiveness of the proposed framework for various types of parametric PDEs. Furthermore, our systematic comparisons provide practical guidelines for selecting a proper extrapolation method depending on the available information, desired accuracy, and required inference speed.

42 ENGINEERING↗
DOE OSTI2022

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

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

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Understanding the Role of Hydroxyl Functionalization in Linear Poly(Ethylenimine) for Oxidation‐Resistant Direct Air Capture of CO 2

Aminopolymer-based adsorbents are a prominent class of materials being used for direct air capture of CO 2 at the industrial scale. However, improving their working lifetime, specifically by increasing their resilience to oxidative degradation, remains an ongoing challenge. Toward this end, functionalization of aminopolymers with non-amine functionalities such as hydroxyls has emerged in recent years as a promising strategy toward improving adsorbent lifetime. Although there is a growing body of work demonstrating the effectiveness of this approach and investigating the origin of this improved stability, studies to date have primarily focused on branched aminopolymer systems such as branched poly(ethylenimine). In this work, hydroxyl-functionalized linear poly(ethylenimine) is used to continue to probe the underlying protective mechanism of this strategy. A combination of thermogravimetric analysis, NMR relaxometry, differential scanning calorimetry, and computational simulations is used to better understand the relationship between the extent of chemical functionalization, physical properties, and adsorbent performance.

amine-based adsorbent↗
DOE OSTI2021

Identifiability and predictability of integer- and fractional-order epidemiological models using physics-informed neural networks

Here we analyze a plurality of epidemiological models through the lens of physics-informed neural networks (PINNs) that enable us to identify time-dependent parameters and data-driven fractional differential operators. In particular, we consider several variations of the classical susceptible-infectious-removed (SIR) model by introducing more compartments and fractional-order and time-delay models. We report the results for the spread of COVID-19 in New York City, Rhode Island and Michigan states and Italy, by simultaneously inferring the unknown parameters and the unobserved dynamics. For integer-order and time-delay models, we fit the available data by identifying time-dependent parameters, which are represented by neural networks. In contrast, for fractional differential models, we fit the data by determining different time-dependent derivative orders for each compartment, which we represent by neural networks. We investigate the structural and practical identifiability of these unknown functions for different datasets, and quantify the uncertainty associated with neural networks and with control measures in forecasting the pandemic.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Inverse rendering of fusion plasmas: Inferring plasma composition from imaging systems

In this work, we develop a differentiable rendering pipeline for visualising plasma emission within tokamaks, and estimating the gradients of the emission and estimating other physical quantities. Unlike prior work, we are able to leverage arbitrary representations of plasma quantities and easily incorporate them into a non-linear optimisation framework. The efficiency of our method enables not only estimation of a physically plausible image of plasma, but also recovery of the neutral Deuterium distribution from imaging and midplane measurements alone. We demonstrate our method with three different levels of complexity showing first that a poloidal neutrals density distribution can be recovered from imaging alone, second that the distributions of neutral Deuterium, electron density and electron temperature can be recovered jointly, and finally, that this can be done in the presence of realistic imaging systems that incorporate sensor cropping and quantisation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Development of Manganese and Iron Mixed Metal Oxides for Thermochemical Energy Storage

Energy storage plays an essential role for sustainable and renewable energy use in an effort to reduce carbon emissions. Thermochemical energy storage (TCES) is a favorable alternative to fossil fuel energy systems and a promising solution to the intermittency problem of renewable energy sources (such as solar energy). TCES stores energy in the form of chemical bonds using reversible redox reactions (such as manganese metal oxides). This study developed manganese oxide (Mn2O3) and iron oxide (Fe2O3) with an inert dispersant of silicon dioxide (SiO2) (Mn-Fe-Si) for potential large scale energy storage application. These metal oxides are nontoxic, cheap, and relatively abundant. The materials were prepared with varying metal oxide content using physical mixing in a ball mill, pressure pelletization and calcination at different temperatures. Redox cycle tests were conducted using differential scanning calorimetry combined with thermogravimetric analysis (TGA-DSC). The Mn-Fe-Si materials had higher cyclability and energy release compared to the pure metal oxides and manganese oxide with SiO2. The calcination temperature had a significant effect on releasable energy due to the interactions between the oxides. The Mn-Fe-Si mixture calcined at high temperature (1000oC) showed the highest average energy storage density.

Wilkinson, Olivia↗