Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Equation”

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 19 records

Symmetry Determining Equations of the Rankine-Hugoniot Equations for Variable Velocity Shock Waves

The “constant velocity piston” problem (Fig. 1), also known as the “piston problem,” is a standard model for a one dimensional, in our case linear, symmetric shock wave moving through an inviscid, perfect gas. The model can be divided into two regions - a perturbed section on the left and an unperturbed section on the right - by a moving shock wave moving left to right. Both the perturbed and unperturbed sections, i.e. the shocked and unshocked regions, respectively, obey the Eulerian conservation equations; however, at the exact location of the shock, there is a mathematical discontinuity not satisfied by the Euler equations. To ensure continuity and conservation of certain quantities when crossing between the unshocked and shocked regions, we evoke a series of equations derived from the Eulerian conservation equations, called the Rankine-Hugoniot equations, or “jump” equations as it is often referred to in the literature on the topic. The classical constant-velocity piston problem assumes the piston features a constant driving velocity (among many other willing suspensions of belief required in the pursuit of a first principles equation model); consequent to this assumption is a constant-velocity shock and a constant-velocity shocked flow state. However, using Lie Group Theory (LGT), also known as symmetry analysis, we can attempt to reinterpret the model with a shock wave of variable velocity in time and space. An extension of the model in this way opens up the possibility for obtaining new analytical solutions to the piston problem for certain shock velocity models. In this report, we use LGT to derive the symmetry determining equations (SDEs), whose solutions are Lie groups, which permit analytical solutions. In the future, we can then use the SDEs to define constraint equations on the shock velocity model and what the successive solutions to the Euler equations might be based off such constraints. This report is structured as follows: Section 2 provides a brief derivation of the Rankine-Hugoniot (“jump”) equations; Section 3 gives an overview of Lie group theory; Section 4 derives the SDEs of the jump equations; Section 5 derives the Euler conservation equations for fluids; and Section 6 presents concluding remarks and opportunities for future studies.

42 ENGINEERING

Black Box Equations of State: Creating Semi-analytic Solutions to the Noh Problem and Verifying Equation of State Interfaces

The objective of this report is threefold. First, it details a method for deriving a semi-analytic solution to the Noh Problem when using a “black-box” equation of state. Such capability allows us to perform verification on complicated, more realistic equations of state. Examples include Steinberg equations of state for materials and tabulated equations of state. The second objective is to apply the methodology to verify the singularity-eos equation of state library. We do so by solving the Rankine-Hugoinot jump conditions for the Noh Problem, ensuring singularity derives the correct solution and comparing the error to an exact implementation of the equation of state. The third objective is to perform verification of the xRAGE Eulerian hydrodynamics code when interfaced with singularity. We provide the theory, analysis, documentation for a python implementation of the proposed solver, and verification results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC

Symmetry Determining Equations of the Euler Equations

This purpose of this document is to provide a clear derivation of the symmetry determining equations (SDEs) of the Euler equations to supplement the lecture document. The Euler equations are a coupled system of three first-order quasilinear hyperbolic partial differential equations (PDEs). The Euler equations are a system of partial differential equations (PDEs) for fluid velocity v(t, r), fluid density ρ(t, r), and fluid pressure P(t, r), where t is time and r is the spatial coordinate.

97 MATHEMATICS AND COMPUTING

BPZ equations for higher degenerate fields and nonperturbative Dyson-Schwinger equations

In the two-dimensional Liouville conformal field theory, correlation functions involving a degenerate field satisfy partial differential equations due to the decoupling of the null descendant field. On the other hand, the instanton partition function of a four-dimensional N = 2 supersymmetric theory in the Ω -background at a special point of the parameter space also satisfies a partial differential equation resulting from the constraints of the gauge field configurations. This partial differential equation can be proved using the nonperturbative Dyson-Schwinger equations. We show for the next-to-simplest case that the partial differential equations obtained from two different perspectives can be identified, thereby confirming an assertion of the Bogomol’nyi-Prasad-Sommerfield/conformal field theory correspondence. Published by the American Physical Society 2024

Astronomy & Astrophysics

High-order limiting methods using maximum principle bounds derived from the Boltzmann equation I: Euler equations

The use of limiting methods for high-order numerical approximations of hyperbolic conservation laws generally requires defining an admissible region/bounds for the solution. In this work, we present a novel approach for computing solution bounds and limiting for the Euler equations through the kinetic representation provided by the Boltzmann equation, which allows for extending limiters designed for linear advection directly to the Euler equations. Given an arbitrary set of solution values to compute bounds over (e.g., numerical stencil) and a desired linear advection limiter, the proposed approach yields an analytic expression for the admissible region of particle distribution function values, which may be numerically integrated to yield a set of bounds for the density, momentum, and total energy. Further, these solution bounds are shown to preserve positivity of density/pressure/internal energy and, when paired with a limiting technique, can robustly resolve strong discontinuities while recovering high-order accuracy in smooth regions without any ad hoc corrections (e.g., relaxing the bounds). This approach is demonstrated in the context of an explicit unstructured high-order discontinuous Galerkin/flux reconstruction scheme for a variety of difficult problems in gas dynamics, including cases with extreme shocks and shock-vortex interactions. Furthermore, this work presents a foundation for limiting techniques for more complex macroscopic governing equations that can be derived from an underlying kinetic representation for which admissible solution bounds are not well-understood.

42 ENGINEERING

Energy-momentum-conserving stochastic differential equations and algorithms for the nonlinear Landau-Fokker-Planck equation

Coulomb collision is a fundamental diffusion process in plasmas that can be described by the Landau-Fokker-Planck (LFP) equation or the stochastic differential equation (SDE). While energy and momentum are conserved exactly in the LFP equation, they are conserved only on average by the conventional corresponding SDEs, suggesting that the underlying stochastic process may not be well defined by such SDEs. Here, in this study, we derive new SDEs with exact energy-momentum conservation for the Coulomb collision by factorizing the collective effect of field particles into individual particles and enforcing Newton's third law. These SDEs, when interpreted in the Stratonovich sense, have a particularly simple form that represents pure diffusion between particles without drag. To demonstrate that the new SDEs correspond to the LFP equation, we develop numerical algorithms that converge to the SDEs and preserve discrete conservation laws. Simulation results are presented in a benchmark of various relaxation processes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Towards a Quantum Algorithm for the Incompressible Nonlinear Navier-Stokes Equations

In this work, we present novel concepts for quantum algorithms to solve transient, nonlinear partial differential equations (PDEs). The challenge lies in how to effectively represent, encode, process, and evolve the nonlinear system of PDEs on quantum computers. We will discuss the new techniques using the incompressible Navier-Stokes equations as an example, because it represents the fundamental nonlinear feature and yet removes certain complexity in physics, allowing us to focus on the design of quantum algorithms. Previous attempts solving nonlinear PDEs in quantum computation have often involved storing multiple copies of solutions or employing linearizations. Neither is practical due to exponential scaling with evolution time or insufficient solution accuracy. We propose a new framework based on matrix product states (MPSs) and matrix product operators (MPOs), in addition to the Krylov subspace methods. For example, the solution variables of the Navier-Stokes equations are represented by MPSs, and the linear and nonlinear terms are processed by MPOs. The time evolution of the operators is attained by a fast-forwarding algorithm using Krylov subspace methods. Furthermore, we discuss various techniques for efficient encoding of MPSs, measurement reduction for MPOs, and use of tensor operations to treat multi-variate, multi-physics characteristics of Navier-Stokes.

Gopalakrishnan Meena, Murali [ORNL] (ORCID:0000000

Novel scalings of neutron star properties from analyzing dimensionless Tolman–Oppenheimer–Volkoff equations

The Tolman–Oppenheimer–Volkoff (TOV) equations govern the radial evolution of pressure and energy density in static neutron stars (NSs) in hydrodynamical equilibrium. Using the reduced pressure and energy density with respect to the NS central energy density, the original TOV equations can be recast into dimensionless forms. While the traditionally used integral approach for solving the original TOV equations require an input nuclear Equation of State (EOS), the dimensionless TOV equations can be anatomized by using the reduced pressure and energy density as polynomials of the reduced radial coordinate without using any input nuclear EOS. It has been shown in several of our recent works that interesting and novel perspectives about NS core EOS can be extracted directly from NS observables by using the latter approach. Our approach is based on intrinsic and perturbative analyses of the dimensionless (IPAD) TOV equations (IPAD-TOV). In this review article, we first discuss the length and energy density scales of NSs as well as the dimensionless TOV equations for scaled variables and their perturbative solutions near NS cores. We then review several new insights into NS physics gained from solving perturbatively the scaled TOV equations. Whenever appropriate, comparisons with the traditional approach from solving the original TOV equations will be made. In particular, we first show that the nonlinearity of the TOV equations basically excludes a linear EOS for dense matter in NS cores. We then show that perturbative analyses of the scaled TOV equations enable us to reveal novel scalings of the NS mass, radius and the compactness with certain combinations of the NS central pressure and energy density. Thus, observational data on either mass, radius or compactness can be used to constrain directly the core EOS of NS matter independent of the still very uncertain nuclear EOS models. As examples, the EOS of the densest visible matter in our Universe before the most massive neutron stars collapse into black holes (BHs) as well as the central EOS of a canonical or a 2.1 solar mass NS are extracted without using any nuclear EOS model. In addition, we show that causality in NSs sets an upper bound of about 0.374 for the ratio of pressure over energy density and correspondingly a lower limit for trace anomaly in supra-dense matter. We also demonstrate that the strong-field gravity plays a fundamental role in extruding a peak in the density/radius profile of the speed of sound squared (SSS) in massive NS cores independent of the nuclear EOS. Finally, some future perspectives of NS research using the new approach reviewed here by solving perturbatively the dimensionless TOV equations are outlined.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS

Gravitational waves from binary neutron star mergers with a spectral equation of state

In numerical simulations of binary neutron star systems, the equation of state of the dense neutron star matter is an important factor in determining both the physical realism and the numerical accuracy of the simulations. Some equations of state used in simulations are C 2 or smoother in the pressure/density relationship function, such as a polytropic equation of state, but may not have the flexibility to model stars or remnants of different masses while keeping their radii within known astrophysical constraints. Other equations of state, such as tabular or piece-wise polytropic, may be flexible enough to model additional physics and multiple stars' masses and radii within known constraints, but are not as smooth, resulting in additional numerical error. We will study in this paper a recently developed family of equation of state, using a spectral expansion with sufficient free parameters to allow for a larger flexibility than current polytropic equations of state, and with sufficient smoothness to reduce numerical errors compared to tabulated or piece-wise polytropic equations of state. We perform simulations at three mass ratios with a common chirp mass, using two distinct spectral equations of state, and at multiple numerical resolutions. We evaluate the gravitational waves produced from these simulations, comparing the phase error between resolutions and equations of state, as well as with respect to analytical models. From our simulations we estimate that the phase difference at merger for binaries with a dimensionless weighted tidal deformability difference greater than Δ$\tilde{Λ}$=55 can be captured by the SpEC code for these equations of state.

79 ASTRONOMY AND ASTROPHYSICS

Score-Based Physics-Informed Neural Networks for High-Dimensional Fokker–Planck Equations

The Fokker-Planck (FP) equation is a foundational partial differential equation (PDE) in stochastic processes involving Brownian motions. However, the curse of dimensionality (CoD) poses a formidable challenge when dealing with high-dimensional FP equations. Although Monte Carlo simulation and (vanilla) Physics-Informed Neural Networks (PINNs) have shown the potential to tackle CoD, both methods exhibit significant numerical errors in high dimensions when dealing with the probability density function (PDF) associated with Brownian motion. The point-wise PDF values tend to decrease exponentially as dimensionality increases, surpassing the precision of numerical simulations and resulting in substantial errors. In addition, due to its massive sampling, Monte Carlo fails to offer fast sampling. Modeling the logarithm likelihood (LL) via vanilla PINNs transforms the FP equation into a notoriously difficult Hamilton-Jacobi-Bellman (HJB) equation, which is impractical for PINN learning, whose error grows rapidly with dimension. To this end, we propose a novel approach utilizing a score-based solver to fit the score function in stochastic differential equations (SDEs). The score function, defined as the gradient of the LL, plays a fundamental role in inferring LL and PDF and enables fast SDE sampling, offering an effective means to overcome the CoD. Three fitting methods, Score Matching (SM), Sliced Score Matching (SSM), and Score-PINN, are introduced, each contributing unique advantages in computational complexity, accuracy, and generality. The proposed score-based SDE solver operates in two stages: first, employing score matching or Score-PINN to acquire the score function; and second, solving the LL via an ordinary differential equation (ODE) using the obtained score function. Comparative evaluations across these methods showcase varying trade-offs. The proposed methodology is evaluated across diverse SDEs, including anisotropic Ornstein-Uhlenbeck processes, geometric Brownian motion, and Brownian motion with varying eigenspace. We also test various distributions, including Gaussian, Log-normal, Laplace, and Cauchy distributions. The numerical results demonstrate the score-based SDE solver’s stability, speed, and performance across different experimental settings, solidifying its potential as a solution to CoD for high-dimensional FP equations.

97 MATHEMATICS AND COMPUTING

The Freundlich isotherm equation best represents phosphate sorption across soil orders and land use types in tropical soils of Puerto Rico

Biomass production in the lowland wet tropical forest is greater than in any other biome, and it is typically limited by soil phosphorus (P) availability. However, the mechanisms involved in the P cycle remain poorly represented in Earth System Models (ESMs). Soil P sorption processes are key in the P cycle and for understanding the extent of P limitation for plant productivity. Currently, a few ESMs include isotherm equations to model these processes. Although the Langmuir equation is widely cited, other isotherm equations may better describe sorption in tropical soils. Here, we use a diverse range of soil samples from Puerto Rico to test the validity of the Langmuir, Freundlich, and Temkin equation. We found that across four soil orders (Inceptisols, Mollisols, Oxisols, Ultisols), and forested and cultivated land use types, the Freundlich equation best represented soil P sorption. Furthermore, the Langmuir and the Temkin equations poorly represent soil P adsorption, especially at low P concentrations. Specifically, the Langmuir equation underestimated soil P adsorption by 40% and the Temkin equation overestimated adsorption by 76%. We also found, as expected, that soil clay content and pH were the most important parameters explaining the variability of the Freundlich (K f ) constant. Greater clay content and lower pH, common in highly weathered Ultisols and Oxisols which are abundant in the tropics, led to greater K f values. Overall, our results suggest that a diversity of soils can prompt underestimation of P sorption when using the Langmuir isotherm, which leads to an overestimation of available P that can have repercussions on ESM predictions of the P cycle and tropical forest productivity.

Earth system models

Singularity-EOS: Performance Portable Equations of State and Mixed Cell Closures

We present Singularity-EOS, a new performance-portable library for equations of state and related capabilities. Singularity-EOS provides a large set of analytic equations of state, such as the Gruneisen equation of state, and tabulated equation of state data under a unified interface. It also provides support capabilities around these equations of state, such as Python wrappers, solvers for finding pressure-temperature equilibrium between multiple equations of state, and a unique modifier framework, allowing the user to transform a base equation of state, for example by shifting or scaling the specific internal energy. All capabilities are performance portable, meaning they compile and run on both CPU and GPU for a wide variety of architectures.

97 MATHEMATICS AND COMPUTING

Identifying Differential Equations in Fourier Domain (FourierIdent)

We investigate identifying differential equations in the frequency domain. Fourier analysis is an important tool in theoretical analysis and numerical solvers of differential equations, yet there is limited work in exploring this connection in the identification of differential equations. This paper aims to identify the underlying differential equation in the frequency domain, from a given single realization of the differential equation perturbed by noise. Such setting imposes difficulties which are different from other identification methods where computation is carried out in the physical domain. We propose several ways to mitigate the challenges arising from noise in data and large differences in the magnitudes of frequency responses. The main takeaways are that identifying differential equations solely in the frequency domain is challenging, the method we propose is based on a form of domain partitions in the frequency domain, and this method shows benefits for complex data even with high level of noise. We introduce a Fourier feature denoising, and define the meaningful data region and the core regions of features to reduce the effect of noise in the frequency domain and to enhance the accuracy in coefficient identification. The proposed method is tested on various differential equations with linear, nonlinear, and high-order derivative feature terms, and shows advantages on complex data with many frequency modes, even under high level of noise.

97 MATHEMATICS AND COMPUTING

Directional finite difference method for directly solving 3D gyrokinetic field equations with enhanced accuracy

The gyrokinetic (GK) field equation is a three-dimensional (3D) elliptic equation, but it is often simplified to a set of two-dimensional (2D) equations by assuming that the field does not vary along a specific direction. However, this simplification can introduce inevitable 0th-order numerical errors, as nonlinear mode coupling in toroidal geometry can produce undesirable harmonic modes that violate the assumption. In this work, we propose a novel directional finite difference method (FDM) with a local coordinate transformation to better resolve the target field of interest. The directional FDM can accurately solve 3D GK field equations without simplifications, which can overcome the limitations of conventional methods. The accuracy and efficiency of different FDMs are analyzed in great detail for a variety of geometries, from simple 2D Cartesian coordinates to realistic 3D curvilinear coordinates. The 0th-order numerical errors of simplified 2D GK equations were found to be more problematic for low-harmonic modes and low aspect ratio geometries such as spherical tokamaks. On the other hand, the directional 3D FDM can accurately resolve a much wider range of harmonic modes aligned to the direction of interest, including the low-harmonic modes. In conclusion, we demonstrate that the directional 3D FDM is a highly effective algorithm for solving the 3D GK field equations, achieving accuracy improvements of 10 to 100 times or more, particularly for low-harmonic modes in spherical tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Numerical schemes for 3-wave kinetic equations: A complete treatment of the collision operator

In our previous work Walton and Tran (2023), numerical schemes for a simplified version of 3-wave kinetic equations, in which only the simple forward-cascade terms of the collision operators are kept, have been successfully designed, especially to capture the long time dynamics of the equation given the multiple blow-up time phenomenon. In this second work in the series, we propose numerical treatments for the complete 3-wave kinetic equations, in which the complete, much more complicated collision operators are fully considered based on a novel conservative form of the equation. Here we then derive an implicit finite volume scheme to solve the equation. The new discretization uses an adaptive time-stepping method which allows for the simulations to be carried to very long times. Our computed solutions are compared with previously derived long-time asymptotic estimates for the decay rate of total energy of time-dependent solutions of 3-wave kinetic equations and found to be in excellent agreement.

97 MATHEMATICS AND COMPUTING

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

Validating corrosion models: A comparison of governing equations

Experimental validation of Finite Element Method (FEM) models varying electrochemical governing equations, inclusion of chemical reactions, and time on the resultant damage profile for two galvanic couples is explored. Two anode materials (Magnesium AZ31 and Carbon Steel) in contact with a cathode (Stainless Steel 304 L) were modeled in/exposed to NaCl (1 and 0.1 M respectively for the anode materials) for up to one week. The physics approach, inclusion of chemical reactions, and the boundary conditions required to accurately represent the damage profile in FEM models depended on the galvanic couple materials and, ultimately, the corrosion rate. For high rates of corrosion (i.e., magnesium anode), the Nernst-Planck equation with Electroneutrality was sufficient to describe the damage, while, for low rates of corrosion (i.e., carbon steel anode), the Laplace equation was sufficient. In all cases, the most complete governing equation (Nernst-Planck-Poisson Equation) was not necessary to accurately describe the damage. Precipitation reactions in solution also played a critical role in the predicted damage profile, especially for high corrosion rate systems. Finally, for short time periods (< 6 h), the choice of governing equations does not significantly influence damage profile results. Overall, the choice of physics to reduce error in simulations relies on the boundary conditions, geometry, conductivity of the solution, electrochemical potential differences, and time of exposure. The above results are discussed with regard to accuracy and computational savings.

Carbon steel

Euler equations and the Sod shock tube problem

The Euler equations are a subset of the magnetohydrodynamic (MHD) equations in the infinitely collisional, unmagnetized limit. MHD modeling is central to many areas of plasma physics, ranging from low-temperature glow discharges to inertial confinement fusion. An important aspect of the Euler equations is their ability to describe states with discontinuities, such as shock waves. A standard benchmark test for numerical implementation of the Euler equations is the Sod shock tube. In this test, the system is initialized at rest with a pressure and density discontinuity, which results in a shock wave traveling into the low-pressure region and a rarefaction wave traveling into the high-pressure region. Starting with the presentation of the Euler equations, a numerical algorithm is presented here to solve these equations in one dimension. This is followed by an overview of the Sod shock tube problem that includes the precise initial setup and the analytic solution. Finally, the analytic solution is compared with results from numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY