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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Landau discriminants

Scattering amplitudes in quantum field theories have intricate analytic properties as functions of the energies and momenta of the scattered particles. In perturbation theory, their singularities are governed by a set of nonlinear polynomial equations, known as Landau equations, for each individual Feynman diagram. The singularity locus of the associated Feynman integral is made precise with the notion of the Landau discriminant, which characterizes when the Landau equations admit a solution. In order to compute this discriminant, we present approaches from classical elimination theory, as well as a numerical algorithm based on homotopy continuation. These methods allow us to compute Landau discriminants of various Feynman diagrams up to 3 loops, which were previously out of reach. For instance, the Landau discriminant of the envelope diagram is a reducible surface of degree 45 in the three-dimensional space of kinematic invariants. We investigate geometric properties of the Landau discriminant, such as irreducibility, dimension and degree. In particular, we find simple examples in which the Landau discriminant has codimension greater than one. Furthermore, we describe a numerical procedure for determining which parts of the Landau discriminant lie in the physical regions. In order to study degenerate limits of Landau equations and bounds on the degree of the Landau discriminant, we introduce Landau polytopes and study their facet structure. Finally, we provide an efficient numerical algorithm for the computation of the number of master integrals based on the connection to algebraic statistics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Causal diamonds, cluster polytopes and scattering amplitudes

The “amplituhedron” for tree-level scattering amplitudes in the bi-adjoint φ 3 theory is given by the ABHY associahedron in kinematic space, which has been generalized to give a realization for all finite-type cluster algebra polytopes, labelled by Dynkin diagrams. In this letter we identify a simple physical origin for these polytopes, associated with an interesting (1 + 1)-dimensional causal structure in kinematic space, along with solutions to the wave equation in this kinematic “spacetime” with a natural positivity property. The notion of time evolution in this kinematic spacetime can be abstracted away to a certain “walk”, associated with any acyclic quiver, remarkably yielding a finite cluster polytope for the case of Dynkin quivers. The A n–3 , B n–1 /C n–1 and D n polytopes are the amplituhedra for n-point tree amplitudes, one-loop tadpole diagrams, and full integrand of one-loop amplitudes. We also introduce a polytope D¯ n , which chops the D n polytope in half along a symmetry plane, capturing one-loop amplitudes in a more efficient way.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Open associahedra and scattering forms

We continue the study of open associahedra associated with bi-color scattering amplitudes initiated in ref. [1]. We focus on the facet geometries of the open associahedra, uncovering many new phenomena such as fiber-product geometries. We then provide novel recursion procedures for calculating the canonical form of open associahedra, generalizing recursion relations for bounded polytopes to unbounded polytopes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Renormalizing two-fermion operators in the SMEFT via supergeometry

We extend the geometric framework of field-space covariance for loop computations, thereby unifying the treatment of scalars, fermions, and gauge bosons in effective field theories. This allows us to derive a manifestly covariant formula for one-loop UV divergences that includes contributions from mixed boson-fermion graphs. The result is expressed in terms of geometric invariants of the field-space supermanifold. As a demonstration of this formula, we compute the renormalization group equations for two-fermion operators at the dimension-eight level in the Standard Model Effective Field Theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Preserving nonlinear constraints in variational flow filtering data assimilation

Data assimilation aims to estimate the states of a dynamical system by optimally combining sparse and noisy observations of the physical system with uncertain forecasts produced by a computational model. The states of many dynamical systems of interest obey nonlinear physical constraints, and the corresponding dynamics is confined to a certain sub-manifold of the state space. Standard data assimilation techniques applied to such systems yield posterior states lying outside the manifold, violating the physical constraints. This work focuses on particle flow filters which use stochastic differential equations to evolve state samples from a prior distribution to samples from an observation-informed posterior distribution. The variational Fokker-Planck (VFP)—a generic particle flow filtering framework—is extended to incorporate non-linear, equality state constraints in the analysis. To this end, two algorithmic approaches that modify the VFP stochastic differential equation are discussed: (i) VFPSTAB, to inexactly preserve constraints with the addition of a stabilizing drift term, and (ii) VFPDAE, to exactly preserve constraints by treating the VFP dynamics as a stochastic differential-algebraic equation (SDAE). Additionally, an implicit-explicit time integrator is developed to evolve the VFPDAE dynamics. The strength of the proposed approach for constraint preservation in data assimilation is demonstrated on three test problems: the double pendulum, Korteweg-de-Vries, and the incompressible Navier-Stokes equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Effects of the CPT-even and Lorentz violation on the Bhabha scattering at finite temperature

In this paper a Lorentz-violating CPT-even non-minimal coupling term is considered. A new interaction term between fermions and photons emerges. In this context, the differential cross-section for Bhabha scattering at finite temperature is calculated. The temperature effects are introduced using the Thermo-Field Dynamics (TFD) formalism. It is shown that the differential cross-section is changed due to both effects, Lorentz violation and finite temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Deuteration Effects on the Physical and Optoelectronic Properties of Donor–Acceptor Conjugated Polymers

The significant differences in scattering cross sections between deuterium and protium are unique to neutron scattering techniques and have been a long-standing area of interest within the neutron scattering community. Researchers have explored selective deuteration to manipulate scattering contrast in soft matter systems, leading to the widespread use of deuterium labeling in materials development. As deuteration changes the atomic mass, it alters physical properties such as molecular volume, polarizability, and polarity, which in turn may affect noncovalent interactions and crystal ordering. Despite previous studies, there remains a limited understanding of how deuteration impacts donor–acceptor (DA) conjugated polymers. To address this, we synthesized deuterated DPP polymers and systematically investigated the effects of side-chain deuteration on their thermal stability, crystal packing, morphology, and optoelectronic properties. We found that deuteration increased the melting and crystallization temperatures of DPP polymers, although it did not significantly alter their morphology, molecular packing, or charge mobility. These properties were assessed by using atomic force microscopy (AFM), X-ray scattering, and thin-film transistor device measurements, respectively, for DPP polymers. Our work shows that deuterium labeling could be a powerful method for controlling scattering length density, enabling neutrons to study the structure and dynamics of conjugated polymers without impacting their electronic performance.

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

Learning Constitutive Relations From Soil Moisture Data via Physically Constrained Neural Networks

Abstract The constitutive relations of the Richardson‐Richards equation encode the macroscopic properties of soil water retention and conductivity. These soil hydraulic functions are commonly represented by models with a handful of parameters. The limited degrees of freedom of such soil hydraulic models constrain our ability to extract soil hydraulic properties from soil moisture data via inverse modeling. We present a new free‐form approach to learning the constitutive relations using physically constrained neural networks. We implemented the inverse modeling framework in a differentiable modeling framework, JAX, to ensure scalability and extensibility. For efficient gradient computations, we implemented implicit differentiation through a nonlinear solver for the Richardson‐Richards equation. We tested the framework against synthetic noisy data and demonstrated its robustness against varying magnitudes of noise and degrees of freedom of the neural networks. We applied the framework to soil moisture data from an upward infiltration experiment and demonstrated that the neural network‐based approach was better fitted to the experimental data than a parametric model and that the framework can learn the constitutive relations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The influence of hydrostatic pressure and annealing conditions on the magnetostructural transitions in MnCoGe

In this work, the phase transitions of stoichiometric MnCoGe alloys were studied by systematically varying the annealing conditions and applying hydrostatic pressure. First-order martensitic structural transitions from the Ni2In-type hexagonal austenite phase to the TiNiSi-type orthorhombic martensite phase spanned a wide temperature window (>200 K) as a result of quenching the samples at temperatures ranging from the solid phase at 700 °C to the liquid phase at 1150 °C. Despite the large variation in their structural transition temperatures, the changes in cell parameters across the structural transitions and the Curie–Weiss temperatures of the martensite/austenite phase were relatively small. For the sample quenched from 800 °C, coupled magnetostructural transitions were observed, and the largest maximum magnetic entropy change was found to be -ΔS max = 33:6 J/kg K for a 7-T field change. The coupled magnetostructural transitions and the corresponding magnetic entropy enhancements were found to also be achievable by applying hydrostatic pressures. Meanwhile, as the quenching temperatures or hydrostatic pressures increased, the first-order martensitic structural transition shifted toward lower temperature until it was ultimately absent, in which case only the crystal structure and magnetic transition of the Ni 2 In-type hexagonal austenite phase were present.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Comparison of private flux region instability in conventional and super-X divertor configurations

Understanding turbulence in the divertor leg of tokamaks is essential to predicting the heat deposition profile on the divertor plate. This in turn is important for evaluating advanced divertor configurations, such as the super-X divertor. Within the divertor region, the private flux region is of interest because it is relatively unaffected by turbulence extending from the outboard midplane, so instabilities in this region could have a particularly pronounced effect on transport. These instabilities are modeled using the Arbitrary Topology Equation Reader (ArbiTER) eigenvalue code. Eigenmodes are examined further by comparing physics models to determine the fundamental mechanisms behind their formation, and quantifying the effect of individual terms. This analysis is conducted on both conventional and super-X divertors to compare these effects. Here, the resulting analysis reveals the presence of a geodesic curvature driven instability that is significantly more pronounced in the super-X configuration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effects of diffusion barriers on reaction wave stability in Co/Al reactive multilayers

Bimetallic, reactive multilayers are uniformly structured materials composed of alternating sputter-deposited layers that may be ignited to produce self-propagating mixing and formation reactions. These nanolaminates are most commonly used as rapid-release heat sources. The specific chemical composition at each metal/metal interface determines the rate of mass transport in a mixing and formation reaction. The inclusion of engineered diffusion barriers at each interface will not only inhibit solid-state mixing but also may impede the self-propagating reactions by introducing instabilities to wavefront morphology. Here, this work examines the effect of adding diffusion barriers on the propagation of reaction waves in Co/Al multilayers. The Co/Al system has been shown to exhibit a reaction propagation instability that is dependent on the bilayer thickness, which allows for the occurrence of unstable modes in otherwise stable designs from the inclusion of diffusion barriers. Based on the known stability criteria in the Co/Al multilayer system, the way in which the inclusion of diffusion barriers changes a multilayer's heat of reaction, thermal conductivity, and material mixing mechanisms can be determined. These factors, in aggregate, lead to changes in the wavefront velocity and stability.

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

Measuring differential flow angle fluctuations in relativistic nuclear collisions

Recently studies of the differential nature of the flow angle fluctuations, known as event plane angular decorrelation, indicated that measurements that assume a common symmetry plane may need to consider the flow angle fluctuations effect. Using the HIJING and AMPT models, it is shown that the flow angle fluctuations measurements, obtained with the two- and three-subevents correlation method, could have significant nonflow effects associated with long- and short-range nonflow correlations. The current study demonstrates the four-subevents cumulant methods' ability to reduce the nonflow effects. Finally, it is further argued that the measurements using the four-subevents correlation method can be used to provide accurate quantification of the differential flow angle fluctuations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Data-Driven Invertible Neural Surrogates of Atmospheric Transmission

We present Data-Driven Invertible Neural Surrogates of Atmospheric transmission, or DINSAT. DINSAT is a novel framework for inferring an atmospheric transmission profile from a spectral scene. This framework leverages a lightweight, physics-based simulator that is automatically tuned -- by virtue of autodifferentiation and differentiable programming -- to construct a surrogate atmospheric profile to model the observed data. The framework has utility in (i) performing atmospheric correction, (ii) recasting spectral data between various modalities (e.g. radiance and reflectance at the surface and at the sensor), and (iii) inferring atmospheric transmission profiles, such as absorbing bands and their relative magnitudes. We demonstrate the utility of these methods by performing a canonical atmospheric correction task for the purposes of further analysis - in this case, target detection within a scene.

Koch, James V.↗
DOE OSTI2023

Interaction of Concentrated (68%) Aqueous Nitric Acid Solutions with Cheesecloth: A TGA-DSC-MS Study

The exposure of cotton cheesecloth to concentrated aqueous solutions of nitric acid reduced the onset of combustion by 30-40°C (from ~350°C to ~310°C) relative to that of both unexposed cheesecloth and cheesecloth that was rinsed immediately following exposure to the same solution. Cheesecloth exposed to concentrated (68%) nitric acid underwent three exothermic pre-combustion reactions which onset at 125-135°C, 160- 170°C and 200-210°C. These reactions give off mostly water vapor, with lesser amounts of NO and carbon dioxide. Exposure to concentrated nitric acid damages the cellulose polymers, oxidizing them, reducing them in length, and increasing their reactivity in response to thermal stress. Cheesecloth which was rinsed after exposure behaved similarly to unexposed cheesecloth, indicating that NC was not formed on exposure to the concentrated aqueous nitric acid solution. The addition of a 239Pu – 241Am spike in an aqueous nitric – hydrofluoric solution to the cheesecloth reduced the onset temperature of the first exothermic reaction to 126°C and lowered the onset temperature of combustion from ~310°C to ~275°C. When dried on to cellulose substrates in the presence of excess nitrate, metal ions may act as oxidizers, or may catalyze the oxidation of the carbonaceous material, lowering ignition temperatures. No evidence for the formation of nitrocellulose following the exposure of cheesecloth to nitric acid was observed during this study

36 MATERIALS SCIENCE↗
DOE OSTI2023

Domain Aware Deep-learning Algorithms Integrated with Scientific-computing Technologies (DADAIST)

This technical report summarized the contribution of the DADAIST project funded by the Data Model Convergence Initiative via the Laboratory Directed Research and Development (LDRD) investments at Pacific Northwest National Laboratory (PNNL). Specifically, we report the development of the NeuroMANCER (Neural Modules with Adaptive Nonlinear Constraints and Efficient Regularizations), a new open-source Scientific Machine Learning library for formulating and solving parametric constrained optimization problems, physics-informed system identification, and parametric optimal control problems. NeuroMANCER is using differentiable programming to combine modern data-driven models and optimization modeling language into a coherent algorithmic and software framework. NeuroMANCER is a Pytorch-based framework and adopts much of its philosophy focused on research and development, rapid prototyping, and streamlined deployment. Strong emphasis is given to extensibility, interoperability with the PyTorch ecosystem, and quick adaptability to custom domain problems. Neuromancer repository contains a comprehensive library of differentiable modules, including custom activation functions, matrix factorizations, deep learning architectures, neural differential equations, differential equation solvers, implicit layers such as iterative solvers, high-level API for symbolic expressions, API for modeling and control of dynamical systems, and extensive set of tutorial code examples in the form of python scripts and jupyter notebooks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Machine-learning based model reduction for partial differential equations

We develop a novel synergistic approach between model reduction and machine learning. The specific goal of this project is to aid in the construction of reduced order models for basis functions that are custom-made to represent the solution of partial differential equations. Partial differential equations (PDEs) are one of the main mathematical tools for describing physical phenomena. However, due to either efficiency or necessity, for many real-world problems, we are interested in constructing reduced order models (ROMs) which focus only on the explicit computation of subsets of the active spatio-temporal scales in the problem, while treating the interaction with the rest of the scales approximately. The task of accurate representation of such interactions (usually called memory terms) constitutes a vast area of research known as model reduction. PI Stinis has significant expertise in the construction of ROMs for complex systems. In addition, in recent work with the project key participant Qadeer, they have utilized machine learning to acquire custom-made basis functions (CBFs) to expand the solutions of PDEs. In the proposed work, we will merge the two concepts by constructing ROMs for subsets of the CBFs needed to represent the solution of a PDE. Specifically, we will use the Mori-Zwanzig model reduction formalism to construct ROMs for subsets of CBFs for nonlinear PDEs of various complexity, as well as investigate the usage of CBFs in the spectral vanishing viscosity method for problems that can form shocks in finite time. The outcome of the research is aimed to be proof-of-concept about a novel synergistic approach between model reduction and machine learning, thus advancing the field of scientific machine learning. Such a capability will benefit the efficient modeling of physical systems appearing in various areas of interest to the DOE.

97 MATHEMATICS AND COMPUTING↗