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.