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 1,405 records · Page 78

Coupled Gravity and Elevation Measurement of Ice Sheet Mass Change

During June 2003, we measured surface gravity at six locations about a glaciological measurement site located on the South-central Greenland Ice. We operated a GPS unit for 90 minutes at each site -the unit was operated simultaneously with a base station unit in Sondrestrom Fjord so as to enable differential, post-processing of the data. We installed an aluminum, accumulation-rate-pole at each site. The base section of the pole also served as the mount for the GPS antenna. Two gravimeters were used simultaneously at each site. Measurements were repeated at each site with at time lapse of at least 50 minutes. We measured snow physical properties in two shallow pits The same measurement sites were occupied in 1981 and all were part of a hexagonal network of geodetic and glaciological measurements established by The Ohio State University in 1980. Additional gravity observations were acquired at three of the sites in 1993 and 1995. Gravity data were collected in conjunction with Doppler satellite measurements of position and elevation in 1981 and global positioning system measurements subsequently. The use of satellite navigation techniques permitted reoccupation of the same sites in each year to within a few 10 s of meters or better. After detrending the gravity data, making adjustments for tides and removing the residual effects of local spatial gradients in gravity, we observe an average secular decrease in gravity of about 0.01 milligal/year, but with tenths of milligal variations about the mean trend. The trend is consistent with a nearly linear increase in surface elevation of between 7 to 10 c d y r (depending on location) as measured by repeated airborne laser altimeter, surface Doppler satellite and GPS elevation measurements. Differences between the residual gravity anomalies after free air correction may be attributable to local mass changes. This project is a collaboration between the Byrd Polar Research Center of the Ohio State University and the Arctic Technology Center of the Danish Technical University.

Jezek, K. C.↗

Digital Twin Framework for PIP-II Linac: AI-Driven Multi-Scale Modeling from Ion Source to 800 MeV

The PIP-II linac will enable >1.2 MW beam power for DUNE, requiring unprecedented operational reliability across its warm front-end (RFQ, MEBT) and five distinct SRF sections operating at 162.5/325/650 MHz. We present a comprehensive digital twin framework uniquely combining a fully differentiable fast beam transport code with neural network surrogates trained on high-fidelity PIC simulations, capturing space charge and nonlinear dynamics beyond traditional envelope codes while achieving 10⁴× speedup at <1% accuracy. End-to-end differentiability enables gradient-based optimization across 500+ parameters simultaneously—previously impossible with conventional tools—while the model incorporates static/dynamic errors and serves as a virtual commissioning platform for diverse hardware integration. The framework facilitates reinforcement learning for pulsed/CW mode transitions, predictive maintenance through anomaly detection, and autonomous tuning algorithm development with real-time execution capability. Validation against physics simulations shows excellent agreement for the front-end, with initial results demonstrating potential for 30% commissioning time reduction and proactive fault mitigation, providing a scalable blueprint for operating next-generation high-intensity accelerators.

Pathak, Abhishek [Fermilab] (ORCID:000000021704208↗

Constraining the destruction rate of K 40 in stellar nucleosynthesis through the study of the Ar 40 ( p , n ) K 40 reaction

Background: We present that K 40 plays a significant role in the radiogenic heating of Earth-like exoplanets, which can affect the development of a habitable environment on their surfaces. The initial amount of K 40 in the interior of these planets depends on the composition of the interstellar clouds from which they formed. Within this context, nuclear reactions that regulate the production of K 40 during stellar evolution can play a critical role. Purpose: In this study, we constrain for the first time the astrophysical reaction rate of K 40 ( n , p ) Ar 40 , which is responsible for the destruction of K 40 during stellar nucleosynthesis. We provide to the nuclear physics community high-resolution data on the cross section and angular distribution of the Ar 40 ( p , n ) K 40 reaction. These are important to various applications involving Ar 40 . The associated reaction rate of the Ar 40 ( p , n ) K 40 process addresses a reaction rate gap in the Joint Institute for Nuclear Astrophysics REACLIB database in the region of intermediate-mass isotopes. Methods: We performed differential cross-section measurements on the Ar 40 ( p , n ) K 40 reaction, for six energies in the center-of-mass system between 3.2 and 4.0 MeV and various angles between 0 ° and 135 ° . The experiment took place at the Edwards Accelerator Laboratory at Ohio University using the beam swinger target location and a standard neutron time-of-flight technique. We extracted total and partial cross sections by integrating the double differential cross sections we measured. Results: The total and partial cross sections varied with energy due to the contribution from isobaric analog states and Ericson type fluctuations. The energy-averaged neutron angular distributions were symmetrical relative to 90 ° . Based on the experimental data, local transmission coefficients were extracted and were used to calculate the astrophysical reaction rates of Ar 40 ( p , n ) K 40 and K 40 ( n , p ) Ar 40 reactions. The new rates were found to vary significantly from the theoretical rates in the REACLIB library. We implemented the new rates in network calculations to study nucleosynthesis via the slow neutron capture process, and we found that the produced abundance of K 40 is reduced by up to 10% compared to calculations with the library rates. At the same time, the above result removes a significant portion of the previous theoretical uncertainty on the K 40 yields from stellar evolution calculations. Conclusions: Our results support a destruction rate of K 40 in massive stars via the K 40 ( n , p ) Ar 40 reaction that is larger compared to previous estimates. The rate of K 40 destruction via the K 40 ( n , p ) Ar 40 reaction now has a dramatically reduced uncertainty based on our measurement. Lastly, this result directly affects the predicted stellar yields of K 40 from nucleosynthesis, which is a critical input parameter for the galactic chemical evolution models that are currently employed for the study of significant properties of exoplanets.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Hierarchical Model Reduction Driven by Machine Learning for Parametric Advection-Diffusion-Reaction Problems in the Presence of Noisy Data

Abstract We propose a new approach to generate a reliable reduced model for a parametric elliptic problem, in the presence of noisy data. The reference model reduction procedure is the directional HiPOD method, which combines Hierarchical Model reduction with a standard Proper Orthogonal Decomposition, according to an offline/online paradigm. In this paper we show that directional HiPOD looses in terms of accuracy when problem data are affected by noise. This is due to the interpolation driving the online phase, since it replicates, by definition, the noise trend. To overcome this limit, we replace interpolation with Machine Learning fitting models which better discriminate relevant physical features in the data from irrelevant unstructured noise. The numerical assessment, although preliminary, confirms the potentialities of the new approach.

97 MATHEMATICS AND COMPUTING↗

Revealing the effect of local stresses on twin growth mechanisms in titanium using synchrotron X-ray diffraction

Deformation twinning has a significant impact on the evolution of microstructure and mechanical response in hexagonal close packed (hcp) metals. Understanding the physical mechanisms associated with twin nucleation and growth processes is important to enable the broad use of hcp metals. While both nucleation and growth are conditioned by the local stress state, few studies have quantified it in bulk deformed samples. In this study, we reveal the effects of local stresses on twin thickening using in-situ synchrotron experiments with differential aperture X-ray microscopy. High purity Ti is deformed under four-point bending to activate {1012} tensile twins. Further, 3D stress fields with a spatial resolution of 0.5μm are mapped in the vicinity of low and high macroscopic Schmid factor (MSF) twins growing inside a grain. Reconciling this experimental analysis with recent twin growth models reveals that the growth of low MSF twins is limited by the nucleation of twin growth defects in the high-purity Ti sample. While growth-inducing defect nucleation may occur in many places on the interface of the high resolved shear stress (RSS)/MSF twin, the nucleation rate of twinning disconnections for the low MSF twin is high only at the end of the twin near a stress concentration. Once nucleated, these defects easily propagate into lower RSS regions of the grain resulting in twin growth. This work highlights the importance of stress concentrations not only for twin (embryo) nucleation, but also to the growth process of twins, particularly in low MSF twins.

36 MATERIALS SCIENCE↗

Effect of Added Gadolinium Oxide on the Thermal Air Oxidation of Uranium Dioxide

To develop a more reliable and stable UO2-based nuclear fuel, the Pacific Northwest National Laboratory (PNNL) investigated modifying fuel with several soluble lanthanides and zirconium. This article provides the results of these studies investigating gadolinium doping at levels up to 10 mass%. The authors characterized and compared commercially- and PNNL-prepared gadolinium-doped UO2 to determine the oxygen-to-metal ratio, elemental distribution, chemical composition, physical appearance, lattice parameters, and grain structure using atomic force microscopy, scanning electron microscopy coupled with energy dispersive spectroscopy, and X-ray diffractometry. After confirming PNNL-prepared UO2 and (Gd,U)O2 were similar to commercially prepared UO2 and (Gd,U)O2, we measured the thermal behavior of these gadolinium-doped UO2 materials to air oxidation using differential scanning calorimetry and thermogravimetric analysis. Addition of gadolinium stabilized the first oxidation product U4¬O9/U3O7 and slowed the subsequent oxidation to U3O8. Comparison of our measured two-step oxidation of UO2 to U3O8 at 325°C to common gas/solid kinetic reaction models found that each oxidation step is best described as a convolution of kinetic behaviors; the gadolinium insertion into the UO2 lattice enforces a significant alteration in mechanism and oxidation rate. Noticeable changes in mechanism become apparent between 1% and 3% gadolinium content.

Uranium Dioxide, Gadolinium Oxide, Thermal Air Oxi↗

Nematic fluctuations in an orbital selective superconductor Fe 1+ y Te 1– x Se x

Fe 1+y Te 1–x Se x is characterized by its complex magnetic phase diagram and highly orbital-dependent band renormalization. Despite this, the behavior of nematicity and nematic fluctuations, especially for high tellurium concentrations, remains largely unknown. Here we present a study of both B 1 g and B 2 g nematic fluctuations in Fe 1+y Te 1–x Se x (0 ≤ x ≤ 0.53) using the technique of elastoresistivity measurement. We discovered that the nematic fluctuations in two symmetry channels are closely linked to the corresponding spin fluctuations, confirming the intertwined nature of these two degrees of freedom. We also revealed an unusual temperature dependence of the nematic susceptibility, which we attributed to a loss of coherence of the d xy orbital. Our results highlight the importance of orbital differentiation on the nematic properties of iron-based materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Probing putative orbital differentiation effects via Eu 2+ spin dynamics in Sr 1-x Eu x Fe 2 As 2

Here, in this work, we report x-ray powder diffraction, elemental analysis, electrical resistivity, magnetic susceptibility, specific heat, and electron spin resonance (ESR) in single crystals of Sr 1-x Eu x Fe 2 As 2 . We observed a breakdown of the previously reported scaling between the Eu 2+ Korringa relaxation rate obtained from ESR and the spin density wave temperature evolution for Sr-rich samples. This result suggests a distinct evolution of the orbital differentiation of the Fe 3d bands along the Sr-based series when compared to the Ba counterpart. We argue that this difference is related to a larger splitting between the structural (tetragonal-to-orthorhombic) and the Fe-driven spin density wave transitions induced by Eu doping in this series. In fact, our results indicate that the two transitions follow an opposite x-Eu dependence for Sr-concentrated samples. Our work shows that Sr 1-x Eu x Fe 2 As 2 series and the comparison with their Ba-based counterparts are exciting platforms to be explored for understanding the interplay among orbital differentiation, magnetism, and structural distortions in the iron pnictides

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Mathematical Foundations for Nonlocal Interface Problems: Multiscale Simulations of Heterogeneous Materials (Final LDRD Report)

Nonlocal models provide a much-needed predictive capability for important Sandia mission applications, ranging from fracture mechanics for nuclear components to subsurface flow for nuclear waste disposal, where traditional partial differential equations (PDEs) models fail to capture effects due to long-range forces at the microscale and mesoscale. However, utilization of this capability is seriously compromised by the lack of a rigorous nonlocal interface theory, required for both application and efficient solution of nonlocal models. To unlock the full potential of nonlocal modeling we developed a mathematically rigorous and physically consistent interface theory and demonstrate its scope in mission-relevant exemplar problems.

97 MATHEMATICS AND COMPUTING↗

Nuclear Data: What is the big deal? [Slides]

This presentation begins with a look at the practical application: differential and integral data. It also includes a nuclear data evaluation: cross section formalisms, data processing and cross section library generation, computer tools, and criticality safety assessment: nuclear data. The presentation ends with some concluding remarks.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Examining Graphite Degradation in Molten Salt Environments: A Chemical, Physical, and Material Analysis

Molten-salt reactors (MSRs) are Generation IV nuclear reactors that use liquid salt as a coolant and/or fuel. In several MSR designs, graphite serves as a moderator and/or reflector. However, due to limited experimental data and operational experience, our understanding of graphite behavior in molten salt environments remains incomplete. This report aims to identify the degradation mechanisms of nuclear graphite in MSRs, detail the mechanisms of each factor, and provide an initial assessment of their impact on the structural integrity of graphite components. This assessment is based on an extensive literature review and insights from subject matter experts. Furthermore, given the limited data, a modeling strategy using existing Grizzly software is proposed for a more thorough analysis where appropriate. Additionally, it presents mitigation strategies where applicable. The report covers physical degradation mechanisms such as infiltration, erosion, and abrasion, as well as chemical degradation mechanisms including fluorination, intercalation, corrosion, and oxidation. Molten salt can infiltrate the porous structure of graphite, leading to several detrimental effects. Entrapment of fissile products within the graphite pores can cause radiation damage and could pose challenges in the handling and disposal of contaminated components. The differential thermal expansion between the infiltrated salt and graphite, along with internal stress from pressurized molten salt and volumetric heating, can compromise the structural integrity of graphite. To mitigate these effects, employing ultra-fine graphite grades and applying sealants and coatings are effective strategies. A computational model based on coupled solid mechanics and heat transfer phenomena could be used to predict the internal stresses using Grizzly software. In pebble-bed MSRs, graphite fuel pebbles can cause abrasion against reactor components due to friction and wear. The severity of wear is influenced by various factors such as temperature, environment, and the presence of lubricants. Tribological studies reveal that higher temperatures and molten salt environments, such as FLiBe, significantly reduce wear rates compared to dry conditions. Additionally, the chemical composition of the salt can further optimize graphite's tribological performance. Long-term wear effects can be modeled by incorporating surface defects into the geometry and predict stresses under thermal and radiation effects using Grizzly software. Chemical degradation of graphite in a molten salt environment can occur through fluorination and intercalation. Fluorination can occur via replacement of hydrogen or oxygen atoms, or at the active sites, but does not cause structural degradation. Intercalation, on the other hand, can lead to exfoliation, where layers of graphite separate and peel away, damaging the graphite. Protective coatings can enhance graphite's resistance to intercalation. Graphite generally exhibits good chemical stability in molten salt environments, though it can corrode under specific conditions, particularly in the presence of impurities or oxidants. Studies have shown that protective coatings, such as plasma-sprayed partially stabilized zirconia (PSZ), can effectively prevent such degradation. Corrosion behavior varies significantly with different graphite grades and coating applications, underscoring the need for detailed studies on uncoated and coated graphite to understand and mitigate corrosion mechanisms in MSRs. Research indicates that the presence of oxidants and impurities can accelerate graphite degradation in molten salts, making it essential to explore acceptable impurity limits. Oxidation is another critical degradation mechanism, leading to weight loss and structural damage due to the formation of CO and CO 2 from the reaction of carbon atoms with oxygen. This process creates new porosity and compromises graphite's integrity. While extensive research on graphite oxidation has been conducted for gas-cooled reactors, studies specific to MSRs are limited. Findings from the coal industry suggest that molten alkali metal salts can significantly accelerate graphite oxidation, a hypothesis worth exploring for fluoride salts in MSRs. Understanding oxidation behavior in MSRs is vital for developing protective measures. The analysis of post-irradiated graphite from the MSRE experiment demonstrated exceptional chemical compatibility with molten fluoride salt, suggesting that the extent of chemical attack on graphite largely depends on the salt's infiltration capability. Therefore, the use of ultra-fine grade graphite could help mitigate chemical degradation effects. Existing oxidation modeling capabilities in Grizzly, which use reaction-diffusion equations to model graphite-air interactions, could be adapted to simulate the chemical degradation effects of graphite in molten salt environments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Calculation of the motion of an airplane under the influence of irregular disturbances

The application of mathematical advances made in electricity and other branches to problems of airplane dynamics is demonstrated. The Heaviside-Bromwich methods of solution of linear differential equations are described and it is shown how these methods avoid the consideration of boundary conditions and of particular or complementary integrals. It is pointed out that if the solution of the differential equation is obtained for the case of a unit disturbance, the effect of varying disturbances may be found therefrom by Carson's theorem. A graphical solution of Carson's integral for irregular disturbances is given. The procedure of obtaining unit solutions of the equations is then taken up and the analogy between Heaviside's symbolic series solution and a physical procedure of approximation is shown. It is suggested that a fictitious impulsive disturbance be used in the treatment of initial motions.

Jones, R. T.↗

Maximum likelihood identification of aircraft parameters with unsteady aerodynamic modelling

A simplified aerodynamic force model based on the physical principle of Prandtl's lifting line theory and trailing vortex concept has been developed to account for unsteady aerodynamic effects in aircraft dynamics. Longitudinal equations of motion have been modified to include these effects. The presence of convolution integrals in the modified equations of motion led to a frequency domain analysis utilizing Fourier transforms. This reduces the integro-differential equations to relatively simple algebraic equations, thereby reducing computation time significantly. A parameter extraction program based on the maximum likelihood estimation technique is developed in the frequency domain. The extraction algorithm contains a new scheme for obtaining sensitivity functions by using numerical differentiation. The paper concludes with examples using computer generated and real flight data

Keskar, D. A.↗

Coupled Gravity and Elevation Measurements of Ice Sheet Mass Change

We measured surface gravity and position at ten locations about two glaciological measurement networks located on the South-central Greenland Ice during June 2004. Six of the individual sites of the first network were occupied the previous year. At the repeat sites we were able to measure annual accumulation rate and surface displacement by referencing measurements to aluminum poles left in the firn the previous year. We occupied 4 additional sites at a second measurement network for the first time since initial observations were last made at the network in 1981. At each individual site, we operated a GPS unit for 90 minutes - the unit was operated simultaneously with a base station unit in Sondrestrom Fjord so as to enable differential, post-processing of the data. We installed an aluminum, accumulation-rate-pole at each site. The base section of the pole also served as the mount for the GPS antenna. A new, Scintrex gravimeter was used at each site and relative gravity measurements were tied to the network of absolute gravity stations in Sondrestrom. We measured snow physical properties in two shallow pits. This report summarizes our observations and data analysis.

Jezek, K. C.↗

Extending Science from Lunar Laser Ranging

The Lunar Laser Ranging (LLR) experiment has accumulated 50 years of range data of improving accuracy from ground stations to the laser retroreflector arrays (LRAs) on the lunar surface. The upcoming decade offers several opportunities to break new ground in data precision through the deployment of the next generation of single corner-cube lunar retroreflectors and active laser transponders. This is likely to expand the LLR station network. Lunar dynamical models and analysis tools have the potential to improve and fully exploit the long temporal baseline and precision allowed by millimetric LLR data. Some of the model limitations are outlined for future efforts. Differential observation techniques will help mitigate some of the primary limiting factors and reach unprecedented accuracy. Such observations and techniques may enable the detection of several subtle signatures required to understand the dynamics of the Earth- Moon system and the deep lunar interior. LLR model improvements would impact multi- disciplinary fields that include lunar and planetary science, Earth science, fundamental physics, celestial mechanics and ephemerides.

Vishnu Viswanathan↗

$\mathrm{SageNet}$: Fast Neural Network Emulation of the Stiff-amplified Gravitational Waves from Inflation

Accurate modeling of the inflationary gravitational waves (GWs) requires time-consuming, iterative numerical integrations of differential equations to take into account their backreaction on the expansion history. To improve computational efficiency while preserving accuracy, we present the Stiff-amplified Gravitational-wave Emulator Network (SageNet), a deep learning framework designed to replace conventional numerical solvers (code available at https://github.com/YifangLuo/SageNet). SageNet employs a long short-term memory architecture to emulate the present-day energy density spectrum of the inflationary GWs with possible stiff amplification, Ω GW (f). Trained on a data set of 25,689 numerically generated solutions, SageNet allows accurate reconstructions of Ω GW (f) and generalizes well to a wide range of cosmological parameters; 90.9% of the test emulations with randomly distributed parameters exhibit errors of under 4%. In addition, SageNet demonstrates its ability to learn and reproduce the artificial, adaptive sampling patterns in numerical calculations, which implement denser sampling of frequencies around changes in spectral indices in Ω GW (f). The dual capability of learning both physical and artificial features of the numerical GW spectra establishes SageNet as a robust alternative to exact numerical methods. Finally, our benchmark tests show that SageNet reduces the computation time from tens of seconds to milliseconds, achieving a speedup of ∼10 4 times over standard CPU-based numerical solvers with the potential for further acceleration on GPU hardware. These capabilities make SageNet a powerful tool for accelerating Bayesian inference procedures for extended cosmological models. In a broad sense, the SageNet framework offers a fast, accurate, and generalizable solution to modeling cosmological observables whose theoretical predictions demand costly differential equation solvers.

Astronomy data modeling↗

NEML2: An efficient and modular multiphysics constitutive modeling library for hybrid computing environments

This paper presents NEML2, an open-source, high-performance library developed for constitutive material modeling, designed to support the flexible and modular development of models for complex material behavior. Building on the foundational structure of its predecessor, NEML, the NEML2 library introduces significant improvements, including enhanced vectorization, automatic differentiation, and seamless integration with PyTorch, facilitating the application of machine learning techniques in material simulations. NEML2 provides a C++ backend with Python bindings, enabling users to create custom material models that can be executed efficiently on both CPU and GPU platforms. The library also supports coupling with Multiphysics simulation frameworks like MOOSE, making it suitable for realistic simulations involving coupled physical processes. Rigorous quality assurance through unit and regression testing ensures the reliability of results, while the extensible, user-friendly design encourages collaboration and reproducibility across the scientific community. This paper provides an overview of NEML2’s architecture, core features, and applications, highlighting its impact on accelerating material qualification and advancing computational methods in materials science.

GPU↗