Wide-ranged multiphase equation of state for iron and model variations addressing uncertainties in high-pressure melting
Not Available
36 MATERIALS SCIENCE↗
SEARCH · Engineering Papers
Results for “shock wave”
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.
Structural study of hcp and liquid iron under shock compression up to 275 GPa
Here we combine nanosecond laser shock compression with in situ picosecond x-ray diffraction to provide structural data on iron up to 275 GPa. We constrain the extent of hcp-liquid coexistence, the onset of total melt, and the structure within the liquid phase. Our results indicate that iron, under shock compression, melts completely by 258(8) GPa. A coordination number analysis indicates that iron is a simple liquid at these pressure-temperature conditions. We also perform texture analysis between the ambient body-centered-cubic (bcc) α, and the hexagonal-closed-packed (hcp) high-pressure ε-phase. We rule out the Rong-Dunlop orientation relationship (OR) between the α and ε - phase s. However, we cannot distinguish between three other closely related ORs: Burger's, Mao-Bassett-Takahashi, and Potter's OR. The solid-liquid coexistence region is constrained from a melt onset pressure of 225(3) GPa from previously published sound speed measurements and full melt [246.5(1.8)–258(8) GPa] from x-ray diffraction measurements, with an associated maximum latent heat of melting of 623 J/g. This value is lower than recently reported theoretical estimates and suggests that the contribution to the earth's geodynamo energy budget from heat release due to freezing of the inner core is smaller than previously thought. Melt pressures for these nanosecond shock experiments are consistent with gas gun shock experiments that last for microseconds, indicating that the melt transition occurs rapidly.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Evidence for dissociation in shock-compressed methane
Theory and experiments show that with increasing pressure, the chemical bonds of methane rearrange, leading to the formation of complex polymers and then to dissociation. However, there is disagreement on the exact conditions where these changes take place. In this study, methane samples were precompressed in diamond-anvil cells and then shock compressed to pressures reaching 400 GPa, the highest pressures yet explored in methane. Furthermore, the results reveal a qualitative change in the Hugoniot curve at 80-150 GPa, which is interpreted as a signature of dissociation based on thermodynamic calculations and theoretical predictions.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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↗
Melting anisotropy in crystalline solids
Despite the long and extensive history of melting studies, experiments to determine the dependence of melting on crystal orientation are lacking. Using longitudinal sound speed measurements in aluminum single crystals shock compressed along < 100 > and < 110 > to 168 GPa, we address this need and show that the melting transition (both onset and completion) is strongly anisotropic. Furthermore, the present findings demonstrate the need to consider the role of deformation induced lattice defects on the melting transition and provide a general approach to examine melting anisotropy in crystalline solids.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Diffuse scattering from dynamically compressed single-crystal zirconium following the pressure-induced α → ω phase transition
The prototypical α → ω phase transition in zirconium is an ideal test bed for our understanding of polymorphism under extreme loading conditions. After half a century of study, a consensus had emerged that the transition is realized via one of two distinct displacive mechanisms, depending on the nature of the compression path. However, recent dynamic-compression experiments equipped with diffraction diagnostics performed in the past few years have revealed new transition mechanisms, demonstrating that our understanding of the underlying atomistic dynamics and transition kinetics is in fact far from complete. We present classical molecular dynamics simulations of the α → ω phase transition in single-crystal zirconium shock compressed along the [0001] axis using a machine-learning-class potential. The transition is predicted to proceed primarily via a modified version of the two-stage Usikov-Zilberstein mechanism, whereby the high-pressure ω phase heterogeneously nucleates at boundaries between grains of an intermediate β phase. We further observe the fomentation of atomistic disorder at the junctions between β grains, leading to the formation of highly defective interstitial material between the ω grains. We directly compare synthetic x-ray diffraction patterns generated from our simulations with those obtained using femtosecond diffraction in recent dynamic-compression experiments, and show that the simulations produce the same unique, anisotropic diffuse scattering signal unlike any previously seen from an elemental metal. Our simulations suggest that the diffuse signal arises from a combination of thermal diffuse scattering, nanoparticlelike scattering from residual kinetically stabilized α and β grains, and scattering from interstitial defective structures. Published by the American Physical Society 2024
36 MATERIALS SCIENCE↗
Refractive index of lithium fluoride at high dynamic stresses
Alkali halides are prototypical ionic solids whose refractive index is a fundamental property related to their lattice structure and ionic polarizability. In particular, lithium fluoride (LiF) has the largest band gap of any known transparent material and maintains transparency at high pressures, making it well suited for use as an optical window in dynamic compression experiments. While empirical fits to the measured density dependence of the refractive index have been provided, a model that is valid over the entire pressure range of experiments and based on a polarization-based theoretical description is lacking. We present a refractive index model for dynamically compressed LiF based on the Lorentz-Lorenz equation, where the molecular polarizability is determined using a single oscillator model with a strain polarizability parameter Λ=0.73. Here, we show that our modeling approach provides an excellent match to the LiF refractive index data for both shock and ramp compression experiments to 900 GPa (density compression greater than threefold). Additionally, updated fits are provided to determine the refractive index correction for LiF windows used in dynamic compression experiments.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Melting and transport properties of Al 2 O 3 at extreme conditions
The high-pressure, temperature phase diagram and transport properties of materials are of broad interest to planetary sciences and high-energy-density sciences and applications. Alumina (Al 2 O 3 ) or its various forms (e.g., solutions with other oxides or silicates) are important constituents in Earth’s and super-Earths’ mantles, common window materials in dynamic compression experiments, and standard pressure calibrators in diamond-anvil-cell experiments. Its structures and transport properties are of particular importance but have not been well studied at above 100 GPa pressures. Based on extensive first-principles molecular dynamics calculations, we obtain atomic level insights on structural differences and accurate results on electrical conductivity (σ dc ), thermal conductivity (κ), and reflectivity of Al 2 O 3 in several phases at pressures of ~130–1300 GPa and temperatures of 4000–20,000 K. We find the solid-to-liquid state changes of Al 2 O 3 are accompanied by an insulator-to-semimetal transition with enhanced σ dc and κ, which is similar to other silicates and oxides and can facilitate magnetic field generation in deep interiors of rocky planets. We also find the solid-to-solid transition (from Rh 2 O 3 (II) to CaIrO 3 ) is accompanied by negligible changes in σ dc and κ at 160 GPa, contrary to previous expectations, whereas another higher pressure phase transition (CaIrO 3 to U 2 S 3 at 450 GPa) is associated with increases in σ dc and κ. Furthermore, we show the transport properties can be enhanced by oxygen vacancy in Al 2 O 3 , but are not significantly affected by forming solutions with MgSiO 3 , particularly when in the solid states.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Shock-driven amorphization and melting in Fe 2 O 3
We present measurements on Fe 2 O 3 amorphization and melt under laser-driven shock compression up to 209(10) GPa via time-resolved in situ x-ray diffraction. At 122(3) GPa, a diffuse signal is observed indicating the presence of a noncrystalline phase. Structure factors have been extracted up to 182(6) GPa showing the presence of two well-defined peaks. A rapid change in the intensity ratio of the two peaks is identified between 145(12) and 151(12) GPa, indicative of a phase change. The noncrystalline diffuse scattering is consistent with shock amorphization of Fe 2 O 3 between 122(3) and 145(12) GPa, followed by an amorphous-to-liquid transition above 151(12) GPa. Upon release, a noncrystalline phase is observed alongside crystalline α–Fe 2 O 3 . The extracted structure factor and pair distribution function of this release phase resemble those reported for Fe 2 O 3 melt at ambient pressure.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Platinum liquid-vapor phase boundary mapped by fluid flyer experiments
We report a direct measurement of the temperature and density of a metal along its liquid-vapor coexistence (L-V) curve. By shocking platinum to a high-pressure liquid, we imparted sufficient heat for subsequent isentropic release to place it in a state on the boundary between the liquid and vapor phases. Released material in the liquid phase acted as a high velocity flyer pinned to the L-V curve. We measured velocity and radiant emission of the flyer as well as the interface motion and transiting shock states induced in a downstream window material by its impact. We used these measurements to calculate temperature and density of the L-V curve state which we compare to density functional theory predictions.
Equations of state↗
Benchmarking boron carbide equation of state using computation and experiment
Boron carbide ( B 4 C ) is of both fundamental scientific and practical interest due to its structural complexity and how it changes upon compression, as well as its many industrial uses and potential for use in inertial confinement fusion (ICF) and high-energy density physics experiments. Here, we report the results of a comprehensive computational study of the equation of state (EOS) of B 4 C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up to 61 megabar from planar shock experiments performed at the National Ignition Facility (NIF). Our computational methods include path integral Monte Carlo, activity expansion, as well as all-electron Green's function Korringa-Kohn-Rostoker and molecular dynamics that are both based on density functional theory. We calculate the pressure-internal energy EOS of B 4 C over a broad range of temperatures ( ~ 6 × 10 3 – 5 × 10 8 K) and densities (0.025–50 g / cm 3 ). We assess that the largest discrepancies between theoretical predictions are ≲ 5 % near the compression maximum at 1– 2 × 10 6 K. This is the warm-dense state in which the K shell significantly ionizes and has posed grand challenges to theory and experiment. By comparing with different EOS models, we find a Purgatorio model (LEOS 2122) that agrees with our calculations. The maximum discrepancies in pressure between our first-principles predictions and LEOS 2122 are ~ 18 % and occur at temperatures between 6 × 10 3 – 2 × 10 5 K, which we believe originate from differences in the ion thermal term and the cold curve that are modeled in LEOS 2122 in comparison with our first-principles calculations. To account for potential differences in the ion thermal term, we have developed three new equation-of-state models that are consistent with theoretical calculations and experiment. We apply these new models to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Constraining physical models at gigabar pressures
High-energy-density (HED) experiments in convergent geometry are able to test physical models at pressures beyond hundreds of millions of atmospheres. The measurements from these experiments are generally highly integrated and require unique analysis techniques to procure quantitative infor mation. This work describes a methodology to constrain the physics in convergent HED experiments by adapting the methods common to many other fields of physics. As an example, a mechanical model of an imploding shell is constrained by data from a thin-shelled direct-drive exploding-pusher experiment on the OMEGA Laser System using Bayesian inference, resulting in the reconstruction of the shell dynamics and energy transfer during the implosion. The model is tested by analyzing synthetic data from a 1-D hydrodynamics code and is sampled using a Markov chain Monte Carlo to generate the posterior distributions of the model parameters. As a result, the goal of this work is to demonstrate a general methodology that can be used to draw conclusions from a wide variety of HED experiments.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
First-principles equation of state database for warm dense matter computation
We put together a first-principles equation of state (FPEOS) database for matter at extreme conditions by combining results from path integral Monte Carlo and density functional molecular dynamics simulations of the elements H, He, B, C, N, O, Ne, Na, Mg, Al, and Si as well as the compounds LiF , B 4 C , BN , CH 4 , CH 2 , C 2 H 3 , CH , C 2 H , MgO , and MgSiO 3 . For all these materials, we provide the pressure and internal energy over a density-temperature range from ~ 0.5 to 50 g cm - 3 and from ~ 10 4 to 10 9 K, which are based on ~ 5000 different first-principles simulations. We compute isobars, adiabats, and shock Hugoniot curves in the regime of L - and K -shell ionization. Invoking the linear mixing approximation, we study the properties of mixtures at high density and temperature. Furthermore, we derive the Hugoniot curves for water and alumina as well as for carbon-oxygen, helium-neon, and CH-silicon mixtures. We predict the maximal shock compression ratios of H 2 O , H 2 O 2 , Al 2 O 3 , CO , and CO 2 to be 4.61, 4.64, 4.64, 4.89, and 4.83, respectively. Finally we use the FPEOS database to determine the points of maximum shock compression for all available binary mixtures. We identify mixtures that reach higher shock compression ratios than their end members. We discuss trends common to all mixtures in pressure-temperature and particle-shock velocity spaces. In the Supplemental Material, we provide all FPEOS tables as well as computer codes for interpolation, Hugoniot calculations, and plots of various thermodynamic functions.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Improved ion heating in fast ignition by pulse shaping
The fast ignition paradigm for inertial fusion offers increased gain and tolerance of asymmetry by compressing fuel at low entropy and then quickly igniting a small region. Because this hotspot rapidly disassembles, the ions must be heated to ignition temperature as quickly as possible, but most ignitor designs directly heat electrons. A constant-power ignitor pulse, which is generally assumed, is suboptimal for coupling energy from electrons to ions. Using a simple model of a hotspot in isochoric plasma, a novel pulse shape to maximize ion heating is presented in analytical form. Bounds are derived on the maximum ion temperature attainable by electron heating only. Moreover, arranging for faster ion heating allows a smaller hotspot, improving fusion gain. As a result, under representative conditions, the optimized pulse can reduce ignition energy by over 20%.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Microphysics of shock-grain interaction for inertial confinement fusion ablators in a fluid approach
Ablator materials used for inertial confinement fusion, such as high-density carbon (HDC) and beryllium, have grain structure which may lead to small-scale density nonuniformity and the generation of perturbations when the materials are shocked and compressed. Here, we use a combination of a linear theory of shock interaction with density nonuniformity [Velikovich et al., Phys. Plasmas 14, 072706 (2007)] and numerical simulations to study shock interaction with a model representation of HDC grains. While the shock-grain interaction is nonlinear, the linear theory shows some key features of the shock-grain interaction, which also hold for the (nonlinear) simulations. The postshock perturbations are made up of sonic reflections off of grain boundaries and vorticity deposition along them, with the latter dominating the perturbed energy content. The mean (per mass) postshock perturbed kinetic energy decreases with increasing grain size, but energy will be deposited at increasing spatial scale. From the perspective of the postshock perturbed energy, the detailed linear theory largely supports a proposed method [S. Davidovits et al., Phys. Plasmas 29, 112708 (2022)] for deresolving the grains (in a similar grains model) that treats the grains statistically. Finally, our simulation results highlight the influence of thermal conduction on the perturbation dynamics at grain scales.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Sensitivity of the dynamic-shell target to laser drive nonuniformities
The dynamic-shell concept for inertial confinement fusion (ICF) uses an initially homogeneous target and a carefully shaped laser pulse to form a shell and implode it. The laser pulse consists of a series of pickets that drive shocks into the target. The first few shocks converge inwards and rebound from the center of the target, creating an expanding, low-density plasma. Subsequent shocks are launched into the expanding plasma and eventually coalesce to form a shell, which is then imploded with a traditional ICF laser pulse. This study describes radiation-hydrodynamic simulations that investigate the sensitivity of dynamic-shell targets to imperfections in the laser drive. A one-dimensional (1D) study looks at mistiming and power variations in the pickets and a two-dimensional (2D) study examines irradiation perturbations imposed by the laser-beam geometry. Simulations show that less than ~2% power imbalance or 200 ps timing variation in the pickets is sufficient to keep the yield above 90% of the maximum. Additionally, the 2D simulations show that 72 or more beams are required to keep irradiation nonuniformities low enough to obtain fusion yields close to that of 1D simulations.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Calibration of reactive burn and Jones-Wilkins-Lee parameters for simulations of a detonation-driven flow experiment with uncertainty quantification
Here, uncertainties in the explosive-specific parameters of the Jones-Wilkins-Lee (JWL) equation of state (EOS) are carefully considered in hydrodynamic simulations of an explosive experiment to minimize the error in the flow prediction. Experimental data of the leading shock position in the transverse direction over time serves as the prediction metric for quantifying simulation prediction error. The uncertainty quantification technique, global sensitivity analysis, is utilized to determine the JWL parameters to which the transverse shock propagation is most sensitive. A polynomial response surface (PRS) is constructed in the space of the most influential JWL parameters, and the point of minimum error between the experimental data and the PRS yields calibrated JWL parameters for the experimental flow. The simulation results following the parameter calibration show good agreement with the experimental data. It was found that two significant parameters, the heat release per unit mass of reactant Q and JWL model exponent R 1 are strongly related, which makes it difficult to identify accurate values.
36 MATERIALS SCIENCE↗
Visualizations of direct fuel injection effects in a supersonic cavity flameholder
Not Available
MATHEMATICS AND COMPUTING↗