The high-pressure dynamic response of titanium dioxide (TiO 2) is not only of interest because of its numerous industrial applications but also because of its structural similarities to silica (SiO 2). We performed plate impact experiments in a two-stage light gas gun, at peak stresses from 64 to 221 GPa to determine the TiO 2 response along the Hugoniot. The lower stress experiment at 64 GPa shows an elastic behavior followed by an elastic-plastic transition, whereas the high stress experiments above 64 GPa show a single wave structure. Previous shock studies have shown the presence of high-pressure phases (HPP) I (26 GPa) and HPP II (100 GPa); however, our data suggest that the HPP I phase is stable up to 150 GPa. Using a combination of data from our current study and our previous Z-data, we determine that TiO 2 likely melts on the Hugoniot at 157 GPa. Furthermore, our data confirm that TiO 2 is not highly incompressible as shown by a previous study.
The strength of brittle porous media is of concern in numerous applications, for example, earth penetration, crater formation, and blast loading. Thus, it is of importance to possess techniques that allow for constitutive model calibration within the laboratory setting. The goal of the current work is to demonstrate an experimental technique allowing for strength assessment of porous media subjected to shock loading, which can be implemented into pressure-dependent yield surfaces within numerical simulation schemes. As a case study, the deviatoric response of distended α-SiO2 has been captured in a tamped Richtmyer-Meshkov instability (RMI) environment at a pressure regime of 4-10 GPa. Hydrocode simulations were used to interpret RMI experimental data, and a resulting pressure-dependent yield surface akin to the often employed modified Drucker-Prager model was calibrated. Simulations indicate that the resulting jet length generated by the RMI is sensitive to the porous media strength, thereby providing a feasible experimental platform capable of capturing the pressurized granular deviatoric response. Furthermore, in efforts to validate the RMI-calibrated strength model, a set of Mach-lens experiments was performed and simulated with the calibrated pressure-dependent yield surface. Excellent agreement between the resulting Mach-lens length in experiment and simulation provides additional confidence to the RMI yield-surface calibration scheme.
We discuss major challenges in modeling giant impacts between planetary bodies, focusing on the equations of state (EOS). During the giant impact stage of planet formation, rocky planets are melted and partially vaporized. However, most EOS models fail to reproduce experimental constraints on the thermodynamic properties of the major minerals over the required phase space. Here, we present an updated version of the widely-used ANEOS model that includes a user-defined heat capacity limit in the thermal free energy term. Our revised model for forsterite (Mg2SiO4), a common proxy for the mantles of rocky planets, provides a better fit to material data over most of the phase space of giant impacts. We discuss the limitations of this model and the Tillotson equation of state, a commonly used alternative model.
We report the atomic- and nanosecond-scale quantification of kinetics of a shock-driven phase transition in Zr metal. We uniquely make use of a multiple shock-and-release loading pathway to shock Zr into the β phase and to create a quasisteady pressure and temperature state shortly after. Coupling shock loading with in situ time-resolved synchrotron x-ray diffraction, we probe the structural transformation of Zr in the steady state. Our results provide a quantified expression of kinetics of formation of β-Zr phase under shock loading: transition incubation time, completion time, and crystallization rate.
Understanding and quantifying the uncertainties in experimental results are crucial to properly interpreting simulations based on those results. While methods are reasonably well established for estimating those uncertainties in high-pressure shock experiments on homogeneous materials, it is much more difficult to treat relatively low-pressure experiments where shock rise times are significant and material strength is not negligible. Sample heterogeneity further complicates the issue, especially when that heterogeneity is not characterized in each sample. Here, we extend the Monte Carlo impedance matching approach used in high-pressure Z experiments to low-pressure experiments on heterogeneous porous materials. The approach incorporates uncertainties not only in the equation of state of the impedance matching standard but also those associated with its strength. In addition, we also examine approaches for determining material heterogeneity and evaluate its effect on the experimental results.
Townsend, Joshua P.; Root, Seth R.; Davies, Erik D.; Lemke, Raymond W.; Bliss, David E.; Fratanduono, Dayne E.; Kraus, Richard K.; Millot, Marius M.; Spaulding, Dylan S.; Stewart, Sarah S.; Jacobsen, Stein J.; Shulenburger, Luke S.
Root, Seth R.; Townsend, Joshua P.; Davies, Erik D.; Lemke, Raymond.(.; Bliss, David E.; Fratanduono, Dayne E.; Kraus, Richard K.; Millot, Marius M.; Spaulding, Dylan S.; Shulenburger, Luke S.; Stewart, Sarah S.; Jacobsen, Stein J.
The shock Hugoniot for full-density and porous CeO2 was investigated in the liquid regime using ab initio molecular dynamics (AIMD) simulations with Erpenbeck's approach based on the Rankine-Hugoniot jump conditions. The phase space was sampled by carrying out NVT simulations for isotherms between 6000 and 100 000 K and densities ranging from ρ=2.5 to 20g/cm3. The impact of on-site Coulomb interaction corrections +U on the equation of state (EOS) obtained from AIMD simulations was assessed by direct comparison with results from standard density functional theory simulations. Classical molecular dynamics (CMD) simulations were also performed to model atomic-scale shock compression of larger porous CeO2 models. Results from AIMD and CMD compression simulations compare favorably with Z-machine shock data to 525 GPa and gas-gun data to 109 GPa for porous CeO2 samples. Using results from AIMD simulations, an accurate liquid-regime Mie-Grüneisen EOS was built for CeO2. In addition, a revised multiphase SESAME-Type EOS was constrained using AIMD results and experimental data generated in this work. This study demonstrates the necessity of acquiring data in the porous regime to increase the reliability of existing analytical EOS models.
Here, we report real-time observations of a phase transition in the ionic solid CaF2, a model AB2 structure in high-pressure physics. Synchrotron x-ray diffraction coupled with dynamic loading to 27.7 GPa, and separately with static compression, follows, in situ, the fluorite to cotunnite structural phase transition, both on nanosecond and on minute time scales. Using Rietveld refinement techniques, we examine the kinetics and hysteresis of the transition. Our results give insight into the kinetic time scale of the fluorite-cotunnite phase transition under shock compression, which is relevant to a number of isomorphic compounds.
Measurement of the Hugoniot and sound velocity provides information on the bulk modulus and Grüneisen parameter of a material at extreme conditions. The capability to launch multilayered (copper/aluminum) flyer plates at velocities in excess of 20 km/s with the Sandia Z accelerator has enabled high-precision sound-velocity measurements at previously inaccessible pressures. For these experiments, the sound velocity of the copper flyer must be accurately known in the multi-Mbar regime. Here we describe the development of copper as an absolutely calibrated sound-velocity standard for high-precision measurements at pressures in excess of 400 GPa. Using multilayered flyer plates, we performed absolute measurements of the Hugoniot and sound velocity of copper for pressures from 500 to 1200 GPa. These measurements enabled the determination of the Grüneisen parameter for dense liquid copper, clearly showing a density dependence above the melt transition. Combined with earlier data at lower pressures, these results constrain the sound velocity as a function of pressure, enabling the use of copper as a Hugoniot and sound-velocity standard for pressures up to 1200 GPa.
In this study, we have developed a conceptual design of a next-generation pulsed-power accelerator that is optmized for driving megajoule-class dynamic-material-physics experiments at pressures as high as 1 TPa. The design is based on an accelerator architecture that is founded on three concepts: single-stage electrical-pulse compression, impedance matching, and transit-time-isolated drive circuits. Since much of the accelerator is water insulated, we refer to this machine as Neptune. The prime power source of Neptune consists of 600 independent impedance-matched Marx generators. As much as 0.8 MJ and 20 MA can be delivered in a 300-ns pulse to a 16-mΩ physics load; hence Neptune is a megajoule-class 20-MA arbitrary waveform generator. Neptune will allow the international scientific community to conduct dynamic equation-of-state, phase-transition, mechanical-property, and other material-physics experiments with a wide variety of well-defined drive-pressure time histories. Because Neptune can deliver on the order of a megajoule to a load, such experiments can be conducted on centimeter-scale samples at terapascal pressures with time histories as long as 1 μs.
Mixtures of light elements with heavy elements are important in inertial confinement fusion. We explore the physics of molecular scale mixing through a validation study of equation of state (EOS) properties. Density functional theory molecular dynamics (DFT-MD) at elevated temperature and pressure is used to obtain the thermodynamic state properties of pure xenon, ethane, and various compressed mixture compositions along their principal Hugoniots. To validate these simulations, we have performed shock compression experiments using the Sandia Z-Machine. A bond tracking analysis correlates the sharp rise in the Hugoniot curve with the completion of dissociation in ethane. The DFT-based simulation results compare well with the experimental data along the principal Hugoniots and are used to provide insight into the dissociation and temperature along the Hugoniots as a function of mixture composition. Interestingly, we find that the compression ratio for complete dissociation is similar for several compositions suggesting a limiting compression for C-C bonded systems.
We use Sandia's Z machine and magnetically accelerated flyer plates to shock compress liquid krypton to 850 GPa and compare with results from density-functional theory (DFT) based simulations using the AM05 functional. We also employ quantum Monte Carlo calculations to motivate the choice of AM05. We conclude that the DFT results are sensitive to the quality of the pseudopotential in terms of scattering properties at high energy/temperature. A new Kr projector augmented wave potential was constructed with improved scattering properties which resulted in excellent agreement with the experimental results to 850 GPa and temperatures above 10 eV (110 kK). Finally, we present comparisons of our data from the Z experiments and DFT calculations to current equation of state models of krypton to determine the best model for high energy-density applications.
Planar shock experiments were conducted on granular tungsten carbide (WC) and tantalum oxide (Ta{sub 2}O{sub 5}) using the Z machine and a 2-stage gas gun. Additional shock experiments were also conducted on a nearly fully dense form of Ta{sub 2}O{sub 5}. The experiments on WC yield some of the highest pressure results for granular materials obtained to date. Because of the high distention of Ta{sub 2}O{sub 5}, the pressures obtained were significantly lower, but the very high temperatures generated led to large contributions of thermal energy to the material response. These experiments demonstrate that the Z machine can be used to obtain accurate shock data on granular materials. The data on Ta{sub 2}O{sub 5} were utilized in making improvements to the P-{lambda} model for high pressures; the model is found to capture the results not only of the Z and gas gun experiments but also those from laser experiments on low density aerogels. The results are also used to illustrate an approach for generating an equation of state using only the limited data coming from nanoindentation. Although the EOS generated in this manner is rather simplistic, for this material it gives reasonably good results.
Hydrocarbon foams are commonly used in HEDP experiments, and are subject to shock compression from tens to hundreds of GPa. Modeling foams is challenging due to the heterogeneous character of the foam. A quantitative understanding of foams under strong dynamic compression is sought. We use Sandia's ALEGRA-MHD code to simulate 3D mesoscale models of pure poly(4-methyl-1-petene) (PMP) foams. We employ two models of the initial polymer-void structure of the foam and analyze the statistical properties of the initial and shocked states. We compare the simulations to multi-Mbar shock experiments at various initial foam densities and flyer impact velocities. Scatter in the experimental data may be a consequence of the initial foam inhomogeneity. We compare the statistical properties the simulations with the scatter in the experimental data.
The noble gas xenon is a particularly interesting element. At standard pressure xenon is an fcc solid which melts at 161 K and then boils at 165 K, thus displaying a rather narrow liquid range on the phase diagram. On the other hand, under pressure the melting point is significantly higher: 3000 K at 30 GPa. Under shock compression, electronic excitations become important at 40 GPa. Finally, xenon forms stable molecules with fluorine (XeF{sub 2}) suggesting that the electronic structure is significantly more complex than expected for a noble gas. With these reasons in mind, we studied the xenon Hugoniot using DFT/QMD and validated the simulations with multi-Mbar shock compression experiments. The results show that existing equation of state models lack fidelity and so we developed a wide-range free-energy based equation of state using experimental data and results from first-principles simulations.