Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equation 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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Investigation of warm dense matter by combining static and dynamic compression

Shock-compressing material samples that are precompressed with diamond-anvil cells (DACs) enables benchmarking equation of state measurements over a wide range of high-pressure conditions, in particular, warm dense matter relevant to planetary science, astrophysics, and inertial confinement fusion. These targets provide optical access to the shocked sample, so optical velocity interferometry and pyrometry are typical diagnostics. We present a scaling for the coupling of the laser drive to the shock pressure in DACs based on over 150 prior experiments that can be used to design future experiments. For initial static pressures ≤ 5 GPa, using a sapphire window mitigates the effect of photoionization compared to diamond windows and allows the measurement of shock pressures up to 1 TPa with optical diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Bubble Point Measurements of cis-1,1,1,4,4,4-Hexafluorobutene [R-1336mzz(Z)] + trans-1,2-Dichloroethene [R-1130(E)] mixtures

Saturation pressures of pure R-1336mzz(Z) and R-1130(E) and bubble point pressures of three R-1336mzz(Z)/1130(E) blends were measured from 265 K to 360 K. For each pure refrigerant or refrigerant blend, a total of twenty unique saturation pressures or bubble points were measured. In total 100 unique state points were obtained. Presently, no Helmholtz-energy-explicit type equation of state (EoS) is available for R-1130(E). While an extended corresponding states EoS for R-1130(E) is available to estimate the properties of the R-1336mzz(Z)/1130(E) blend, this model does not resolve the azeotropic behavior of the mixture. Therefore, the perturbed-chain statistical-associating fluid theory (PC-SAFT) EoS is used to model the vapor–liquid equilibria of the R-1336mzz(Z)/1130(E) blend. PC-SAFT model parameters are reported, and the overall performance of the model is characterized by deviations from the experimental data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Enhancement of PyARC for Westinghouse Electric Company’s Lead Fast Reactor Design and Modeling (Final TCF Report)

Westinghouse Electric Company is a nuclear reactor vendor headquartered in the U.S. that is developing advanced reactor technology for the U.S. and global markets. Westinghouse has been relying on the neutronics Argonne Reactor Codes (ARC) executed through the NEAMS Workbench and its PyARC module that are developed under the DOE-NE Nuclear Energy Advanced Modeling and Simulation (NEAMS) and Advanced Reactor Technology (ART) – Fast Reactor programs. Through this user experience, Westinghouse identified several enhancements that would benefit the ARC codes’ usability by the US industry and therefore its commercialization potential. The enhancements were proposed to deliver both improvements in workflow and analysis capabilities to better support effective fast reactor core design and analysis to the nuclear industry. The PyARC workflow was extended in this project by integrating non-neutronic ARC codes DASSH and NUBOW-3D. The Ducted Assembly Steady-State Heat equation (DASSH) code is developed at ANL to perform steady-state thermal hydraulic sub-channel analysis in liquid metal fast reactor assemblies to determine optimized coolant flow and temperature distributions, which in this project was updated and validated for lead fast reactor (LFR) applications. The interface between REBUS and NUBOW-3D were improved in this project to assess the impact of the core restraint design and thermal induced expansion effects on the reactivity of the core, and to model the deformations of the fuel assemblies induced by temperature and irradiation. Finally, the ARC models that were extensively verified and validated through various SFR-based modeling benchmarks are extended in this project through code-to-code comparison on relevant LFR-specific neutronics benchmarks against Monte-Carlo neutronic solutions. Overall, this work enables verification of the capability of the ARC codes for a wide range of Generation-IV reactor designs. The outcome of this project is the release of a comprehensive modeling toolkit of validated, robust and efficient codes, as well as their user interface, that enables industry to perform a wide range of fast reactor analyses for design and licensing of their concepts.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

On the expansion of non-ideal copper plasma into vacuum

The expansion of non-ideal copper plasma into vacuum is analyzed for the conditions typical for explosive electron emission in vacuum arcs. The computational model solves the Euler equations with an equation of state (EoS) for weakly non-ideal plasma, taking into account the ionization energy correction due to electron–ion Coulomb coupling. We have ascertained that the EoS has insignificant influence on the plasma expansion while the plasma properties and composition change drastically. Based on our simulation results, we discuss the validity of the “frozen” state theory often used in the vacuum arc plasma diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Chiral Effective Field Theory’s Impact on Advancing Quantum Monte Carlo Methods

Thirty years ago, Steven Weinberg published his seminal paper on “Nuclear Forces from chiral Lagrangians” which has revolutionized the field of theoretical nuclear physics. Nowadays, interactions derived from chiral effective field theory are routinely used to describe nuclear systems ranging from atomic nuclei to the dense matter explored in the core of neutron stars with theoretical uncertainty estimates. Here, in our contribution to the special issue “Celebrating 30 years of Steven Weinberg’s paper Nuclear Forces from Chiral Lagrangians”, we focus on the impact that chiral effective field theory interactions have played in advancing microscopic studies of atomic nuclei and the nuclear-matter equation of state using quantum Monte Carlo methods.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Reshock analysis for PMMA driven above the threshold for chemical decomposition

Polymethyl methacrylate (PMMA) is a transparent thermoplastic often used in shock compression science. Previous proceedings (AIP Conf. Proc. 845, 131 (2006)) described plate impact experiments in which PMMA was reshocked to pressures of up to ~130 GPa, well above that at which it decomposes on its principal Hugoniot. Because some of the conclusions of the original analysis were surprising, the results were reanalyzed in a later proceeding (ibid. 1426, 771 (2012)). In addition to our own reanalysis of the data, we performed hydrodynamic simulations of the experiments based on a new equation of state for PMMA shocked above its threshold for decomposition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Hydrodynamic computations of high-power laser drives generating metal ejecta jets from surface grooves

Understanding dynamic fragmentation in shock-loaded metals and predicting properties of the resulting ejecta are of considerable importance for both basic and applied science. The nature of material ejection has been shown to change drastically when the free surface melts on compression or release. In this work, we present hydrodynamic simulations of laser-driven microjetting from micron-scale grooves on a tin surface. We study microjet formation across a range of shock strengths from drives that leave the target solid after release to drives that induce shock melting in the target. The shock-state particle velocity (Up) varies from 0.3 to 3 km/s and the shock breakout pressure is 3–120 GPa. The microjet tip velocity is 1–8 km/s and the free-surface velocity varies from 0.1 to 5 km/s. Two tin equations of state are examined: a “soft” model (LEOS 501) where the target melts for U p > 1 km/s and a more detailed multiphase model (SESAME 2161) that melts for U p > 1.4 km/s. We use these two models to examine the influence of phase change and the choice of the material model on microjet formation and evolution. We observe in our computational results that jet formation can be classified into three regimes: a low-energy regime where material strength affects jet formation, a moderate-energy regime dominated by the changing phase of tin material, and a high-energy regime where results are insensitive to the material model and jet formation is described by an idealized steady-jet theory. Using an ensemble of 2D simulations, we show that these trends hold across a wide range of drive energies and groove angles.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Indirect methods in nuclear astrophysics with relativistic radioactive beams

Reactions with radioactive nuclear beams at relativistic energies have opened new doors to clarify the mechanisms of stellar evolution and cataclysmic events involving stars and during the big bang epoch. Numerous nuclear reactions of astrophysical interest cannot be assessed directly in laboratory experiments. Ironically, some of the information needed to describe such reactions, at extremely low energies (e.g., keVs), can only be studied on Earth by using relativistic collisions between heavy ions at GeV energies. In this contribution, we make a short review of experiments with relativistic radioactive beams and of the theoretical methods needed to understand the physics of stars, adding to the knowledge inferred from astronomical observations. We continue by introducing a more detailed description of how the use of relativistic radioactive beams can help to solve astrophysical puzzles and several successful experimental methods. State-of-the-art theories are discussed at some length with the purpose of helping us understand the experimental results reported. The review is not complete and we have focused most of it to traditional methods aiming at the determination of the equation of state of symmetric and asymmetric nuclear matter and the role of the symmetry energy. Whenever possible, under the limitations of our present understanding of experimental data and theory, we try to pinpoint the information still missing to further understand how stars evolve, explode, and how their internal structure might be. We try to convey the idea that in order to improve microscopic theories for many-body calculations, nuclear structure, nuclear reactions, and astrophysics, and in order to constrain and allow for convergence of our understanding of stars, we still need considerable improvements in terms of accuracy of experiments and the development of new and dedicated nuclear facilities to study relativistic reactions with radioactive beams.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Symmetry and scaling in one-dimensional compressible two-phase flow

Investigations of shock compression of heterogeneous materials often focus on the shock front width and overall profile. The number of experiments required to fully characterize the dynamic response of a material often belie the structure–property relationships governing these aspects of a shock wave. Recent observations measured a pronounced shock-front width on the order of 10 s of ns in particulate composites. We focus on particulate composites with disparate densities and investigate whether the mechanical interactions between the phases are adequate to describe this emergent behavior. The analysis proceeds with a general Mie–Grüneisen equation of state for the matrix material, a general drag force law with general power-law scaling for the particle-matrix coupling of the phases, and a volume fraction-dependent viscosity. Lie group analysis is applied to one-dimensional hydrodynamic flow equations for the self-consistent interaction of particles embedded in a matrix material. The particle phase is characterized by a particle size and volume fraction. The Lie group analysis results in self-similar solutions reflecting the symmetries of the flow. The symmetries lead to well-defined scaling laws, which may be used to characterize the propagation of shock waves in particle composites. An example of the derived scaling laws for shock attenuation and rise time is shown for experimental data on shock-driven tungsten-loaded polymers. A key result of the Lie analysis is that there is a relationship between the exponents characterizing the form of the drag force and the exponent characterizing the shock velocity and its attenuation in a particulate composite. Comparison to recent experiments results in a single exponent that corresponds to a conventional drag force.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

A description of structured waves in shock compressed particulate composites

Dynamic compression of composite materials is of scientific interest because the mechanical mismatch between internal phases challenges continuum theories. Typical assumptions about steady wave propagation and quasi-instantaneous state changes require reexamination along with the need for time-dependent models. To that end, data and models are presented here for the shock compression of an idealized particulate composite. To serve as a generic representative of this material class, a polymer matrix was filled with tungsten particles, ranging from 1 to 50 vol. %. This creates a simple microstructure containing randomly scattered particles with an extreme impedance mismatch to the binding matrix. These materials were parallel plate impact loaded by Al flyers traveling at 1.8–5.0 km/s. Velocimetry provided records of the equilibrium state and the compression wave structure for each case with trends quantified by an empirical fit. The same quantities were also studied as a function of the wave's propagation distance. A homogenized viscoelastic model then made it possible to progress from cataloging trends to making predictions. Starting from a Mie–Gruneisen equation of state, additional time varying terms were added to capture the transient response. After calibration, accurate predictions of the steady wave structure were possible.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

The high explosives & affected targets (HEAT) dataset

Artificial Intelligence (AI) surrogate models offer a computationally efficient alternative to full-physics simulations, yet no existing datasets are publicly available for training, testing, and validation of machine learning models of the dynamics of high-explosive driven shocks through multiple materials. Shock propagation through materials is a computationally challenging problem because simulations must include material-specific equations of state (EOS) along with descriptions of other physical processes such as plastic deformation, phase change, damage processes, fluid instabilities, and multi-material interactions. Shocks are typically initiated by high-velocity impacts or explosive loading. The latter case necessitates the addition of models of reactive materials to represent high-explosive (HE) detonation. Here, to address the lack of an expansive dataset for multi-material shock propagation in the AI/ML community, we present the High-Explosives and Affected Targets (HEAT) Dataset. HEAT is a physics-rich collection of two-dimensional, cylindrically symmetric, simulations generated using an Eulerian, multi-material, shock-propagation code developed at Los Alamos National Laboratory. The dataset includes two partitions: (1) the expanding shock-cylinder (CYL) simulations, Figs. 1, and (2) the Perturbed Layered Interface (PLI) simulations, Fig. 2. Entries in both partitions consist of time series of arrays of thermodynamic fields (pressure, density, and temperature), kinematic fields (position and velocity), and additional fields that depend on thermodynamic and/or kinematic fields (e.g., material stress). Materials in the CYL partition include solids (aluminium, copper, depleted uranium, stainless steel, tantalum, and a generic polymer), a liquid (water), gases (air, nitrogen), and a generic detonating material (high explosive, HE). The PLI partition spans a highly varying geometry but consists of fixed materials across entries: Copper, aluminium, stainless steel, generic polymer, and generic HE. HEAT captures critical phenomena such as momentum transfer, shock propagation, plastic deformation, and thermal effects, making HEAT a valuable benchmark for development of AI/ML emulation of multi-material shock propagation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Bayesian inference of nuclear incompressibility from collective flow in mid-central Au+Au collisions at 400–1500 MeV/nucleon

The incompressibility K of symmetric nuclear matter (SNM) is determined through a Bayesian analysis of collective flow data from Au + Au collisions at beam energies $E = 400 -1500$ MeV/nucleon. This analysis utilizes a Gaussian process (GP) emulator applied to the isospin-dependent quantum molecular dynamics (IQMD) model for heavy-ion collisions, both with and without incorporating the momentum dependence of the single-nucleon potentials. Specifically, at the 68% confidence level, using rapidity and transverse velocity dependence of proton elliptic flow data with and without consideration of the momentum dependence, the inferred incompressibility values are $K=188.9^{+2.9}_{-4.5}$ MeV and $256.1^{+8.2}_{-8.7}$ MeV at $E = 400$ MeV/nucleon, respectively. When the transverse momentum dependence of proton-like directed flow data is included, the inferred incompressibility values become $K=222.3^{+9.0}_{-9.9}$ MeV and $K=285.5^{+6.7}_{-7.3}$ MeV, respectively. Furthermore, we found that the value of K derived from observables of proton elliptic flow increases with beam energy. Finally, this indicates that the equation of state (EoS) of nuclear matter hardens at higher densities and temperatures in reactions with higher beam energies.

Bayesian inference↗
DOE OSTI2022

Constraining the release of Sn to the ambient melting point following shock loading using time-resolved x-ray diffraction

The formation of liquid following release from a shocked state governs the transition from spall to cavitation and the formation of ejecta in metals. In order to build physics-based models of these processes, it is necessary to critically evaluate the relative importance of kinetics and entropy generation during the release along with the accuracy of multiphase equations of state. Tin (Sn) has served as a testbed for a variety of experiments examining strength and ejecta due to its accessible melt boundary and solid–solid phase transitions. This work presents experiments examining the phase evolution of high purity Sn following the shock and release to ambient stress near the melting point. Sn is found to release to states between its ambient solidus and liquidus from approximately 19 to 33 GPa under uniaxial loading, with the two-phase region being characterized by a reduction in the intensity of the (220), (211) β -Sn doublet. Jetting experiments performed at 27–28 GPa exhibit comparable diffraction patterns with what is observed following the uniaxial release. The solid fractions of β-Sn in the ambient mixed phase region are found to decrease linearly with increasing shock stress as increasing liquid Sn is formed. Here, the results provide much needed information for interpreting measurements of dynamic strength at a high strain rate and experiments examining cavitation and shallow bubble collapse in Sn.

, spall↗
DOE OSTI2025

Group Conductivity and Nonadiabatic Born Effective Charges of Disordered Metals, Warm Dense Matter, and Hot Dense Plasma

The average ionization state is a critical parameter in plasma models for charged particle transport, equations of state, and optical response. The dynamical or nonadiabatic Born effective charge (NBEC), calculated via first principles time-dependent density-functional theory, provides exact ionic partitioning of bulk electron response for both metallic and insulating materials. The NBEC can be transformed into a “group conductivity,” i.e., the electron conductivity ascribed to a subset of ions. We show that for disordered metallic systems, such as warm dense matter (WDM) and hot dense plasma, the static limit of the NBEC is different from the average ionization states, but that the ionization state can be extracted from the group conductivity even in mixed systems. Here we demonstrate this approach using a set of archetypical examples, including cold and warm aluminium, low- and high-density WDM carbon, and a WDM carbon-beryllium-hydrogen mixture.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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

New insights into supradense matter from dissecting scaled stellar structure equations

The strong-field gravity in general relativity (GR) realized in neutron stars (NSs) renders the equation of state (EOS) P(ε) of supradense neutron star matter to be essentially nonlinear and refines the upper bound for Φ ≡ P/ε to be much smaller than the special relativity (SR) requirement with linear EOSs, where P and ε are respectively the pressure and energy density of the system considered. Specifically, a tight bound Φ ≲ 0.374 is obtained by perturbatively anatomizing the intrinsic structures of the scaled Tolman–Oppenheimer–Volkoff (TOV) equations without using any input nuclear EOS. New insights gained from this novel analysis provide EOS-model-independent constraints on the properties (e.g., density profiles of the sound speed squared s 2 = dP/dε and trace anomaly Δ = 1/3 – Φ) of cold supradense matter in NS cores. Using the gravity-matter duality in theories describing NSs, we investigate the impact of gravity on supradense matter EOS in NSs. In particular, we show that the NS mass M NS , radius R, and compactness ξ ≡ M NS /R scale with certain combinations of its central pressure and energy density (encapsulating its central EOS). Thus, observational data on these properties of NSs can straightforwardly constrain NS central EOSs without relying on any specific nuclear EOS model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Isotope effects in supercooled H 2 O and D 2 O and a corresponding-states-like rescaling of the temperature and pressure

Water shows anomalous properties that are enhanced upon supercooling. The unusual behavior is observed in both H 2 O and D 2 O, however, with different temperature dependences for the two isotopes. It is often noted that comparing the properties of the isotopes at two different temperatures (i.e., a temperature shift) approximately accounts for many of the observations—with a temperature shift of 7.2 K in the temperature of maximum density being the most well-known example. However, the physical justification for such a shift is unclear. Motivated by recent work demonstrating a “corresponding-states-like” rescaling for water properties in three classical water models that all exhibit a liquid–liquid transition and critical point, the applicability of this approach for reconciling the differences in the temperature- and pressure-dependent thermodynamic properties of H 2 O and D 2 O is investigated here. Utilizing previously published data and equations-of-state for H 2 O and D 2 O, we show that the available data and models for these isotopes are consistent with such a low temperature correspondence. These observations provide support for the hypothesis that a liquid–liquid critical point, which is predicted to occur at low temperatures and high pressures, is the origin of many of water’s anomalies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Simultaneous compression and opacity data from time-series radiography with a Lagrangian marker

Time-resolved radiography can be used to obtain absolute shock Hugoniot states by simultaneously measuring at least two mechanical parameters of the shock, and this technique is particularly suitable for one-dimensional converging shocks where a single experiment probes a range of pressures as the converging shock strengthens. However, at sufficiently high pressures, the shocked material becomes hot enough that the x-ray opacity falls significantly. Additionally, if the system includes a Lagrangian marker such that the mass within the marker is known, this additional information can be used to constrain the opacity as well as the Hugoniot state. In the limit that the opacity changes only on shock heating, and not significantly on subsequent isentropic compression, the opacity of the shocked material can be determined uniquely. More generally, it is necessary to assume the form of the variation of opacity with isentropic compression or to introduce multiple marker layers. Alternatively, assuming either the equation of state or the opacity, the presence of a marker layer in such experiments enables the non-assumed property to be deduced more accurately than from the radiographic density reconstruction alone. An example analysis is shown for measurements of a converging shock wave in polystyrene at the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗