SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

DOE OSTI2025

JCZ3: An Equation of State Based on the Exponential-Six Intermolecular Potential – Update: Formulation

This report presents a comprehensive analysis of the JCZ3 equation of state (EOS) as implemented in the TIGER code, focusing on its thermodynamic framework, mathematical formulations, and implications for modeling gas phase behavior under varying conditions. Overall, the report affirms the TIGER code's foundational robustness and its potential as a valuable tool for simulating high-pressure, high-temperature gas mixtures, while also highlighting areas for further refinement and validation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

thornado-hydro: A Discontinuous Galerkin Method for Supernova Hydrodynamics with Nuclear Equations of State

This paper describes algorithms for non-relativistic hydrodynamics in the toolkit for high-order neutrino radiation hydrodynamics (thornado), which is being developed for multiphysics simulations of core-collapse supernovae (CCSNe) and related problems with Runge–Kutta discontinuous Galerkin (RKDG) methods. More specifically, thornado employs a spectral type nodal collocation approximation, and we have extended limiters — a slope limiter to prevent non-physical oscillations and a bound-enforcing limiter to prevent non-physical states — from the standard RKDG framework to be able to accommodate a tabulated nuclear equation of state (EoS). To demonstrate the efficacy of the algorithms with a nuclear EoS, we first present numerical results from basic test problems in idealized settings in one and two spatial dimensions, employing Cartesian, spherical-polar, and cylindrical coordinates. Then, we apply the RKDG method to the problem of adiabatic collapse, shock formation, and shock propagation in spherical symmetry, initiated with a 15 M ⊙ progenitor. Herein, we find that the extended limiters improve the fidelity and robustness of the RKDG method in idealized settings. The bound-enforcing limiter improves robustness of the RKDG method in the adiabatic collapse application, while we find that slope limiting in characteristic fields is vulnerable to structures in the EoS — more specifically, in the phase transition from nuclei and nucleons to bulk nuclear matter. The success of these applications marks an important step toward applying RKDG methods to more realistic CCSN simulations with thornado in the future.

Nuclear astrophysics (1129)↗

NASA NTRS1992

Solar oscillations and the equation of state

The dependence of frequencies of solar oscillations on the thermodynamic state of the sun is considered. On the basis of an analysis of the equations of stellar structure, as well as the relevant aspects of the properties of the oscillations, it is argued that in the convection zone, information about the equation of state which is relatively unaffected by other uncertainties in the physics of the solar interior can be isolated. The different treatments that have been employed to describe the thermodynamics of stellar plasmas are reviewed. Through application of several of these treatments to the computation of models of the solar envelope, it is demonstrated that the sensitivity of the observed frequencies is in fact sufficient to distinguish even quite subtle features of the physics of solar matter.

Christensen-Dalsgaard, Jorgen↗

NASA NTRS2004

State-Dependent Riccati Equation Regulation of Systems with State and Control Nonlinearities

The state-dependent Riccati equations (SDRE) is the basis of a technique for suboptimal feedback control of a nonlinear quadratic regulator (NQR) problem. It is an extension of the Riccati equation used for feedback control of linear problems, with the addition of nonlinearities in the state dynamics of the system resulting in a state-dependent gain matrix as the solution of the equation. In this paper several variations on the SDRE-based method will be considered for the feedback control problem with control nonlinearities. The control nonlinearities may result in complications in the numerical implementation of the control, which the different versions of the SDRE method must try to overcome. The control methods will be applied to three test problems and their resulting performance analyzed.

Beeler, Scott C.↗

NASA NTRS1973

An equation of state for oxygen and nitrogen

Recent measurements of thermodynamic properties of oxygen and nitrogen have provided data necessary for development of a single equation of state for both fluids. Data are available in summary report and two-part detailed study on thermodynamic properties of oxygen and nitrogen. Same data are used to develop vapor-pressure equation and heat-capacity equation.

Jacobsen, R. T.↗

DOE OSTI2026

Microscopic constraints for the equation of state and structure of neutron stars: A Bayesian model mixing framework

Bayesian model mixing (BMM) is a statistical technique that can combine constraints from different regions of an input space in a principled way. Here we extend our BMM framework for the equation of state (EOS) of strongly interacting matter from symmetric nuclear matter to asymmetric matter, specifically focusing on zero-temperature, charge-neutral, 𝛽-equilibrated matter. We use Gaussian processes (GPs) to infer constraints on the neutron-star matter EOS at intermediate densities from two different microscopic theories: chiral effective-field theory (𝜒⁢EFT) at baryon densities around nuclear saturation, 𝑛 𝐵 ∼ 𝑛 0 , and perturbative QCD at asymptotically high baryon densities, 𝑛 𝐵 ⩾ 20⁢𝑛 0 . The uncertainties of the 𝜒⁢EFT and pQCD EOSs are obtained using the BUQEYE truncation error model. We demonstrate the flexibility of our framework through the use of two categories of GP kernels: conventional stationary kernels and a nonstationary changepoint kernel. We use the latter to explore potential constraints on the dense matter EOS by including exogenous data representing theory predictions and heavy-ion collision measurements at densities ⩾ 2⁢𝑛 0 . We also use our EOSs to obtain neutron-star mass-radius relations and their uncertainties. Finally, our framework, whose implementation will be available through a GitHub repository, provides a prior distribution for the EOS that can be used in large-scale neutron-star inference frameworks.

Bayesian methods↗

DOE OSTI2020

How equation of state selection impacts accuracy near the critical point: Forced convection supercritical CO 2 flow over a cylinder

Increasing application of supercritical fluids leads to increased use of numerical simulations to guide designs. However, the impact of Equations of State (EoS) selection on computational cost and accuracy is not well understood. Laminar flow around a heated cylinder is used as a test case. A robust matrix of 99-simulations is used to analyze the flow of carbon dioxide just above the critical point at the isobaric condition of 8 MPa and temperatures between 305 and 390 K. Four EoS were studied: Ideal Gas Law, cubic EoS’s Peng-Robinson and Soave-Redlich-Kwong, and the Span-Wagner EoS. Results were identical for isothermal cases, but for heated or cooled cases, the high-fidelity Span-Wagner EoS simulations presented drastically different wake structures. Compared to the high-fidelity SW EoS, cEoS’s underestimated Nu values by up to 38% and cEoS’s underestimated CD values by 18–42%. In conclusion, the Span-Wagner EoS required only 17% more computation time then the cubic EoS simulations, making the Span-Wagner EoS recommended for engineering design.

42 ENGINEERING↗

NASA NTRS1985

Parametrized equation of state for electron liquids in the Singwi-Tosi-Land-Sjolander aproximation

Results are reported of a theoretical study of an equation of state for electron liquids (one-component plasmas of electrons embedded in a uniform neutralizing background of positive charges), where there is an interplay between the strong Coulomb-coupling effect and the degrees of Fermi degeneracy. Calculations are based on the Singwi-Tosi-Land-Sjolander approximation (1968). The calculated results are parametrized in the form of analytic formulas for the interaction and excess free energies, which are applicable over a wide range of parameters as long as the electrons are in a paramagnetic fluid state. Unlike the present study, earlier studies concentrated on plasmas where the Fermi degeneracy parameter tended either to zero or to infinity, thus excluding many actual plasmas (stellar interiors, heavy planets such as Jupiter, plasmas in projected inertial confinement fusion experiments, and the liquid metals).

Tanaka, S.↗

DOE OSTI2025

Unified nonparametric equation-of-state inference from the neutron-star crust to perturbative-QCD densities

Perturbative quantum chromodynamics (pQCD), while valid only at densities exceeding those found in the cores of neutron stars, could provide constraints on the dense-matter equation of state (EOS). Here, in this work, we examine the impact of pQCD information on the inference of the EOS using a nonparametric framework based on Gaussian processes (GPs). We examine the application of pQCD constraints through a ``pQCD likelihood,'' and verify the findings of previous works; namely, a softening of the EOS at the central densities of the most massive neutron stars and a reduction in the maximum neutron-star mass. Although the pQCD likelihood can be easily integrated into existing EOS inference frameworks, this approach requires an arbitrary selection of the density at which the constraints are applied. The EOS behavior is also treated differently on either side of the chosen density. To mitigate these issues, we extend the EOS model to higher densities, thereby constructing a ``unified'' description of the EOS from the neutron-star crust to densities relevant for pQCD. In this approach the pQCD constraints effectively become part of the prior. Since the EOS is unconstrained by any calculation or data between the densities applicable to neutron stars and pQCD, we argue for maximum modeling flexibility in that regime. We compare the unified EOS with the traditional pQCD likelihood, and although we confirm the EOS softening, we do not see a reduction in the maximum neutron-star mass or any impact on macroscopic observables. Though residual model dependence cannot be ruled out, we find that pQCD suggests the speed of sound in the densest neutron-star cores has already started decreasing toward the asymptotic limit; we find that the speed of sound squared at the center of the most massive neutron star has an upper bound of $\sim 0.5$ at the 90% level.

equations of state of nuclear matter↗

DOE OSTI2025

General features of the stellar matter equation of state from microscopic theory, new maximum-mass constraints, and causality

The profile of a neutron star probes a very large range of densities, from the density of iron up to several times the density of saturated nuclear matter, and thus no theory of hadrons can be considered reliable if extended to those regions. We emphasize the importance of taking contemporary ab initio theories of nuclear and neutron matter as the baseline for any extension method, which will unavoidably involve some degree of phenomenology. We discuss how microscopic theory, on the one end, with causality and maximum-mass constraints, on the other, set strong boundaries to the high-density equation of state. We present our latest neutron star predictions where we combine polytropic extensions and parametrizations guided by speed of sound considerations. The predictions we show include our baseline neutron star cooling curves.

chiral effective field theory↗

DOE OSTI2020

Role of opacity at the 9 keV back lighter energy used in measuring the equation of state of boron at pressures up to a Gbar

Many experiments have been conducted at the National Ignition Facility to measure the equation of state Hugoniot of plastic, boron, and diamond at extreme pressures up to a Gbar. The “Gbar” design employs a strong spherically converging shock launched through a solid ball of material using a hohlraum radiation drive. The shock front conditions are characterized using X-ray radiography, typically at energies near 9 keV. Here we examine how the opacity of boron at 9 keV changes at high pressures and temperatures. Understanding this is vital to unfolding the density in the shock front as pressures exceed 100 Mbar. We compare opacity calculations from a number of methods including the legacy XSN opacity tables, super transition array (STA), average atom (AVAT), and detailed configuration accounting (DCA) methods. We examine how the changing opacity is correlated with the K-shell occupation calculated using the opacity methods as well as other electronic structure methods such as MECCA, Purgatorio, and an all-electron Kohn-Sham density functional theory calculation that uses optimized all-electron norm-conserving Vanderbilt pseudopotentials (ONCV). We also examine the relative contribution of free-free opacity to the overall opacity, which is dominated by bound-free absorption.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Ti-6Al-4V to over 1.2 TPa: Shock Hugoniot experiments, ab initio calculations, and a broad-range multiphase equation of state

Titanium alloys are used in a large array of applications. In this work we focus our attention on the most used alloy, Ti-6Al-4V (Ti64), which has excellent mechanical and biocompatibility properties with applications in aerospace, defense, biomedical, and other fields. Here we present high-fidelity experimental shock compression data measured on Sandia’s Z machine. We extend the principal shock Hugoniot for Ti64 to more than threefold compression, up to over 1.2 TPa. We use the data to validate our ab initio molecular dynamics simulations and to develop a highly reliable, multiphase equation of state (EOS) for Ti64, spanning a broad range of temperature and pressures. The first-principles simulations show very good agreement with Z data and with previous three-stage gas gun data from Sandia’s STAR facility. The resulting principal Hugoniot and the broad-range EOS and phase diagram up to 10 TPa and 10 5 K are suitable for use in shock experiments and in hydrodynamic simulations. The high-precision experimental results and high-fidelity simulations demonstrate that the Hugoniot of the Ti64 alloy is stiffer than that of pure Ti and reveal that Ti64 melts on the Hugoniot at a significantly lower pressure and temperature than previously modeled.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Reconstructing the neutron star equation of state from observational data via automatic differentiation

Neutron star observables like masses, radii, and tidal deformability are direct probes to the dense matter equation of state (EoS). Here, a novel deep learning method that optimizes an EoS in the automatic differentiation framework of solving inverse problems is presented. The trained neural network EoS yields narrow bands for the relationship between the pressure and speed of sound as a function of the mass density. The results are consistent with those obtained from conventional approaches and the observational bound on the tidal deformability inferred from the gravitational wave event, GW170817.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Using density functional theory to construct multiphase equations of state: tin as an example

We present a general methodology for using density functional theory (DFT) calculations as a basis for constructing multiphase equations of state (EOS). Focusing on tin as an example, we discuss the full process of generating a tabular EOS, starting from first-principles calculations and arriving at a full EOS in agreement with existing experimental data. We begin by describing DFT calculations for the five solid phases and liquid phase of tin, including cold curve, phonon, and DFT-based molecular dynamics calculations. We then discuss methods for incorporating DFT results into materials models used in OpenSesame, the program used to create a full tabular multiphase EOS. Next, we outline a general strategy for adjusting model parameters in OpenSesame to fit to a variety of experimental data, including isobaric data, isothermal data, shock data, solid-solid phase boundary measurements, and measurements of the melt curve. We end with a discussion of the advantages of using DFT data as a foundation for EOS generation and discuss potential strategies for automated EOS generation based on the fitting process highlighted in this work.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1980

Anorthite - Thermal equation of state to high pressures

New shock-wave data are presented for anorthite from which a full high-temperature, high-pressure equation of state is derived. Whereas anorthite has relatively low values of thermal expansion and Grueneisen parameter at zero pressure, it is found that these attain relatively high values in the high density state corresponding to the high-pressure phase Hugoniot but decrease upon compression as expected. It is noted that higher order anharmonic contributions decrease more rapidly with pressure and that the thermal expansion therefore saturates to a high temperature value at pressures above about 100 GPa. Reduction of the Hugoniot data permits shock temperatures to be calculated; it also yields a principal adiabat for the high pressure branch of the Hugoniot. The initial bulk modulus of this adiabat is essentially identical to that of anorthite, whereas the initial density is about 3.40 Mg/cu m.

Jeanloz, R.↗

NASA NTRS1986

Temperature effects on the universal equation of state of solids

Recently it has been argued based on theoretical calculations and experimental data that there is a universal form for the equation of state of solids. This observation was restricted to the range of temperatures and pressures such that there are no phase transitions. The use of this universal relation to estimate pressure-volume relations (i.e., isotherms) required three input parameters at each fixed temperature. It is shown that for many solids the input data needed to predict high temperature thermodynamical properties can be dramatically reduced. In particular, only four numbers are needed: (1) the zero pressure (P=0) isothermal bulk modulus; (2)it P=0 pressure derivative; (3) the P=0 volume; and (4) the P=0 thermal expansion; all evaluated at a single (reference) temperature. Explicit predictions are made for the high temperature isotherms, the thermal expansion as a function of temperature, and the temperature variation of the isothermal bulk modulus and its pressure derivative. These predictions are tested using experimental data for three representative solids: gold, sodium chloride, and xenon. Good agreement between theory and experiment is found.

Vinet, P.↗

NASA NTRS1987

Temperature effects on the universal equation of state of solids

Recently it has been argued based on theoretical calculations and experimental data that there is a universal form for the equation of state of solids. This observation was restricted to the range of temperatures and pressures such that there are no phase transitions. The use of this universal relation to estimate pressure-volume relations (i.e., isotherms) required three input parameters at each fixed temperature. It is shown that for many solids the input data needed to predict high temperature thermodynamical properties can be dramatically reduced. In particular, only four numbers are needed: (1) the zero pressure (P = 0) isothermal bulk modulus; (2) its P = 0 pressure derivative; (3) the P = 0 volume; and (4) the P = 0 thermal expansion; all evaluated at a single (reference) temperature. Explicit predictions are made for the high temperature isotherms, the thermal expansion as a function of temperature, and the temperature variation of the isothermal bulk modulus and its pressure derivative. These predictions are tested using experimental data for three representative solids: gold, sodium chloride, and xenon. Good agreement between theory and experiment is found.

Vinet, Pascal↗

DOE OSTI2026

Calculating shock Hugoniot and isentropes using multiphase equation of state tables and application to shock and release of diamond ablators in inertial confinement fusion implosions

Advances in shock and ramp compression techniques now allow experimental access to unprecedented extreme conditions of pressure and temperature, providing a means to test theoretical models. Here, we describe a simple methodology to compute multi-phase shock Hugoniot and isentropes using multiphase equation of state tables. We treat explicitly the phase coexistence along the phase boundary to reveal the evolution of the sample as it undergoes the phase transformation in adiabatic conditions. We illustrate the method by calculating the predicted shock and shock-and-release behavior of diamond at conditions relevant for the initial stage of inertial confinement fusion implosions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗