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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

The use of centrifugation to study early Drosophila embryogenesis

By the end of 10th nuclear cycle, the somatic nuclei of the Drosophila embryo have migrated to the periphery of the egg. Centrifugation of embryos did not result in the displacement of these nuclei, since cytoskeletal elements anchor them to the cortex. But, mild centrifugal forces displace the centrally located, nascent yolk nuclei. If this increased sensitivity to hypergravity occurs before the beginning of nuclear differentiation during cycle 8, when the nascent yolk and somatic nuclei physically separate, then it would mark the earliest functional difference between these two lineages.

Non-NASA Center↗
NASA NTRS1992

Extracellular matrix as a solid-state regulator in angiogenesis: identification of new targets for anti-cancer therapy

Angiogenesis, the growth of blood capillaries, is regulated by soluble growth factors and insoluble extracellular matrix (ECM) molecules. Soluble angiogenic mitogens act over large distances to initiate capillary growth whereas changes in ECM govern whether individual cells will grow, differentiate, or involute in response to these stimuli in the local tissue microenvironment. Analysis of this local control mechanism has revealed that ECM molecules switch capillary endothelial cells between differentiation and growth by both binding specific transmembrane integrin receptors and physically resisting cell-generated mechanical loads that are applied to these receptors. Control of capillary endothelial cell form and function therefore may be exerted by altering the mechanical properties of the ECM as well as its chemical composition. Understanding of this mechanochemical control mechanism has led to the development of new angiogenesis inhibitors that may be useful for the treatment of cancer.

Review↗
DOE OSTI2020

Quantum Physics without the Physics

This report explains the basic theory and common terminology of quantum physics without assuming any knowledge of physics. It was written by a group of applied mathematicians while they were reading up on the subject. The intended audience consists of applied mathematicians, computer scientists, or anyone else who wants to improve their understanding of quantum physics. We assume that the reader is familiar with fundamental concepts of linear algebra, differential equations, and to some extent the theory of Hilbert spaces. Most of the material can be found in the book by Nielsen and Chuang and in the lecture notes on open quantum systems by Lidar. Another excellent online source of information is Wikipedia, even though most of its articles on quantum physics assume a solid understanding of physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Single-stage gradient-based stellarator coil design: Optimization for near-axis quasi-symmetry

Here we present a new coil design paradigm for magnetic confinement in stellarators. Our approach directly optimizes coil shapes and coil currents to produce a vacuum quasi-symmetric magnetic field with a target rotational transform on the magnetic axis. This approach differs from the traditional two-stage approach in which first a magnetic configuration with desirable physics properties is found, and then coils to approximately realize this magnetic configuration are designed. The proposed single-stage approach allows us to find a compromise between confinement and engineering requirements, i.e., find easy-to-build coils with good confinement properties. Using forward and adjoint sensitivities, we derive derivatives of the physical quantities in the objective, which is constrained by a nonlinear periodic differential equation. In two numerical examples, we compare different gradient-based descent algorithms and find that incorporating approximate second-order derivative information through a quasi-Newton method is crucial for convergence. We also explore the optimization landscape in the neighborhood of a minimizer and find many directions in which the objective is mostly flat, indicating ample freedom to find simple and thus easy-to-build coils.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Physics-aware adaptive checkpointing with shadow systems for nonlinear PDE simulations

Large-scale simulations of nonlinear partial differential equations (PDEs) that exhibit strongly transient behavior and pattern-forming dynamics produce enormous amounts of data, which, even with modern storage systems, cannot be stored for later curation. Current I/O strategies either write dense time series of snapshots, which is often prohibitive in I/O and storage, or store a few checkpoints that enable restart but incur expensive recomputation cost and provide no control over post-restart error growth, especially when lossy compression is used. Moreover, most, if not all, existing strategies take no account of the actual physical state of the system. Here, we present a simple physics-aware I/O framework in which a low-cost shadow system adaptively triggers lossy checkpoints when the shadow system deviates from the fine-scale simulation. The shadow system can be a coarsened replica of the fine-scale simulation that evolves concurrently. This means that checkpoints are taken based on the physical state of the system: fewer checkpoints are triggered when the system is quiescent while more are taken when the system undergoes a rapid change. This type of behavior is observed in many systems such as Brusselator and FitzHugh–Nagumo. We illustrate that our framework maintains stable restarts, keeps fine-scale restart errors bounded by shadow errors, and reconstructs the time history with significantly lower error and storage than interpolating fixed-interval snapshots, with low-cost shadow replay and modest online synchronization overhead.

Gong, Qian [ORNL] (ORCID:0000000235704142)↗
DOE OSTI2023

Loop-by-loop differential equations for dual (elliptic) Feynman integrals

We present a loop-by-loop method for computing the differential equations of Feynman integrals using the recently developed dual form formalism. We give explicit prescriptions for the loop-by-loop fibration of multi-loop dual forms. Then, we test our formalism on a simple, but non-trivial, example: the two-loop three-mass elliptic sunrise family of integrals. We obtain an ε-form differential equation within the correct function space in a sequence of relatively simple algebraic steps. In particular, none of these steps relies on the analysis of q-series. Then, we discuss interesting properties satisfied by our dual basis as well as its simple relation to the known ε-form basis of Feynman integrands. The underlying K3-geometry of the three-loop four-mass sunrise integral is also discussed. Finally, we speculate on how to construct a “good” loop-by-loop basis at three-loop.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

SPIKANs: separable physics-informed Kolmogorov–Arnold networks

Physics-Informed Neural Networks (PINNs) have emerged as a promising method for solving partial differential equations (PDEs) in scientific computing. While PINNs typically use multilayer perceptrons (MLPs) as their underlying architecture, recent advancements have explored alternative neural network structures. One such innovation is the Kolmogorov–Arnold Network (KAN), which has demonstrated benefits over traditional MLPs, including faster neural scaling and better interpretability. The application of KANs to physics-informed learning has led to the development of Physics-Informed KANs (PIKANs), enabling the use of KANs to solve PDEs. However, despite their advantages, KANs often suffer from slower training speeds, particularly in higher-dimensional problems where the number of collocation points grows exponentially with the dimensionality of the system. To address this challenge, we introduce Separable Physics-Informed Kolmogorov–Arnold Networks (SPIKANs). This novel architecture applies the principle of separation of variables to PIKANs, decomposing the problem such that each dimension is handled by an individual KAN. This approach drastically reduces the computational complexity of training without sacrificing accuracy, facilitating their application to higher-dimensional PDEs. Through a series of benchmark problems, we demonstrate the effectiveness of SPIKANs, showcasing their superior scalability and performance compared to PIKANs and highlighting their potential for solving complex, high-dimensional PDEs in scientific computing.

Kolmogorov-Arnold networks↗
DOE OSTI2025

Hidden zeros of the cosmological wavefunction

Motivated by the recent discovery of hidden zeros in particle and string amplitudes, we characterize zeros of individual graph contributions to the cosmological wavefunction of a scalar field theory. We demonstrate that these contributions factorize near these zeros for all tree graphs and provide evidence that this extends to loop graphs as well. We explicitly construct polytopal realizations of the relevant graph associahedra and show that the cosmological zeros have natural geometric and physical interpretations. As a byproduct, we establish an equivalence between the wavefunction coefficients of chain graphs and flat-space Tr(ϕ 3 ) amplitudes, enabling us to leverage the cosmological zeros to uncover the recently discovered hidden zeros of colored amplitudes.

Cosmological models↗
NASA NTRS2013

Extended Analytic Device Optimization Employing Asymptotic Expansion

Analytic optimization of a thermoelectric junction often introduces several simplifying assumptionsincluding constant material properties, fixed known hot and cold shoe temperatures, and thermallyinsulated leg sides. In fact all of these simplifications will have an effect on device performance,ranging from negligible to significant depending on conditions. Numerical methods, such as FiniteElement Analysis or iterative techniques, are often used to perform more detailed analysis andaccount for these simplifications. While numerical methods may stand as a suitable solution scheme,they are weak in gaining physical understanding and only serve to optimize through iterativesearching techniques. Analytic and asymptotic expansion techniques can be used to solve thegoverning system of thermoelectric differential equations with fewer or less severe assumptionsthan the classic case. Analytic methods can provide meaningful closed form solutions and generatebetter physical understanding of the conditions for when simplifying assumptions may be valid.In obtaining the analytic solutions a set of dimensionless parameters, which characterize allthermoelectric couples, is formulated and provide the limiting cases for validating assumptions.Presentation includes optimization of both classic rectangular couples as well as practically andtheoretically interesting cylindrical couples using optimization parameters physically meaningful toa cylindrical couple. Solutions incorporate the physical behavior for i) thermal resistance of hot andcold shoes, ii) variable material properties with temperature, and iii) lateral heat transfer through legsides.

Thermoelectric Generators↗
DOE OSTI2021

Decomposition of Feynman integrals by multivariate intersection numbers

We present a detailed description of the recent idea for a direct decomposition of Feynman integrals onto a basis of master integrals by projections, as well as a direct derivation of the differential equations satisfied by the master integrals, employing multivariate intersection numbers. We discuss a recursive algorithm for the computation of multivariate intersection numbers, and provide three different approaches for a direct decomposition of Feynman integrals, which we dub the straight decomposition, the bottom-up decomposition, and the top-down decomposition. These algorithms exploit the unitarity structure of Feynman integrals by computing intersection numbers supported on cuts, in various orders, thus showing the synthesis of the intersection-theory concepts with unitarity-based methods and integrand decomposition. We perform explicit computations to exemplify all of these approaches applied to Feynman integrals, paving a way towards potential applications to generic multi-loop integrals.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Simultaneous Momentum and Position Measurement and the Instrumental Weyl-Heisenberg Group

The canonical commutation relation, [Q,P]=iℏ, stands at the foundation of quantum theory and the original Hilbert space. The interpretation of P and Q as observables has always relied on the analogies that exist between the unitary transformations of Hilbert space and the canonical (also known as contact) transformations of classical phase space. Now that the theory of quantum measurement is essentially complete (this took a while), it is possible to revisit the canonical commutation relation in a way that sets the foundation of quantum theory not on unitary transformations but on positive transformations. This paper shows how the concept of simultaneous measurement leads to a fundamental differential geometric problem whose solution shows us the following. The simultaneous P and Q measurement (SPQM) defines a universal measuring instrument, which takes the shape of a seven-dimensional manifold, a universal covering group we call the instrumental Weyl-Heisenberg (IWH) group. The group IWH connects the identity to classical phase space in unexpected ways that are significant enough that the positive-operator-valued measure (POVM) offers a complete alternative to energy quantization. Five of the dimensions define processes that can be easily recognized and understood. The other two dimensions, the normalization and phase in the center of the IWH group, are less familiar. The normalization, in particular, requires special handling in order to describe and understand the SPQM instrument.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

An attention-based neural ordinary differential equation framework for modeling inelastic processes

To preserve strictly conservative behavior as well as model the variety of dissipative behavior displayed by solid materials, we propose a significant enhancement to the internal state variable-neural ordinary differential equation (ISV-NODE) framework. In this data-driven, physics-constrained modeling framework internal states are inferred rather than prescribed. The ISV-NODE consists of: (a) a stress model dependent on observable deformation and inferred internal state, and (b) a model of the evolution of the internal states. The enhancements to ISV-NODE proposed in this work are multifold: (a) a partially input convex neural network stress potential provides polyconvexity in terms of observed strain while leaving the inferred state unconstrained, and (b) an internal state flow model uses common latent features to inform novel attention-based gating and drives the flow of internal state only in dissipative regimes. We demonstrated that this architecture can accurately model dissipative and conservative behavior across an isotropic, isothermal elastic-viscoelastic-elastoplastic spectrum with three exemplars, while maintaining fundamental principles by design.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Development of a Single-Phase, Transient, Subchannel Code, within the MOOSE Multi-Physics Computational Framework

Subchannel codes have been widely used for thermal-hydraulics analyses in nuclear reactors. This paper details the development of a novel subchannel code within the Idaho National Laboratory’s (INL) Multi-physics Object Oriented Simulation Environment (MOOSE). MOOSE is a parallel computational framework targeted at the solution of systems of coupled, nonlinear partial differential equations, that often arise in the simulation of nuclear processes. As such, it includes codes/modules able to solve the multiple linear and nonlinear physics that describe a nuclear reactor, under normal operation conditions or accidents. This includes thermal-hydraulics, fuel performance, and neutronics codes, between others. A MOOSE-based subchannel code is a new addition to the fleet of INL-developed codes, based on the MOOSE framework. In this work, we present the derivation of the subchannel equations for a single-phase fluid, we proceed with the description of the algorithm that is used to solve these equations and describe how this algorithm was implemented within MOOSE. We also present how this code can be coupled to the BISON fuel performance code. Next, we verify the friction model and the turbulent mixing model. We calibrate the turbulent modeling parameters for momentum mixing and enthalpy mixing, C T , β. We validate the code using experimental results and last demonstrate the coupling capabilities using a simple example.

42 ENGINEERING↗
NASA NTRS1985

The cosmic ray differential diurnal variation dependences on the zenith angle and the geomagnetic disturbance

Simultaneous and continuous muon measurements in two opposite azimuthal directions under equal zenith angles demonstrated the importance of this method for cosmic ray diurnal variation investigations. Lately these measurements were extended by means of improved telescopes. The obtained cosmic ray diurnal variations were presented as intensity differential curves. Theoretical investigations connected the properties of these curves with some interplanetary spece parameters. The harmonics of these curves were interpreted physically. Some order difference curves were introduced. In earlier works some dependences between the parameters characterizing the first and the second harmonics of the differential intensity curves and the geomagnetic activity were found. Then all measurements were carried out under only one zenith angle. The results of investigations of similar dependences using data of simultaneous measurements under three different zenith angles are presented.

Kavlakov, S.↗
DOE OSTI2025

Preliminary study of auto-differentiation algorithm in beam dynamics with stochastic process

Modern particle accelerator optimization requires sophisticated computational methods to address the inherently stochastic nature of beam dynamics. This research develops a framework applying AD to SDEs that specifically addresses beam dynamics challenges in particle accelerators, focusing on accurately modeling and optimizing beam behavior in regimes dominated by stochastic processes. By incorporating key physical phenomena such as synchrotron radiation, wakefield effects, and quantum excitation, the framework aims to provide auto differentiation on the figure of merit of the phase space evolution and beam dynamics. The methodology will enable effective optimization method in a dynamic system with stochastic process.

Accelerator Physics↗
DOE OSTI2022

A High-Throughput MEMS-Based Differential Scanning Calorimeter for Direct Thermal Characterization of Antibodies

Calorimeters, which can be used for rapid thermal characterization of biomolecules, are getting intense attention in drug development. This paper presents a novel MEMS-based differential scanning calorimeter (DSC) for direct thermal characterization of protein samples. The DSC consisted of a pair of temperature sensors made by vanadium oxide (VOx) film with a temperature coefficient of resistivity of –0.025/K at 300 K, a microfluidic device with high thermal insulation (2.8 K/mW), and a Peltier heater for linear temperature scanning. The DSC exhibited high sensitivity (6.1 µV/µW), low noise (0.4 µW), high scanning rate (45 K/min), and low sample consumption volume (0.63 µL). The MEMS DSC was verified by measuring the temperature-induced denaturation of lysozyme at different pH, and then used to study the thermal stability of a monoclonal antibody (mAb), an antigen-binding fragment (Fab), and a dual variable domain immunoglobulin (DVD-Ig) at pH = 6. The results showed that lysozyme is a stable protein in the pH range of 4.0–8.0. The protein stability study revealed that the transition temperatures of the intact Fab fragment, mAb, and DVD proteins were comparable with conformational stability results obtained using conventional commercial DSC. These studies demonstrated that the MEMS DSC is an effective tool for directly understanding the thermal stability of antibodies in a high-throughput and low-cost manner compared to conventional calorimeters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Designing observables for measurements with deep learning

Many analyses in particle and nuclear physics use simulations to infer fundamental, effective, or phenomenological parameters of the underlying physics models. When the inference is performed with unfolded cross sections, the observables are designed using physics intuition and heuristics. We propose to design targeted observables with machine learning. Unfolded, differential cross sections in a neural network output contain the most information about parameters of interest and can be well-measured by construction. The networks are trained using a custom loss function that rewards outputs that are sensitive to the parameter(s) of interest while simultaneously penalizing outputs that are different between particle-level and detector-level (to minimize detector distortions). We demonstrate this idea in simulation using two physics models for inclusive measurements in deep inelastic scattering. We find that the new approach is more sensitive than classical observables at distinguishing the two models and also has a reduced unfolding uncertainty due to the reduced detector distortions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Time integrator agnostic charge conserving finite element PIC

Developing particle-in-cell (PIC) methods using finite element basis sets, and without auxiliary divergence cleaning methods, was a longstanding problem until recently. It was shown that if consistent spatial basis functions are used, one can indeed create a methodology that was charge conserving, albeit using a leapfrog time stepping method. While this is a significant advance, leapfrog schemes are only conditionally stable and time step sizes are closely tied to the underlying mesh. Ideally, to take full advantage of advances in finite element methods (FEMs), one needs a charge conserving PIC methodology that is agnostic to the time stepping method. This is the principal contribution of this paper. In what follows, we shall develop this methodology, prove that both charge and Gauss’ laws are discretely satisfied at every time step, provide the necessary details to implement this methodology for both the wave equation FEM and Maxwell solver FEM, and finally demonstrate its efficacy on a suite of test problems. The method will be demonstrated by single particle evolution, non-neutral beams with space-charge, and adiabatic expansion of a neutral plasma, where the Debye length has been resolved, and real mass ratios are used.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗