A numerical study on mesoscale simulation of quartzite and sandstone under shock loading

In this paper, we present two numerical models for the mesoscale (grain scale) simulation of planar shock waves in quartzite and sandstone using the in-house hydrocode SOPHIA. The models are compared in terms of their capability to represent physical mechanisms, such as phase transitions in quartz and pore collapse in sandstone, and they are validated by comparison to literature data. The study is part of the MEMIN (Multidisciplinary Experimental and Modeling Impact Research Network) project, which is devoted to the experimental and numerical investigation of the effects of meteorite impact on geological materials from laboratory scale to natural scale. The first model is based on the Smoothed Particle Hydrodynamics (SPH) method. Simulations with rather simplified structures in planar symmetry are presented. The model is used to investigate basic effects of porosity, pore geometry and water saturation. The second model presented is a more detailed, three-dimensional Finite Element (FE) model. With this model, the effects of grain anisotropy and different types of shear strength modeling are studied. In a parameter study, we investigate the influence of these parameters on shock Hugoniot relations, such as shock velocity (U-s) vs. particle velocity (U-p) and compressive longitudinal stress (sigma(L)) vs. U-p. Finally, the models are compared and the specific advantages and disadvantages of the different modeling variants are outlined. (C) 2017 Elsevier Ltd. All rights reserved.