Search NASA⌕ Search

SEARCH · Search NASA

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The root cause of disruptive NTMs and paths to stable operation in DIII-D ITER baseline scenario plasmas

Analyses of the DIII-D ITER Baseline Scenario database support that the disruptive m,n=2,1 magnetic islands are pressure gradient driven, non-linear instabilities seeded in a sequence of stochastic transient magnetic perturbations, and that the current profile relaxation does not affect the m,n=2,1 island onset rate. At low torque, these Neoclassical Tearing Modes are most commonly seeded by non-linear 3-wave coupling when the differential rotation between the q=1 & q=2 rational surfaces approaches zero. Lack of statistically significant difference between the current profiles of stable and unstable states, as well as lack of correlation between the tearing mode onset rate and the current profile relaxation both reject causality between the current profile evolution and the 2,1 magnetic island onsets in these plasmas. These support that preserving the differential rotation between the q=1 and q=2 rational surfaces is key to long pulse stable operation in the plasma scenario planned for ITER, while optimization of the current profile within the explored parameter space may lead to much weaker improvements than sustaining the differential rotation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Glauber-theory analysis of nuclear reactions on a 12 C target with variational Monte Carlo wave functions

The application of Glauber theory has been playing an increasingly important role with the study of unstable or exotic nuclei. Its adaptation to medium and high-energy nucleus-nucleus collisions is severely limited because one has to evaluate the matrix elements of multiple-scattering operators. The extraction of physical observables has been done using ‘approximate’ Glauber theory whose validity is hard to evaluate. Here, we perform a full calculation of the matrix elements using Monte Carlo integration and analyze the elastic differential cross sections and the total reaction cross sections for p+¹²C, ⁴,⁶He+¹²C, and ¹²C+¹²C collisions. We use the variational Monte Carlo wave functions for ⁴,⁶He and ¹²C obtained by using realistic two- and three-nucleon potentials. We demonstrate the performance of the Glauber-theory calculations by comparing with available experimental data. We further discuss the accuracy of the conventional approximate methods in the light of the cumulant expansion for Glauber’s phase-shift function.

Horiuchi, W. [Osaka Metropolitan University (Japan↗
DOE OSTI2024

Optimization of a Welding Procedure for Making Critical Aluminum Welds on the LBNF Absorber Core Block

he LBNF Absorber consists of thirteen 6061-T6 aluminum core blocks. The core blocks are water cooled with de-ionized (DI) water which becomes radioactive during beam operations. The cooling water flows through gun-drilled channels in the core blocks. A weld quality optimization was performed to produce National Aeronautical Standard (NAS) 1514 Class I quality welds on the aluminum core blocks. This was not successful in all cases. An existing Gas Tungsten Arc Welding (GTAW) Welding Procedure Specification (WPS) was fine tuned to minimize, in most cases, and eliminate detectable tungsten inclusions in the welds. All the weld coupons however, passed welding inspection as per the piping code: ASME B31.3 Normal Fluid Service. Tungsten electrode diameter, type, and manufacturer were varied. Some of the samples were pre-heated and others were not. It was observed that larger diameter electrodes, 5/32 in., with pre-heated joints resulted in welds with the least number of tungsten inclusions. It is hypothesized that thinner electrodes breakdown easily and get lodged into the weld pool during the welding process. This breakdown is further enhanced by the large temperature differential between the un-preheated sample and the hot electrode.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Explicit block encodings of boundary value problems for many-body elliptic operators

Simulation of physical systems is one of the most promising use cases of future digital quantum computers. In this work we systematically analyze the quantum circuit complexities of block encoding the discretized elliptic operators that arise extensively in numerical simulations for partial differential equations, including high-dimensional instances for many-body simulations. When restricted to rectangular domains with separable boundary conditions, we provide explicit circuits to block encode the many-body Laplacian with separable periodic, Dirichlet, Neumann, and Robin boundary conditions, using standard discretization techniques from low-order finite difference methods. To obtain high-precision, we introduce a scheme based on periodic extensions to solve Dirichlet and Neumann boundary value problems using a high-order finite difference method, with only a constant increase in total circuit depth and subnormalization factor. We then present a scheme to implement block encodings of differential operators acting on more arbitrary domains, inspired by Cartesian immersed boundary methods. We then block encode the many-body convective operator, which describes interacting particles experiencing a force generated by a pair-wise potential given as an inverse power law of the interparticle distance. This work provides concrete recipes that are readily translated into quantum circuits, with depth logarithmic in the total Hilbert space dimension, that block encode operators arising broadly in applications involving the quantum simulation of quantum and classical many-body mechanics.

Kharazi, Tyler [University of California, Berkeley↗
DOE OSTI2024

The Role of Cu Species in the Regeneration of a Coked Cu/BEA Zeolite Catalyst

Active site occlusion by carbon species, often referred to as coke, is a common deactivation mechanism for heterogeneous catalysts. While it is known that transition metals can lower the temperature required for oxidative coke removal, the roles of ionic species and metal nanoparticles in coke removal are not well understood. This work aims to differentiate how ionic and nanoparticle Cu sites catalyze oxidative regeneration of a coked beta (BEA) zeolite catalyst. This was accomplished by synthesizing catalysts containing exclusively CuOx nanoparticles (NPs) or Cu2+ ions supported on BEA, physically mixing Cu/BEA with a coked BEA catalyst, and monitoring coke removal using in situ spectroscopy. Our results point to an improved combustion activity for CuOx NPs relative to Cu2+ ions, especially in the combustion of graphitic-type coke species. This improved understanding of coke combustion in transition metal containing catalysts informs strategies to improve catalyst regeneration, thereby increasing operational lifetimes.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗
DOE OSTI2025

Parametric matrix models

We present a general class of machine learning algorithms called parametric matrix models. In contrast with most existing machine learning models that imitate the biology of neurons, parametric matrix models use matrix equations that emulate physical systems. Similar to how physics problems are usually solved, parametric matrix models learn the governing equations that lead to the desired outputs. Parametric matrix models can be efficiently trained from empirical data, and the equations may use algebraic, differential, or integral relations. While originally designed for scientific computing, we prove that parametric matrix models are universal function approximators that can be applied to general machine learning problems. After introducing the underlying theory, we apply parametric matrix models to a series of different challenges that show their performance for a wide range of problems. For all the challenges tested here, parametric matrix models produce accurate results within an efficient and interpretable computational framework that allows for input feature extrapolation.

Computational science↗
DOE OSTI2024

Systems Level Fuel Cycle Modeling in TMAP8 - A Demonstration

The tritium migration analysis program (TMAP) has been used for tritium inventory tracking and analysis for several years, and the Multiphysics Object Oriented Simulation Environment (MOOSE)-based TMAP8 has several improvements over TMAP4 and TMAP7, such as support for multiple dimensions and non-cartesian coordinate systems, as well as interoperability with sub-apps at higher and lower length scales. We demonstrate that TMAP8 has the additional capacity to solve systems-level problems using zero-dimensional ordinary differential equations by reproducing a literature model which describes a systems-level fuel-cycle of tritium inventory in a hypothetical fusion power plant. The capacity to run several coupled multi-scale physics calculations as part of a single package will be necessary for accurate blanket and fuel-cycle design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Glauber-Theory Calculations of High-Energy Nuclear Scattering Observables Using Variational Monte Carlo Wave Functions

Experiments using intermediate- to high-energy radioactive nuclear beams present numerous findings. Extracting important properties of physical observables relies on a firm theoretical analysis. Though Glauber theory is believed to work well, no convincing calculation has so far been done. Here, we perform ab initio Glauber theory calculations of both elastic differential cross sections and total reaction cross sections for p+ 12 C, 12 C+ 12 C, and 6 He+ 12 C systems. The wave functions of both 6 He and 12 C are generated by variational Monte Carlo calculations with spatial and spin-isospin correlations induced by realistic two- and three-nucleon potentials. Glauber’s phase-shift function is computed by Monte Carlo integration up to all orders of nucleon-nucleon multiple scatterings. We show an excellent performance of the Glauber description to the selected data on the above systems. We also find that the cumulant expansion of the phase-shift function converges rapidly up to the second order for the above systems. This finding will open up interesting applications for the analysis of high-energy nuclear experiments.

Horiuchi, W. [Osaka Metropolitan University (Japan↗
DOE OSTI2026

Observable CMB B -modes from cosmological phase transitions

A B -mode polarization signal in the cosmic microwave background (CMB) is widely regarded as smoking gun evidence for gravitational waves produced during inflation. Here, we demonstrate that tensor perturbations sourced during noninflationary epochs can yield non-negligible B -mode signals, which can in principle complicate the interpretation of future observational data. As a case study, we consider tensor perturbations sourced in the bubble collision stage of a first-order cosmological phase transition occurring in a secluded dark sector. Although phase transitions arise from causal subhorizon physics, they nevertheless exhibit a white noise power spectrum on superhorizon scales. Power is suppressed on the large scales relevant for CMB B -mode polarization, but it is not necessarily negligible. We show that for appropriately chosen phase transition parameters, the maximal B -mode amplitude can compete with inflationary predictions that can be tested with current and future experiments. These scenarios can be differentiated by performing measurements on multiple angular scales, since the phase transition signal predicts peak power on smaller scales.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Isochronous and period-doubling diagrams for symplectic maps of the plane

Symplectic mappings of the plane serve as key models for exploring the fundamental nature of complex behavior in nonlinear systems. Central to this exploration is the effective visualization of stability regimes, which enables the interpretation of how systems evolve under varying conditions. While the area-preserving quadratic Hénon map has received significant theoretical attention, a comprehensive description of its mixed parameter-space dynamics remain lacking. This limitation arises from early attempts to reduce the full two-dimensional phase space to a one-dimensional projection, a simplification that resulted in the loss of important dynamical features. Consequently, there is a clear need for a more thorough understanding of the underlying qualitative aspects. This paper aims to address this gap by revisiting the foundational concepts of reversibility and associated symmetries, first explored in the early works of G.D. Birkhoff. We extend the original framework proposed by Hénon by adding a period-doubling diagram to his isochronous diagram, which allows to represents the system’s bifurcations and the groups of symmetric periodic orbits that emerge in typical bifurcations of the fixed point. A qualitative and quantitative explanation of the main features of the region of parameters with bounded motion is provided, along with the application of this technique to other symplectic mappings, including cases of multiple reversibility. Modern chaos indicators, such as the Reversibility Error Method (REM) and the Generalized Alignment Index (GALI), are employed to distinguish between various dynamical regimes in the mixed space of variables and parameters. These tools prove effective in differentiating regular and chaotic dynamics, as well as in identifying twistless orbits and their associated bifurcations. Additionally, we discuss the application of these methods to real-world problems, such as visualizing dynamic aperture in accelerator physics, where our findings have direct relevance.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Periodic Korteweg–de Vries soliton potentials generate quasisymmetric magnetic field strength in a finite plasma- β equilibrium

Quasisymmetry (QS) is a hidden symmetry of the magnetic field strength, B , that enables effective confinement of charged particles in a fully three-dimensional (3D) toroidal plasma equilibrium. Such equilibria are typically modeled by the ideal magnetohydrostatic (MHS) equations. The nonlinear, overdetermined nature of the QS MHS equations severely complicates our understanding of the interplay between 3D shaping, equilibrium properties such as pressure and rotational transform, and B . Progress has been made through expansions near the magnetic axis; however, a more comprehensive theory is desirable. Using a combination of analysis and regression on a large dataset of numerically optimized quasisymmetric stellarators, we demonstrate that there is a hidden lower dimensionality of B on a magnetic flux surface with connections to the theory of periodic solitons. We show that B on a flux surface is determined by three or at most four flux functions, each of which determines a critical value of the derivative of B along the field line. While being consistent with the near-axis models, our results are global and hold even on the last closed flux surface.

Differential geometry↗
DOE OSTI2024

Study of electro-physical properties of CsPbBr3 perovskite single crystals

Methods have been developed for the synthesis and growth of lead halide perovskites CsPbBr3, which can be used as detectors of optical, X-ray and gamma-radiation. Single crystals of these compounds were grown by the Bridgman method. The melt temperature was 580 °C, the gradient at the crystallization front was 5 K/cm, and the ampoule lowering speed was 3 mm/hour. In this work, the electro-physical properties of the inorganic perovskite CsPbBr3 were investigated. Crystals for research were cut from the beginning and end of the ingot. The optical absorption spectra were measured depending on the crystal thickness and temperature, the differential resistivity dif and the product  were determined from the I-V characteristics. It was established that dif practically does not change for samples made from the beginning and the end of the ingot, and is dif  21010 Ωcm at 100 V and 293 K.

Sklyarchuk, V.↗
DOE OSTI2025

Bacterial-Retted Hemp Fiber/PLA Composites

The push for sustainability in all facets of manufacturing has led to an increased interest in biomass as an alternative to non-renewable materials. Hemp bast fiber mats were produced from a bacterial retting process, named BFM, as the fiber reinforcement. The objective of this study was to evaluate the feasibility of laminating BFM with polylactic acid (PLA) for a composite panel product. Since both BFM and PLA are biodegradable, the resulting BFM-PLA composites will be 100% biodegradable. PLA pallets were processed into thin polymer sheets which served as the matrix. The BFM and PLA plates were laminated in five layers and compression-molded into composite panels. Experiments were conducted on the three BFM-to-PLA ratios (35/65, 45/55, and 50/50). Mechanical properties (tensile and bending properties) and physical properties (thickness swell and water absorption) were tested and compared to the currently commercial sheet molding compound (SMC) from fiber glass. The thermal behavior of the BFM/PLA composites was characterized using dynamic mechanical analysis (DMA) and differential scanning calorimetry (DSC). The developed BFM/PLA composite product is a sustainable alternative to existing synthetical fiber-reinforced polymer (FRP) that is biodegradable in landfill at the end of life.

Smith, Lee M. (ORCID:0000000264390357)↗
DOE OSTI2025

Machine Learning Techniques for Data Reduction of Climate Applications

Scientists conduct large-scale simulations to compute derived quantities-of-interest (QoI) from primary data. Often, QoI are linked to specific features, regions, or time intervals, such that data can be adaptively reduced without compromising the integrity of QoI. For many spatiotemporal applications, these QoI are binary in nature and represent presence or absence of a physical phenomenon. We present a pipelined compression approach that first uses neural-network-based techniques to derive regions where QoI are highly likely to be present. Then, we employ a Guaranteed Autoencoder (GAE) to compress data with differential error bounds. GAE uses QoI information to apply low-error compression to only these regions. This results in overall high compression ratios while still achieving downstream goals of simulation or data collections. Experimental results are presented for climate data generated from the E3SM Simulation model for downstream quantities such as tropical cyclone and atmospheric river detection and tracking. These results show that our approach is superior to comparable methods in the literature.

Li, Xiao [University of Florida]↗
DOE OSTI2025

Field-level reconstruction from foreground-contaminated 21-cm maps

Current and upcoming 21-cm experiments will soon be able to map 21-cm spatial fluctuations in three dimensions for a wide range of redshifts. However, bright foreground contamination and the nature of radio interferometry create significant challenges, making it difficult to access rich cosmological information from the Fourier modes that lie within the “foreground wedge”. Here, in this work, we introduce two approaches aiming to reconstruct the full 21-cm density field, including the missing modes in the wedge: (a) a field-level inference under an effective field theory (EFT) framework; (b) a diffusion-based deep generative model trained on simulations. Under the EFT framework, we implement a fully differentiable forward model that maps the initial conditions of matter fluctuations to the observed, foreground-filtered 21-cm maps. This enables a gradient-based sampler to simultaneously sample the initial conditions and bias parameters, allowing a physically motivated mode reconstruction. Alternatively, we apply a variational diffusion model to perform 21-cm density reconstruction at the map level. Our model is trained on semi-numerical simulations over a wide range of astrophysical parameters. Our results from both approaches should provide improved cosmological constraints from the field level and also enable cross-correlation between experiments that have little or no overlapping modes.

cosmological perturbation theory↗
DOE OSTI2024

Fiducial and differential cross-section measurements of electroweak $W\gamma jj$ production in $pp$ collisions at $\sqrt{s} = 13$ TeV with the ATLAS detector

The observation of the electroweak production of a W boson and a photon in association with two jets, using pp collision data at the Large Hadron Collider at a centre of mass energy of $\sqrt{s} = 13$ TeV, is reported. The data were recorded by the ATLAS experiment from 2015 to 2018 and correspond to an integrated luminosity of 140 fb -1 . This process is sensitive to the quartic gauge boson couplings via the vector boson scattering mechanism and provides a stringent test of the electroweak sector of the Standard Model. Events are selected if they contain one electron or muon, missing transverse momentum, at least one photon, and two jets. Multivariate techniques are used to distinguish the electroweak $W\gamma jj$ process from irreducible background processes. The observed significance of the electroweak $W\gamma jj$ process is well above six standard deviations, compared to an expected significance of 6.3 standard deviations. Fiducial and differential cross sections are measured in a fiducial phase space close to the detector acceptance, which are in reasonable agreement with leading order Standard Model predictions from MADGRAPH5+PYTHIA8 and SHERPA. The results are used to constrain new physics effects in the context of an effective field theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

A structural underpinning of the lower critical solution temperature (LCST) behavior behind temperature-switchable liquids

In this work, we use state-of-the-art X-ray scattering and molecular dynamics to analyze amine-water mixtures that show the unusual lower critical solution temperature (LCST) behavior. The goal is to provide direct experimental evidence for the entropy-lowering molecular cluster formation hypothesized as necessary for LCST behavior. Differential wide-angle X-ray scattering and pair distribution analysis and small-angle X-ray scattering measurements were combined with molecular modeling and liquid-liquid equilibrium measurements, revealing direct experimental evidence for the hypothesis. However, the response of the amine phase to accommodating water is even more subtle than the simple hypothesis suggests, with the formation of robust nanoscale reverse micelles. The techniques developed in this paper can be expected to yield insights in the use of temperature-switchable liquids in solvent extraction and other separations, and the stabilization of organelles in living cells that do not have physical membranes but do require compositional gradients to operate.

36 MATERIALS SCIENCE↗