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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Dynamic properties of FeCrMnNi, a high entropy alloy

The goal of this paper was to assess the dynamic properties of a new class of materials, High entropy alloys (HEAs). Specifically, gas gun experiments coupled with recovery were performed on a specific HEA, FeCrMnNi, to measure its equation of state (EOS) and its spall strength, which is a measure of the stress required to nucleate voids under shock loading conditions. While there has been a plethora of work performed to assess the mechanical properties of these HEAs under uniaxial stress conditions as a function of strain rate, investigation of its properties in dynamic extremes remains rare. The current work fills this gap in knowledge for this novel class of materials. Our results show that the measured EOS for this material at one velocity was in reasonable agreement with an estimated Hugoniot. Furthermore, while the spall strength for this material was found to be ~1.9 GPa, with some variation based on sample location, the main failure mode was brittle. This brittle failure leading to formation of cracks in the material is different than the ductile failure observed in pure iron and its alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Feasibility and performance of the staged Z-pinch: A one-dimensional study with FLASH and MACH2

Z-pinch platforms constitute a promising pathway to fusion energy research. Here, we present a one-dimensional numerical study of the staged Z-pinch (SZP) concept using the FLASH and MACH2 codes. We discuss the verification of the codes using two analytical benchmarks that include Z-pinch-relevant physics, building confidence on the codes' ability to model such experiments. Then, FLASH is used to simulate two different SZP configurations: a xenon gas-puff liner (SZP1*) and a silver solid liner (SZP2). The SZP2 results are compared against previously published MACH2 results, and a new code-to-code comparison on SZP1* is presented. Using an ideal equation of state and analytical transport coefficients, FLASH yields a fuel convergence ratio (CR) of approximately 39 and a mass-averaged fuel ion temperature slightly below 1 keV for the SZP2 scheme, significantly lower than the full-physics MACH2 prediction. For the new SZP1* configuration, full-physics FLASH simulations furnish large and inherently unstable CRs (>300) but achieve fuel ion temperatures of many kilo-electron volts. While MACH2 also predicts high temperatures, the fuel stagnates at a smaller CR. The integrated code-to-code comparison reveals how magnetic insulation, heat conduction, and radiation transport affect platform performance and the feasibility of the SZP concept.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The Radius of PSR J0740+6620 from NICER and XMM-Newton Data

PSR J0740+6620 has a gravitational mass of 2.08 ± 0.07 M ⊙ , which is the highest reliably determined mass of any neutron star. As a result, a measurement of its radius will provide unique insight into the properties of neutron star core matter at high densities. Here we report a radius measurement based on fits of rotating hot spot patterns to Neutron Star Interior Composition Explorer (NICER) and X-ray Multi-Mirror (XMM-Newton) X-ray observations. We find that the equatorial circumferential radius of PSR J0740+6620 is 13.7$^{+2.6}_{–1.5}$ km (68%). We apply our measurement, combined with the previous NICER mass and radius measurement of PSR J0030+0451, the masses of two other ~2 M ⊙ pulsars, and the tidal deformability constraints from two gravitational wave events, to three different frameworks for equation-of-state modeling, and find consistent results at ~1.5-5 times nuclear saturation density. For a given framework, when all measurements are included, the radius of a 1.4 M ⊙ neutron star is known to ±4% (68% credibility) and the radius of a 2.08 M ⊙ neutron star is known to ±5%. The full radius range that spans the ±1σ credible intervals of all the radius estimates in the three frameworks is 12.45 ± 0.65 km for a 1.4 M ⊙ neutron star and 12.35 ± 0.75 km for a 2.08 M ⊙ neutron star.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

High-Strain Rate Spall Strength Measurement for CoCrFeMnNi High-Entropy Alloy

In this study, we experimentally investigate the high stain rate and spall behavior of Cantor high-entropy alloy (HEA), CoCrFeMnNi. First, the Hugoniot equations of state (EOS) for the samples are determined using laser-driven CoCrFeMnNi flyers launched into known Lithium Fluoride (LiF) windows. Photon Doppler Velocimetry (PDV) recordings of the velocity profiles find the EOS coefficients using an impedance mismatch technique. Following this set of measurements, laser-driven aluminum flyer plates are accelerated to velocities of 0.5–1.0 km/s using a high-energy pulse laser. Upon impact with CoCrFeMnNi samples, the shock response is found through PDV measurements of the free surface velocities. From this second set of measurements, the spall strength of the alloy is found for pressures up to 5 GPa and strain rates in excess of 106 s−1. Further analysis of the failure mechanisms behind the spallation is conducted using fractography revealing the occurrence of ductile fracture at voids presumed to be caused by chromium oxide deposits created during the manufacturing process.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Experiments and simulations of isochorically heated warm dense carbon foam at the Texas Petawatt Laser

An experimental and simulation study of warm dense carbon foams at ambient density (n e ~10 21 cm –3 ) is presented. This study of isochorically heated foams is motivated by their potential application in carbon-atmosphere white-dwarf envelopes, where there are modeling uncertainties due to the equation of state. The foams are heated on an approximately picosecond time scale with a laser-accelerated proton beam. The cooling and expansion of the heated foams can be modeled with appropriately initialized radiation-hydrodynamics codes; xRAGE code is used in this work. The primary experimental diagnostic is the streaked optical pyrometer, which images a narrow band of radiation from the rear surface of the heated material. Presented are xRAGE modeling results for both solid aluminum targets and carbonized resorcinol-formaldehyde foam targets, showing that the foam appears to cool slowly on the pyrometer because of partial transparency. So that simulations of cooling foam are processed properly, it is necessary to account for finite optical depth in the photosphere calculation, and the methods for performing that calculation are presented in depth.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Accuracy, transferability, and computational efficiency of interatomic potentials for simulations of carbon under extreme conditions

Large-scale atomistic molecular dynamics (MD) simulations provide an exceptional opportunity to advance the fundamental understanding of carbon under extreme conditions of high pressures and temperatures. However, the fidelity of these simulations depends heavily on the accuracy of classical interatomic potentials governing the dynamics of many-atom systems. Here, this study critically assesses several popular empirical potentials for carbon, as well as machine learning interatomic potentials (MLIPs), in their ability to simulate a range of physical properties at high pressures and temperatures, including the diamond equation of state, its melting line, shock Hugoniot, uniaxial compressions, and the structure of liquid carbon. Empirical potentials fail to accurately predict the behavior of carbon under high pressure–temperature conditions. In contrast, MLIPs demonstrate quantum accuracy, with Spectral Neighbor Analysis Potential (SNAP) and atomic cluster expansion (ACE) being the most accurate in reproducing the density functional theory results. ACE displays remarkable transferability despite not being specifically trained for extreme conditions. Furthermore, ACE and SNAP exhibit superior computational performance on graphics processing unit-based systems in billion atom MD simulations, with SNAP emerging as the fastest. In addition to offering practical guidance in selecting an interatomic potential with a fine balance of accuracy, transferability, and computational efficiency, this work also highlights transformative opportunities for groundbreaking scientific discoveries facilitated by quantum-accurate MD simulations with MLIPs on emerging exascale supercomputers.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Thermoelastic properties of tungsten at simultaneous high pressure and temperature

In this work, the compressional ( P ) and shear wave velocities ( S ) and unit cell volumes (densities) of polycrystalline tungsten (W) have been measured simultaneously up to 10.5 GPa and 1073 K using ultrasonic interferometry in conjunction with x-ray diffraction and x-radiography techniques. Thermoelastic properties of W were derived using different methods. We obtained the isothermal bulk modulus K T 0 = 310.3(1.5) GPa, its pressure derivative K ' T 0 = 4.4(3), its temperature derivative at constant pressure ( ∂ K T / ∂ T ) P = – 0.0138 ( 1 ) GPa K – 1 and at constant volume ( ∂ K T / ∂ T ) V = – 0.0050 GPa K – 1 , the thermal expansion α ( 0, T ) = 1.02(27) × 10 –5 + 7.39(3.2) × 10 –9 T (K –1 ), as well as the pressure derivative of thermal expansion ( ∂ α / ∂ P ) T = – 1.44 ( 1 ) × 10 – 7 K – 1 GP a – 1 based on the high-temperature Birch–Murnaghan equation of state (EOS), the Vinet EOS, and thermal pressure approach. Finite strain analysis allowed us to derive the elastic properties and their pressure/temperature derivatives independent of the choice of pressure scale. A least-squares fitting yielded K S 0 = 314.5(2.5) GPa, K S 0 ' = 4.45(9), (∂ K S /∂ T ) P = – 0.0076(6) GPa K –1 , G 0 = 162.4(9) GPa, G 0 ' = 1.8(1), (∂ G /∂ T ) P = – 0.0175(9) GPa K –1 , and α 298 K = 1.23 × 10 – 5 K – 1 . Fitting current data to the Mie–Grüneisen–Debye EOS with derived θ 0 = 383.4 K yielded γ 0 = 1.81 ( 6 ) and q = 0.3 . The thermoelastic parameters obtained from various approaches are consistent with one another and comparable with previous results within uncertainties. Our current study provides a complete and self-consistent dataset for the thermoelastic properties of tungsten at high P – T conditions, which is important to improve the theoretical modeling of these materials under dynamic conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The phase diagram of Ti-6Al-4V at high-pressures and high-temperatures

We report results from a series of diamond-anvil-cell synchrotron x-ray diffraction and large-volume-press experiments, and calculations, to investigate the phase diagram of commercial polycrystalline high-strength Ti-6Al-4V alloy in pressure–temperature space. Up to ~30 GPa and 886 K, Ti-6Al-4V is found to be stable in the hexagonal-close-packed, or α phase. The effect of temperature on the volume expansion and compressibility of α–Ti-6Al-4V is modest. The martensitic α → ω (hexagonal) transition occurs at ~30 GPa, with both phases coexisting until at ~38–40 GPa the transition to the ω phase is completed. Between 300 K and 844 K the α → ω transition appears to be independent of temperature. ω–Ti-6Al-4V is stable to ~91 GPa and 844 K, the highest combined pressure and temperature reached in these experiments. Pressure–volume–temperature equations-of-state for the α and ω phases of Ti-6Al-4V are generated and found to be similar to pure Ti. A pronounced hysteresis is observed in the ω–Ti-6Al-4V on decompression, with the hexagonal structure reverting back to the α phase at pressures below ~9 GPa at room temperature, and at a higher pressure at elevated temperatures. Based on our data, we estimate the Ti-6Al-4V α–β–ω triple point to occur at ~900 K and 30 GPa, in good agreement with our calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Symmetric Nuclear Matter from the Strong Interaction

We study the equation of state of symmetric nuclear matter at zero temperature over a wide range of densities using two complementary theoretical approaches. At low densities, up to twice nuclear saturation density, we compute the energy per particle based on modern nucleon-nucleon and three-nucleon interactions derived within chiral effective field theory. For higher densities, we derive for the first time constraints in a Fierz-complete setting directly based on quantum chromodynamics using functional renormalization group techniques. Here, we find remarkable consistency of the results obtained from both approaches as they come together in density and the natural emergence of a maximum in the speed of sound 𝑐𝑆 at supranuclear densities. The presence of this maximum appears tightly connected to the formation of a diquark gap. Notably, this maximum is observed to exceed the asymptotic value 𝑐$^{2}_{𝑆}$ =1/3 while its exact position in terms of the density cannot yet be determined conclusively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Bayesian Exploration of Phenomenological EoS of Neutron/Hybrid Stars with Recent Observations

The description of the stellar interior of compact stars remains as a big challenge for the nuclear astrophysics community. The consolidated knowledge is restricted to density regions around the saturation of hadronic matter ρ 0 =2.8 × 10 14 g cm -3 , regimes where our nuclear models are successfully applied. As one moves towards higher densities and extreme conditions up to the quark/gluons deconfinement, little can be said about the microphysics of the equation of state (EoS). Here, we employ a Markov Chain Monte Carlo (MCMC) strategy to access the variability at high density regions of polytropic piecewise models for neutron star (NS) EoS or possible hybrid stars, i.e., a NS with a small quark-matter core. With a fixed description of the hadronic matter for low density, below the nuclear saturation density, we explore a variety of models for the high density regimes leading to stellar masses near to 2.5 M ⊙ , in accordance with the observations of massive pulsars. The models are constrained, including the observation of the merger of neutrons stars from VIRGO-LIGO and with the pulsar observed by NICER. In addition, we also discuss the possibility of the use of a Bayesian power regression model with heteroscedastic error. The set of EoS from the Laser Interferometer Gravitational-Wave Observatory (LIGO) was used as input and treated as the data set for the testing case.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Scale translation yields insights into gas adsorption under nanoconfinement

This work describes a scale-translating simulation framework to investigate gas adsorption behavior in nanoconfined pores. The framework combines molecular simulations (MSs), equation of state (EoS), and lattice Boltzmann (LB) simulations. MSs reveal the physics of methane adsorption in nano-sized pores, where input values of fugacity coefficients are optimized based on EoS predictions. Then, an LB free-energy model, which incorporates a viral EoS, upscales intermolecular forces and estimates adsorption behavior via a proposed fluid–wall interaction model. Armed with the values of the LB interaction parameter as a function of pressure, the LB model is used to predict fluid behavior in irregular nanopores, and the results are validated against reference MS data. The LB model is then used to study adsorption behavior at a continuum scale in representative organic shale nanopores based on finely characterized Vaca Muerta shale samples. Furthermore, the results show that methane adsorption could significantly increase contained fluids by 10%–25% in pores smaller than 20 nm. However, in larger pores (40 nm to 90 nm), adsorption's impact diminishes to 2%–3%, suggesting sorption's negligible role beyond a 40 nm pore size.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Correlation between time and angular alignment in molecular dynamics simulations of heavy ion collisions

Neutron-proton equilibration is a process which has been used to study the density dependence of the symmetry energy term in the nuclear equation-of-state. Here, this study utilizes constrained molecular dynamics (CoMD) simulations of 70 Zn + 70 Zn with collision energies of 35 and 45 MeV/nucleon. An algorithm is used which searches through CoMD events and identifies the PLF* after it separates from the target and determines its lifetime, Δt. It also determines the fragments that the PLF* breaks apart into and determines their angular alignment. This technique gives an opportunity to explore how the average alignment of dynamically produced fragments, $\langle$α $\rangle$ dyn , evolves with PLF* lifetime. An approximately linear relationship was determined with d$\langle$α $\rangle$ dyn /dΔt = 0.98 ± 0.08 rad / zs and 1.06 ± 0.09 rad / zs for the 35 and 45 MeV / nucleon, respectively, indicating a correlation with magnitude consistent with classically determined values which were used for prior experimental studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Laser-driven ramp-compression experiments on the national ignition facility

This report details the analyses and related uncertainties in measuring longitudinal-stress–density paths in indirect laser-driven ramp equation-of-state (EOS) experiments. Experiments were conducted at the National Ignition Facility (NIF) located at the Lawrence Livermore National Laboratory. The NIF can deliver up to 2 MJ of laser energy over 30 ns and provide the necessary laser power and control to ramp compress materials to TPa pressures (1 TPa = 10 × 10 6 atmospheres). These data provide low-temperature solid-state EOS data relevant to the extreme conditions found in the deep interiors of giant planets. In these experiments, multi-stepped samples with thicknesses in the range of 40–120 µm experience an initial shock compression followed by a time-dependent ramp compression to peak pressure. Interface velocity measurements from each thickness combine to place a constraint on the Lagrangian sound speed as a function of particle velocity, which in turn allows for the determination of a continuous stress–density path to high levels of compressibility. In this report, we present a detailed description of the experimental techniques and measurement uncertainties and describe how these uncertainties combine to place a final uncertainty in both stress and density. Here, we address the effects of time-dependent deformation and the sensitivity of ramp EOS techniques to the onset of phase transformations.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Classification of compact objects and model comparison using EOS knowledge

Nuclear theory and experiments, alongside astrophysical observations, constrain the equation of state (EOS) of supranuclear-dense matter. Conversely, knowledge of the EOS allows an improved interpretation of nuclear or astrophysical data. In this article, we use several established constraints on the EOS and the new NICER measurement of PSR J0437-4715 to comment on the nature of the primary companion in GW230529 and the companion of PSR J0514-4002E. We find that, with a probability of ≳84% and ≳68%, respectively, both objects are black holes. These likelihoods increase to above 95% when one uses GW170817’s remnant as an upper limit on the TOV mass. We also demonstrate that the current knowledge of the EOS substantially disfavors high masses and radii for PSR J⁢0030+0451, inferred recently when combining NICER with XMM-Newton background data and using particular hot-spot models. Lastly, we also use our obtained EOS knowledge to comment on measurements of the nuclear symmetry energy, finding that the large value predicted by the PREX-II measurement displays some mild tension with other constraints on the EOS.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Deformation of binary and boron-doped Ni 3 Al alloys at high pressures studied with synchrotron x-ray diffraction

In situ x-ray synchrotron diffraction experiments were carried out on nickel-based high-strength superalloys under pressure to understand their deformation mechanism using a diamond anvil cell (DAC). Radial x-ray diffraction determines the room-temperature equations of state and yield strengths of binary Ni 3 Al alloy and 500 ppm boron-doped Ni 3 Al to pressures of 20 and 46 GPa, respectively. Crystallographic preferred orientations observed in these superalloys due to anisotropic stress field in DAC indicate the onset of plastic deformation. Inverse pole figure analysis reveals that the underlying deformation mechanisms change from an octahedral slip to a simultaneous activation of octahedral and cube slips upon doping with boron. The yield-strength values were found to increase with pressure and are comparable to those determined from axial diffraction experiments. The results indicate that the yield strength of Ni 3 Al:B is about 0.5 GPa higher (at pressures below 20 GPa) due to grain boundary strengthening by boron. It is shown that due to high elastic anisotropy of Ni 3 Al alloy, the yield-strength estimations from diffraction experiments strongly depend on the micromechanical model used to convert the measured elastic strains into stresses.

42 ENGINEERING↗
DOE OSTI2023

Uncertainty quantification for neutrino opacities in core-collapse supernovae and neutron star mergers

We perform an extensive study of the correlations between the neutrino-nucleon inverse mean free paths (IMFPs) and the underlying equation of states (EoSs). Strong interaction uncertainties in the neutrino mean free path are investigated in different density regimes. The nucleon effective mass, the nucleon chemical potentials, and the residual interactions in the medium play an important role in determining neutrino-nucleon interactions in a density-dependent manner. Here we study how the above quantities are constrained by an EoS consistent with (i) nuclear mass measurements, (ii) proton-proton scattering phase shifts, and (iii) neutron star observations. We then study the uncertainties of both the charged current and the neutral current neutrino-nucleon inverse mean free paths due to the variation of these quantities, using the Hartree-Fock+random phase approximation method. Finally, we calculate the Pearson correlation coefficients between (i) the EoS-based quantities and the EoS-based quantities; (ii) the EoS-based quantities and the IMFPs; (iii) the IMFPs and the IMFPs. We find a strong impact of residual interactions on neutrino opacity in the spin and spin-isospin channels, which are not well constrained by current nuclear modelings.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Raman scattering of rhenium for secondary pressure calibration

With the increasing number of 100 s GPa experiments in the diamond anvil cell (DAC), improved accuracy in secondary pressure calibrations to extreme pressures is essential. The rhenium equation of state has been proposed as a pressure calibrant via x-ray diffraction with potentially broad applications as it is commonly used as a gasket material in DAC experiments. In this work, we conducted Raman spectroscopy experiments on rhenium in the DAC and report the pressure shift of the E2g mode, a refined high-pressure C44 and mode-Grüneisen parameter above 200 GPa. We used flat, beveled, and toroidal diamond anvils under quasi-hydrostatic and non-hydrostatic conditions. By measuring the E2g mode from the culet edge to the center, we analyzed pressure distribution based on culet type and distance from the anvil center. The shift in the E2g mode can be expressed as a function of pressure, and diamond edge measurements appear reliable across all anvil types. Comparing the center and edge pressures reveals anvil cupping, offering insights into predicting or preventing anvil failure during materials properties measurements at extreme conditions.

Diamond anvil cells↗
DOE OSTI2021

The S π RIT time projection chamber

The SAMURAI Pion Reconstruction and Ion-Tracker Time Projection Chamber (SπRIT TPC) was designed to enable measurements of heavy ion collisions with the SAMURAI spectrometer at the RIKEN Radioactive Isotope Beam Factory and provide constraints on the Equation of State of neutron-rich nuclear matter. The SπRIT TPC has a 50.5 cm drift length and an 86.4 cm × 134.4 cm 2 pad plane with 12,096 pads that are equipped with the Generic Electronics for TPCs. The SπRIT TPC allows excellent reconstruction of particles and provides isotopic resolution for pions and other light charged particles across a wide range of energy losses and momenta. Here, details of the SπRIT TPC are presented, along with discussion of the TPC performance based on cosmic ray and charged particles emitted in heavy ion collisions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗