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.
Demaske, B.; Hudspeth, M.; Mandal, A.; Jensen, B.; Crum, R.; Vogler, Tracy V.
Mesoscale simulations of a LiF impactor colliding with a PMMA capsule containing aluminum powder (r00 = 1.5 g/cc) have been performed to investigate shock-induced melting in porous metals. Impact velocities of 1-2.5 km/s are chosen to coincide with in situ X-ray diffraction experiments, which provide direct evidence of shock-induced melting in aluminum powders. Mesoscale simulations show shock heating within the powder is highly nonuniform and melting remains incomplete over hundreds of nanoseconds behind the shock front despite equilibrium pressure-temperature states from continuum simulations lying above the experimental melt line. Such incomplete melting behavior is consistent with X-ray diffraction data obtained in experiment. For an impact velocity of ~1 km/s, mesoscale simulations predict re-solidification behind the shock front as high-temperature regions are cooled below the melt line. Reducing the grain size of the powder by a factor of two leads to a reduction in the time required to reach complete melt such that total melting of the powder may be observed experimentally for an impact velocity of 2.42 km/s.
A novel experimental methodology is presented to study the deviatoric response of powders in shock regimes. The powders are confined to a cylindrical wedge volume, and a projectile-driven shock wave with a sinusoidally varying front propagates through the powder. The perturbed shock wave exhibits a damping behavior due to irreversible processes of viscosity and strength (deviatoric) of the powder with propagation through increasing powder thicknesses. The inclined surface of the wedge is polished and coated to establish a diffuse surface suitable for reflecting incident laser light into a high-speed camera imaging at 5 MHz. Images of the contrast loss upon shock wave arrival at the observation surface are post-processed for qualitative and quantitative information. New data of shock damping behavior with parameters of perturbation wavelength and initial shock strength are presented for powders of copper, tantalum, and tungsten carbide as well as their mixtures. We present the first full-field images showing additional spatial disturbances on the perturbed shock front that appear dependent on particle material and morphology.
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.
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.
CTH is an Eulerian hydrocode developed by Sandia National Laboratories (SNL) to solve a wide range of shock wave propagation and material deformation problems. Adaptive mesh refinement is also used to improve efficiency for problems with a wide range of spatial scales. The code has a history of running on a variety of computing platforms ranging from desktops to massively parallel distributed-data systems. For the Trinity Phase 2 Open Science campaign, CTH was used to study mesoscale simulations of the hypervelocity penetration of granular SiC powders. The simulations were compared to experimental data. A scaling study of CTH up to 8192 KNL nodes was also performed, and several improvements were made to the code to improve the scalability.