On resolving meso-scale calculations of pore-collapse-generated hotspots in energetic crystals for consistency with atomistic models

Meso-scale calculations of pore collapse and hotspot formation in energetic crystals provide closure models to macro-scale hydrocodes for predicting the shock sensitivity of energetic materials. To this end, previous works obtained atomistics-consistent material models for two common energetic crystals, HMX and RDX, such that pore collapse calculations adhered closely to molecular dynamics (MD) results on key features of energy localization, particularly the profiles of the collapsing pores, appearance of shear bands, and the transition from viscoplastic to hydrodynamic collapse. However, some important aspects such as the temperature distributions in the hotspot were not as well captured. One potential issue was noted but not resolved adequately in those works, namely the grid resolution that should be employed in the meso-scale calculations for various pore sizes and shock strengths. Conventional computational mechanics guidelines for selecting meshes as fine as possible, balancing computational effort, accuracy and grid independence, were shown not to produce physically consistent features associated with shear localization. Here, we examine the physics of pore collapse, shear band evolution and structure, and hotspot formation, for both HMX and RDX; we then evaluate under what conditions atomistics-consistent models yield physically correct (considering MD as ground truth) hotspots for a range of pore diameters, from nm to microns, and for a wide range of shock strengths. The study provides insights into the effects of pore size and shock strength on pore collapse and hotspots, identifying aspects such as size-independent behaviors, and proportion of energy contained in shear as opposed to jet impact-heated regions of the hotspot. Areas for further improvement of atomistics-consistent material models are also indicated.