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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Comparing matching prescriptions between pre-equilibrium and hydrodynamic models in high-energy nuclear collisions

Abstract State-of-the-art simulations of high-energy nuclear collisions rely on hybrid setups, involving in particular a pre-equilibrium stage to let the system evolve from a far-from-equilibrium initial condition towards a near-equilibrated state after which fluid dynamics can be applied meaningfully. A known issue is the mismatch between the equation of state in the fluid-dynamical evolution and the effective one in the previous stage, which leads to discontinuities at the interface between the two models. Here we introduce a new matching prescription at this interface, based on the entropy, and we compare it with the standard one relying on local energy conservation. We study the behavior of various quantities at the switching time between the models and investigate a number of final-state hadronic observables. For the latter, we show that they are not modified significantly by the choice of matching prescription, provided an appropriate normalization is chosen for the initial state. In turn, our approach reduces sizeably the ratio of bulk over thermodynamic pressure at the beginning of the fluid-dynamical stage.

Borghini, Nicolas (ORCID:0000000239063258)↗
DOE OSTI2025

Relativistic core–valence-separated equation-of-motion coupled-cluster singles and doubles method: Efficient implementation and benchmark calculations

An efficient implementation for the relativistic exact two-component core–valence-separated equation-of-motion coupled-cluster singles and doubles (X2C-CVS-EOM-CCSD) method is reported. The explicit exclusion of pure valence excitations in the EOM-CCSD excited-state eigenvalue equations significantly improves the efficiency for calculations of core-excited states. Benchmark relativistic CVS-EOM-CC calculations with systematic inclusion of relativistic, correlation, and basis-set effects are shown to provide highly accurate results for core ionized and excited states involving heavy atoms.

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

The importance of few-nucleon forces in chiral effective field theory

We study the importance of few-nucleon forces in chiral effective field theory for describing many-nucleon systems. A combinatorial argument suggests that three-nucleon forces-which are conventionally regarded as next-to-next-to-leading order-should accompany the two-nucleon force already at leading order (LO) starting with mass number $A$ ≃10–20. We find that this promotion enables the first realistic description of the 16 O ground state based on a renormalization-group-invariant LO interaction. We also performed coupled-cluster calculations of the equation of state for symmetric nuclear matter and our results indicate that LO four-nucleon forces could play a crucial role for describing heavy-mass nuclei. In conclusion, the enhancement mechanism we found is very general and could be important also in other many-body problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

R -mode Stability of GW190814's Secondary Component as a Supermassive and Superfast Pulsar

The nature of GW190814's secondary component m 2 of mass (2.50–2.67)M ⊙ in the mass gap between the currently known maximum mass of neutron stars and the minimum mass of black holes is currently under hot debate. Among the many possibilities proposed in the literature, the m 2 was suggested as a superfast pulsar while its r-mode stability against the run-away gravitational radiation through the Chandrasekhar-Friedman-Schutz mechanism is still unknown. Previously, Fortin et al. constructed a sample of 33 unified equations of state using the same nuclear interactions from the crust to the core consistently; from that sample we use those equations that fulfill all currently known astrophysical and nuclear physics constraints to compare the minimum frequency required for m 2 to rotationally sustain a mass greater than 2.50 M ⊙ with the critical frequency above which the r-mode instability occurs. Here, we use two extreme damping models assuming that the crust is either perfectly rigid or elastic. Using the stability of 19 observed low-mass X-ray binaries as an indication that the rigid crust damping of the r-mode dominates within the models studied, we find that m 2 is r-mode-stable while rotating with a frequency higher than 870.2 Hz (0.744 times its Kepler frequency of 1169.6 Hz) as long as its temperature is lower than about 3.9 × 10 7 K, further supporting the proposal that GW190814's secondary component is a supermassive and superfast pulsar.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Classical and quantum simulations of 1+1-dimensional ${\mathbb{Z}}_{2}$ gauge theory at finite temperature and density

Simulating strongly coupled gauge theories at finite temperature and density is a longstanding challenge in nuclear and high-energy physics with fundamental implications for condensed matter physics. Here, we simulate such systems using minimally entangled typical thermal state (METTS) approaches, which combine classical random sampling with imaginary-time evolution, implementable on either classical or quantum computers, to estimate thermal averages of observables. We study 1+1-dimensional ${\mathbb{Z}}_{2}$ gauge theory coupled to spinless fermionic matter, which maps onto a local quantum spin chain. We benchmark both a classical matrix-product-state implementation of METTS and a recently proposed adaptive variational approach for near-term quantum devices, focusing on the equation of state and measures of fermion confinement. Of particular importance is the choice of basis for METTS sampling, which impacts both the sampling overhead and quantum circuit complexity. Our work sets the stage for future studies of strongly coupled gauge theories using classical and quantum hardware.

Chen, I-Chi [Iowa State Univ., Ames, IA (United St↗
DOE OSTI2023

Spark Channel Dynamics of Electrostatic Discharges

When two differently-charged objects are brought in close proximity to each other, the resulting high electric fields can cause electron avalanche breakdown of the air gap separating the objects, a process known as electrostatic discharge (ESD). If enough initial charge is stored on the objects, the electrical breakdown can proceed to ionize the air to such a degree that a highly conductive filament of plasma forms in the gap, known as a spark channel. The spark electrically bridges the air gap, resulting in a rapid pulse of current that neutralizes the charge difference. The current pulse produces significant heating of the gas in the spark, resulting in dissociation, ionization, thermal radiation, and hydrodynamic expansion. ESD presents a hazard to electrically-sensitive devices, with consequences such as economic losses (e.g. damaged electronics) or unsafe response (e.g. unintended ignition of flammable gas mixtures, initiation of detonators, etc.). For this thesis, the ESD spark is taken to occur between two conducting electrodes, with the spark channel being axisymmetric in a cylindrical coordinate system centered on the channel. An RLC-type circuit is used for the discharge model of the ESD event. The spark is treated as a time-dependent resistance that is in series with a capacitance, an inductance, and (optionally) a load resistance representing a “victim” component under threat from the ESD event. The primary motivation of this work is to use a numerical hydrodynamic model to understand the energy dissipation and transport processes in the spark. The model consists of the compressible Euler equations of mass, momentum, and energy conservation together with an Eddington/P1 approximation for thermal radiation transport. To close the hydrodynamic system, an equation of state (EOS) was fitted from tabular data for air that accounts for the dissociation and ionization of air species. The hydrodynamic equations are solved using a conservative Lagrangian finite volume method. These partial differential equations are coupled to the circuit equations by calculation of the spark resistance via numerical integration of the electrical conductivity of the channel. Computational results are compared against experimental measurements of discharge current and radial density of the spark channel.

42 ENGINEERING↗
DOE OSTI2021

Excited-state potentials for modelling dense plasmas from first principles

The modelling of dense plasmas using finite-temperature density functional theory has proven very successful in determining transport properties and the equation of state of systems where quantum many-body effects and correlations play a key role in their structure. Here we show how excited-state projector augmented- wave potentials can be used to extend these calculations to explicitly model core-hole states, allowing for the calculation of the electronic structure of a range of integer charge configurations embedded in a dense plasma environment. Our excited-state potentials show good agreement with all-electron calculations at finite-temperatures, motivating their use as an efficient approach in modelling from first principles both the structure of strongly-coupled non-equilibrium plasmas and their interaction with intense X-rays.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Unraveling trace anomaly of supradense matter via neutron star compactness scaling

The trace anomaly Δ ≡ 1/3 −𝑃/𝜖 =1/3 −𝜙 quantifies the possibly broken conformal symmetry in supradense matter under pressure 𝑃 at energy density 𝜖. Perturbative QCD (pQCD) predicts a vanishing Δ at extremely high energy or baryon densities when the conformal symmetry is realized but its behavior at intermediate densities reachable in neutron stars (NSs) is still very uncertain. The extraction of Δ from NS observations strongly depends on the employed model for nuclear equation of state (EOS). Using the IPAD-TOV method based on an intrinsic and perturbative analysis of the dimensionless (IPAD) Tolman-Oppenheimer-Volkoff (TOV) equations that are further verified numerically by using 10 5 EOSs generated randomly with a metamodel in a very broad EOS parameter space constrained by terrestrial nuclear experiments and astrophysical observations, here we first show that the compactness 𝜉 ≡ 𝐺⁡𝑀 NS /𝑅⁢𝑐 2 ≡ 𝑀 NS /𝑅 of a NS with mass 𝑀 NS and radius 𝑅 scales very accurately with $\bar{Π}$ c ≡ $Π$ c · (1 +18⁢X/25) ≡ X/(1 +3⁢X 2 +4⁢X) · (1 +18⁢X/25) where X ≡ 𝜙 c = 𝑃 c /𝜖 c is the ratio of pressure over energy density at NS centers. The scaling of NS compactness thus enables one to readily read off the central trace anomaly Δ c = 1/3 −X directly from the observational data of either the mass-radius or red-shift measurements. Finally, we then demonstrate indeed that the available NS data themselves from recent X-ray and gravitational wave observations can determine model insensitively the trace anomaly as a function of energy density in NS cores, providing a stringent test of existing NS models and a clear guidance in a new direction for further understanding the nature and EOS of supradense matter.

nuclear astrophysics↗
DOE OSTI2024

Multidimensional Rietveld refinement of high-pressure neutron diffraction data of PbNCN

High-pressure neutron powder diffraction data from PbNCN were collected on the high-pressure diffraction beamline SNAP located at the Spallation Neutron Source (SNS) of Oak Ridge National Laboratory (Tennessee, USA). The diffraction data were analyzed using the novel method of multidimensional (two dimensions for now, potentially more in the future) Rietveld refinement and, for comparison, employing the conventional Rietveld method. To achieve two-dimensional analysis, a detailed description of the SNAP instrument characteristics was created, serving as an instrument parameter file, and then yielding both cell and spatial parameters as refined under pressure for the first time for solid-state cyanamides/carbodiimides. The bulk modulus B 0 = 25.1 (15) GPa and its derivative B ′ 0 = 11.1 (8) were extracted for PbNCN following the Vinet equation of state. Surprisingly, an internal transition was observed beyond 2.0 (2) GPa, resulting from switching the bond multiplicities (and bending direction) of the NCN 2− complex anion. The results were corroborated using electronic structure calculation from first principles, highlighting both local structural and chemical bonding details.

Meinerzhagen, Yannick↗
DOE OSTI2022

Design of high-pressure iron Rayleigh–Taylor strength experiments for the National Ignition Facility

Iron is an important metal, scientifically and technologically. It is a common metal on Earth, forming the main constituent of the planet's inner core, where it is believed to be in solid state at high pressure and high temperature. It is also the main component of many important structural materials used in quasistatic and dynamic conditions. Laser-driven Rayleigh–Taylor instability provides a means of probing material strength at high pressure and high temperature. The unavoidable phase transition in iron at relatively low pressure induces microstructural changes that ultimately affect its strength in this extreme regime. This inevitable progression can make it difficult to design experiments and understand their results. Here, we address this challenge with the introduction of a new approach: a direct-drive design for Rayleigh–Taylor strength experiments capable of reaching up to 400 GPa over a broad range of temperatures. We use 1D and 2D hydrodynamic simulations to optimize target components and laser pulse shape to induce the phase transition and compress the iron to high pressure and high temperature. At the simulated pressure–temperature state of 350 GPa and 4000 K, we predict a ripple growth factor of 3–10 depending on the strength with minimal sensitivity to the equation of state model used. The growth factor is the primary observable, and the measured value will be compared to simulations to enable the extraction of the strength under these conditions. These experiments conducted at high-energy laser facilities will provide a unique way to study an important metal.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Machine learning the metastable phase diagram of covalently bonded carbon

Abstract Conventional phase diagram generation involves experimentation to provide an initial estimate of the set of thermodynamically accessible phases and their boundaries, followed by use of phenomenological models to interpolate between the available experimental data points and extrapolate to experimentally inaccessible regions. Such an approach, combined with high throughput first-principles calculations and data-mining techniques, has led to exhaustive thermodynamic databases (e.g. compatible with the CALPHAD method), albeit focused on the reduced set of phases observed at distinct thermodynamic equilibria. In contrast, materials during their synthesis, operation, or processing, may not reach their thermodynamic equilibrium state but, instead, remain trapped in a local (metastable) free energy minimum, which may exhibit desirable properties. Here, we introduce an automated workflow that integrates first-principles physics and atomistic simulations with machine learning (ML), and high-performance computing to allow rapid exploration of the metastable phases to construct “metastable” phase diagrams for materials far-from-equilibrium. Using carbon as a prototypical system, we demonstrate automated metastable phase diagram construction to map hundreds of metastable states ranging from near equilibrium to far-from-equilibrium (400 meV/atom). We incorporate the free energy calculations into a neural-network-based learning of the equations of state that allows for efficient construction of metastable phase diagrams. We use the metastable phase diagram and identify domains of relative stability and synthesizability of metastable materials. High temperature high pressure experiments using a diamond anvil cell on graphite sample coupled with high-resolution transmission electron microscopy (HRTEM) confirm our metastable phase predictions. In particular, we identify the previously ambiguous structure of n -diamond as a cubic-analog of diaphite-like lonsdaelite phase.

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

Simultaneous x-ray diffraction measurements of nine pressure calibrants to 140 GPa

Accurate pressure calibration is fundamental to quantitative high-pressure science, yet inconsistencies among widely used secondary pressure scales persist. In this work, we report simultaneous synchrotron x-ray diffraction meatogether in diamond anvil cells to 140 GPa. This simultaneous-measurement strategy minimizes transitive errors and run-to-run inconsistencies inherent to paired-measurement approaches, yielding direct experimental volume-volume (V-V) relations among the calibrants. These V-V relations link the compression state of each phase to that of every other phase measured in the same experiment. By anchoring these V-V relations to the reduced 300 K equation of state (EOS) of copper derived from ramp-compression measurements, we derive an internally consistent set of Cu-referenced pressure scales via Vinet EOS fits. Because the primary experimental result is the V-V dataset, it provides an EOS-independent constraint that can be rereferenced if primary standards are revised. The dataset further enables a direct experimental cross-check of consistency between the adopted reduced 300 K Cu EOS and the corresponding reduced 300 K ramp-derived EOSs for Pt, Au, Ta, and hcp-Fe under a fixed V-V constraint. The residuals quantify small, material-dependent mismatches between each phase EOS and the Cu reference EOS.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Identifying the nature of the QCD transition in relativistic collision of heavy nuclei with deep learning

Using deep convolutional neural network (CNN), the nature of the QCD transition can be identified from the final-state pion spectra from hybrid model simulations of heavy-ion collisions that combines a viscous hydrodynamic model with a hadronic cascade “after-burner”. Two different types of equations of state (EoS) of the medium are used in the hydrodynamic evolution. The resulting spectra in transverse momentum and azimuthal angle are used as the input data to train the neural network to distinguish different EoS. Different scenarios for the input data are studied and compared in a systematic way. A clear hierarchy is observed in the prediction accuracy when using the event-by-event, cascade-coarse-grained and event-fine-averaged spectra as input for the network, which are about 80%, 90% and 99%, respectively. A comparison with the prediction performance by deep neural network (DNN) with only the normalized pion transverse momentum spectra is also made. High-level features of pion spectra captured by a carefully-trained neural network were found to be able to distinguish the nature of the QCD transition even in a simulation scenario which is close to the experiments.

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

Learning robust parameter inference and density reconstruction in flyer plate impact experiments

Estimating physical parameters or material properties from experimental observations is a common objective in many areas of physics and material science. In many experiments, especially in shock physics, radiography is the primary means of observing the system of interest. However, radiography does not provide direct access to key state variables, such as density, which prevents the application of traditional parameter estimation approaches. Here we focus on flyer plate impact experiments on porous materials, and resolving the underlying parameterized equation of state (EoS) and crush porosity model parameters given radiographic observation(s). We use machine learning as a tool to demonstrate with high confidence that using only high impact velocity data does not provide sufficient information to accurately infer both EoS and crush model parameters, even with fully resolved density fields or a dynamic sequence of images. We thus propose an observable data set consisting of low and high impact velocity experiments/simulations that capture different regimes of compaction and shock propagation, and proceed to introduce a generative machine learning approach which produces a posterior distribution of physical parameters directly from radiographs. We demonstrate the effectiveness of the approach in estimating parameters from simulated flyer plate impact experiments, and show that the obtained estimates of EoS and crush model parameters can then be used in hydrodynamic simulations to obtain accurate and physically admissible density reconstructions. Finally, we examine the robustness of the approach to model mismatches, and find that the learned approach can provide useful parameter estimates in the presence of out-of-distribution radiographic noise and previously unseen physics, thereby promoting a potential breakthrough in estimating material properties from experimental radiographic images.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Perturbative treatment of nonlocal chiral interactions in auxiliary-field diffusion Monte Carlo calculations

Nuclear many-body systems, ranging from nuclei to neutron stars, are some of the most interesting physical phenomena in our universe, and quantum Monte Carlo (QMC) approaches are among the most accurate many-body methods currently available to study them. In recent decades, interactions derived from chiral effective field theory (EFT) have been widely adopted in the study of nuclear many-body systems. One drawback of the QMC approach is the requirement that the nuclear interactions need to be local, whereas chiral EFT interactions usually contain nonlocalities. In this work, we leverage the capability of computing second-order perturbative corrections to the ground-state energy in order to develop a self-consistent approach to including nonlocal operators in QMC calculations. In conclusion, we investigate both the deuteron and the neutron-matter equation of state in order to show the robustness of our technique and pave the way for future QMC calculations at higher orders in the EFT, where nonlocal operators cannot be avoided.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Fast and flexible neutrino decoupling. Part I. The Standard Model

Cosmological determinations of the number of relativistic neutrino species, N eff , are becoming increasingly accurate, and further improvements are expected both from CMB and BBN data. Given this context, we update the evaluation of N eff and the current entropy density via the momentum-averaged approach. This allows for a numerically fast description of neutrino decoupling, easily portable to an array of new physics scenarios. We revisit all aspects of this approach, including collision terms with full electron mass dependence, finite temperature QED corrections to the equation of state, neutrino oscillations, and the modelling of neutrino ensembles with effective chemical potentials. For integrated observables, our results differ by less than 0.04% from the solution of the momentum-dependent evolution equation. We outline how to extend the approach to BSM settings, and will highlight its power in Part II. To facilitate the practical implementation, we release a Mathematica and python code within nudec_BSM_v2, easily linkable to BBN codes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Transport physics dependence of Bohm speed in presheath–sheath transition

We report the ion exit flow speed at the sheath entrance is constrained by the Bohm criterion, which is used as a boundary condition for simulations that do not resolve the sheath region. Traditional Bohm criterion analysis invokes the equation of state and, thus, ignores transport physics in the sheath transition problem. An expression for the Bohm speed away from the asymptotic limit is derived from a set of anisotropic plasma transport equations. The thermal force, collisional temperature isotropization, and heat flux enter into the evaluation of the Bohm speed. By comparison with kinetic simulation results, this expression is shown to be accurate in the presheath–sheath transition region rather than a single point at the sheath entrance over a broad range of collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Nonequilibrium statistical thermodynamics of multicomponent interfaces

Nonequilibrium interfacial thermodynamics has important implications for crucial biological, physical, and industrial-scale transport processes. Here, we discuss a theory of local equilibrium for multiphase multicomponent interfaces that builds upon the “sharp” interface concept first introduced by Gibbs, allowing for a description of nonequilibrium interfacial processes such as those arising in evaporation, condensation, adsorption, etc. By requiring that the thermodynamics be insensitive to the precise location of the dividing surface, one can identify conditions for local equilibrium and develop methods for measuring the values of intensive variables at the interface. We then use extensive, high-precision nonequilibrium molecular dynamics (NEMD) simulations to verify the theory and establish the validity of the local equilibrium hypothesis. In particular, we demonstrate that equilibrium equations of state are also valid out of equilibrium, and can be used to determine interfacial temperature and chemical potential(s) that are consistent with nonequilibrium generalizations of the Clapeyron and Gibbs adsorption equations. We also show, for example, that, far from equilibrium, temperature or chemical potential differences need not be uniform across an interface and may instead exhibit pronounced discontinuities. However, even in these circumstances, we demonstrate that the local equilibrium hypothesis and its implications remain valid. These results provide a thermodynamic foundation and computational tools for studying or revisiting a wide variety of interfacial transport phenomena.

36 MATERIALS SCIENCE↗