SEARCH · Engineering Papers

Results for “EOS”

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 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Operando TEM study of a working copper catalyst during ethylene oxidation

Active catalysts are typically metastable, and their surface state depends on the gas-phase chemical potential and reaction kinetics. To gain relevant insights into structure-performance relationships, it is essential to investigate catalysts under their operational conditions. Here, we use operando TEM combining real-time observations with online mass spectrometry (MS) to study a Cu catalyst during ethylene oxidation. We identify three distinct regimes characterized by varying structures and states that show different selectivities with temperature, and elucidate the reaction pathways with the aid of theoretical calculations. Our findings reveal that quasi-static Cu 2 O at low temperatures is selective towards ethylene oxide (EO) and acetaldehyde (AcH) via an oxometallacycle (OMC) pathway. In the dynamic Cu 0 /Cu 2 O oscillation regime at medium temperatures, partially reduced and strained oxides decrease the activation energies associated with partial oxidation. At high temperatures, the catalyst is predominantly Cu 0 , partially covered by a monolayer Cu 2 O. While Cu 0 is extremely efficient in dehydrogenation and eventual combustion, the monolayer oxide favors direct EO formation. These results challenge conclusions drawn from ultra-high vacuum studies that suggested metallic copper would be a selective epoxidation catalyst and highlight the need for operando study under realistic conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Bayesian inference of multi-messenger astrophysical data: Joint and coherent inference of gravitational waves and kilonovae

Multi-messenger observations of binary neutron star mergers can provide information on the neutron star’s equation of state (EOS) above the nuclear saturation density by directly constraining the mass-radius diagram. We present a Bayesian framework for joint and coherent analyses of multi-messenger binary neutron star signals. As a first application, we analyze the gravitational-wave GW170817 and the kilonova (kN) AT2017gfo data. These results are then combined with the most recent X-ray pulsar analyses of PSR J0030+0451 and PSR J0740+6620 to obtain new EOS constraints.We extend the bajes infrastructure with a joint likelihood for multiple datasets, support for various semi-analytical kN models, and numerical-relativity (NR)-informed relations for the mass ejecta, as well as a technique to include and marginalize over modeling uncertainties. The analysis of GW170817 used the TEOBResumS effective-one-body waveform template to model the gravitational-wave signal. The analysis of AT2017gfo used a baseline multicomponent spherically symmetric model for the kN light curves. Various constraints on the mass-radius diagram and neutron star properties were then obtained by resampling over a set of ten million parameterized EOSs, which was built under minimal assumptions (general relativity and causality).

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

QCD material parameters at zero and non-zero chemical potential from the lattice

Using an eighth-order Taylor expansion in baryon-chemical potential, we recently obtained the (2+1) flavor QCD equation of state (EoS) at non zero conserved charge chemical potentials from the lattice. We focused on strangeness-neutral, isospin-symmetric QCD matter, which closely resembles the situation encountered in heavy-ion collision experiments. Using this EoS, we present here results on various QCD material parameters; in particular we compute the specific heat, speed of sound, and compressibility along appropriate lines of constant physics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Equilibrium expectations for non-Gaussian fluctuations near a QCD critical point

With the highly anticipated results from the Beam Energy Scan II program at RHIC being recently revealed, an understanding of particle-number fluctuations and their significance as a potential signature of a possible QCD critical point is crucial. Early works that embarked on this endeavor sought to estimate the fluctuations due to the presence of a critical point assuming they stay in equilibrium. From these results came the proposal to focus efforts on higher, non-Gaussian, moments of the event-by-event distributions, in particular of the number of protons. These non-Gaussian moments are especially sensitive to critical fluctuations, as their magnitudes are proportional to high powers of the critical correlation length. As the equation of state provides key input for hydrodynamical simulations of heavy-ion collisions, we estimate equilibrium fluctuations from the BEST equation of state (EoS) that includes critical features from the 3D Ising Model. In particular, the proton factorial cumulants and their dependence on non-universal mapping parameters is investigated within the BEST EoS. Furthermore, the correlation length, as a central quantity for the assessment of fluctuations in the vicinity of a critical point, is also calculated in a consistent manner with the scaling equation of state. An understanding of the equilibrium estimates of proton factorial cumulants will be useful for further comparison to estimates of out-of-equilibrium fluctuations in order to determine the magnitude of the observable fluctuations to be expected in heavyion collision experiments, in which the time spent near a critical point is short.

Karthein, Jamie M. [Massachusetts Institute of Tec↗

DOE OSTI2024

Empirical correlations between the function of entropy ( Z S ) and net artificial viscous work in a shock physics hydrocode

Entropy is a state variable that may be obtained from any thermodynamically complete equation of state (EOS). However, hydrocode calculations that output the entropy often contain numerical errors; this is not because of the EOS, but rather the solution techniques that are used in hydrocodes (especially Eulerian) such as convection, remapping, and artificial viscosity. Here, in this work, empirical correlations are investigated to reduce the errors in entropy without altering the solution techniques for the conservation of mass, momentum, and energy. Specifically, these correlations are developed for the function of entropy Z S , and they depend upon the net artificial viscous work, as determined via Sandia National Laboratories’ shock physics hydrocode CTH. These results are a continuation of a prior effort to implement the entropy-based CREST reactive burn model in CTH, and they are presented here to stimulate further interest from the shock physics community. Future work is planned to study higher-dimensional shock waves, shock wave interactions, and possible ties between the empirical correlations and a physical law.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Stellarator fusion systems enabled by arrays of planar coils

We present an overview of a novel electromagnetic coil configuration for stellarators and its application to two near-term fusion systems. The novel coil configuration is the planar coil stellarator, able to implement precisely-quasisymmetric 3D magnetic fields using a set of planar, plasma-encircling coils and a set of planar, field-shaping coils situated on a surface surrounding the plasma. This configuration combines the stellarator’s advantages of steady-state operation, stability, low recirculating power fraction, and a mature physics basis, with the benefits of using simpler, planar coils which allow for a maintenance scheme leveraging large ports, and the ability to control magnets individually. The initial near-term use case considered is a steady state deuterium–deuterium stellarator neutron source, called Eos. The second near-term use case considered is a deuterium–tritium stellarator fusion pilot plant, called Helios, that would be approximately twice the linear dimension of the Eos design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Fast ion confinement in quasi-axisymmetric stellarator equilibria

This report presents an initial analysis of the fast ion confinement and losses within quasi-axisymmetric stellarator equilibria in consideration by Thea Energy. The equilibria have not yet been explicitly optimized for fast particle confinement and require validation. Modeling with the ASCOT5 code is used to directly examine the fast ion transport. The particle tracking simulations are purely (neo)classical in nature and simply contain the supplied equilibrium and collisions (pitch-angle, energy slowing, and velocity diffusion) from supplied thermal profiles. Uniform marker deposition is used to probe the general confinement properties of the equilibria while a realistic beam-born population is provided from the BEAMS3D code and an alpha particle population is calculated from a fusion source integrator. A first wall is included and defines the loss boundary. Analysis for NBI ions within Thea Energy’s conceptual Eos neutron source are presented along with alpha particles in an enlarged DT-plasma. The fast ions are assessed in regards to their confinement time, pitch, energy, and spatial coordinates. For each population, the impact of collisions and orbit drifts are discussed. It is found that NBI born ions in Eos are strongly confined until slowing-down, owing largely to the tangential injection geometry, while 22% of DT-born alpha energy is lost in the scaled device, indicating that any fusion pilot plant design optimization should include metrics for fast ion confinement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Exploring the role of 𝑑* hexaquarks on quark deconfinement and hybrid stars

We investigate the impact of the d*( 2380) hexaquark on the equation of state (EOS) of dense matter within hybrid stars (HSs) using the chiral mean-field model (CMF). The hexaquark is included as a new degree of freedom in the hadronic phase, and its influence on the deconfinement transition to quark matter is explored. We reparametrize the CMF model to ensure compatibility with recent astrophysical constraints, including the observation of massive pulsars and gravitational wave events. Our results show that the presence of d* significantly modifies the EOS, leading to a softening at high densities and a consequent reduction in the predicted maximum stellar masses. Furthermore, we examine the possibility of a first-order deconfinement phase transition within the context of the extended stability branch of slow stable HSs (SSHSs). We find that the presence of hexaquarks can delay the deconfinement phase transition and reduce the associated energy density gap, affecting the structure and stability of HSs. Our results suggest that, as the hexaquark appearance tends to destabilize stellar configurations, fine-tuning of model parameters is required to obtain both the presence of hexaquarks and quark deconfinement in these systems. In this scenario, the SSHS branch plays a crucial role in obtaining HSs with hexaquarks that satisfy current astrophysical constraints. Our work provides new insights into the role of exotic particles like d* in dense matter and the complex interplay between hadronic and quark degrees of freedom inside compact stellar objects.

Nuclear astrophysics↗

DOE OSTI2024

Sound speed and Grüneisen parameter up to three terapascal in shock-compressed iron

This paper presents the first sound speed and Grüneisen parameter data for fluid iron compressed to 3 TPa (30 million atmospheres) and 20 g/cm 3 on the Hugoniot. Both the sound speed and Grüneisen parameter are derivatives of the equation of state (EOS), and thus tightly constrain the contours of the EOS surface. The sound speed data are systematically lower than expected from a simple extrapolation of previous data. The Grüneisen parameter shows a 30% drop at pressures and temperatures above the melt transition. Furthermore, while some models compare well with either the sound speed or Grüneisen parameter, none of today’s state-of-the-art models can explain both sets of data. Furthermore these new data will provide pivotal benchmarks for both future theoretical EOSs of warm dense iron and modeling planetary states and processes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Bayesian mixture model approach to quantifying the empirical nuclear saturation point

The equation of state (EOS) in the limit of infinite symmetric nuclear matter exhibits an equilibrium density, $n_0 \approx 0.16 \, \mathrm{fm}^{-3}$, at which the pressure vanishes and the energy per particle attains its minimum, $E_0 \approx -16 \, \mathrm{MeV}$. Although not directly measurable, the nuclear saturation point $(n_0,E_0)$ can be extrapolated by density functional theory (DFT), providing tight constraints for microscopic interactions derived from chiral effective field theory (EFT). However, when considering several DFT predictions for $(n_0,E_0)$ from Skyrme and Relativistic Mean Field (RMF) models together, a discrepancy between these model classes emerges at high confidence levels that each model prediction's uncertainty cannot explain. How can we leverage these DFT constraints to rigorously benchmark nuclear saturation properties of chiral interactions? To address this question, we present a Bayesian mixture model that combines multiple DFT predictions for $(n_0,E_0)$ using an efficient conjugate prior approach. The inferred posterior distribution for the saturation point's mean and covariance matrix follows a Normal-inverse-Wishart class, resulting in posterior predictives in the form of correlated, bivariate $t$-distributions. The DFT uncertainty reports are then used to mix these posteriors using an ordinary Monte Carlo approach. At the 95\% credibility level, we estimate $n_0 \approx 0.157 \pm 0.010 \, \mathrm{fm}^{-3}$ and $E_0 \approx -15.97 \pm 0.40 \, \mathrm{MeV}$ for the marginal (univariate) $t$-distributions. Combined with chiral EFT calculations of the pure neutron matter EOS, we obtain bivariate normal distributions for the nuclear symmetry energy and its slope parameter evaluated at $n_0$: $S_v \approx 32.0 \pm 1.1 \, \mathrm{MeV}$ and $L\approx 52.6\pm 8.1 \, \mathrm{MeV}$ (95\%), respectively. Furthermore, our Bayesian framework is publicly available, so practitioners can readily use and extend our results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Four-dimensional QCD equation of state with multiple chemical potentials

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

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Evolution of in-medium baryon-baryon scattering cross sections and stiffness of dense nuclear matter from Bayesian analyses of FOPI proton-flow excitation functions

Within a Bayesian statistical framework using a Gaussian process (GP) emulator for an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model simulator of heavy-ion reactions with momentum-independent Skyrme interactions, we infer from the proton directed and elliptical flow in mid-central Au + Au reactions, measured at beam energies from 150 to 1200 MeV/nucleon by the FOPI Collaboration, the posterior probability distribution functions (PDFs) of the in-medium baryon-baryon scattering cross section (BBSCS) modification factor 𝑋 (with respect to their free-space values) and the stiffness parameter 𝐾 of dense nuclear matter. Here, we find that the most probable value of 𝑋 evolves from around 0.7 to 1.0 as the beam energy 𝐸 beam /𝐴 increases. On the other hand, the posterior PDF⁡(𝐾) may have dual peaks having roughly the same height or extended shoulders at high 𝐾 values. More quantitatively, the posterior PDF⁡(𝐾) changes from having a major peak around 220 MeV, characterizing a soft EOS in the reaction at 𝐸 beam /𝐴 = 150 MeV, to one that peaks around 320 MeV, indicating a stiff EOS in the reactions at 𝐸 beam /𝐴 higher than about 600 MeV. The transition from soft to stiff happens in mid-central Au + Au reactions at beam energies around 250 MeV/nucleon in which 𝐾 = 220 MeV and 𝐾 = 320 MeV are approximately equally probable. Altogether, the FOPI proton-flow excitation function data indicate a gradual hardening of hot and dense nuclear matter as its density and temperature increase in reactions with higher beam energies.

nuclear matter↗

DOE OSTI2024

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

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

Astronomy & Astrophysics↗

DOE OSTI2024

Impact of the newly revised gravitational redshift of x-ray burster GS 1826-24 on the equation of state of supradense neutron-rich matter

Thanks to the recent advancement in producing rare isotopes and measuring their masses with unprecedented precision, the updated nuclear masses around the waiting-point nucleus 64 Ge in the rapid-proton capture process have led to a significant revision of the surface gravitational redshift of the neutron star (NS) in GS 1826-24 by refitting its x-ray burst light curve using Modules for Experiments in Stellar Astrophysics (MESA). The resulting NS compactness ξ is between 0.183 and 0.259 at 95% confidence level, and its upper boundary is significantly smaller than the maximum ξ previously known. Incorporating these new data within a comprehensive Bayesian statistical framework, we investigate its impact on the Equation of State (EOS) of supradense neutron-rich matter and the required spin frequency for GW190814’s minor m 2 with mass 2.59 ± 0.05⁢M ⊙ to be a rotationally stable pulsar. We found that the EOS of high-density symmetric nuclear matter (SNM) has to be softened significantly while the symmetry energy at supersaturation densities stiffened compared to our prior knowledge from earlier analyses using data from both astrophysical observations and terrestrial nuclear experiments. In particular, the skewness J 0 characterizing the stiffness of high-density symmetric nuclear matter (SNM) decreases significantly, while the slope L, curvature K sym , and skewness Jsym of nuclear symmetry energy all increase appreciably compared to their fiducial values. Here, we also found that the most probable spin rate for the m 2 to be a stable pulsar is very close to its mass-shedding limit once the revised redshift data from GS 1826-24 is considered, making the m 2 unlikely the most massive NS observed so far.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Neutron matter from local chiral effective field theory interactions at large cutoffs

Neutron matter is an important many-body system that provides valuable constraints for the equation of state (EOS) of neutron stars. Neutron-matter calculations employing chiral effective field theory (EFT) interactions have been extensively used for this purpose. Among the various many-body methods, quantum Monte Carlo (QMC) methods stand out due to their nonperturbative nature and the achievable precision. However, QMC methods require local interactions as input, which leads to the appearance of stronger regulator artifacts compared to nonlocal interactions. To circumvent this, we employ large-cutoff interactions derived within chiral EFT (400 MeV ≤ Λ 𝑐 ≤ 700MeV) for studies of pure neutron matter. These interactions have been adjusted to nucleon-nucleon scattering phase shifts, the triton binding energy, as well as the triton 𝛽-decay half-life. We find that regulator artifacts significantly decrease with increasing cutoff, leading to a significant reduction of uncertainties in the neutron-matter EOS. We discuss implications for the symmetry energy and demonstrate how our new calculations lead to a reduction in the theoretical uncertainty of predicted neutron-star radii by up to 30% for low-mass stars.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Systematics of the chemical freeze-out line in the high baryon density regime explored at SIS100

The systematic uncertainties of chemical freeze-out fits at SIS100 energies (Au+Au reactions at $\sqrt{s_{NN}}$ = 3 – 5 GeV) are studied using UrQMD simulations. Although hadron production in UrQMD does not occur on a sharp chemical freeze-out hyper-surface, the extracted fit quality is shown to be very good. The extracted chemical parameters depend on the selected hadron species as well as the underlying equation of state (EoS) of the matter. Including light nuclei and anti-protons in the fit increases the expected freeze-out temperature, while a stiffer EoS increases the obtained chemical potential. Similarly, the baryon densities extracted by the thermal fits depend on the choice of hadrons as well as the underlying equation of state. These results are important for the upcoming CBM@FAIR physics program and highlight that a degree of caution is advised when one relates the chemical freeze-out curve to features on the QCD phase diagram like the critical endpoint or a possible phase transition.

Hofmann, Emma Lilith [GSI Helmholtzzentrum für Sch↗

DOE OSTI2025

Single-shot electro-optic sampling with arbitrary terahertz polarization

With the recent development of diversity electro-optic sampling (DEOS), significant progress has been made in the range of applicability of single-shot EOS measurements, allowing broadband THz waveforms to be captured in a single shot over large temporal windows. In addition to the decrease in acquisition time compared to standard multishot data acquisition, this technique allows measurements on systems far from equilibrium with large shot-to-shot noise or with irreversible or poorly repeatable dynamics. Although DEOS has been demonstrated and verified for linearly polarized THz waveforms, we investigate the effects resulting from the presence of a secondary polarization component. This imposes new challenges for accurate waveform reconstruction, and opens the opportunity to measure out complex polarization states such as arbitrary elliptically polarized THz field. We demonstrate a single-shot diversity-electro-optic-sampling-based approach to capture both x- and y-THz fields simultaneously with a single (110)-cut EO crystal for THz polarimetry and ellipsometry over a wide range of frequencies.

Lenz, Maximilian [University of California, Los An↗

DOE OSTI2024

IDAES-PSE 2.7.0 Release

The Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform is a versatile computational environment offering extensive process systems engineering (PSE) capabilities for optimizing the design and operation of complex, interacting technologies and systems. IDAES enables users to efficiently search vast, complex design spaces to discover the lowest cost solutions while supporting the full process modeling lifecycle, from conceptual design to dynamic optimization and control. The extensible, open platform empowers users to create models of novel processes and rapidly develop custom analyses, workflows, and end-user applications. IDAES-PSE 2.7.0 Release Highlights New features: AutoScaler and CustomScalerBase classes: Such tools are the core of the new scaling framework being implemented in IDAES. Wider adoption of scaling tools among users will result in quicker and more robust model solutions. Scaler for equilibrium reactor and saponification properties: These scaler models are examples to follow for how to use the new scaling tools. ONNX Surrogate support from Optimization & Machine Learning Toolkit (OMLT): ONNX is an open standard format to save and load ML/AI models that is widely supported by all major frameworks. This capability makes it easier for IDAES users to create surrogate models and use them without having to support each framework individually. 1D Membrane Model for CO2 Capture and Utilization: Supports ongoing efforts for modeling and optimizing polymer membrane processes for CO2 capture and conversion into formic acid. StreamScaler unit model: Unrelated to the CustomScalerBase, this unit model allows a stream’s extensive variables to be scaled by a fixed factor. This allows streams being processed by multiple units in parallel to be scaled down to unit scale and scaled back up to process scale. Bug fixes or improvements: Scaling, EoS, Diagnostics tool, Modular Properties, tests & documentation Deprecations: Old Cubic EoS

AS↗