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 577 records · Page 32

Symplectic machine learning model for fast simulation of space-charge effects

Symplectic simulation of space-charge effects is crucial for the design and operation of high-intensity particle accelerators. Traditional methods for simulating these effects are often computationally expensive, resulting in significant overhead. In this work, we introduce a generative model based on a U-Net architecture within a generative adversarial network framework to efficiently simulate space-charge effects. The model is trained to predict the transverse multiparticle space-charge Hamiltonian, which can be physically computed using a gridless spectral method. The one-step symplectic transverse transfer map for the particles is then obtained by differentiating the predicted Hamiltonian. Benchmarking results demonstrate that this generative model achieves an order of magnitude higher computational efficiency compared to the spectral method, providing a highly efficient alternative for simulating space-charge effects with a large number of particles. By maintaining symplecticity, the model effectively preserves the phase-space structure and mitigates nonphysical errors in long-term simulations. This model has been integrated into jutrack, a novel autodifferentiable accelerator modeling code developed in the julia programming language.

Beam code development & simulation techniques↗

Photoproduction of η c , b near threshold

We analyze the photoproduction of η c , b off a proton in the threshold region, in terms of C-odd gluonic correlations in the off-forward proton matrix element. Near threshold, the skewness is large, leading to a production amplitude that is dominated by four C-odd twist-3 gluon GPDs. We use the QCD instanton vacuum to estimate these C-odd contributions in the proton. The results are used to estimate the differential cross sections for coherent photoproduction of η c , b in the threshold region, at current electron facilities. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

A bacterial sensor taxonomy across earth ecosystems for machine learning applications

Microbial communities have evolved to colonize all ecosystems of the planet, from the deep sea to the human gut. Microbes survive by sensing, responding, and adapting to immediate environmental cues. This process is driven by signal transduction proteins such as histidine kinases, which use their sensing domains to bind or otherwise detect environmental cues and “transduce” signals to adjust internal processes. We hypothesized that an ecosystem’s unique stimuli leave a sensor “fingerprint,” able to identify and shed insight on ecosystem conditions. To test this, we collected 20,712 publicly available metagenomes from Host-associated, Environmental, and Engineered ecosystems across the globe. We extracted and clustered the collection’s nearly 18M unique sensory domains into 113,712 similar groupings with MMseqs2. We built gradient-boosted decision tree machine learning models and found we could classify the ecosystem type (accuracy: 87%) and predict the levels of different physical parameters (R2 score: 83%) using the sensor cluster abundance as features. Feature importance enables identification of the most predictive sensors to differentiate between ecosystems which can lead to mechanistic interpretations if the sensor domains are well annotated. To demonstrate this, a machine learning model was trained to predict patient’s disease state and used to identify domains related to oxygen sensing present in a healthy gut but missing in patients with abnormal conditions. Moreover, since 98.7% of identified sensor domains are uncharacterized, importance ranking can be used to prioritize sensors to determine what ecosystem function they may be sensing. Furthermore, these new predictive sensors can function as targets for novel sensor engineering with applications in biotechnology, ecosystem maintenance, and medicine.

97 MATHEMATICS AND COMPUTING↗

Collaborative Research: Understanding Nanoparticle-Plasma Interactions in Dusty Non-Thermal Plasmas by Nanoparticle Probes and Online Aerosol Characterization (Final Technical Report)

This report describes results from the University of Minnesota’s contribution to studies performed under the award DE-SC0018202 (principal Investigator: Christopher J. Hogan Jr.. This award is a collaborative research grant with support equally provided to Case Western Reserve University and to the University of Illinois (with principal investigator R. Mohan Sankaran). To date ten peer-reviewed published papers have resulted from the studies performed with University of Minnesota personnel contributing; this report describes the results of these studies as well as several unpublished studies, grouped into seven separate study types. The overarching goal of the carried out studies was to understand nanoparticle behavior and dynamics in non-thermal plasma synthesis reactors, specifically adapting aerosol measurement technology to provide online measurements of particle size, charge, and structure in plasma reactor systems.

42 ENGINEERING↗

DS-TIDE: Harnessing Dynamical Systems for Efficient Time-Independent Differential Equation Solving

Time-Independent Differential Equations (TIDEs) are central to modeling equilibrium behavior across a wide range of scientific and engineering domains, from electrostatics to porous media flow. Conventional numerical solvers offer reliable solutions but incur significant computational costs due to fine-grained discretization and iterative procedures. Machine learning-based approaches address this by replacing iterative solving processes with one-time inference; however, their sophisticated models require extensive training resources that often exceed those of traditional solvers. Consequently, designing a TIDE solver that achieves high accuracy, broad applicability, and exceptional computational efficiency remains a fundamental challenge. In this paper, we propose DS-TIDE, a novel hardware solver that is inspired by, and subsequently leverages, the intrinsic connection between Dynamical Systems (DS) and Differential Equations (DEs) to efficiently and accurately solve TIDEs. DS-TIDE employs a CMOS-compatible DS-based processor, whose physical states evolve under carefully designed DE-driven dynamics and naturally converge to equilibrium -- the solution of the target TIDE -- within ~1µs on a ~1-watt DS-TIDE processor. To enhance expressivity, DS-TIDE incorporates Heterogeneous Dynamics with Temporal Layering (HDTL), which solves TIDEs through a three-stage DS evolution -- conditioning, solving, and decoding -- each governed by specialized dynamics. The entire evolution process is analogous to an infinitely deep neural network temporally unrolled, offering the system the capability of representing complex equations. Furthermore, DS-TIDE is equipped with an on-device DS-DE Auto-Alignment mechanism that dynamically adapts intrinsic hardware dynamics within milliseconds, effectively aligning the system’s dynamics to diverse target DEs. Experimental results across TIDEs from a wide range of scientific and engineering domains demonstrate that DS-TIDE achieves ~10^3× speedup, ~10^5× energy savings, and competitive or superior accuracy compared to state-of-the-art numerical and ML-based solvers.

Liu, Chuan↗

Local lattice distortions and electronic orders in strongly correlated systems by resonant total x-ray scattering: A case study of A Pt 2 X 2 intermetallics ( A = U , Ce, or La and X = Si or Ge)

Revealing local lattice distortions in physical systems is nontrivial, largely because such structural features often are not necessarily amenable to traditional crystallographic investigations. This poses limits to our understanding of the underlying physic, particularly in the case of strongly correlated systems where structural, electronic and magnetic degrees of freedom are intertwined. Using resonant total x-ray scattering coupled to differential pair distribution function analysis, we reveal the presence of pronounced local lattice distortions in the square net Pt planes in the prototypical strongly correlated APt 2 X 2 intermetallics (A = U, Ce or La and X = Si or Ge). The distortions are present before charge density wave, magnetic, Kondo lattice coherence and/or superconducting orders emerge in these materials and, as density functional calculations suggest, are likely to affect them. Finally, our study sheds light on the poorly known interactions between electronic orders and structural disorder in strongly correlated systems and also demonstrates an advanced experimental approach to determine it that is relevant to any physical system showing deviations from perfect crystallinity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

From light neutrino decay phenomenology to muon neutrino cross-section measurement at MINERvA experiment

This thesis tackles two essential topics of neutrino physics: neutrino decay and neutrino cross-section measurement. First, the invisible and visible neutrino decay is analyzed through a phenomenological approach, considering future long-baseline neutrino experiments such as DUNE and a hypothetical neutrino beam toward the ANDES laboratory. The study takes into account the ν μ and ν e disappearance and appearance, respectively, for both FHC and RHC flux modes. The results showed a negligible matter effect for DUNE but significantly more notable at ANDES. At 90% C.L., the sensitivity to the decay parameter α3 can be as small as 2 × 10 -6 eV 2 for a chosen coupling. The impact of neutrino decay in the determination of θ 23 and δ CP were also shown. Second, the double-differential cross-section measurement for ν μ -carbon interactions with three-momentum transfer |$\vec{q}$| < 1.2 GeV obtained with medium energy exposures in the NuMI beam at MINERvA experiment are reported. The measurement is presented as a function of |$\vec{q}$| and E avail and reviews different interaction models and nuclear effects along quasi-elastic to resonance processes to define a new model for a better agreement. The double differential cross sections are compared to the MnvTunes, GENIE, and NuWro predictions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Strictly Enforcing Invertibility and Conservation in CNN-Based Super Resolution for Scientific Datasets

Abstract Recently, deep convolutional neural networks (CNNs) have revolutionized image “super resolution” (SR), dramatically outperforming past methods for enhancing image resolution. They could be a boon for the many scientific fields that involve imaging or any regularly gridded datasets: satellite remote sensing, radar meteorology, medical imaging, numerical modeling, and so on. Unfortunately, while SR-CNNs produce visually compelling results, they do not necessarily conserve physical quantities between their low-resolution inputs and high-resolution outputs when applied to scientific datasets. Here, a method for “downsampling enforcement” in SR-CNNs is proposed. A differentiable operator is derived that, when applied as the final transfer function of a CNN, ensures the high-resolution outputs exactly reproduce the low-resolution inputs under 2D-average downsampling while improving performance of the SR schemes. The method is demonstrated across seven modern CNN-based SR schemes on several benchmark image datasets, and applications to weather radar, satellite imager, and climate model data are shown. The approach improves training time and performance while ensuring physical consistency between the super-resolved and low-resolution data. Significance Statement Recent advancements in using deep learning to increase the resolution of images have substantial potential across the many scientific fields that use images and image-like data. Most image super-resolution research has focused on the visual quality of outputs, however, and is not necessarily well suited for use with scientific data where known physics constraints may need to be enforced. Here, we introduce a method to modify existing deep neural network architectures so that they strictly conserve physical quantities in the input field when “super resolving” scientific data and find that the method can improve performance across a wide range of datasets and neural networks. Integration of known physics and adherence to established physical constraints into deep neural networks will be a critical step before their potential can be fully realized in the physical sciences.

54 ENVIRONMENTAL SCIENCES↗

Parametrizations of the spin density matrix

Highlights: • Proposed two methods for parametrizing the SDM of spinning, unstable particles. • All non-negativity conditions fulfilled; adaptability to various constraints. • First method: based on a theorem; applicable in a simple and flexible way. • Second method: particularly suitable for especial conditions. • Applications: two of them interesting in the search for new physics. We propose for the spin density matrix two parametrizations which automatically fulfill the non-negativity conditions, without setting any bound on the parameters. The first one relies on a theorem, that we prove, and it is rather simple and easily adaptable to some specific reactions, where, for example, parity is conserved or angular momentum conservation entails selection rules. Moreover, in the case when the rank is less than the order of the density matrix, we show how to improve the fits to the data, either by implementing previous suggestions, or by elaborating an alternative method, for which we prove a second theorem. Our second parametrization is a variant of previous treatments, it appears suitable for some particular processes. Moreover, we discuss about the possibility of inferring the elements of the density matrix from the differential decay width. Last, we illustrate various examples of current interest, both in strong and weak decays; some of them may be helpful in the investigation of physics beyond the standard model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Evolution of Electric Sector Water Use Under Alternative Electricity Futures

Increasing competition for water resources in the United States could create future challenges for allocating and using thermal cooling water in the U.S. electric power sector. While thermal power plant retirements and the growth of wind and solar technologies can reduce national aggregate power sector cooling water use, local water constraints and growing demand for agricultural or municipal supply could create greater needs for higher-cost alternative water supplies such as groundwater or recycled wastewater. For some regions, these incentives could change future electricity planning and operational decisions. These relationships and impacts are studied here using the National Renewable Energy Laboratory (NREL) Regional Energy Deployment System (ReEDS), a national electric sector planning model that has recently been upgraded to include a highly detailed representation of thermal cooling water demand and supply. Thermal power technologies are differentiated by both cooling technology and water source type to track and constrain thermal cooling water use in a way that incorporates both physical and legal water considerations in the United States. This capability is exercised under a range of electricity sector futures with alternative technology costs, fuel prices, and water constraints to illustrate ways that U.S. electric sector water use could evolve under uncertain future electric sector drivers. In exploring changes to regional generation and transmission planning and operation, water requirements, and cost, we highlight the environmental and economic impacts of future power sector water decisions.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Dynamic modelling of slip in a wind turbine spherical roller main bearing

This paper considers the problem of the dynamic modelling of macro slip in spherical roller bearings. By revisiting the fundamental physics which drive these systems, potential issues in existing models have been identified. Furthermore, in pure rolling conditions it was found that governing differential equations become "stiff", requiring the use of implicit methods of time integration. The problem of individual roller macro slip in a wind turbine main bearing is then investigated using a simplified representation of system dynamics. Model results indicate clear links between slip/friction and the operational strategy of the wind turbine, as well as significantly higher frictional effects in the downwind main bearing row. Due to modelling simplifications, these results should not yet be considered conclusive, with further work required.

17 WIND ENERGY↗

Experimental investigation and thermodynamic modeling of the binary Pt–Sn system

Here, the thermodynamics of the Pt–Sn system is revisited with attention to phase equilibria involving PtSn 4 , a topological quantum material with interesting physical properties. The phase relationship in the Sn-rich region of the Pt–Sn binary system was investigated by a combination of differential scanning calorimetry and equilibration experiments. Based on the present experiments and a critical evaluation of the available thermodynamic and phase diagram data reported in the literature, thermodynamic assessments of the entire Pt–Sn system were performed using the CALPHAD approach. Two sets of thermodynamic parameters were optimized, where the liquid phase was described using the Bragg-Williams random mixing model and the associate solution model, respectively. The modeling results are compared with the reliable experimental data and are shown to provide good representations of the phase equilibria and thermodynamic properties in this system.

36 MATERIALS SCIENCE↗

Characterization of partial wetting by CMAS droplets using multiphase many-body dissipative particle dynamics and data-driven discovery based on PINNs

The molten sand that is a mixture of calcia, magnesia, alumina and silicate, known as CMAS, is characterized by its high viscosity, density and surface tension. The unique properties of CMAS make it a challenging material to deal with in high-temperature applications, requiring innovative solutions and materials to prevent its buildup and damage to critical equipment. Here, we use multiphase many-body dissipative particle dynamics simulations to study the wetting dynamics of highly viscous molten CMAS droplets. The simulations are performed in three dimensions, with varying initial droplet sizes and equilibrium contact angles. We propose a parametric ordinary differential equation (ODE) that captures the spreading radius behaviour of the CMAS droplets. The ODE parameters are then identified based on the physics-informed neural network (PINN) framework. Subsequently, the closed-form dependency of parameter values found by the PINN on the initial radii and contact angles are given using symbolic regression. Finally, we employ Bayesian PINNs (B-PINNs) to assess and quantify the uncertainty associated with the discovered parameters. In brief, this study provides insight into spreading dynamics of CMAS droplets by fusing simple parametric ODE modelling and state-of-the-art machine-learning techniques.

Mechanics↗

The nonexistence of a paddlewheel effect in superionic conductors

Since the 1980s, the paddlewheel effect has been suggested as a mechanism to boost lithium-ion diffusion in inorganic materials via the rotation of rotor-like anion groups. However, it remains unclear whether the paddlewheel effect, defined as large-angle anion group rotations assisting Li hopping, indeed exists; furthermore, the physical mechanism by which the anion-group dynamics affect lithium-ion diffusion has not yet been established. In this work, we differentiate various types of rotational motions of anion groups and develop quaternion-based algorithms to detect, quantify, and relate them to lithium-ion motion in ab initio molecular dynamics simulations. Our analysis demonstrates that, in fact, the paddlewheel effect, where an anion group makes a large angle rotation to assist a lithium-ion hop, does not exist and thus is not responsible for the fast lithium-ion diffusion in superionic conductors, as historically claimed. Instead, we find that materials with topologically isolated anion groups can enhance lithium-ion diffusivity via a more classic nondynamic soft-cradle mechanism, where the anion groups tilt to provide optimal coordination to a lithium ion throughout the hopping process to lower the migration barrier. This anion-group disorder is static in nature, rather than dynamic and can explain most of the experimental observations. Our work substantiates the nonexistence of the long-debated paddlewheel effect and clarifies any correlation that may exist between anion-group rotations and fast ionic diffusion in inorganic materials.

36 MATERIALS SCIENCE↗

AMReX v2024

The software framework, AMReX, supports the development of block-structured adaptive mesh refinement (AMR) algorithms for solving systems of partial differential equations. AMR reduces the computational cost and memory footprint compared to a uniform mesh while preserving the essential local descriptions of different physical processes in complex multiphysics algorithms. AMR uses a hierarchical representation of the solution at multiple levels of resolution where the solution on each level is defined on the union of data containers at that resolution. These data containers, which represent the solution over a logically rectangular subregion of the domain, can contain field data defined on a mesh, Lagrangian particles or combinations of both. In addition to these basic data types, AMReX supports a multilevel embedded boundary representation of complex geometry; linear solvers for cell-centered and nodal data; asynchronous I/O in a native format readable by ParaView, VisIt and yt; and interfaces to hypre and PETSc solvers. AMReX enables applications to run on distributed memory architectures with multicore CPUs and with GPU accelerators. AMReX uses a lightweight abstraction layer that effectively hides the details of the architecture from the application. The framework currently supports CUDA, HIP and SYCL for GPU acceleration and OpenMP for multi-core CPU architectures.

Almgren, Ann↗