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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Advancing material modeling in hydrocodes using a concurrent finite-element and molecular dynamics multiscale framework

We present a multiscale simulation framework that couples the finite-element method with molecular dynamics. Bypassing traditional equations of state (EOS) by using in-line atomistic simulations, the method offers the advantage of incorporating detailed microscale physics not easily represented with coarse-grained models. Coupling consistency with the continuum code is ensured through the use of lifting and restriction operators, in line with heterogeneous multiscale methods. The concurrent continuum-atomistic framework is validated through comparison with experimental results and conventional EOS models, and demonstrated in a shock-driven hydrodynamic flow simulation under extreme conditions. We further evaluate the framework's usability by comparing it to state-of-the-art EOS models of deuterium. A computational performance study reveals that the atomistic EOS evaluation is a feasible alternative to conventional approaches, and demonstrates a weak scaling of 99% efficiency. These results highlight the framework's potential for large-scale multiscale modeling across a broad range of materials and conditions.

Computer science↗
DOE OSTI2020

High-pressure polymorphism in pyridine

Single crystals of the high-pressure phases II and III of pyridine have been obtained by in situ crystallization at 1.09 and 1.69 GPa, revealing the crystal structure of phase III for the first time using X-ray diffraction. Phase II crystallizes in P 2 1 2 1 2 1 with Z ′ = 1 and phase III in P 4 1 2 1 2 with Z ′ = ½. Neutron powder diffraction experiments using pyridine-d 5 establish approximate equations of state of both phases. The space group and unit-cell dimensions of phase III are similar to the structures of other simple compounds with C 2v molecular symmetry, and the phase becomes stable at high pressure because it is topologically close-packed, resulting in a lower molar volume than the topologically body-centred cubic phase II. Phases II and III have been observed previously by Raman spectroscopy, but have been mis-identified or inconsistently named. Raman spectra collected on the same samples as used in the X-ray experiments establish the vibrational characteristics of both phases unambiguously. The pyridine molecules interact in both phases through CH...π and CH...N interactions. The nature of individual contacts is preserved through the phase transition between phases III and II, which occurs on decompression. A combination of rigid-body symmetry mode analysis and density functional theory calculations enables the soft vibrational lattice mode which governs the transformation to be identified.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Negative linear compressibility in Se at ultra-high pressure above 120 GPa

A series of in situ synchrotron X-ray diffraction (XRD) measurements were carried out, combined with first-principles calculations, to study structural phase transitions of selenium at high pressures and room temperature. Several phase transitions were observed, among which an isostructural phase transition was found at around 120 GPa for the first time. Evolved from the rhombohedral (space group R$\bar{3}$m) structure (Se-V), the new phase (Se-V') exhibited an interesting increase of lattice parameter a at pressures from 120 to 148 GPa, known as negative linear compressibility (NLC). The discovery of NLC behavior observed in this work is mainly attributed to the accuracy and fine steps controlled by the membrane system for in situ XRD data collected with an exposure time of 0.5 s. After 140 GPa, a body-centered cubic (b.c.c.) structure Se-VI (space group Im$\bar{3}$m) was formed, which remains stable up to 210 GPa, the highest pressure achieved in this study. The bulk moduli of phases Se-V, Se-V' and Se-VI were estimated to be 83 ± 2, 321 ± 2 and 266 ± 7 GPa, respectively, according to the P–V curve fit by the third-order Birch–Murnaghan equation of state. The Se-V' phase shows a bulk modulus almost 4 times larger than that of the Se-V phase, which is mainly due to the effect of its NLC. NLC in a higher pressure range is always more significant in terms of fundamental mechanism and new materials discovery, yet it has barely been reported at pressures above 100 GPa. This will hopefully inspire future studies on potential NLC behaviors in other materials at ultra-high pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Greybox Thermal Parameter Identification of Electric Machine Stators

The parameters of electric machine thermal equivalent circuit networks are difficult to predict due to material and manufacturing uncertainties. In this paper, a Greybox system identification approach is used to identify parameters of electric machine stator lumped parameter thermal networks (LPTNs). LPTNs provide a low order, computationally efficient, dynamic model of temperatures at specific locations. Second and third order LPTN model structures are defined as state space equations with stator thermal parameters to be identified. To test the Greybox electric machine stator thermal system identification, five stator motorette prototypes were constructed with controlled variations in slot fill and slot liner thickness. The variation in the motorette thermal parameters and thermal time constants are detected using the Greybox identification. Special attention is given to the impact of sampling rate and Greybox data record length on parameter estimation accuracy.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

A High-Fidelity Modeling Tool to Support the Design of Oxy-Combustors for Direct-Fired Supercritical CO2 Cycles

The challenge in the design of oxy-combustors for direct-fired supercritical CO 2 (sCO 2 ) cycles is in addressing disparate performance metrics and objectives. Key design parameters to consider include, among others, injector design for mixing and flame stability, split of recycled CO 2 diluent between injectors and cooling films, target flame temperatures to control noncondensable products, and strategies to inject the diluent CO 2 for film cooling and thermal control. In order to support novel oxy-combustor designs, a high-fidelity yet numerically efficient modeling framework based on the CRUNCH CFD® flow solver is presented, featuring key physics-based submodels relevant in this regime. For computational efficiency in modeling large kinetic sets, a flamelet/progress variable (FPV) based tabulated-chemistry approach is utilized featuring a three-stream extension to allow for the simulation of the CO 2 film cooling stream in addition to the fuel and oxidizer streams. Finite rate chemistry effects are modeled in terms of multiple progress variables for the primary flame as well as for slower-evolving chemical species such as NO x and SO x contaminants. Real fluid effects are modeled using advanced equations of states. The predictive capabilities of this computationally tractable design support tool are demonstrated on a conceptual injector design for an oxy-combustor operating near 30 MPa. Simulations results provide quantitative feedback on the effectiveness of the film cooling as well as the level of contaminants (CO, NO, and N) in the exhaust due to impurities entering from the injectors. These results indicate that this framework would be a useful tool for refining and optimizing the oxy-combustor designs as well as risk mitigation analyses.

Engineering↗
DOE OSTI2020

Dealing with Measurement Uncertainties as Nuisance Parameters in Bayesian Model Calibration

In the presence of model discrepancy, the calibration of physics-based models for physical parameter inference is a challenging problem. Lack of identifiability between calibration parameters and model discrepancy requires additional identifiability constraints to be placed on the model discrepancy to obtain unique physical parameter estimates. If these assumptions are violated, the inference for the calibration parameters can be systematically biased. In many applications, such as in dynamic material property experiments, many of the calibration inputs refer to measurement uncertainties. Here, we develop a metric for identifying overfitting of these measurement uncertainties, propose a prior capable of reducing this overfitting, and show how this leads to a diagnostic tool for validation of physical parameter inference. The approach is demonstrated for a benchmark example and applied for a material property application to perform inference on the equation of state parameters of tantalum.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Elastic Bayesian Model Calibration

Functional data are ubiquitous in scientific modeling. For instance, quantities of interest are modeled as functions of time, space, energy, density, etc. Uncertainty quantification methods for computer models with functional response have resulted in tools for emulation, sensitivity analysis, and calibration that are widely used. However, many of these tools do not perform well when the computer model’s parameters control both the amplitude variation of the functional output and its alignment (or phase variation). This paper introduces a framework for Bayesian model calibration when the model responses are misaligned functional data. The approach generates two types of data out of the misaligned functional responses: (1) aligned functions so that the amplitude variation is isolated and (2) warping functions that isolate the phase variation. These two types of data are created for the computer simulation data (both of which may be emulated) and the experimental data. The calibration approach uses both types so that it seeks to match both the amplitude and phase of the experimental data. The framework is careful to respect constraints that arise, especially when modeling phase variation, and is framed in a way that it can be done with readily available calibration software. In conclusion, we demonstrate the techniques on two simulated data examples and on two dynamic material science problems: a strength model calibration using flyer plate experiments and an equation of state model calibration using experiments performed on the Sandia National Laboratories’ Z-machine.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Thermodynamics conditions of matter in the neutrino decoupling region during neutron star mergers

In this work we investigate the thermodynamics conditions at which neutrinos decouple from matter in neutron star merger remnants by post-processing results of merger simulations. We find that the matter density and the neutrino energies are the most relevant quantities in determining the decoupling surface location. For mean energy neutrinos (~9, 15 and 24 MeV for ν e , $\bar{ν}$ e and ν μ,τ , respectively) the transition between diffusion and free-streaming conditions occurs around 10 11 g cm -3 for all neutrino species. Weak and thermal equilibrium freeze-out occurs deeper (several 10 12 g cm -3 ) for heavy-flavor neutrinos than for $\bar{ν}$ e and ν e (≳10 11 g cm -3 ). Decoupling temperatures are broadly in agreement with the average neutrino energies, with softer equations of state characterized by ~1 MeV larger decoupling temperatures. Neutrinos streaming at infinity with different energies come from different remnant parts. While low-energy neutrinos (~3 MeV) decouple at ρ~10 13 g cm -3 , T~10 MeV and Y e ≲0.1 close to weak equilibrium, high-energy ones (~50 MeV ) decouple from the disk at ρ~10 9 g cm -3 , T~2 MeV and Y e ≳0.25. The presence of a massive NS or a BH influences the neutrino thermalization. While in the former case decoupling surfaces are present for all relevant energies, the lower maximum density ( ≲10 12 g cm -3 ) in BH-torus systems does not allow softer neutrinos to thermalize and diffuse.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Neutrino emission from binary neutron star mergers: characterising light curves and mean energies

Abstract Neutrinos are copiously emitted by neutron star mergers, due to the high temperatures reached by dense matter during the merger and its aftermath. Neutrinos influence the merger dynamics and shape the properties of the ejecta, including the resulting r -process nucleosynthesis and kilonova emission. In this work, we analyse neutrino emission from a large sample of binary neutron star merger simulations in Numerical Relativity, covering a broad range of initial masses, nuclear equation of state and viscosity treatments. We extract neutrino luminosities and mean energies, and compute quantities of interest such as the peak values, peak broadnesses, time averages and decrease time scales. We provide a systematic description of such quantities, including their dependence on the initial parameters of the system. We find that for equal-mass systems the total neutrino luminosity (several $$10^{53}{\hbox {erg}~{\hbox {s}}^{-1}}$$ 10 53 erg s - 1 ) decreases as the reduced tidal deformability increases, as a consequence of the less violent merger dynamics. Similarly, tidal disruption in asymmetric mergers leads to systematically smaller luminosities. Peak luminosities can be twice as large as the average ones. Electron antineutrino luminosities dominate (initially by a factor of 2-3) over electron neutrino ones, while electron neutrinos and heavy flavour neutrinos have similar luminosities. Mean energies are nearly constant in time and independent on the binary parameters. Their values reflect the different decoupling temperature inside the merger remnant. Despite present uncertainties in neutrino modelling, our results provide a broad and physically grounded characterisation of neutrino emission, and they can serve as a reference point to develop more sophisticated neutrino transport schemes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Impact of symmetry energy on sound speed and spinodal decomposition in dense neutron-rich matter

Using a meta model for nuclear Equation of State (EOS) with its parameters constrained by astrophysical observations and terrestrial nuclear experiments, we examine effects of nuclear EOS especially its symmetry energy $E$ sym ($ρ$) term on the speed of sound squared $C$ $^2_s$($ρ$) and the critical density ρt where $C$ $^2_s$($ρ_t$) vanishes (indicating the onset of spinodal decomposition) in both dense neutron-rich nucleonic matter relevant for relativistic heavy-ion collisions and the cold $n$ + $p$ + $e$ + $μ$ matter in neutron stars at $β$- equilibrium. Unlike in nucleonic matter with fixed values of the isospin asymmetry $δ$, in neutron stars with a density dependent isospin profile $δ$($ρ$) determined consistently by the $β$-equilibrium and charge neutrality conditions, the $C$ $^2_s$($ρ$) almost always show a peak and then vanishes at $ρ_t$. Furthermore, the latter strongly depends on the high-density behavior of $E$ sym ($ρ$) if the skewness parameter $J$ 0 characterizing the stiffness of high-density symmetric nuclear matter (SNM) EOS is not too far above its currently known most probable value of about -190 MeV inferred from recent Bayesian analyses of neutron star observables. Moreover, in the case of having a super-soft $E$ sym ($ρ$) that is decreasing with increasing density above about twice the saturation density of nuclear matter, the $ρ$ $t$ is significantly lower than the density where the $E$ sym ($ρ$) vanishes (indicating the onset of isospin-separation instability and pure neutron matter formation) in neutron star cores.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

QLBT: a linear Boltzmann transport model for heavy quarks in a quark-gluon plasma of quasi-particles

Abstract We develop a new heavy quark transport model, QLBT, to simulate the dynamical propagation of heavy quarks inside the quark-gluon plasma (QGP) created in relativistic heavy-ion collisions. Our QLBT model is based on the linear Boltzmann transport (LBT) model with the ideal QGP replaced by a collection of quasi-particles to account for the non-perturbative interactions among quarks and gluons of the hot QGP. The thermal masses of quasi-particles are fitted to the equation of state from lattice QCD simulations using the Bayesian statistical analysis method. Combining QLBT with our advanced hybrid fragmentation-coalescence hadronization approach, we calculate the nuclear modification factor $$R_\mathrm {AA}$$ R AA and the elliptic flow $$v_2$$ v 2 of D mesons at the Relativistic Heavy-Ion Collider and the Large Hadron Collider. By comparing our QLBT calculation to the experimental data on the D meson $$R_\mathrm {AA}$$ R AA and $$v_2$$ v 2 , we extract the heavy quark transport parameter $$\hat{q}$$ q ^ and diffusion coefficient $$D_\mathrm {s}$$ D s in the temperature range of $$1-4~T_\mathrm {c}$$ 1 - 4 T c , and compare them with the lattice QCD results and other phenomenological studies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Quark stars with 2.6 $M_\odot $ in a non-minimal geometry-matter coupling theory of gravity

This work analyses the hydrostatic equilibrium configurations of strange stars in a non-minimal geometry-matter coupling (GMC) theory of gravity. Those stars are made of strange quark matter, whose distribution is governed by the MIT equation of state. The non-minimal GMC theory is described by the following gravitational action: f(R,L) = R/2 + L + σRL, where R represents the curvature scalar, L is the matter Lagrangian density, and σ is the coupling parameter. When considering this theory, the strange stars become larger and more massive. In particular, when σ = 50 km 2 , the theory can achieve the 2.6 M ⊙ , which is suitable for describing the pulsars PSR J2215+5135 and PSR J1614-2230, and the mass of the secondary object in the GW190814 event. The 2.6 M ⊙ is a value hardly achievable in General Relativity, even considering fast rotation effects, and is also compatible with the mass of PSR J0952-0607 (M = 2.35 ± 0.17 M ⊙ ), the heaviest and fastest pulsar in the disk of the Milky Way, recently measured, supporting the possible existence of strange quark matter in its composition. The non-minimal GMC theory can also give feasible results to describe the macroscopical features of strange star candidates.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Deep learning assisted jet tomography for the study of Mach cones in QGP

Abstract Mach cones are expected to form in the expanding quark-gluon plasma (QGP) when energetic quarks and gluons traverse the hot medium at a velocity faster than the speed of sound in high-energy heavy-ion collisions. The shape of the Mach cone and the associated diffusion wake are sensitive to the initial jet production location and the propagation direction of the parton shower relative to the radial flow because of the distortion caused by the collective expansion of the QGP and the large density gradient. The shape of jet-induced Mach cones and their distortions in heavy-ion collisions provide a unique and direct probe of the dynamical evolution and the equation of state of QGP. However, it is difficult to identify the Mach cone and the diffusion wake in current experimental measurements of final hadron distributions because they are averaged over all possible initial jet production locations and parton-shower propagation directions. To overcome this difficulty, we develop a deep learning assisted jet tomography which uses the full information of the final hadrons from jets to localize the initial jet production positions. This method can help to constrain the initial regions of jet production in heavy-ion collisions and enable a differential study of Mach-cones with different path lengths and orientations relative to the radial flow of the QGP in heavy-ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Correlations between maximum mass of neutron stars and the nuclear matter properties and the constraints from PSR J0740+6620 and GW190814

The correlations between the maximum mass of neutron stars (NSs) with the properties of nuclear matter at saturation density ρ 0 , i.e., the incompressibility K 0 , the symmetry energy E sym (ρ 0 ), and its slope L(ρ 0 ), the isoscalar effective mass m$^{*}_{s0}$, and the parameter f I , have been investigated using several sets of Skyrme interactions based on the momentum-dependent (MD) Skyrme interactions SAMi-J27 (soft), SAMi-J31 (stiff), and SAMi-J35 (super-stiff), as well as the momentum-independent (MI) Skyrme interaction SkT5. Only one physical quantity is varied at a time for a given set of Skyrme interactions. It is found that the maximum mass of NSs has (1) a strong positive linear correlation with K 0 for stiff and super-stiff equation of state (EOS) and no correlation with K 0 for soft EOS; (2) a negative logarithm correlation with the symmetry energy E sym (ρ 0 ) for all the EOSs; (3) a positive logarithm correlation with the slope L(ρ 0 ) for all the EOSs; (4) a positive logarithm correlation with the parameter f I for MD EOS; (5) a negative power-law correlation with m$^{*}_{s0}$ for MD EOS. Based on these established correlations, the constraints to the nuclear matter properties from the mass of PSR J0740+6620 and GW190814 are investigated. Furthermore, we find that while rotation enhances the maximum mass of NSs, it does not effect the underlying correlations between the maximum mass and the aforementioned nuclear matter parameters.

Zhou, M. [Shaanxi Normal University, Xi’an (China)↗
DOE OSTI2026

$\Sigma ^{+}$ production in pp collisions at $\sqrt{\textit{s}}=13$ TeV

The measurement of $\Sigma ^{+}$ production in pp collisions at $\sqrt{s}=13$ TeV is presented. The measurement is performed at midrapidity in both minimum-bias and high-multiplicity pp collisions at $\sqrt{s}=13$ TeV. The $\Sigma ^{+}$ is reconstructed via its weak-decay topology in the decay channel $\Sigma ^{+} \rightarrow \mathrm{{p}} + \pi ^{0}$ with $\pi ^{0} \rightarrow \gamma + \gamma$. In a novel approach, the neutral pion is reconstructed by combining photons that convert in the detector material with photons measured in the calorimeters. The transverse-momentum $(p_{\textrm{T}})$ distributions of the $\Sigma ^{+}$ and its rapidity densities dN/dy in both event classes are reported. The $p_{\textrm{T}}$ spectrum in minimum-bias collisions is compared to QCD-inspired event generators. The ratio of $\Sigma ^{+}$ to previously measured Λ baryons is in good agreement with calculations from the Statistical Hadronization Model. The high efficiency and purity of the novel reconstruction method for $\Sigma ^{+}$ presented here will enable future studies of the interaction of $\Sigma ^{+}$ with protons in the context of femtoscopic measurements, which could be crucial for understanding the equation of state of neutron stars.

Abualrob, I. J. [University of Houston] (ORCID:000↗
DOE OSTI2026

Beyond $\rho ^{2/3}$ scaling: microscopic origins and multimessengers of high-density nuclear symmetry energy

The nature and Equation of State (EOS) of dense neutron-rich matter are still very poorly known, while they have broad impacts on many interesting issues in both astrophysics and nuclear physics. In particular, the nuclear symmetry energy $E_{\textrm{sym}}(\rho )$ encoding the cost to make nuclear matter more neutron-rich has been the most uncertain component of the EOS of dense neutron-rich nucleonic matter. It significantly affects the radii, tidal deformations, cooling rates, and frequencies of various oscillation modes of isolated neutron stars as well as the strain amplitude and frequencies of gravitational waves from their mergers, besides its many effects on the structures of nuclei as well as the dynamics and observables of their collisions. Siemens (1970s) observed that $E_{\textrm{sym}}(\rho )$ scales as $(\rho /\rho _0)^{2/3}$ near the saturation density $\rho _0$ of nuclear matter, since both the kinetic part and the potential contribution (quadratic in momentum) exhibit this dependence. The scaling holds if: (1) the nucleon isoscalar potential is quadratic in momentum, and (2) the isovector interaction is weakly density-dependent. After examining many empirical evidences and understanding theoretical findings in the literature, we conclude that: (1) Siemens’ $\rho ^{2/3}$ scaling is robust and serves as a valuable benchmark for both nuclear theories and experiments up to $2\rho _0$ but breaks down at higher densities, (2) Experimental and theoretical findings about $E_{\textrm{sym}}(\rho )$ up to $2\rho _0$ are broadly consistent, but uncertainties remain large for its curvature $K_{\textrm{sym}}(\rho )$ and higher-order parameters, (3) Above $2\rho _0$ , uncertainties grow due to poorly constrained spin-isospin-dependent tensor and three-body forces as well as the resulting nucleon short-range correlations. Looking forward, combining signals from both observations of neutron stars and terrestrial heavy-ion reaction experiments is the most promising path to finally constraining the high-density $E_{\textrm{sym}}(\rho )$ and the EOS of supradense neutron-rich matter. Multiple examples of community efforts to further constrain the high-density $E_{\textrm{sym}}(\rho )$ using both real and mocked data of present and future high-precision observations of neutron stars, as well as heavy-ion collisions involving high-energy rare isotopes, are briefly reviewed.

Li, Bao-An [East Texas A & M Univ., Commerce, TX (↗
DOE OSTI2020

Review of uncertainties in the cosmic supernova relic neutrino background

We review the computation of and associated uncertainties in the current understanding of the relic neutrino background due to core-collapse supernovae, black hole formation and neutron star merger events. We consider the current status of uncertainties due to the nuclear equation of state (EoS), the progenitor masses, the source supernova neutrino spectrum, the cosmological star formation rate, the stellar initial mass function, neutrino oscillations, and neutrino self-interactions. We summarize the current viability of future neutrino detectors to distinguish the nuclear EoS and the temperature of supernova neutrinos via the detected relic supernova neutrino spectrum.

Astronomy & Astrophysics↗