Search NASA⌕ Search

SEARCH · Search NASA

Results for “EQUATIONS - MOTION”

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 415 records · Page 23

Deep transfer operator learning for partial differential equations under conditional shift

Transfer learning enables the transfer of knowledge gained while learning to perform one task (source) to a related but different task (target), hence addressing the expense of data acquisition and labelling, potential computational power limitations and dataset distribution mismatches. Here, we propose a new transfer learning framework for task-specific learning (functional regression in partial differential equations) under conditional shift based on the deep operator network (DeepONet). Task-specific operator learning is accomplished by fine-tuning task-specific layers of the target DeepONet using a hybrid loss function that allows for the matching of individual target samples while also preserving the global properties of the conditional distribution of the target data. Inspired by conditional embedding operator theory, we minimize the statistical distance between labelled target data and the surrogate prediction on unlabelled target data by embedding conditional distributions onto a reproducing kernel Hilbert space. We demonstrate the advantages of our approach for various transfer learning scenarios involving nonlinear partial differential equations under diverse conditions due to shifts in the geometric domain and model dynamics. Our transfer learning framework enables fast and efficient learning of heterogeneous tasks despite considerable differences between the source and target domains.

42 ENGINEERING↗

Master equation study of hydrogen abstraction from HCHO by OH via a chemically activated intermediate

The abstraction reaction of hydrogen from formaldehyde by OH radical plays an important role in formaldehyde oxidation. The reaction involves a bimolecular association to form a chemically activated hydrogen-bonded reaction complex followed by a unimolecular reaction of the complex to generate the products. The reaction rate is usually considered to be pressure-independent by assuming equilibrium between the reactants and the complex. However, our nonequilibrium calculations based on the chemically significant eigenmode of the master equation, carried out with our recently developed TUMME program, indicate that the reaction complex makes the rate constant dependent on pressure at low temperatures (T < 200 K). The calculations include anharmonicity, variational effects, and multi-dimensional tunneling. We find that the reaction rate constant reaches a low-pressure limit at pressures below 10 Torr over the whole investigated temperature range (20–1800 K), which explains why the available low-temperature experiments, which are for pressures below 2 Torr, did not observe the pressure dependence. A new extension of the TUMME master-equation program is used to explore the time evolutions of the concentrations of the OH radical and the complex under pseudo-first-order conditions. The time-dependent evolution of the concentrations of the complex at a low temperature provide direct evidence for the stabilization of the reaction complex at high pressures, and it shows the negligible role of the stabilized reaction complex at low pressures. Finally, the picture that emerges is qualitatively consistent with our previous study of the reaction of methanol with OH in that the tunneling in the unimolecular step from the complex to the products affects the phenomenological reaction rate constants differently at high and low pressures and leads to a significant pressure effect.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Active Learning of Microgrid Frequency Dynamics Using Neural Ordinary Differential Equations

Accurate frequency modelling of inverter‐based resource (IBR)‐dominated power systems is crucial for ensuring stable, reliable and resilient operations, particularly given their inherent low‐inertia characteristics and fast dynamics that traditional swing equation‐based models inadequately capture. This paper explores neural ordinary differential equations (Neural ODEs) as a computationally efficient, data‐driven framework for modelling power system frequency dynamics, specifically within microgrids integrating high penetrations of distributed energy resources (DERs). The developed neural ODEs framework incorporates a neural network architecture designed to capture input dynamics. By actively perturbing the system with a known signal, the Python‐based neural ODEs framework was trained using measured system states and inputs, without the need for detailed system information. The framework, tested on a model of the Cordova, AK, microgrid, achieved a goodness of fit ranging from 60% to 99% across different state variables and maintained a mean square error in the 10 -6 p.u. range under square and step excitation signals. The proposed approach demonstrated robustness to measurement noise and initial condition variations while maintaining low computational complexity suitable for real‐time power system control applications. Furthermore, transfer learning enabled the neural ODEs model to adapt to the following changes in system topology or generator dispatch, highlighting its effectiveness for dynamic microgrids with frequently evolving configurations and diverse DERs.

Aryal, Tara [South Dakota State Univ., Brookings, ↗

Proton and Λ flow and the equation of state at high density

Results on proton and Λ flow, calculated with the UrQMD model that incorporates different realistic density dependent equations of state, are presented. It is shown that the proton and hyperon flow shows sensitivity to the equation of state and especially to the appearance of a phase transition at densities below 4n0. Even though qualitatively hyperons and protons exhibit the same beam energy dependence of the flow, the quantitative results are different. In this context it is suggested that the hyperon measurements can be used to study the density dependence of the hyperon interaction in high density QCD matter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Verification of an improved equation-free projective integration method for neoclassical plasma-profile evolution in tokamak geometry

A brute-force, long-time gyrokinetic simulation of plasma profile evolution in magnetic fusion devices is not desirable due to large computational resource requirements and a possible accumulation of numerical error. The equation-free projective integration method of Keverekidis et al. [Commun. Math. Sci. 1(4), 715–762 (2003)] is one of the outstanding candidates in projecting micro-scale simulations to a longer timescale. However, its application to tokamak plasma has not been fruitful due to the appearance of spurious transient oscillations in the lifting process, which are present when the kinetic simulations are initialized with a simplified model distribution function and which make the kinetic simulations to deviate from the desired paths. In this work, a kinetically informed lifting algorithm is added to the equation-free projective integration method, which is then verified in the electrostatic gyrokinetic particle-in-cell code XGCa [R. Hager and C. S. Chang, Phys. Plasmas 23, 042503 (2016)] for a neoclassical ion heat transport problem with adiabatic electrons. This new lifting operator is demonstrated to control spurious transients, enabling an over four-times reduction in the overall computing time in the time-evolution of the ion temperature profile in an axisymmetric toroidal plasma. Further reduction in the computing time is found to be limited due to the stability properties of the linear least squares projective integrator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonuniqueness of generalized quantum master equations for a single observable

When deriving exact generalized master equations for the evolution of a reduced set of degrees of freedom, one is free to choose what quantities are relevant by specifying projection operators. However, obtaining a reduced description does not always need to be achieved through projections-one can also use conservation laws for this purpose. Such an operation should be considered as distinct from any kind of projection; that is, projection onto a single observable yields a different form of master equation compared to that resulting from a projection followed by the application of a constraint. Here, we give a simple example to show this point and give relationships that the different memory kernels must satisfy to yield the same dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On the Cauchy problem for Boltzmann equation modeling a polyatomic gas

In the present article, we consider the Boltzmann equation that models a polyatomic gas by introducing one additional continuous variable, referred to as microscopic internal energy. We establish existence and uniqueness theory in the space homogeneous setting for the full non-linear case, under an extended Grad-type assumption on transition probability rates, which comprises hard potentials for both the relative speed and internal energy with the rate in the interval 0,2, multiplied by an integrable angular part and integrable partition functions. The Cauchy problem is resolved by means of an abstract ordinary differential equation (ODE) theory in Banach spaces for the initial data with finite and strictly positive gas mass and energy, finite momentum, and additionally finite K* polynomial moment, with K* depending on the rate of the transition probability and the structure of a polyatomic molecule or its internal degrees of freedom. Moreover, we prove that polynomially and exponentially weighted Banach space norms associated with the solution are both generated and propagated uniformly in time.

Physics↗

An analytic and complete equation of state for condensed phase materials

Analytic equations of state (EOS) are intended to reproduce theoretical and experimental data in a single phase portion of the thermodynamic space. We devise a complete and thermodynamically consistent model with four distinct features: (1) a reference isotherm that remains thermodynamically stable, (2) a flexible specific heat model based on a fourth-order rational polynomial, (3) a Grüneisen parameter that depends on specific volume and temperature, and (4) pressure and internal energy functions that can be inverted analytically in temperature. The model aims to improve the accuracy of existing equations of state while remaining computationally efficient. To demonstrate its features, we include calibrations for single-crystal pentaerythritol tetranitrate (PETN), liquid nitromethane (NM), and hexagonal close-packed beryllium (Be) metal. The parameter optimization uses the specific heat capacity, Grüneisen parameter, and static compression curves obtained from density functional theory for the crystalline solids and molecular dynamics simulations for liquid NM. We also present a velocity autocorrelation function that yields accurate phonon densities of states for the EOS calibration from the molecular dynamics trajectories. Each of the three calibrations is constrained to enforce the ambient state from experimental measurements and validated against experimental Hugoniot data from multiple sources. We also include one-dimensional hydrodynamic simulations of the isentropic compression experiments for beryllium conducted at the Z facility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Tutorial: The equations of electron emission and their evaluation

Electron emission and transport through and over potential barriers is an essential process requiring modeling and simulation to meet the design needs and characterization of an exceedingly broad range of technologically important devices and processes. The simulation and description of thermal, field, and photoemission, and the related concerns of space–charge affected electron flow, often make use of specialized formulations developed in the early days of quantum mechanics. Advancements in the utilization of electron sources and particularly the simulation of devices and applications using advanced particle-in-cell and trajectory methods for beam optics codes create a strong need for a pedagogical account of the emission models to ensure correct numerical evaluation of their equations. This Tutorial starts from simple phenomenological accounts and progressively builds to comprehensive models emphasizing straightforward and often rapid calculation. It recommends formulations to supplant the canonical Richardson–Laue–Dushman (thermal), Fowler–Nordheim (field), Fowler–DuBridge (photo), and Baroody (secondary) equations and provides a useful formulation of space–charge affected flow commonly described by the Child–Langmuir relation that takes into account cathode dependence on surface field.

Physics↗

A meshless stochastic method for Poisson–Nernst–Planck equations

A plethora of biological, physical, and chemical phenomena involve transport of charged particles (ions). Its continuum-scale description relies on the Poisson–Nernst–Planck (PNP) system, which encapsulates the conservation of mass and charge. The numerical solution of these coupled partial differential equations is challenging and suffers from both the curse of dimensionality and difficulty in efficiently parallelizing. We present a novel particle-based framework to solve the full PNP system by simulating a drift–diffusion process with time- and space-varying drift. We leverage Green’s functions, kernel-independent fast multipole methods, and kernel density estimation to solve the PNP system in a meshless manner, capable of handling discontinuous initial states. The method is embarrassingly parallel, and the computational cost scales linearly with the number of particles and dimension. We use a series of numerical experiments to demonstrate both the method’s convergence with respect to the number of particles and computational cost vis-à-vis a traditional partial differential equation solver.

Chemistry↗

Radiation drive designed to extend the pressure ranges measured in Gbar equation of state experiments at the National Ignition Facility

We present the design and demonstration of a Shock-Strengthening hohlraum radiation temperature drive in the Gbar experimental platform at the National Ignition Facility intended to increase the pressure range measured in a single experiment. Previously published experiments by Döppner et al. measured the equation of state in polystyrene from 25 to 60 Mbar. Recent experimental data of the Shock-Strengthening drive initially demonstrated a much larger pressure range from 15 to 110 Mbar using the same peak radiation temperature and experimental platform. The Shock-Strengthening drive starts with a low temperature foot that launches a weak shock into the sample and is followed by a continuous increase in radiation temperature to strengthen the leading shock. The additional strengthening increases the pressure within the sample beyond what is achievable by convergence alone. Design features of the Shock-Strengthening drive and accompanying radiation hydrodynamics simulations are used to illustrate the method by which the pressure range is increased from previous experiments. This method of modifying the radiation temperature drive can be used on the Gbar platform to significantly increase the range for equation of state data collected in a single experiment for many materials.

Physics - Plasma physics↗

On “incorrectnes” of the rate equation for photon density of semiconductor lasers

The authors of works [1, 2] published in Quantum Electronics assert that the rate equation for photon density with a term taking into account the contribution of spontaneous emission to the laser mode is erroneous (i.e., the spontaneous term should be excluded from consideration) and criticise works [3, 4]. In the present work, we analyse previously published papers and show that the necessity of taking into account the contribution of spontaneous emission to the laser mode follows from the quantum-mechanical rate equation for semiconductor lasers. The results obtained in [1, 2] are analysed. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Neural network reconstruction of the dense matter equation of state from neutron star observables

The Equation of State (EoS) of strongly interacting cold and hot ultra-dense QCD matter remains a major challenge in the field of nuclear astrophysics. With the advancements in measurements of neutron star masses, radii, and tidal deformabilities, from electromagnetic and gravitational wave observations, neutron stars play an important role in constraining the ultra-dense QCD matter EoS. Here, in this work, we present a novel method that exploits deep learning techniques to reconstruct the neutron star EoS from mass-radius (M-R) observations. We employ neural networks (NNs) to represent the EoS in a model-independent way, within the range of ~1-7 times the nuclear saturation density. The unsupervised Automatic Differentiation (AD) framework is implemented to optimize the EoS, so as to yield through TOV equations, an M-R curve that best fits the observations. We demonstrate that this method works by rebuilding the EoS on mock data, i.e., mass-radius pairs derived from a randomly generated polytropic EoS. The reconstructed EoS fits the mock data with reasonable accuracy, using just 11 mock M-R pairs observations, close to the current number of actual observations.

79 ASTRONOMY AND ASTROPHYSICS↗

A topological way of finding solutions to the Yang–Mills equation

We propose a systematic way of finding solutions to the classical Yang–Mills equation with nontrivial topology. This approach is based on one of the Wightman axioms for quantum field theory, which is referred to as the form invariance condition in this paper. For a given gauge group and a spacetime with certain isometries, thanks to this axiom that imposes strong constraints on the general ansatz, a systematic way of solving the Yang–Mills equation can be obtained in both flat and curved spacetimes. In order to demonstrate this method, we recover various known solutions as special cases, as well as producing new solutions not previously reported in the literature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Parametrically driven nonlinear Dirac equation with arbitrary nonlinearity

Here, the damped and parametrically driven nonlinear Dirac equation with arbitrary nonlinearity parameter $\kappa$ is analyzed, when the external force is periodic in space and given by $\mathcal{f(x) = r}$ cos(Kx), both numerically and in a variational approximation using five collective coordinates (time dependent shape parameters of the wave function). Our variational approximation satisfies exactly the low-order moment equations. Because of competition between the spatial period of the external force λ = 2π/K , and the soliton width $l_s$, which is a function of the nonlinearity $\kappa$ as well as the initial frequency ω 0 of the solitary wave, there is a transition (at fixed ω 0 ) from trapped to unbound behavior of the soliton, which depends on the parameters $\mathcal{r}$ and K of the external force and the nonlinearity parameter $\kappa$. We previously studied this phenomena when $\kappa$ = 1 where we showed that for λ $\gg l_s$ the soliton oscillates in an effective potential, while forλ $\ll l_s$ it moves uniformly as a free particle. In this paper we focus on the $\kappa$ dependence of the transition from oscillatory to particle behavior and explicitly compare the curves of the transition regime found in the collective coordinate approximation as a function of $\mathcal{r}$ and K when $\kappa$ = 1/2,1,2 at fixed value of the frequency ω 0 . Since the solitary wave gets narrower for fixed ω 0 as a function of $\kappa$, we expect and indeed find that the regime where the solitary wave is trapped is extended as we increase $\kappa$.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Towards dynamical low-rank approximation for neutrino kinetic equations. Part I: Analysis of an idealized relaxation model

Dynamical low-rank approximation (DLRA) is an emerging tool for reducing computational costs and provides memory savings when solving high-dimensional problems. Here, in this work, we propose and analyze a semi-implicit dynamical low-rank discontinuous Galerkin (DLR-DG) method for the space homogeneous kinetic equation with a relaxation operator, modeling the emission and absorption of particles by a background medium. Both DLRA and the discontinuous Galerkin (DG) scheme can be formulated as Galerkin equations. To ensure their consistency, a weighted DLRA is introduced so that the resulting DLR-DG solution is a solution to the fully discrete DG scheme in a subspace of the standard DG solution space. Similar to the standard DG method, we show that the proposed DLR-DG method is well-posed. We also identify conditions such that the DLR-DG solution converges to the equilibrium. Numerical results are presented to demonstrate the theoretical findings.

97 MATHEMATICS AND COMPUTING↗

Numerical methods and hypoexponential approximations for gamma distributed delay differential equations

Abstract Gamma distributed delay differential equations (DDEs) arise naturally in many modelling applications. However, appropriate numerical methods for generic gamma distributed DDEs have not previously been implemented. Modellers have therefore resorted to approximating the gamma distribution with an Erlang distribution and using the linear chain technique to derive an equivalent system of ordinary differential equations (ODEs). In this work, we address the lack of appropriate numerical tools for gamma distributed DDEs in two ways. First, we develop a functional continuous Runge–Kutta (FCRK) method to numerically integrate the gamma distributed DDE without resorting to Erlang approximation. We prove the fourth-order convergence of the FCRK method and perform numerical tests to demonstrate the accuracy of the new numerical method. Nevertheless, FCRK methods for infinite delay DDEs are not widely available in existing scientific software packages. As an alternative approach to solving gamma distributed DDEs, we also derive a hypoexponential approximation of the gamma distributed DDE. This hypoexponential approach is a more accurate approximation of the true gamma distributed DDE than the common Erlang approximation but, like the Erlang approximation, can be formulated as a system of ODEs and solved numerically using standard ODE software. Using our FCRK method to provide reference solutions, we show that the common Erlang approximation may produce solutions that are qualitatively different from the underlying gamma distributed DDE. However, the proposed hypoexponential approximations do not have this limitation. Finally, we apply our hypoexponential approximations to perform statistical inference on synthetic epidemiological data to illustrate the utility of the hypoexponential approximation.

97 MATHEMATICS AND COMPUTING↗

Convex relaxation for Fokker–Planck equation

We propose an approach to directly estimate the moments or marginals for a high-dimensional equilibrium distribution in statistical mechanics by solving the high-dimensional Fokker–Planck equation in terms of low-order cluster moments or marginals. With this approach, we bypass the exponential complexity of estimating the full high-dimensional distribution and directly solve the simplified partial differential equations for low-order moments/marginals. Moreover, the proposed moment/marginal relaxation is fully convex and can be solved via off-the-shelf solvers. We further propose a time-dependent version of the convex programs to study non-equilibrium dynamics. In a specific setting, we show the proposed method can recover a mean-field-type equilibrium density. Numerical results are provided to demonstrate the performance of the proposed algorithm for high-dimensional systems.

Chen, Yian↗