Mesoscale simulations of shock compaction of a granular ceramic: effects of mesostructure and mixed-cell strength treatment

Abstract

The shock response of granular materials is important in a variety of contexts but the precise dynamics of grains during compaction is poorly understood. Here we use 2D mesoscale numerical simulations of the shock compaction of granular tungsten carbide to investigate the effect of internal structure within the particle bed and ‘stiction’ between grains on the shock response. An increase in the average number of contacts with other particles, per particle, tends to shift the Hugoniot to higher shock velocities, lower particle velocities and lower densities. This shift is sensitive to inter-particle shear resistance. Eulerian shock physics codes approximate friction between, and interlocking of, grains with their treatment of mixed cell strength (stiction) and here we show that this has a significant effect on the shock response. When studying the compaction of particle beds it is not common to quantify the pre-compaction internal structure, yet our results suggest that such differences should be taken into account, either by using identical beds or by averaging results over multiple experiments.