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,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Quantification of carbon phases in detonation soot

Soot is a powdery substance composed of solid residues that results from the incomplete combustion of hydrocarbons. Soot holds valuable information because its chemical makeup can provide answers for condensed carbon kinetics in HE detonations by quantitively identifying what carbon species remain. For carbon-rich explosives, this is particularly important because carbon kinetics is a major pathway for energy release. Current efforts have been able to qualitatively identify the carbons and confirm their presence, but we haven’t been able to successfully quantify them. For example, small angle x-ray diffraction has been fielded to qualitatively identify carbon phases as a function of time during a detonation (work by Trevor Willey and Mike Nielson). Other work (ice experiments) has focused on measuring the “effective” static cold curve equation of state (EOS) of the recovered carbon soots.

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

Cluster Dynamics Simulations of Intra-Granular Fission Gas Bubble Size and Pressure Evolution in UO 2

Fission gases such as xenon (Xe) play a critical role in determining the behavior and response of nuclear fuel. Given that Xe has little solubility in UO 2 , it accumulates and forms bubbles, which significantly impact fuel performance. Intra- and inter-granular bubble nucleation and growth can lead to fuel swelling, and once bubbles interconnect at grain boundaries, fission gas can be released into the plenum. At low temperatures, limited uranium vacancy mobility can restrict swelling, therefore causing the bubbles to become highly pressurized. Consequently, this can induce micro-cracking, promote fission gas release (increasing the likelihood of cladding failure), and even lead to fuel pulverization under accident conditions such as a loss of coolant accident. As bubble evolution is strongly influenced by local temperature and fission rate, markedly different behavior occurs across the radial profile of the fuel pellet. Capturing the mechanisms that underpin bubble evolution is therefore important to predict these behaviors in the fuel. Previous models describing important mechanisms informed by lower length scale simulations have been developed under the NEAMS program. These can describe the evolution of a single bubble type (i.e., single value for radius and pressure) at each position in the pellet, for instance using the Centipede cluster dynamic code. However, in reality, a full distribution in bubble sizes and pressures exists within the microstructure at a given position in the pellet. To address this the cluster dynamics code Xolotl, which can predict Xe and vacancy phase space (i.e., bubble distributions) for intra-granular bubbles, has been used before. Prior work benchmarked the Xolotl code against the Centipede cluster dynamics code to ensure compatibility and to verify that mobile defect properties are adequately transferred between the two codes, along with some physics improvements. In this work, we go further by introducing a physics-based set of improvements that will allow us to accurately predict bubble size distributions and internal bubble pressures under representative UO 2 irradiation conditions. The improvements include (i) coupling bubble-defect reaction energies to a virial equation of state (EOS), (ii) including a bubble surface tension contribution, (iii) incorporating radiation-induced re-solution of Xe and vacancies, (iv) enabling pressure-driven dislocation loop punching through an effective emission of interstitial clusters informed by interstitial loop energetics, (v) accounting for radiation induced athermal diffusion of Xe, and (vi) implementing a Booth-type grain boundary sink representation for all mobile defects and defect clusters. After these modifications, we observe good agreement of Xolotl fission gas bubble size and concentration predictions with legacy experimental measurements. Additionally, it allows the distribution of Xe bubble pressures and radius to also be predicted and compared to data produced through the Advanced Fuels Campaign (AFC) program. Here, we have done this by running simulations under conditions similar to the AFC post-irradiation examination (PIE) samples irradiated at North Anna 2 light water reactor (LWR). Our results shows excellent agreement with these experimental measurements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Effect of Heat Treatment on Spallation of Additively manufactured AlSi10Mg

The eutectic alloy AlSi10Mg is commonly used for additive manufacturing (AM), as its high Si content improves printability. This abundance of Si forms Si-rich precipitates whose structure is heavily dependent on the thermal history and have been shown to strongly influence its mechanical properties. In this presentation we present the results of experiments to study the effect of the Si-structure on the equation of state and spallation strength of AM AlSi10Mg.

Elmslie, Timothy Alexander [Sandia National Labora↗
DOE OSTI2025

Dual-unitary shadow tomography

We introduce a classical shadow tomography scheme based on dual-unitary brick-wall circuits termed "dual-unitary shadow tomography" (DUST). For this we study operator spreading and Pauli weight dynamics in one-dimensional qubit systems, evolved by random two-local dual-unitary gates arranged in a brick-wall structure, ending with a final measurement layer. We do this by deriving general constraints on the Pauli weight transfer matrix and specializing to the case of dual-unitarity. We first show that dual-unitaries must have a minimal amount of entropy production. Remarkably, we find that operator spreading in these circuits have a rich structure resembling that of relativistic quantum field theories, with massless chiral excitations that can decay or fuse into each other, which we call left- or right-movers. We develop a mean-field description of the Pauli weight in terms of $\rho(x,t)$, which represents the probability of having nontrivial support at site $x$ and depth $t$ starting from a fixed weight distribution. We develop an equation of state for $\rho(x,t)$, and simulate it numerically using Monte Carlo simulations. Lastly, we demonstrate that the fast-thermalizing properties of dual-unitary circuits make them better at predicting large operators than shallow brick-wall Clifford circuits. Our results are robust to finite-size effects due to the chirality of dual-unitary brick-wall circuits.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Advancements in High-Q development for novel medium-velocity 644 MHz 5-cell elliptical superconducting RF cavities for continuous-wave operation in heavy-ion linacs

The Facility for Rare Isotope Beams is a world-leading center for experimental nuclear physics research, and relies on a first-in-kind superconducting RF driver linac to supply 400 kW beams of a uniquely large range of particles, from protons through the heaviest uranium ions. The FRIB400 project proposes to double the end energy of the current superconducting FRIB linac for the heaviest uranium ions from 200 MeV per nucleon (MeV/u) to 400 MeV/u, which equates to approximately 1 GeV for protons. The increased rare isotope production from higher-energy drive beams would deliver new and exciting capabilities to FRIB users, which include extending the facility's reach along the neutron drip-line, greatly increasing the facility's rare isotope yield, and significantly improving the precision of parameter measurements for the nuclear matter equation of state, which is a particularly timely complement to the recent advent of multi-messenger astronomy, and the study of neutron star mergers. The unique set of design parameters defined by the 400 MeV/u (uranium) energy goal, the physical size, continuous-wave operation, and cryogenic capacity of the extant facilities at FRIB, requires development of a first-in-kind beta = v/c = 0.65, 644 MHz 5-cell elliptical superconducting RF (SRF) cavity, with field-leading efficiency, measured by the cavity's intrinsic quality factor, Q_0. The field of SRF has recently experienced a very active period of development, which demonstrated the potential to raise the Q_0 of 1.3 GHz cavities by factors of two to three or more. In addition to uncovering fascinating new properties of the fundamental physics of superconducting RF phenomena, these improvements suggest substantial improvements could be made over the minimum Q_0 of the proposed FRIB400 upgrade SRF cavities, if these new techniques can be adapted to the beta = 0.65, 644 MHz regime. This work presents the first validation of the the novel beta = 0.65, 644 MHz 5-cell elliptical superconducting RF cavity design, and the first results of both conventional and advanced RF surface preparation techniques. In particular, this work focuses on two recently-developed high-Q_0 recipes: N-doping, and furnace baking, setting the current world-record Q_0 for this type of resonant cavity at 3.8 x 10^10 at the FRIB400 operating gradient of 17.5 MV/m. We also find evidence that these results may be further improved upon, and make the case for further RF surface preparation recipe optimization. We then pair these high-power RF investigations with techniques borrowed from material science to try to better understand the ways in which specific properties of the niobium material, such as grain structure and superconducting flux pinning force, can be related to its superconducting RF performance. These experiments better elucidate the effects of high-temperature annealing on niobium samples, which suggest methods of further optimizing cavity flux expulsion and thus SRF performance.Read

McGee, Kellen E.↗
DOE OSTI2025

Bulk viscosity of two-flavor color superconducting quark matter in neutron star mergers

This work investigates the bulk viscosity of warm, dense, neutrino-transparent, color-superconducting quark matter, where damping of density oscillations in the kHz frequency range arises from weak-interaction-driven direct Urca processes involving quarks. We study the two-flavor red-green paired color-superconducting (2SC) phase, while allowing for the presence of unpaired strange quarks and blue color light quarks of all flavors. Our calculations are based on the SU(3) Nambu–Jona–Lasinio model, extended to include both vector interactions and the ‘t Hooft determinant term. The primary focus is on how variations in the NJL Lagrangian parameters—specifically, the diquark and vector coupling strengths—affect both the static properties of quark matter, such as its equation of state and composition, and its dynamical behavior, including bulk viscosity and associated damping timescales. We find that the bulk viscosity and corresponding damping timescale can change by more than an order of magnitude upon varying the vector coupling by a factor of two at high densities and by a lesser degree at lower densities. This sensitivity primarily arises from the susceptibility of 2SC matter, with a smaller contribution from modifications to the weak interaction rates. In comparison, changes in the diquark coupling have a more limited impact. The damping of density oscillations in 2SC matter is similar quantitatively to nucleonic matter and can be a leading mechanism of dissipation in merging hybrid stars containing color superconducting cores.

Gravitational Waves↗
DOE OSTI2020

Experimental and Computational Studies on Superhard Material Rhenium Diboride under Ultrahigh Pressures

An emerging class of superhard materials for extreme environment applications are compounds formed by heavy transition metals with light elements. In this work, ultrahigh pressure experiments on transition metal rhenium diboride (ReB 2 ) were carried out in a diamond anvil cell under isothermal and non-hydrostatic compression. Two independent high-pressure experiments were carried out on ReB 2 for the first time up to a pressure of 241 GPa (volume compression V/V 0 = 0.731 ± 0.004), with platinum as an internal pressure standard in X-ray diffraction studies. The hexagonal phase of ReB 2 was stable under highest pressure, and the anisotropy between the a-axis and c-axis compression increases with pressure to 241 GPa. The measured equation of state (EOS) above the yield stress of ReB 2 is well represented by the bulk modulus K 0 = 364 GPa and its first pressure derivative K 0 ´ = 3.53. Corresponding density-functional-theory (DFT) simulations of the EOS and elastic constants agreed well with the experimental data. DFT results indicated that ReB 2 becomes more ductile with enhanced tendency towards metallic bonding under compression. The DFT results also showed strong crystal anisotropy up to the maximum pressure under study. The pressure-enhanced electron density distribution along the Re and B bond direction renders the material highly incompressible along the c-axis. Our study helps to establish the fundamental basis for anisotropic compression of ReB 2 under ultrahigh pressures

36 MATERIALS SCIENCE↗
DOE OSTI2024

Pressure-Induced Exciton Formation and Superconductivity in Platinum-Based Mineral Sperrylite

We report a comprehensive study of Sperrylite (PtAs 2 ), the main platinum source in natural minerals, as a function of applied pressures up to 150 GPa. While no structural phase transition is detected from pressure-dependent X-ray measurements, the unit cell volume shrinks monotonically with pressure following the third-order Birch–Murnaghan equation of state. The mildly semiconducting behavior found in pure synthesized crystals at ambient pressures becomes more insulating upon increasing the applied pressure before metalizing at higher pressures, giving way to the appearance of an abrupt decrease in resistance near 3 K at pressures above 92 GPa consistent with the onset of a superconducing phase. The pressure evolution of the calculated electronic band structure reveals the same physical trend as our transport measurements, with a non-monotonic evolution explained by a hole band that is pushed below the Fermi energy and an electron band that approaches it as a function of pressure, both reaching a touching point suggestive of an excitonic state. A Lifshitz transition of the electronic structure and an increase in the density of states may naturally explain the onset of superconductivity in this material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Study on the High-Pressure Behavior of Goethite up to 32 GPa Using X-Ray Diffraction, Raman, and Electrical Impedance Spectroscopy

Goethite is a major iron-bearing sedimentary mineral on Earth. In this study, we conducted in situ high-pressure x-ray diffraction, Raman, and electrical impedance spectroscopy measurements of goethite using a diamond anvil cell (DAC) at room temperature and high pressures up to 32 GPa. We observed feature changes in both the Raman spectra and electrical resistance at about 5 and 11 GPa. However, the x-ray diffraction patterns show no structural phase transition in the entire pressure range of the study. The derived pressure-volume (P-V) data show a smooth compression curve with no clear evidence of any second-order phase transition. Fitting the volumetric data to the second-order Birch–Murnaghan equation of state yields V 0 = 138.9 ± 0.5 Å 3 and K 0 = 126 ± 5 GPa.

58 GEOSCIENCES↗
DOE OSTI2020

High-Pressure Sound Velocity Measurements of Liquids Using In Situ Ultrasonic Techniques in a Multianvil Apparatus

Sound velocity and equation of state of liquids provide important constraints on the generation, presence, and transport of silicate and metallic melts in the Earth’s interior. Unlike their solid counterparts, these properties of liquids pose great technical challenges to high-pressure measurements and are poorly constrained. Here we present the technical developments that have been made at the GSECARS beamline 13-ID-D of the Advanced Photon Source for the past several years for determination of sound velocity of liquids using the ultrasonic techniques in a 1000-ton Kawai-type multianvil apparatus. Temperature of the sound velocity measurements has been extended to ~2400 K at 4 GPa and ~2000 K at 8 GPa to enable studies of liquids with very high melting temperatures, such as the silicate liquids.

58 GEOSCIENCES↗
DOE OSTI2020

Structure and Behavior of the Ni End-Member Schreibersite Ni 3 P under Compression to 50 GPa

To better understand the potential presence of light element alloys of Fe and Ni in the Earth’s interior, the crystal structure and compressional behavior of the Ni-P binary compound, schreibersite (Ni 3 P), have been investigated using synchrotron X-ray diffraction experiments. Both powder and two single-crystal samples of synthetic Ni 3 P (in different orientations with respect to the loading axis of the diamond anvil cell) were compressed up to approximately 50 GPa at ambient temperature. The compressional data obtained for Ni3P were fitted with a 3rd order Birch–Murnaghan equation of state. All data indicated that the c/a ratio of unit cell parameters remained approximately constant up to about 30 GPa but then increased progressively with pressure, exhibiting a second slight discontinuity at approximately 40 GPa. The changes in unit cell parameters at ~30 GPa and ~40 GPa suggested discontinuous changes in magnetic ordering. Moreover, the threshold of these subtle discontinuities is sensitive to the stress state and orientation of the crystal in the diamond anvil cell. This study is the first report on the compressional behavior of both powder and single-crystal schreibersite at high-pressure (up to 50 GPa). It offers insights into the effects of Ni 3 P components on the compressional behavior of the Earth’s core.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Structure and Stability of Iron Fluoride at High Pressure–Temperature and Implication for a New Reservoir of Fluorine in the Deep Earth

Fluorine (F) is the most abundant halogen in the bulk silicate Earth. F plays an important role in geochemical and biological systems, but its abundance and distribution in the terrestrial mantle are still unclear. Recent studies suggested that F reservoirs in the deep mantle are potentially hosted in terrestrial oxide minerals, especially in aluminous bridgmanite. However, the knowledge about the formation and stability field of fluoride in the Earth’s interior is rare. In this study, we combine in situ laser-heated diamond anvil cell, synchrotron X-ray diffraction, and first-principles structure search to show that a new tetragonal structure of FeF 3 is stable at pressures of 78–130 GPa and temperatures up to ~1900 K. Simulation predicted the tetragonal phase takes a much denser structure due to the rotation of FeF 6 octahedral units. The equations of states of tetragonal FeF3 are determined by experiment and verified by simulation. Our results indicate that FeF 3 can be a potential key phase for storing F in the Earth’s lower mantle and may explain some mantle-derived magma with high F concentration.

58 GEOSCIENCES↗
DOE OSTI2020

Invisible QCD as Dark Energy

We account for the late time acceleration of the Universe by extending the Quantum Chromodynamics (QCD) color to a SU(3) invisible sector (IQCD). If the Invisible Chiral symmetry is broken in the early universe, a condensate of dark pions (dpions) and dark gluons (dgluons) forms. The condensate naturally forms due to strong dynamics similar to the Nambu–Jona-Lasinio mechanism. As the Universe evolves from early times to present times the interaction energy between the dgluon and dpion condensate dominates with a negative pressure equation of state and causes late time acceleration. We conclude with a stability analysis of the coupled perturbations of the dark pions and dark gluons.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Data Release 1 of the Dark Energy Spectroscopic Instrument

In 2021 May the Dark Energy Spectroscopic Instrument (DESI) collaboration began a 5 yr spectroscopic redshift survey to produce a detailed map of the evolving three-dimensional structure of the Universe between z = 0 and z ≈ 4. DESI’s principal scientific objectives are to place precise constraints on the equation of state of dark energy, the gravitationally driven growth of large-scale structure, and the sum of the neutrino masses, and to explore the observational signatures of primordial inflation. We present DESI DR1, which consists of all data acquired during the first 13 months of the DESI main survey, as well as a uniform reprocessing of the DESI Survey Validation data, which were previously made public in the DESI Early Data Release. The DR1 main survey includes high-confidence redshifts for 18.7M objects, of which 13.1M are spectroscopically classified as galaxies, 1.6M as quasars, and 4M as stars, making DR1 the largest sample of extragalactic redshifts ever assembled. We summarize the DR1 observations, the spectroscopic data-reduction pipeline and data products, large-scale structure catalogs, value-added catalogs, and describe how to access and interact with the data. In addition to fulfilling its core cosmological objectives with unprecedented precision, we expect DR1 to enable a wide range of transformational astrophysical studies and discoveries.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Constraints on the Muon Fraction and Density Profile in Neutron Stars

Muons in neutron stars (NSs) play especially important roles in addressing several interesting new physics questions associated with detecting as well as understanding interactions and astrophysical effects of muonphilic dark matter particles. The key model inputs for studying the latter are the total muon mass M μ , the muon mass fraction M μ /M NS over the NS mass M NS , and the muon radial density profile ρ μ (r) in NSs of varying masses. We investigate these quantities within a minimum model for the core of NSs consisting of neutrons, protons, electrons, and muons using an explicitly isospin-dependent parametric equation of state (EOS) constrained by available nuclear laboratory experiments and the latest astrophysical observations of NS masses, radii, and tidal deformabilities. Here, we found that the absolutely maximum muon mass M μ and its mass fraction M μ /M NS in the most massive NSs allowed by causality are about 0.025 M ⊙ and 1.1%, respectively. For the most massive NS of mass 2.14 M ⊙ observed so far, they reduce to about 0.020 M ⊙ and 0.9%, respectively. We also study respective effects of individual parameters describing the EOS of high-density neutron-rich nucleonic matter on the muon contents in NSs with varying masses. We found that the most important but uncertain nuclear physics ingredient for determining the muon contents in NSs is the high-density nuclear symmetry energy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Bayesian Inference of the Symmetry Energy of Superdense Neutron-rich Matter from Future Radius Measurements of Massive Neutron Stars

In this work, using as references the posterior probability distribution functions of the equation of state (EOS) parameters inferred from the radii of canonical neutron stars (NSs) reported by the LIGO/VIRGO and NICER Collaborations based on their observations of GW170817 and PSR J0030+0451, we investigate how future radius measurements of more massive NSs will improve our current knowledge about the EOS of superdense neutron-rich nuclear matter, especially its symmetry energy term. Within the Bayesian statistical approach using an explicitly isospin-dependent parametric EOS for the core of NSs, we infer the EOS parameters of superdense neutron-rich nuclear matter from three sets of imagined mass–radius correlation data representing typical predictions by various nuclear many-body theories, that is, the radius stays the same, decreases, or increases with increasing NS mass within ±15% between 1.4 and 2.0 M$_\odot$. The corresponding NS average density increases quickly or slowly or slightly decreases as the NS mass increases from 1.4 to 2.0 M$_\odot$. While the EOSs of symmetric nuclear matter (SNM) inferred from the three data sets are approximately the same, the corresponding symmetry energies above about twice the saturation density of nuclear matter are very different, indicating that the radii of massive NSs carry important information about the high-density behavior of nuclear symmetry energy with little influence from the remaining uncertainties of the SNM EOS at suprasaturation densities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

On the Minimum Radius of Very Massive Neutron Stars

Prospects of establishing the radii of massive neutron stars in PSR J1614 – 2230 and PSR J0740 + 6620 from Neutron Star Interior Composition ExploreR and Chandra observatories hold the potential to constrain the equation of state (EoS) of matter to densities well beyond those encountered in canonical stars of mass $\sim 1.4\,{M}_{\odot }$. In this work, we investigate the relation between the radii of very massive neutron stars up to the maximum mass, ${M}_{\max }$, supported by dense matter EoSs. Results from models with hadronic matter are contrasted with those that include a first-order hadron-to-quark phase transition. We find that a lower bound on ${M}_{\max }$ with an upper bound on the radius of massive pulsars serves to rule out quark matter that is too soft, and an upper bound on ${M}_{\max }$ with a lower bound on the radius of massive pulsars strongly disfavors a transition into quark matter that is too stiff appearing at low densities. In conclusion, the complementary role played by radius inferences from future gravitational-wave events of inspiraling binary neutron stars is also briefly discussed.

79 ASTRONOMY AND ASTROPHYSICS↗