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

Inclusive and Inelastic Scattering in Neutrino-Nucleus Interactions

Neutrino-nucleus cross section measurements can provide both insights into nuclear physics and important data that can be used to improve model predictions used for neutrino oscillation physics. Two measurements of neutrino-nucleus cross sections were performed using data from the MINERvA experiment, each probing different classes of neutrino interactions. Double and single-differential flux-integrated measurements of inclusive charged current neutrino-nucleus cross sections at a peak neutrino energy of 3.5 GeV are presented as a function of the longitudinal and transverse momentum of the muon produced in the interaction. Additionally, an analysis of charged-current deep inelastic scattering (DIS) of muon neutrinos was performed in carbon, iron, lead and hydrocarbon in a neutrino beam with a peak energy of 6 GeV. Cross sections were measured in each material, and ratios were also taken between the cross sections of each of the materials and hydrocarbon, allowing for examination of nuclear dependencies of the neutrino cross section. Absolute DIS cross sections as a function of neutrino energy and flux integrated differential cross sections as a function of the Bjorken-x scaling variable were both measured.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Adaptive Space-Time Methods for Large Scale Optimal Design

When modeling complex physical systems with advanced dynamics, such as shocks and singularities, many classic methods for solving partial differential equations can return inaccurate or unusable results. One way to resolve these complex dynamics is through r-adaptive refinement methods, in which a fixed number of mesh points are shifted to areas of high interest. The mesh refinement map can be found through the solution of the Monge-Ampére equation, a highly nonlinear partial differential equation. Due to its nonlinearity, the numerical solution of the Monge-Ampére equation is nontrivial and has previously required computationally expensive methods. In this report, we detail our novel optimization-based, multigrid-enabled solver for a low-order finite element approximation of the Monge-Ampére equation. This fast and scalable solver makes r-adaptive meshing more readily available for problems related to large-scale optimal design. Beyond mesh adaptivity, our report discusses additional applications where our fast solver for the Monge-Ampére equation could be easily applied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Partitioned exponential methods for coupled multiphysics systems

Multiphysics problems involving two or more coupled physical phenomena are ubiquitous in science and engineering. This work develops a new partitioned exponential approach for the time integration of multiphysics problems. After a possible semi-discretization in space, the class of problems under consideration is modeled by a system of ordinary differential equations where the right-hand side is a summation of two component functions, each corresponding to a given set of physical processes. The partitioned-exponential methods proposed herein evolve each component of the system via an exponential integrator, and information between partitions is exchanged via coupling terms. Here, the traditional approach to constructing exponential methods, based on the variation-of-constants formula, is not directly applicable to partitioned systems. Rather, our approach to developing new partitioned-exponential families is based on a general-structure additive formulation of the schemes. Two method formulations are considered, one based on a linear-nonlinear splitting of the right hand component functions, and another based on approximate Jacobians. The paper develops classical (non-stiff) order conditions theory for partitioned exponential schemes based on particular families of T-trees and B-series theory. Several practical methods of third order are constructed that extend the Rosenbrock-type and EPIRK families of exponential integrators. Several implementation optimizations specific to the application of these methods to reaction-diffusion systems are also discussed. Numerical experiments reveal that the new partitioned-exponential methods can perform better than traditional unpartitioned exponential methods on some problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Convergence Criteria for Multiphysics Simulations

The behavior of engineered systems is often influenced by multiple physical phenomena, such as mechanical deformation, heat transfer, and chemical species transport and reactions. There are often strong interactions between these phenomena, and there is increasing interest in applying coupled-physics models to improve understanding of physical behavior under complex environmental conditions. Multiple simulation frameworks that facilitate coupled-physics simulations are in widespread use, and these employ a variety of techniques to account for interactions between those physics. Many frameworks solve the physics models independently and transfer results between them. Alternatively, a single monolithic system of equations for every physics model can be formed and solved. Each of these approaches has its benefits and drawbacks, and the optimal approach varies depending on the nature of the problem. The open-source MOOSE framework was developed targeting solution of large-scale multiphysics problems. Although it provides options for all these coupling approaches, its standard approach for multiphysics solutions is to form and solve a single monolithic system of equations containing the unknowns for all physics models. MOOSE provides a streamlined approach for users to define the solution variables, the terms in the partial differential equations pertaining to each variable, and interactions between solution variables. One aspect of the monolithic solution approach that can be problematic, however, is defining appropriate convergence criteria for the nonlinear system. A standard approach is to determine convergence is to simply take a norm of the residual vector corresponding to the full vector of unknowns. However, if the residual vector contains variables for multiple physics models, the magnitudes of those variables can differ significantly, and the variables can converge at significantly different rates from each other. It is important to ensure that the variables for each of the physics are converged, and also ensure that the convergence criteria are not excessively stringent in cases when there is little change in the solution. This talk presents representative multiphysics problems to highlight these issues, and shows strategies for convergence criteria in MOOSE that are robust for multiphysics models under a variety of conditions.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Rapid Inference of Logic Gate Neural Networks for Anomaly Detection in High Energy Physics

The increasing data rates and complexity of detectors at the Large Hadron Collider (LHC) necessitate fast and efficient machine learning models, particularly for rapid selection of what data to store, known as triggering. Building on recent work in differentiable logic gates, we present a public implementation of a Convolutional Differentiable Logic Gate Neural Network (CLGN). We apply this to detecting anomalies at the Level-1 Trigger at CMS using public data from the CICADA project. We demonstrate that the CLGN achieves physics performance on par with or superior to conventional quantized neural networks. We also synthesize an LGN for a Field-Programmable Gate Array (FPGA) and show highly promising FPGA characteristics, notably zero Digital Signal Processor (DSP) resource usage. This work highlights the potential of logic gate networks for high-speed, on-detector inference in High Energy Physics and beyond.

FOS: Physical sciences↗
DOE OSTI2024

One-sweep moment-based semi-implicit-explicit integration for gray thermal radiation transport

Thermal radiation transport (TRT) is a time dependent, high dimensional partial integro-differential equation. In practical applications such as inertial confinement fusion, TRT is coupled to other physics such as hydrodynamics, plasmas, etc., and the timescales one is interested in capturing are often much slower than the radiation timescale. As a result, TRT is treated implicitly, and due to its stiffness and high dimensionality, is often a dominant computational cost in multiphysics simulations. Here we develop a new approach for implicit-explicit (IMEX) integration of gray TRT in the deterministic SN setting, which requires only one sweep per stage, with the simplest first-order method requiring only one sweep per time step. The partitioning of equations is done via a moment-based high-order low-order formulation of TRT, where the streaming operator and first two moments are used to capture the asymptotic stiff regimes of the streaming limit and diffusion limit. Absorption-reemission is treated explicitly, and although stiff, is sufficiently damped by the implicit solve that we achieve stable accurate time integration without incorporating the coupling of the high order and low order equations implicitly. Due to nonlinear coupling of the high-order and low-order equations through temperature-dependent opacities, to facilitate IMEX partitioning and higher-order methods, we use a semi-implicit integration approach amenable to nonlinear partitions. In conclusion, results are demonstrated on thick Marshak and crooked pipe benchmark problems, demonstrating orders of magnitude improvement in accuracy and wallclock compared with the standard first-order implicit integration typically used.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Properties of different LLM-105 preparations

LLM-105 (2,6-diamino-3,5-dinitropyrazine-1-oxide) has been prepared at several different size quantities (~10–70 Kg) and not subjected to further purification, such as recrystallization. Chemical and physical properties were compared—small-scale safety testing (impact, friction, ESD), thermal (chemical reactivity, differential scanning calorimetry (DSC), thermogravimetric analysis (TGA), one-dimensional time-to-explosion (ODTX)), morphological (particle-size distribution (PSD), scanning electron microscopy (SEM), powder x-ray diffraction (PXRD), skeletal density). Each preparation had similar properties to the other preparations except for PSD and PXRD. For PSD, all preparations exhibited a mono-modal distribution with a variation in median values from ~ 40 μm to over 80 μm. All PXRD values indicated the same morphology with minor variations in crystal orientations. Possible polymorphism was also observed in a few cases. SEM images and PXRD indicate all preparations to be diamond-type and X-type crystals with some twinning. Impact, friction, and ESD and thermal sensitivities values are consistent with these type of crystal habits as seen in vast literature studies.

Explosive characterization↗
DOE OSTI2024

cppTPSA/pyTPSA: a C++/Python package for truncated power series algebra

The truncated power series algebra (TPSA), also referred to as differential algebra (DA), is a well-established and widely used method in particle accelerator physics and astronomy. The most straightforward usage of TPSA/DA is to calculate the Taylor expansion of a given function at a specific point up to order ?. In recent years, as the application of TPSA/TA has been extended to other fields, a reusable implementation of TPSA/DA as a modern C++ library or other high level programming language like Python has become desirable. The cppTPSA package implements TPSA/DA in C++11 and provides developers a convenient library with which to build advanced TPSA/DA-based methods. A Python 3 library, pyTPSA, has also been developed based on the C++ lib.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Updates to the n+ 63,65 Cu Angular Distributions [Slides]

The performance of the benchmark suite is very sensitive to changes in the angular distribution data. The quasi-differential measurements performed by Blain et al at RPI provide a valuable constraint. Furthermore, ENDF/B-VIII.0 disagrees with their measurements consistently at 300 keV (the transition between RRR and high energy). The next step is to conservatively adjust the Legendre coefficients near 300 keV, validating performance against both critical benchmarks and the RPI quasi-differential measurements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Pool boiling critical heat flux studies of accident tolerant fuel cladding materials

Pool boiling experiments have been performed at atmospheric pressure on potential accident tolerant fuel (ATF) cladding materials (Zirlo®, Cr coated Zirlo®, FeCrAl coated Zirlo®, and monolithic SiC) to evaluate the Critical Heat Flux (CHF) points. Corrosion resistant coatings of Cr and FeCrAl alloy were deposited on flat Zirlo® samples using cold spray technology. The as-prepared samples after uniform surface polishing were subjected to autoclave tests at 360 °C water and 18.6 MPa for 360 h to simulate prototypic corrosion of the ATF cladding materials in LWR normal operation. Surface characteristics potentially influencing CHF such as surface morphology, roughness, static contact angle, and surface chemistry were characterized using a suite of characterization methods including scanning electron microscopy, 3D optical profilometry, optical contact angle measurements, and x-ray photoelectron spectroscopy (XPS). Thermo-physical properties of the samples such as density, thermal conductivity, and heat capacity were measured by differential scanning calorimetry and laser flash thermal diffusivity measurement. FeCrAl coatings showed CHF values comparable to bare Zirlo® samples, while slightly lower CHF values were observed for Cr coated samples and bulk SiC flats (produced by chemical vapor deposition, CVD). As-deposited Cr coatings showed 67% higher CHF than the as-polished Cr coatings due to its higher surface roughness level. The autoclave testing for the prepared samples generally increased the CHF except for the surface polished Cr coatings. Similar CHF values of bare Zirlo®, FeCrAl coating, and CVD SiC samples were observed after autoclave test but there were negligible change CHF of Cr coatings. XPS studies indicated that the formation or deposition of a few monolayers of hydrophobic carbonaceous species on surface can affect CHF values. The trends in CHF data are discussed in terms of evolution of the surface characteristics of the ATF materials.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Harnessing Quantum Capacitance in 2D Material/Molecular Layer Junctions for Novel Electronic Device Functionality

Two-dimensional (2D) materials promise advances in electronic devices beyond Moore’s scaling law through extended functionality, such as non-monotonic dependence of device parameters on input parameters. However, the robustness and performance of effects like negative differential resistance (NDR) and anti-ambipolar behavior have been limited in scale and robustness by relying on atomic defects and complex heterojunctions. In this paper, we introduce a novel device concept that utilizes the quantum capacitance of junctions between 2D materials and molecular layers. We realized a variable capacitance 2D molecular junction (vc2Dmj) diode through the scalable integration of graphene and single layers of stearic acid. The vc2Dmj exhibits NDR with a substantial peak-to-valley ratio even at room temperature and an active negative resistance region. The origin of this unique behavior was identified through thermoelectric measurements and ab initio calculations to be a hybridization effect between graphene and the molecular layer. The enhancement of device parameters through morphology optimization highlights the potential of our approach toward new functionalities that advance the landscape of future electronics.

2D materials↗
DOE OSTI2023

Structural inference of networked dynamical systems with universal differential equations

Networked dynamical systems are common throughout science in engineering; e.g., biological networks, reaction networks, power systems, and the like. For many such systems, nonlinearity drives populations of identical (or near-identical) units to exhibit a wide range of nontrivial behaviors, such as the emergence of coherent structures (e.g., waves and patterns) or otherwise notable dynamics (e.g., synchrony and chaos). Here, we seek to infer (i) the intrinsic physics of a base unit of a population, (ii) the underlying graphical structure shared between units, and (iii) the coupling physics of a given networked dynamical system given observations of nodal states. These tasks are formulated around the notion of the Universal Differential Equation, whereby unknown dynamical systems can be approximated with neural networks, mathematical terms known a priori (albeit with unknown parameterizations), or combinations of the two. We demonstrate the value of these inference tasks by investigating not only future state predictions but also the inference of system behavior on varied network topologies. The effectiveness and utility of these methods are shown with their application to canonical networked nonlinear coupled oscillators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

When Do Extended Physics-Informed Neural Networks (XPINNs) Improve Generalization?

Physics-informed neural networks (PINNs) have become a popular choice for solving high-dimensional partial differential equations (PDEs) due to their excellent approximation power and generalization ability. Recently, extended PINNs (XPINNs) based on domain decomposition methods have attracted considerable attention due to their effectiveness in modeling multiscale and multiphysics problems and their parallelization. However, theoretical understanding of their convergence and generalization properties remains unexplored. In this study, we take an initial step towards understanding how and when XPINNs outperform PINNs. Specifically, for general multilayer PINNs and XPINNs, we first provide a prior generalization bound via the complexity of the target functions in the PDE problem and a posterior generalization bound via the posterior matrix norms of the networks after optimization. Moreover, based on our bounds, we analyze the conditions under which XPINNs improve generalization. Concretely, our theory shows that the key building block of XPINN, namely, the domain decomposition, introduces a tradeoff for generalization. On the one hand, XPINNs decompose the complex PDE solution into several simple parts, which decreases the complexity needed to learn each part and boosts generalization. On the other hand, decomposition leads to less training data being available in each subdomain, and hence such a model is typically prone to overfitting and may become less generalizable. Empirically, we choose five PDEs to show when XPINNs perform better than, similar to, or worse than PINNs, hence demonstrating and justifying our new theory.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Neutron elastic and inelastic scattering differential cross sections on carbon

Elastic and inelastic neutron scattering angular distributions were measured on natural carbon samples to confirm existing experimental data and evaluations in the fast neutron region and to guide improvements in resonance parameters, where needed. Sixty-four (n,n') differential cross section measurements were performed at 45 incident neutron energies between 0.5 and 8.0 MeV. Experimental angle-integrated elastic scattering cross sections are consistent with ENDF/B-VIII.0 values with the exception of the region from 3.2 to 4.0 MeV where our results are ~3% higher. In the 3.4 to 3.6 MeV region our differential cross sections are slightly lower at forward angles and somewhat higher at backward angles than the ENDF/B-VIII.0 calculations; however, our results are consistent with measured data from other research groups. The first- through fourth-order elastic scattering Legendre coefficients from fits to these experimental data are consistent with ENDF/B-VIII.0 values across this range of incident neutron energies. Inelastic scattering cross sections were measured at 12 incident neutron energies between 5.6 and 7.8 MeV. Angle-integrated cross sections agree well with the 1978 Perey and ENDF/B-VIII.0 values; however, the angular-distribution shapes deviate somewhat from the ENDF/B-VIII.0 calculations. These data will allow refinement of the resonance parameter description in this difficult energy region.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

3D reconstruction and neural rendering for adversarial machine learning

While evasion attacks on computer vision systems have been widely studied, creating attacks that remain effective under significant changes in viewpoint continues to be challenging. Traditional approaches often rely on affine transformations of images, but these approaches degrade at larger perspective shifts and often produce unrealistic or ineffective perturbations. Recent methods use differentiable renderers to improve viewpoint robustness, but they typically depend on manually constructed 3D models. We introduce a semi-automated pipeline that generates physically printable and perspective-invariant adversarial patches using only a small set of 2D images. Our method integrates 3D reconstruction, neural rendering, adversarial patch optimization, and an object detection victim model into a unified workflow. We use 2D Gaussian Splatting for high fidelity mesh reconstruction and FlexPara for surface parameterization that produces texture maps suitable for patch editing. Together, these components form a fully differentiable pipeline in PyTorch3D that links texture modification to model outputs, enabling efficient optimization of patches that remain effective across many viewpoints. The complete process, from image capture to patch printing and physical evaluation, can be completed within a few hours. We demonstrate the effectiveness of the resulting patches through attacks on the YOLOv8 object detection model and discuss remaining challenges and opportunities for improving robustness and scalability.

Singhvi, Vivaan [ORNL] (ORCID:0009000586288221)↗
DOE OSTI2024

Triple Differential Muon Antineutrino Charge Current Inclusive Cross Section Measurement in NOvA

NOvA is a long-baseline accelerator-based neutrino experiment based in the USA. For its physics goals, NOvA uses two functionally identical detectors. The Near Detector is situated at Fermilab, 1 km from the neutrino target, and the Far Detector is located at Ash River, MN, a distance of 810 km from the neutrino source. The ND sees the high intensity of the neutrino beam due to its close proximity to the neutrino target, giving us a unique opportunity to do high-precision neutrino cross-section measurements. In this talk, we present our latest results of the muon antineutrino charge current inclusive cross-section measurement in the NOvA ND. The new measurement is a triple differential cross section in antimuon kinematic phase-space and in the total energy of all observable final state hadrons. We also compare our data results to various neutrino generator predictions, for example, comparisons to GENIE, NuWro, GiBUU, and NEUT neutrino generators are presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Bootstrap bounds on closed Einstein manifolds

A compact Riemannian manifold is associated with geometric data given by the eigenvalues of various Laplacian operators on the manifold and the triple overlap integrals of the corresponding eigenmodes. This geometric data must satisfy certain consistency conditions that follow from associativity and the completeness of eigenmodes. We show that it is possible to obtain nontrivial bounds on the geometric data of closed Einstein manifolds by using semidefinite programming to study these consistency conditions, in analogy to the conformal bootstrap bounds on conformal field theories. These bootstrap bounds translate to constraints on the tree-level masses and cubic couplings of Kaluza-Klein modes in theories with compact extra dimensions. We show that in some cases the bounds are saturated by known manifolds.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Quantum Algorithm for Linear Non-unitary Dynamics with Near-Optimal Dependence on All Parameters

We introduce a family of identities that express general linear non-unitary evolution operators as a linear combination of unitary evolution operators, each solving a Hamiltonian simulation problem. This formulation can exponentially enhance the accuracy of the recently introduced linear combination of Hamiltonian simulation (LCHS) method [An, Liu, and Lin, Physical Review Letters, 2023]. For the first time, this approach enables quantum algorithms to solve linear differential equations with both optimal state preparation cost and near-optimal scaling in matrix queries on all parameters.

Applied Dynamical Systems↗