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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Constraining the nonanalytic terms in the isospin-asymmetry expansion of the nuclear equation of state

In this report we examine the properties of the isospin-asymmetry expansion of the nuclear equation of state from chiral two- and three-body forces. We focus on extracting the high-order symmetry energy coefficients that consist of both normal terms (occurring with even powers of the isospin asymmetry) as well as terms involving the logarithm of the isospin asymmetry that are formally nonanalytic around the expansion point of isospin-symmetric nuclear matter. These coefficients are extracted from numerically precise perturbation theory calculations of the equation of state coupled with a new set of finite difference formulas that achieve stability by explicitly removing the effects of higher-order terms in the expansion. We consider contributions to the symmetry energy coefficients from both two- and three-body interactions. It is found that the coefficients of the logarithmic terms are generically larger in magnitude than those of the normal terms from second-order perturbation theory diagrams, but overall the normal terms give larger contributions to the ground state energy. The high-order isospin-asymmetry terms are especially relevant at large densities where they affect the proton fraction in β-equilibrium matter, and in particular we find that at twice saturation density they can reduce the proton fraction by up to 0.02.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Neutron stars with a crossover equation of state

The question of whether quark matter exists in neutron stars is a long standing one. Generally one finds that a first order phase transition from baryons to quarks softens the equation of state so much that the star would collapse into a black hole. We consider a crossover equation of state, similar to the crossover that is found in lattice QCD studies at finite temperature and zero or small baryon chemical potentials. Finally, we find that with reasonable parameters it may be possible to support neutron stars up to about 2.2 solar masses. In that case 1% to 10% of the pressure would be contributed by quark matter in the central core of the highest mass stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Radius and equation of state constraints from massive neutron stars and GW190814

Motivated by the unknown nature of the 2.50–2.67M⊙ compact object in the binary merger event GW190814, we study the maximum neutron star mass based on constraints from low-energy nuclear physics, neutron star tidal deformabilities from GW170817, and simultaneous mass-radius measurements of PSR J0030+045 from NICER. Our prior distribution is based on a combination of nuclear modeling valid in the vicinity of normal nuclear densities together with the assumption of a maximally stiff equation of state at high densities, a choice that enables us to probe the connection between observed heavy neutron stars and the transition density at which conventional nuclear physics models must break down. Here we demonstrate that a modification of the highly uncertain suprasaturation density equation of state beyond 2.64 times normal nuclear density is required in order for chiral effective field theory models to be consistent with current neutron star observations and the existence of 2.6M⊙ neutron stars. We also show that the existence of very massive neutron stars strongly impacts the radii of ≈ 2.0M⊙ neutron stars (but not necessarily the radii of 1.4 M⊙ neutron stars), which further motivates future NICER radius measurements of PSR J1614 – 2230 and PSR J0740+6620.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Mapping out the thermodynamic stability of a QCD equation of state with a critical point using active learning

The Beam Energy Scan Theory (BEST) collaboration's equation of state (EoS) incorporates a three-dimensional Ising model critical point into the quantum chromodynamics (QCD) equation of state from lattice simulations. Furthermore, it contains four free parameters related to the size and location of the critical region in the QCD phase diagram. Certain combinations of the free parameters lead to acausal or unstable realizations of the EoS that should not be considered. In this work, we use an active learning framework to rule out pathological EoS efficiently. We find that checking stability and causality for a small portion of the parameters' range is sufficient to construct algorithms that perform with > 96 % accuracy across the entire parameter space. Though in this work we focus on a specific case, our approach can be generalized to any EoS containing a parameter space-class correspondence.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Four-dimensional QCD equation of state with multiple chemical potentials

Here, we construct a four-dimensional version of the equation of state (EoS) model neos, neos-4d, as a function of the temperature and chemical potentials of baryon, electric charge, and strangeness for the hot and dense quantum chromodynamics (QCD) matter created in relativistic nuclear collisions. This EoS enables multiple conserved charge current evolution in a relativistic fluid. Input from lattice QCD simulations and a hadron resonance gas model is considered for constructing the equation of state. We investigate its applicability to the relativistic hydrodynamic description of nuclear collisions and present a method for efficient numerical implementation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Impact of O4 future detections on the determination of the dense matter equations of state

Here, in view of the next LIGO-Virgo-KAGRA Observing period O4 (to start in Spring 2023), we address the question of the ability of the interferometers network to discriminate among different neutron stars equation of states better than what was possible with the observation of the binary neutron stars merger GW170817. We show that the observation of an event similar to GW170817 during O4 would allow us to resolve the dimensionless effective tidal deformability $\tilde{Λ}$ within an uncertainty 7 times better than the one obtained in O2. Thanks to the expected increase in sensitivities, we show that any GW170817-like single-event within a distance of 100 Mpc would imply significantly improved constraints of the neutron stars equations of state. We also illustrate the important impact of the noise in the analysis of the signal, showing how it can impact the effective tidal deformability probability density function for large signal-to-noise ratio.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Central speed of sound, the trace anomaly, and observables of neutron stars from a perturbative analysis of scaled Tolman-Oppenheimer-Volkoff equations

The central speed of sound (SS) measures the stiffness of the equation of state (EOS) of superdense neutron star (NS) matter. Its variations with density and radial coordinate in NSs in conventional analyses often suffer from uncertainties of the specific nuclear EOS used. Using the central SS and NS mass/radius scaling obtained from solving perturbatively the scaled Tolman-Oppenheimer-Volkoff (TOV) equations, we study the variations of SS, trace anomaly and several closely related properties of NSs in an EOSmodel- independent manner. We find that the SS increases with the reduced central pressure Pˆ c ≡ P c =ε c (scaled by the central energy density εc), and the conformal bound for SS tends to break down for NSs with masses higher than about 1.9M ⊙ . The ratio P=ε is upper bounded as P=ε ≲ 0.374 around the centers of stable NSs. We demonstrate that it is an intrinsic property of strong field gravity and is more relevant than the perturbative QCD bound on it. While a sharp phase transition at high densities characterized by a sudden vanishing of SS in cores of massive NSs are basically excluded, the probability for a continuous crossover signaled by a peaked radial profile of SS is found to be enhanced as Pˆ c decreases, implying it likely happens near the centers of massive NSs. Moreover, a new and more stringent causality boundary as R max =km ≳ 4.73$M$$^{max}_{NS}$ =M ⊙ þ 1.14 for the NS mass-radius curve is found to be excellently consistent with observational data on NS masses and radii. Here, new constraints on the ultimate energy density and pressure allowed in NSs before collapsing into black holes are obtained and compared with earlier predictions in the literature

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Finite density QCD equation of state: Critical point and lattice-based 𝑇′ expansion

Here, we present a novel construction of the QCD equation of state (EoS) at finite baryon density. Our work combines a recently proposed resummation scheme for lattice QCD results with the universal critical behavior at the QCD critical point. This allows us to obtain a family of equations of state in the range 0 ≤ 𝜇 𝐵 ≤ 700 MeV and 25 MeV ≤ 𝑇 ≤ 800 MeV, which match lattice QCD results near 𝜇 𝐵 =0 while featuring a critical point in the 3D Ising model universality class. The position of the critical point can be chosen within the range accessible to beam-energy scan heavy-ion collision experiments. The strength of the singularity and the shape of the critical region are parametrized using a standard parameter set. We impose stability and causality constraints and discuss the available ranges of critical point parameter choices, finding that they extend beyond earlier parametric QCD EoS proposals. We present thermodynamic observables, including baryon density, pressure, entropy density, energy density, baryon susceptibility and speed of sound, that cover a wide range in the QCD phase diagram relevant for experimental exploration.

Astronomy & Astrophysics↗
DOE OSTI2020

Benchmarking boron carbide equation of state using computation and experiment

Boron carbide ( B 4 C ) is of both fundamental scientific and practical interest due to its structural complexity and how it changes upon compression, as well as its many industrial uses and potential for use in inertial confinement fusion (ICF) and high-energy density physics experiments. Here, we report the results of a comprehensive computational study of the equation of state (EOS) of B 4 C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up to 61 megabar from planar shock experiments performed at the National Ignition Facility (NIF). Our computational methods include path integral Monte Carlo, activity expansion, as well as all-electron Green's function Korringa-Kohn-Rostoker and molecular dynamics that are both based on density functional theory. We calculate the pressure-internal energy EOS of B 4 C over a broad range of temperatures ( ~ 6 × 10 3 – 5 × 10 8 K) and densities (0.025–50 g / cm 3 ). We assess that the largest discrepancies between theoretical predictions are ≲ 5 % near the compression maximum at 1– 2 × 10 6 K. This is the warm-dense state in which the K shell significantly ionizes and has posed grand challenges to theory and experiment. By comparing with different EOS models, we find a Purgatorio model (LEOS 2122) that agrees with our calculations. The maximum discrepancies in pressure between our first-principles predictions and LEOS 2122 are ~ 18 % and occur at temperatures between 6 × 10 3 – 2 × 10 5 K, which we believe originate from differences in the ion thermal term and the cold curve that are modeled in LEOS 2122 in comparison with our first-principles calculations. To account for potential differences in the ion thermal term, we have developed three new equation-of-state models that are consistent with theoretical calculations and experiment. We apply these new models to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

How Well Do We Know the Neutron-Matter Equation of State at the Densities Inside Neutron Stars? A Bayesian Approach with Correlated Uncertainties

Here, we introduce a new framework for quantifying correlated uncertainties of the infinite-matter equation of state derived from chiral effective field theory (𝜒⁢EFT ). Bayesian machine learning via Gaussian processes with physics-based hyperparameters allows us to efficiently quantify and propagate theoretical uncertainties of the equation of state, such as 𝜒⁢EFT truncation errors, to derived quantities. We apply this framework to state-of-the-art many-body perturbation theory calculations with nucleon-nucleon and three-nucleon interactions up to fourth order in the 𝜒⁢EFT expansion. This produces the first statistically robust uncertainty estimates for key quantities of neutron stars. We give results up to twice nuclear saturation density for the energy per particle, pressure, and speed of sound of neutron matter, as well as for the nuclear symmetry energy and its derivative. At nuclear saturation density, the predicted symmetry energy and its slope are consistent with experimental constraints.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Dynamic Parameter Estimation with Physics-based Neural Ordinary Differential Equations

Accurate estimation of dynamic parameters of gen-erators is crucial to building a reliable model for dynamical studies and reliable operation of the power system. This paper develops a physics-based neural ordinary differential equations (ODE) approach to learn the parameters of generator dynamic model using phasor measurement units (PMU) data. We design a physics-based neural network to represent the swing equations of the power system dynamics. A loss function is defined as the difference between dynamic simulation results from the physics-based neural networks and pseudo PMU measurements. The parameters of generator dynamic model are iteratively updated using the neural ODEs and the adjoint method. By exploiting the mini-batch scheme in neural ODE training, the parameter estimation performance is significantly improved. Numerical study results on a 3-machine 9-bus system show that the proposed algorithm outperforms state-of-the-art baseline method in both computation time and dynamic parameter estimation accuracy.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Multidimensional Tests of a Finite-Volume Solver for MHD With a Real-Gas Equation of State

This work considers two algorithms of a finite-volume solver for the MHD equations with a real-gas equation of state (EOS). Both algorithms use a multistate form of Harten-Lax-Van Leer approximate Riemann solver as formulated for MHD discontinuities. This solver is modified to use the generalized sound speed from the real-gas EOS. Two methods are tested: EOS evaluation at cell centers and at flux interfaces where the former is more computationally efficient. A battery of 1D and 2D tests are employed: convergence of 1D and 2D linearized waves, shock tube Riemann problems, a 2D nonlinear circularly polarized Alfvén wave, and a 2D magneto-Rayleigh-Taylor instability test. The cell-centered EOS evaluation algorithm produces unresolvable thermodynamic inconsistencies in the intermediate states leading to spurious solutions while the flux-interface EOS evaluation algorithm robustly produces the correct solution. The linearized wave tests show this inconsistency is associated with the magnetosonic waves and the magneto-Rayleigh-Taylor instability test demonstrates simulation findings where the spurious solution leads to an unphysical simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Fully quantum algorithm for mesoscale fluid simulations with application to partial differential equations

Fluid flow simulations marshal our most powerful computational resources. In many cases, even this is not enough. Quantum computers provide an opportunity to speed up traditional algorithms for flow simulations. We show that lattice-based mesoscale numerical methods can be executed as efficient quantum algorithms due to their statistical features. This approach revises a quantum algorithm for lattice gas automata to reduce classical computations and state preparation at every time step. For this, the algorithm approximates the qubit relative phases and subtracts them at the end of each time step. Phases are evaluated using the iterative phase estimation algorithm and subtracted using single-qubit rotation phase gates. Further, this method optimizes the quantum resource required and makes it more appropriate for near-term quantum hardware. We also demonstrate how the checkerboard deficiency that the D1Q2 scheme presents can be resolved using the D1Q3 scheme. The algorithm is validated by simulating two canonical partial differential equations: the diffusion and Burgers' equations on different quantum simulators. We find good agreement between quantum simulations and classical solutions for the presented algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Asymptotic Relaxation of Moment Equations for a Multi-species, Homogeneous BGK Model

Multi-species BGK models describe the dynamics of rarefied gases with constituent particles of different elements or compounds with potentially nontrivial velocity distributions. Here, in this paper, moment equations for the bulk velocities, energies, and temperatures of a spatially homogeneous multi-species BGK model are examined. A key challenge in analyzing these equations is the fact that the collision frequencies are allowed to depend on the species temperatures, which allows for more realistic simulations of dilute gas flow. Therefore, a positive lower bound is established for the species temperatures. With this lower bound, a global existence and uniqueness of solutions to the coupled velocity-energy ODE system is established. The lower bound also enables a proof of exponential decay to a unique steady-state solution. Numerical results are presented to demonstrate how the bulk velocities and temperatures relax for large times.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Strangeness-neutral equation of state for QCD with a critical point

We present a strangeness-neutral equation of state for QCD that exhibits critical behavior and matches lattice QCD results for the Taylor-expanded thermodynamic variables up to fourth-order in µ B /T. It is compatible with the SMASH hadronic transport approach and has a range of temperatures and baryonic chemical potentials relevant for phase II of the Beam Energy Scan at RHIC. We provide an updated version of the software BES-EoS, which produces an equation of state for QCD that includes a critical point in the 3D Ising model universality class. This new version also includes isentropic trejectories and the critical contribution to the correlation length. Since heavy-ion collisions have zero global net-strangeness density and a fixed ratio of electric charge to baryon number, the BES-EoS is more suitable to describe this system. Comparison with the previous version of the EoS is thoroughly discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Power System Frequency Dynamics Modeling, State Estimation, and Control using Neural Ordinary Differential Equations (NODEs) and Soft Actor-Critic (SAC) Machine Learning Approaches

With the global energy transition of the electric power system, grid control, supervision, and protection is becoming more challenging. With the increasing integration of renewable energy sources (RES), the system dynamics are changing, causing traditional power system dynamic modeling with swing equation-based modeling approaches to fail. Additionally, the converter-dominated power grid is decreasing the system inertia, making the power system more fragile to the frequency swings. This paper first investigates and compares the application of a model-based Kalman filter state estimation approach with (i) a model-free machine learning approach --- neural ordinary differential equations (NODEs) --- and (ii) a data-driven system identification (SysId) approach to model and infer critical state values of the power system frequency dynamics. Then a model predictive control (MPC) framework is compared to a model-free Soft Actor-Critic (SAC) reinforcement learning (RL) control algorithm in providing efficient fast frequency response (FFR) to the power system frequency dynamics. The approaches are compared in terms of their performance goals as well as their per-timestep computational efficiency. Furthermore, the comparative study for state estimation shows that for the model-free requirement, both NODEs and SysId can provide accurate state estimates; however, with increasing model complexity, NODEs can be a better choice for model identification. Similarly, the results from the FFR comparative study show that the SAC RL-based FFR, once trained, outperforms MPC with better control signals and faster computation time, making the SAC RL-based FFR better option for providing FFR to the power system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A sixth order Mehrstellen scheme with an application to the Method of Local Corrections for the 3D Poisson equation

We present a sixth order finite difference scheme for Poisson’s equation when discretized with the compact 27-point stencil based on Mehrstellen corrections of the forcing function term f. Our approach results in a sixth order accurate solution error as opposed to a fourth-order error imposed by the classical Mehrstellen correction for the 19-point and 27-point stencils. The present study is a continuation of former work of Spotz and Carey (1996) on compact finite difference schemes for Poisson’s equation where sixth order convergence may be obtained under the assumption that the fourth order derivatives of f are determined analytically. Specifically, we show that sixth order convergence can still be attained when only values of f at grid points are available. The sixth order Mehrstellen scheme is further coupled with a Method of Local Corrections (MLC) 3D Poisson solver improving to sixth order accuracy the results reported in Kavouklis and Colella (2019). The MLC test case considered involves an adaptive grid that comprises 7.5 billion cells.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Constraining the Equation of State of Dense Neutron-rich Matter (Final Report)

Using a new Time Projection Chamber (TPC), called the SpiRIT (SAMURAI pion Reconstruction Ion Tracker) constructed and funded by DOE office of science (award # DE-SC0004835), we performed two experiments with the SAMURAI spectrometer at RIKEN, Japan to study the equation of state of neutron rich matter. As a result of the project, the SpiRIT collaboration, an international collaboration consisting of groups from US, Japan, Korea, Poland, China and Germany, has been formed to pursue the science opportunities provided by the SpiRIT TPC. After the experiments, we developed the software to analyze the SpiRIT experiments and extract constraints of symmetry energy at supra-saturation densities. As a result 270 TB of data have been obtained. Over 10 technical papers on the TPC and 4 science papers on the nuclear equation of state have been published. A total of six (2 US, 1 Japanese and 1 Korean) PhD students have graduated based on their research on the EOS research using the SpiRIT data. In addition, the data also stimulated theorists on transport models to form the “Transport Model Evaluation Project” collaboration to improve the transport model to compare calculations with the data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗