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

Hadamard products and BPS networks

We study examples of fourth-order Picard-Fuchs operators that are Hadamard products of two second-order Picard-Fuchs operators. Each second-order Picard-Fuchs operator is associated with a family of elliptic curves, and the Hadamard product computes period integrals on the fibred product of the two elliptic surfaces. We construct 3-cycles on this geometry as the union of 2-cycles in the fibre over contours on the base. We then use the special Lagrangian condition to constrain the contours on the base. This leads to a construction that is reminiscent of spectral networks and exponential networks that have previously appeared in string theory literature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Phospholipase D and phosphatidylinositol-4-phosphate 5-kinase 1 are involved in the regulation of oligodendrocyte morphological differentiation

Highlights: • PLD and PIP5K1 are involved in the regulation of oligodendroglial cell morphological differentiation. • Inhibition of PLD or PIP5K1 decreases differentiation or myelin marker expression. • PLD and PIP5K1 regulate differentiation of oligodendrocyte precursor cells. • Inhibition of PLD or PIP5K1 decreases differentiation marker expression of oligodendrocyte precursor cells. Oligodendroglial cells (oligodendrocytes) differentiate to form the myelin that wraps neuronal axons in the central nervous system (CNS). This myelin sheath supports the propagation of saltatory conduction and protects axons from physical stresses. When oligodendrocytes do not normally differentiate to myelinate axons, their key functions as oligodendrocytes in the CNS are severely impaired. The molecular mechanics that control differentiation still remain to be clarified. Arf6 belongs to the small GTPase family and is known to be a positive regulator of oligodendrocyte differentiation. Here, we show that the phospholipase D (PLD) and phosphatidylinositol-4-phosphate 5-kinase 1 (PIP5K1) molecules, the major effectors of Arf6, are involved in the regulation of oligodendrocyte differentiation. Knockdown of PLD1 or PIP5K type 1γ (PIP5K1C) by their respective specific siRNAs in mouse oligodendroglial FBD-102b cells inhibited morphological differentiation into structures bearing myelin-like processes; this finding is consistent with the concurrent changes in expression of differentiation and myelin marker proteins. Treatment with VU0155069 or UNC3230, specific inhibitors of PLD and PIP5K1, respectively, blunted morphological differentiation and decreased expression of myelin and differentiation marker proteins. Similar results have been obtained in studies using primary oligodendrocytes. These results suggest that the major Arf6 effector molecules PLD and PIP5K1 are among the molecules involved in the regulation of morphological differentiation in oligodendrocytes prior to myelination.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Measurements of differential cross-sections in four-lepton events in 13 TeV proton-proton collisions with the ATLAS detector

Measurements of four-lepton differential and integrated fiducial cross-sections in events with two same-flavour, opposite-charge electron or muon pairs are presented. The data correspond to 139 fb -1 of $ \sqrt{s} $ = 13 TeV proton-proton collisions, collected by the ATLAS detector during Run 2 of the Large Hadron Collider (2015–2018). The final state has contributions from a number of interesting Standard Model processes that dominate in different four-lepton invariant mass regions, including single Z boson production, Higgs boson production and on-shell ZZ production, with a complex mix of interference terms, and possible contributions from physics beyond the Standard Model. The differential cross-sections include the four-lepton invariant mass inclusively, in slices of other kinematic variables, and in different lepton flavour categories. Also measured are dilepton invariant masses, transverse momenta, and angular correlation variables, in four regions of four-lepton invariant mass, each dominated by different processes. The measurements are corrected for detector effects and are compared with state-of-the-art Standard Model calculations, which are found to be consistent with the data. The Z → 4 ℓ branching fraction is extracted, giving a value of (4.41 ± 0.30) × 10 -6 . Constraints on effective field theory parameters and a model based on a spontaneously broken B – L gauge symmetry are also evaluated. Further reinterpretations can be performed with the provided information.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1987

Rotational modulation and flares on RS CVn and BY Dra stars. VI - Physical parameters of the chromospheres/transition regions of V 711 Tau (HR 1099), II Peg and AR Lac during October 1981

Ground-based optical and IUE satellite-ultraviolet observations of three RS CVn stars are combined with density sensitive line ratios and differential emission measure curves to describe the physical conditions in their outer atmospheres. Solar-like densities are found to be representative of average conditions on two of the stars, V 711 Tau and AR Lac. The total radiative losses from these two stars are estimated and found to be larger than the sun by at least two orders of magnitude. Consideration of the volume emitting in two of the principal transition region lines suggests a possible relation between the disk 'filling factor' for these two lines and the dynamo-related Rossby number. Only one hemisphere of the star II Peg, the one showing least evidence of starspots, is similar. On the opposite hemisphere the presence of a discrete emitting region is deduced which is almost coincident in phase with the passage of the dominant optical spot group across the visible disk. The dimensions of this region, and an illustrative interpretation in terms of a large emitting magnetic loop, are discussed. It is compared to large active region loops on the sun.

Byrne, P. B.↗
DOE OSTI2023

Machine learning of hidden variables in multiscale fluid simulation

Solving fluid dynamics equations often requires the use of closure relations that account for missing microphysics. For example, when solving equations related to fluid dynamics for systems with a large Reynolds number, sub-grid effects become important and a turbulence closure is required, and in systems with a large Knudsen number, kinetic effects become important and a kinetic closure is required. By adding an equation governing the growth and transport of the quantity requiring the closure relation, it becomes possible to capture microphysics through the introduction of ‘hidden variables’ that are non-local in space and time. The behavior of the ‘hidden variables’ in response to the fluid conditions can be learned from a higher fidelity or ab-initio model that contains all the microphysics. In our study, a partial differential equation simulator that is end-to-end differentiable is used to train judiciously placed neural networks against ground-truth simulations. We show that this method enables an Euler equation based approach to reproduce non-linear, large Knudsen number plasma physics that can otherwise only be modeled using Boltzmann-like equation simulators such as Vlasov or particle-in-cell modeling.

42 ENGINEERING↗
DOE OSTI2020

hp -VPINNs: Variational physics-informed neural networks with domain decomposition

We formulate a general framework for hp-variational physics-informed neural networks (hp-VPINNs) based on the nonlinear approximation of shallow and deep neural networks and hp-refinement via domain decomposition and projection onto the space of high-order polynomials. The trial space is the space of neural network, which is defined globally over the entire computational domain, while the test space contains piecewise polynomials. Specifically in this study, the hp-refinement corresponds to a global approximation with a local learning algorithm that can efficiently localize the network parameter optimization. Here, we demonstrate the advantages of hp-VPINNs in both accuracy and training cost for several numerical examples of function approximation and in solving differential equations.

42 ENGINEERING↗
DOE OSTI2021

4-manifolds and topological modular forms

We build a connection between topology of smooth 4-manifolds and the theory of topological modular forms by considering topologically twisted compactification of 6d (1, 0) theories on 4-manifolds with flavor symmetry backgrounds. The effective 2d theory has (0, 1) supersymmetry and, possibly, a residual flavor symmetry. The equivariant topological Witten genus of this 2d theory then produces a new invariant of the 4-manifold equipped with a principle bundle, valued in the ring of equivariant weakly holomorphic (topological) modular forms. We describe basic properties of this map and present a few simple examples. As a byproduct, we obtain some new results on ’t Hooft anomalies of 6d (1, 0) theories and a better understanding of the relation between 2d (0, 1) theories and TMF spectra.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Decorated TQFTs and their Hilbert spaces

We discuss topological quantum field theories that compute topological invariants which depend on additional structures (or decorations) on three-manifolds. The q-series invariant $\hat{Z}$(q) proposed by Gukov, Pei, Putrov, and Vafa is an example of such an invariant. We describe how to obtain these decorated invariants by cutting and gluing and make a proposal for Hilbert spaces that are assigned to two-dimensional surfaces in the $\hat{Z}$-TQFT.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Implementing Ordinary Differential Equation Solvers in Rust Programming Language for Modeling Vehicle Powertrain Systems: Preprint

Efficient and accurate ordinary differential equation (ODE) solvers are necessary for powertrain and vehicle dynamics modeling. However, current commercial ODE solvers can be financially prohibitive, leading to a need for accessible, effective, open-source ODE solvers designed for powertrain modeling. Rust is a compiled programming language that has the potential to be used for fast and easy-to-use powertrain models, given its exceptional computational performance, robust package ecosystem, and short time required for modelers to become proficient. However, of the three commonly used (>3,000 downloads) packages in Rust with ODE solver capabilities, only one has more than four numerical methods implemented, and none are designed specifically for modeling physical systems. Therefore, the goal of the Differential Equation System Solver (DESS) was to implement accurate ODE solvers in Rust designed for the component-based problems often seen in powertrain modeling. DESS is a text-based software package that provides a flexible framework for building and solving systems of ODEs. This allows DESS to be included as a dependency for automotive powertrain models that require a variety of solvers and solver configurations. Seven explicit ODE solver methods have been implemented in DESS: Euler’s, Heun’s, midpoint, Ralston’s, classic Runge-Kutta, Bogacki-Shampine, and Cash-Karp. These represent five fixed-step methods and two adaptive-step methods. This paper shows that the solver implementations increase accuracy and computational efficiency compared to Euler's method when modeling a system of three thermal masses in Rust. DESS also includes features designed for modeling component-based physical systems. Users can define relationships between nodes in their system, which the package then translates into a system of equations, leading to simpler and more intuitive code. In the case of a three-thermal-mass system, the user can specify node thermal properties (e.g., thermal capacitance), how nodes are interconnected, and thermal conductance between nodes rather than providing a system of equations. The core contribution from this work is an open-source, text-based Rust package with ODE solvers for automotive powertrain modeling to support cost-free, fast, and accurate simulation.

ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Application of Quantum Machine Learning to High Energy Physics Analysis at LHC using IBM Quantum Computer Simulators and IBM Quantum Computer Hardware

One of the major objectives of the experimental programs at the LHC is the discovery of new physics. This requires the identification of rare signals in immense backgrounds. Using machine learning algorithms greatly enhances our ability to achieve this objective. With the progress of quantum technologies, quantum machine learning could become a powerful tool for data analysis in high energy physics. In this study, using IBM gate-model quantum computing systems, we employ the quantum variational classifier method and the quantum kernel estimator method in two recent LHC flagship physics analyses: $t\bar{t}H$ (Higgs boson production in association with a top quark pair) and $H\rightarrow\mu\mu$ (Higgs boson decays to two muons). We have obtained early results with 10 qubits on the IBM quantum simulator and the IBM quantum hardware. On the quantum simulator, the quantum machine learning methods perform similarly to classical algorithms such as SVM (support vector machine) and BDT (boosted decision tree), which are often employed in LHC physics analyses. On the quantum hardware, the quantum machine learning methods have shown promising discrimination power, comparable to that on the quantum simulator. This study demonstrates that quantum machine learning has the ability to differentiate between signal and background in realistic physics datasets.

Chan, Jay↗
NASA NTRS1987

A framework for qualitative reasoning about solid objects

Predicting the behavior of a qualitatively described system of solid objects requires a combination of geometrical, temporal, and physical reasoning. Methods based upon formulating and solving differential equations are not adequate for robust prediction, since the behavior of a system over extended time may be much simpler than its behavior over local time. A first-order logic, in which one can state simple physical problems and derive their solution deductively, without recourse to solving the differential equations, is discussed. This logic is substantially more expressive and powerful than any previous AI representational system in this domain.

Davis, E.↗
NASA NTRS2007

Application of the Finite Element Method in Atomic and Molecular Physics

The finite element method (FEM) is a numerical algorithm for solving second order differential equations. It has been successfully used to solve many problems in atomic and molecular physics, including bound state and scattering calculations. To illustrate the diversity of the method, we present here details of two applications. First, we calculate the non-adiabatic dipole polarizability of Hi by directly solving the first and second order equations of perturbation theory with FEM. In the second application, we calculate the scattering amplitude for e-H scattering (without partial wave analysis) by reducing the Schrodinger equation to set of integro-differential equations, which are then solved with FEM.

Shertzer, Janine↗
DOE OSTI2021

Symmetry-breaking-induced rare fluctuations in a time-delay dynamic system

Inspired by the experimental and numerical findings, we study the dynamic instabilities of two coupled nonlinear delay differential equations that are used to describe the coherent oscillations between the top and bottom boundary layers in turbulent Rayleigh–Bénard convection. By introducing two sensitivity parameters for the instabilities of the top and bottom boundary layers, we find three different types of solutions, namely in-phase single-period oscillations, multi-period oscillations and chaos. The chaos solution contains rare but large amplitude fluctuations. The statistical properties of these fluctuations are consistent with those observed in the experiment for the massive eruption of thermal plumes, which causes random reversals of the large-scale circulation in turbulent Rayleigh–Bénard convection. Finally, our study thus provides new insights into the origin of rare massive eruptions and sudden changes of large-scale flow pattern that are often observed in convection systems of geophysical and astrophysical scales.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Improved Architectures and Training Algorithms for Deep Operator Networks

Operator learning techniques have recently emerged as a powerful tool for learning maps between infinite-dimensional Banach spaces. Trained under appropriate constraints, they can also be effective in learning the solution operator of partial differential equations (PDEs) in an entirely self-supervised manner. In this work we analyze the training dynamics of deep operator networks (DeepONets) through the lens of Neural Tangent Kernel theory, and reveal a bias that favors the approximation of functions with larger magnitudes. To correct this bias we propose to adaptively re-weight the importance of each training example, and demonstrate how this procedure can effectively balance the magnitude of back-propagated gradients during training via gradient descent. We also propose a novel network architecture that is more resilient to vanishing gradient pathologies. Taken together, our developments provide new insights into the training of DeepONets and consistently improve their predictive accuracy by a factor of 10-50x, demonstrated in the challenging setting of learning PDE solution operators in the absence of paired input-output observations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Differential equations for cosmological correlators

Cosmological fluctuations retain a memory of the physics that generated them in their spatial correlations. The strength of correlations varies smoothly as a function of external kinematics, which is encoded in differential equations satisfied by cosmological correlation functions. In this work, we provide a broader perspective on the origin and structure of these differential equations. As a concrete example, we study conformally coupled scalar fields in a power-law cosmology. The wavefunction coefficients in this model have integral representations, with the integrands being the product of the corresponding flat-space results and “twist factors” that depend on the cosmological evolution. Similar twisted integrals arise for loop amplitudes in dimensional regularization, and their recent study has led to the discovery of rich mathematical structures and powerful new tools for computing multi-loop Feynman integrals in quantum field theory. The integrals of interest in cosmology are also part of a finite-dimensional basis of master integrals, which satisfy a system of first-order differential equations. We develop a formalism to derive these differential equations for arbitrary tree graphs. The results can be represented in graphical form by associating the singularities of the differential equations with a set of graph tubings. Upon differentiation, these tubings grow in a local and predictive fashion. In fact, a few remarkably simple rules allow us to predict — by hand — the equations for all tree graphs. While the rules of this “kinematic flow” are defined purely in terms of data on the boundary of the spacetime, they reflect the physics of bulk time evolution. We also study the analogous structures in tr ϕ 3 theory, and see some glimpses of hidden structure in the sum over planar graphs. This suggests that there is an autonomous combinatorial or geometric construction from which cosmological correlations, and the associated spacetime, emerge.

Cosmological models↗
DOE OSTI2025

Opportunities in multiscale modeling of mosquito-borne flaviviruses

Mosquito-borne flaviviruses, such as Zika, dengue, West Nile, and yellow fever virus, represent a growing public health concern due to their widespread distribution and the severe diseases they cause. These viruses are difficult to control as climate change and urbanization help mosquitoes expand into new areas, increasing the risk of outbreaks. Mathematical models play a key role in understanding their spread, providing insights at every level—from how the virus multiplies inside cells to how it circulates through entire populations. This review examines various approaches used in modeling arboviruses, including microscale models that focus on cellular and molecular dynamics, mesoscale models that address within-host processes, and macroscale models that capture population-level transmission. We briefly summarize the methodology used for models at each scale, which primarily consists of sets of differential equations with parameters that represent physical rates of change for different subprocesses. We particularly highlight how temperature affects virus transmission, which is key to understanding the impact of climate change. We also show how multiscale models can connect viral replication, immune response, and the spread of infection at a larger scale. This is essential for developing better vaccines and treatments, evaluating disease control measures, predicting the impact of climate change, and improving public health responses to outbreaks.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Machine learning feature analysis illuminates disparity between E3SM climate models and observed climate change

In September of 2020, Arctic sea ice extent was the second-lowest on record. State of the art climate prediction uses Earth system models (ESMs), driven by systems of differential equations representing the laws of physics. Previously, these models have tended to underestimate Arctic sea ice loss. The issue is grave because accurate modeling is critical for economic, ecological, and geopolitical planning. We use machine learning techniques, including random forest regression and Gini importance, to show that the Energy Exascale Earth System Model (E3SM) relies too heavily on just one of the ten chosen climatological quantities to predict September sea ice averages. Furthermore, E3SM gives too much importance to six of those quantities when compared to observed data. Finally, identifying the features that climate models incorrectly rely on should allow climatologists to improve prediction accuracy.

54 ENVIRONMENTAL SCIENCES↗