Search NASA⌕ Search

SEARCH · Search NASA

Results for “Equations of state”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

SESAME: The Los Alamos National Laboratory’s Tabulated Equation of State Database Description with Extensions for Multi-phase Representations

Modeling of the thermodynamic equation of state (EOS) of various materials has had a long storied tradition at LANL. As early as 1949 Feynman, Metropolis, and Teller published a paper presenting EOS values for some elements and a methodology for calculating the EOS at high compression [1]. Cowan and Ashkin made notable methodology improvements for compressed materials throughout the 1950’s and beyond [2]. In 1971 Jack Barnes and Jerry Rood created the SESAME database and by 1972 the database became publicly available.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Simultaneous inference of equation of state parameters and unknown data errors with uncertainty quantification via hierarchical Bayesian posterior maximization

Equations of state (EOSs) are a key component in running hydrodynamic simulations as they relate the thermodynamic states for the material. The Davis reactants EOS is commonly used for modeling high explosives (HEs), and the EOS model parameters are calibrated using material specific data. The calibrations are often performed with uncertainty quantification via Bayesian inference to account for uncertainty in the data and generate ensembles of likely parameters. However, there are relatively few HE data sets to use for calibration and many are historical and lack error information. In this work, we simultaneously calibrate the Davis reactants EOS model parameters and unknown data error terms for the high explosive PBX 9501. To quantify the uncertainty in the models and the data, we use a Bayesian framework for the calibration and compute the hierarchical Bayesian posterior distribution with both a posteriori maximization approach and Markov Chain Monte Carlo. In general, we find that, given our assumptions, the two approaches result in similar calibrated parameters, posterior covariance matrices, and insights about the parameters but that the posterior maximization requires far less computational resources.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Astrophysical Equation-of-State Constraints on the Color-Superconducting Gap

We demonstrate that astrophysical constraints on the dense-matter equation of state place an upper bound on the color-superconducting gap in dense matter above the transition from nuclear matter to quark matter. Pairing effects in the color-flavor locked quark matter phase increase the pressure at high density, and if this effect is sufficiently large then the requirements of causality and mechanical stability make it impossible to reach such a pressure in a way that is consistent with what is known at lower densities. The intermediate-density equation of state is inferred by considering extensions of chiral effective field theory to neutron star densities, and conditioning these using current astrophysical observations of neutron star radius, maximum mass, and tidal deformability (PSR J⁢0348+0432, PSR J1624-2230, PSR J⁢0740+6620, GW170817). At baryon number chemical potential μ = 2.6 GeV we find a 95% upper limit on the color-flavor locked pairing gap Δ of 457 MeV using overly conservative assumptions and 216 MeV with more reasonable assumptions. Furthermore, this constraint may be strengthened by future astrophysical measurements as well as by future advances in high-density QCD calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Dynamical ejecta from binary neutron star mergers: Impact of a small residual eccentricity and of the equation of state implementation

Predicting the properties of the matter ejected during and after a neutron star merger is crucial to our ability to use electromagnetic observations of these mergers to constrain the masses of the neutron stars, the equation of state of dense matter, and the role of neutron star mergers in the enrichment of the Universe in heavy elements. Furthermore, our ability to reliably provide such predictions is however limited by a broad range of factors, including the finite resolution of numerical simulations, their treatment of magnetic fields, neutrinos, and neutrino-matter interactions, and the approximate modeling of the equation of state of dense matter. In this manuscript, we study specifically the role that a small residual eccentricity and different implementations of the same equation of state have on the matter ejected during the merger of a 1.3M ⊙ –1.4M ⊙ binary neutron star system. We find that a residual eccentricity e ~ 0.01, as measured ~ 4–6 orbits before merger, causes O(25%–30%) changes in the amount of ejected mass, mainly due to changes in the amount of matter ejected as a result of core bounces during merger. We note that O(1%) residual eccentricities have regularly been used in binary neutron star merger simulations as proxy for circular binaries, potentially creating an additional source of error in predictions for the mass of the dynamical ejecta.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Verification and validation of detonation-shock-dynamics relations for explosives described by general equation of state and chemical reaction models

Detonation shock dynamics is a powerful method to model the behaviour of High Explosives (HE). However in order to use this method, the underlying relationship between the local radius of curvature and the detonation speed must be known. Previous work has developed methods to calculate this effect using simple, single-step Arrhenius and polytropic gas, models for the chemical reaction and the equation of state, respectively. In recent years, more complex models for both reaction rates and equations of state have been developed which show better agreement with experimental data than these simple models, especially when considering condensed phase explosives.. This work presents the governing equations for solving these problems in a way that is generalised to use arbitrary equations of state as well as reaction models which may have more than a single step and multiple product species. This implementation is verified against exact solutions, demonstrating that the equations were implemented properly. The verified algorithm is then validated against experimental data and high fidelity simulations, showing that it is able to make accurate predictions in a regime where the underlying assumptions of the governing equations are valid. Importantly, this approach has many applications: from creating equivalent detonation shock dynamics models for existing reactive burn calibrations for HE; to developing new functional forms and calibrations of reactive burn models for condensed phase high explosives.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1984

Conservative streamtube solution of steady-state Euler equations

This paper presents a new method for solving the steady state Euler equations. The method is similar to streamline curvature methods but has a conserative finite volume formulation which ensures correct shock capturing. Either wall position or wall pressure may be prescribed as boundary conditions, permitting both direct and inverse calculations. In supersonic applications the solution is obtained by space-marching while in subsonic and transonic applications iterative relaxation methods are used. Numerical results are given for: (1) supersonic diffuser with oblique shocks (direct calculation); (2) supersonic jet entering still reservoir (inverse calculation); (3) subsonic bump in a channel with 25 percent blockage (direct and inverse); (4) subsonic high-work turbine cascade (direct); and (5) transonic bump in a channel with 12 percent blockage (direct calculation).

Drela, M.↗
NASA NTRS1991

Shock wave equation of state of serpentine to 150 GPa - Implications for the occurrence of water in the earth's lower mantle

The shock wave equation of state of Mg end-member serpentine was determined to 150 GPa by examining the shock properties of three polycrystalline serpentines: (1) a lizardite serpentine found near Globe (Arizona), (2) an antigorite serpentine from Thurman (New York), and (3) a chrysotile serpentine from Quebec (Canada). The shock wave experiments were carried out using either a two-stage light gas gun or a 40-mm bore propellant. The shock equation of state that was obtained is shown to exhibit four distinct regions: a low-pressure phase, a mixed phase region, a high-pressure phase, and a very high-pressure phase. The high-pressure density and sound speed of an H2O-rich magnesium silicate determined from these experiments indicate that the observed seismic properties of the lower mantle allow the existence of several weight percent of water in the lower mantle.

Tyburczy, James A.↗
NASA NTRS1971

Ride comfort control in large flexible aircraft

The problem of ameliorating the discomfort of passengers on a large air transport subject to flight disturbances is examined. The longitudinal dynamics of the aircraft, including effects of body flexing, are developed in terms of linear, constant coefficient differential equations in state variables. A cost functional, penalizing the rigid body displacements and flexure accelerations over the surface of the aircraft is formulated as a quadratic form. The resulting control problem, to minimize the cost subject to the state equation constraints, is of a class whose solutions are well known. The feedback gains for the optimal controller are calculated digitally, and the resulting autopilot is simulated on an analog computer and its performance evaluated.

Warren, M. E.↗
NASA NTRS1987

Control of SCOLE

A relatively low order model is used to control SCOLE. Drastic truncation of the discretized model is proposed by means of a modal expansion. An open loop eigenvalue problem is illustrated as is truncated modal equations, modal state equations, actual output vector and modal Kalman filter. Also illustrated is independent modal-space control.

Meirovitch, L.↗
NASA NTRS1989

Flux-split algorithms for the multi-dimensional Euler equations with real gases

Upwind algorithms are developed for the numerical solution of the multidimensional Euler equations for real gases. Flux-splitting methods are derived which account for a general equation of state. Approximations to the state equation based on physical arguments result in simplified algorithms which may be implemented into existing perfect-gas codes. Applications of the method to several high-Mach-number high-temperature flows are presented for two and three space dimensions.

Grossman, B.↗
DOE OSTI2025

Direct nonparametric multimessenger constraints on the equation of state of cold dense nuclear matter

We utilize the now substantial amount of astrophysical observations of neutron stars (NSs), along with perturbative quantum chromodynamics (pQCD) calculations at high density, to directly constrain the NS equation of state (EOS). To this end, we construct nonparametric EOS priors by using Gaussian processes trained on 75 EOSs, which include models with either hadrons, hyperons, or quarks at high densities. We create a prior using the full EOS sample (model agnostic), and one prior for each EOS family to test model discrimination. We introduce a novel inference approach, which allows the simultaneous sampling of intrinsic and extrinsic parameters of binary NS mergers, as well as a nonparametric equation of state. We showcase this method in a Bayesian updating scheme by first performing a complete analysis of the binary NS merger event GW170817 with minimal assumptions, and sequentially adding information from x-ray and radio NS observations, along with pQCD calculations. Besides providing standard constraints, such as the pressure at twice nuclear saturation density 𝑝⁡(2⁢𝜌 sat ) = 4.3$^{+0.6}_{−0.6}$ × 10 34 dyne/cm 2 , at 95% confidence level, for the model agnostic prior, our methodology shows how the choice of EOS families used in conditioning changes the inferred astrophysical properties of the EOS, namely tidal deformability and maximum supported NS mass. We find hyperonic priors predicting higher tidal deformabilities for a 1.4⁢𝑀 ⊙ NS, and hadronic priors being preferred by the considered astrophysical data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1977

Equations of state and impact-induced shock-wave attenuation on the moon

Equation of state formulations are considered in a framework that permits comparison with one-dimensional impedance match solutions. The problem considered is the peak pressures attained along the impact symmetry axis when a sphere impacts with a half-space. The regimes of melting and vaporization - in particular, incipiently melted, completely melted, incipiently vaporized, and completely vaporized states - are examined, and the pressures at which critical isentropes intersect the Hugoniots of iron and gabbroic anorthosite are considered. A means of representing the spatial attenuation of shock pressure along the impact axis by two regimes is introduced, and results for the near-field and far-field regime are presented. It is thought that the treatment can be used to obtain quantitative bounds on the impact velocity of the meteorite.

Ahrens, T. J.↗
DOE OSTI2024

Sensitivity of magnetohydrodynamic simulations of Joule-heated conductors to the vaporization curve in equations of state

Magnetohydrodynamic (MHD) simulations of electrically exploded aluminum and copper rods demonstrate a technique to validate equations of state (EOS) for rapidly Joule-heated conductors. The balance of internal and magnetic forces at the conductor-insulator interface drives the metal there along the vaporization phase boundary. Variations between critical points and vaporization curves in existing models predict differing densities and temperatures in MHD simulations for these models. Here, the inclusion of Maxwell constructs in the liquid-vapor biphase region of the EOS caused the rod surface to vaporize earlier in time than unmodified tables with van der Waals loops. Velocimetry of recent experiments is used to validate the location of the vaporization curve in existing EOS models and differentiate between the vapor dome treatments. Dielectric coatings applied to the metal surface restricted the conductor’s expansion and diverted the metal into the warm dense matter regime.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Machine Learning-Driven Conservative-to-Primitive Conversion in Hybrid Piecewise Polytropic and Tabulated Equations of State

We present a novel machine learning (ML)-based method to accelerate conservative-to-primitive inversion, focusing on hybrid piecewise polytropic and tabulated equations of state. Traditional root-finding techniques are computationally expensive, particularly for large-scale relativistic hydrodynamics simulations. To address this, we employ feedforward neural networks (NNC2PS and NNC2PL), trained in PyTorch (2.0+) and optimized for GPU inference using NVIDIA TensorRT (8.4.1), achieving significant speedups with minimal accuracy loss. The NNC2PS model achieves 𝐿 1 and 𝐿 ∞ errors of 4.54 × 10 −7 and 3.44 × 10−6, respectively, while the NNC2PL model exhibits even lower error values. TensorRT optimization with mixed-precision deployment substantially accelerates performance compared to traditional root-finding methods. Specifically, the mixed-precision TensorRT engine for NNC2PS achieves inference speeds approximately 400 times faster than a traditional single-threaded CPU implementation for a dataset size of 1,000,000 points. Ideal parallelization across an entire compute node in the Delta supercomputer (dual AMD 64-core 2.45 GHz Milan processors and 8 NVIDIA A100 GPUs with 40 GB HBM2 RAM and NVLink) predicts a 25-fold speedup for TensorRT over an optimally parallelized numerical method when processing 8 million data points. Moreover, the ML method exhibits sub-linear scaling with increasing dataset sizes. We release the scientific software developed, enabling further validation and extension of our findings. By exploiting the underlying symmetries within the equation of state, these findings highlight the potential of ML, combined with GPU optimization and model quantization, to accelerate conservative-to-primitive inversion in relativistic hydrodynamics simulations.

conservative-to-primitive conversion↗
NASA NTRS1991

Cold equation of state in a strong magnetic field - Effects of inverse beta-decay

The influence of a high magnetic field (B is greater than 10 exp 12 G) on the degenerate matter equation of state appropriate to a neutron star is studied. The regime dominated by relativistic electrons up to the neutron drip density is highlighted. The equilibrium matter composition and equation of state, allowing for inverse beta-decay. Two different equilibrium models are determined: an ideal neutron-proton-electron (npe) gas and the more realistic model of Baym, Pethick, and Sutherland (1971) consisting of a Coulomb lattice of heavy nuclei embedded in an electron gas. For a sufficiently high field strength, the magnetic field has an appreciable effect, changing the adiabatic index of the matter and the nuclear transition densities. The influence of a strong field on some simple nonequilibrium processes, including beta-decay and inverse beta-decay (electron capture) is also considered. The effects produced by the magnetic field are mainly due to the changes in the transverse electron quantum orbits and the allowed electron phase space induced by the field.

Lai, Dong↗
NASA NTRS1976

Equations of state and impact-induced shock-wave attenuation on the moon

Current equation of state formulations, used for finite difference cratering flow calculations, are cast into a framework permitting comparison of peak pressures attained upon impact of a sphere, with a half-space, along the impact symmetry axis, to one dimensional impedance match solutions. On the basis of this formulation and application of thermochemical data, the regimes of melting and vaporization are examined. For the purpose of identifying material which will, upon isentropic release from the impact-induced shock state, result in a solid just brought to its melting point, i.e., incipiently melted (IM); completely melted (CM); just brought to its boiling point, i.e., incipiently vaporized (IV); and completely vaporized (CV) state, the pressures at which the critical isentropes intersect the Hugoniots of iron and gabbroic anorthosite (GA) are examined in detail. The latter rock type is assumed to be representative of the lunar highlands. The Hugoniot pressures, for which IM, CM, IV, and CV will occur upon isentropic expansion, are calculated to be 2.2, 2.6, 4.2, and 16.8 Mbar, respectively.

Ahrens, T. J.↗