Abstract:Hotspot formation due to pore collapse is a key mechanism for initiating detonation of shocked energetic materials. Energy localization at and around the pore collapse site leads to high-temperature hotspots, initiating chemical reactions. Because chemical reaction rates depend sensitively on temperature, predictive continuum models need to get the pore-collapse dynamics and resulting hotspot temperatures right; this imposes stringent demands on the fidelity of thermophysical model forms and parameters, and on the numerical methods employed to perform high-resolution meso-scale calculations. Here, continuum material models for beta-HMX are examined in the context of nanoscale shock-induced pore collapse, treating predictions from molecular dynamics (MD) simulations as ground truth. Using MD-consistent material properties, we show that the currently available strength models for HMX fail to correctly capture pore collapse and hotspot temperatures. Insights from MD are then employed to advance a Modified Johnson-Cook (M-JC) strength model form that captures aspects of shear strain and strain-rate dependency not represented by the standard JC form, but which are shown to be critical for accurately describing the nanoscale physics of shock-induced localization in HMX. The study culminates in a fully MD-determined strength model for beta-HMX that produces continuum pore-collapse results well aligned in all aspects with those predicted by MD, including pore-collapse mechanism and rate, shear-band formation in the collapse zone, and temperature, strain, and stress fields in the hotspot zone and surrounding material. The resulting MD-informed/MD-determined M-JC model should improve the fidelity of simulations to predict the detonation initiation of HMX-based energetic materials in microstructure-aware multi-scale frameworks.

What problem does this paper attempt to address?

This paper mainly discusses the key mechanism of hot spot formation in high-energy crystal β-1,3,5,7-tetranitro-1,3,5,7-tetraazacyclooctane (β-HMX) under impact, which is the high-temperature hot spot generated by pore collapse. The study indicates that accurate prediction of pore collapse dynamics and the resulting hot spot temperature is crucial for understanding and simulating the impact-to-detonation transition of explosive materials. The existing strength models fail to accurately capture the pore collapse and hot spot temperature of β-HMX. Through molecular dynamics (MD) simulations, the researchers found that it is necessary to improve the Johnson-Cook (JC) strength model to include shear strain and strain rate dependency, which are crucial but not fully considered in the standard JC form, in describing the nanoscale impact-induced localization processes in β-HMX. Ultimately, they proposed a modified Johnson-Cook (M-JC) model, which is solely based on MD-determined material properties and highly consistent with the pore collapse results predicted by MD, including collapse mechanisms, rates, shear band formation in the hot spot region, and temperature, strain, and stress fields. The paper emphasizes the importance of accurately predicting hot spot temperature for obtaining chemical reaction rates, which is crucial for predicting the initiation of detonation in HMX-based explosive materials within a multiscale framework with microstructure perception. The study also points out that although single crystal plasticity models can capture the direction-dependence of material strength, there are computational complexity and practicality issues in practical applications. Therefore, the study proposes an MD-based isotropic strength model to simplify the description and capture key aspects of the interaction between impact and HMX microstructure. In summary, the paper aims to address how to accurately predict the dynamics of pore collapse and high-temperature phenomena during hot spot formation in high-energy crystal β-HMX under impact loading by developing continuum models consistent with molecular dynamics, in order to improve the accuracy of predicting detonation initiation.

Atomistics-consistent continuum models and pore-collapse-generated hotspot temperatures in energetic crystals I: beta-1,3,5,7-tetranitro-1,3,5,7-tetrazocane (beta-HMX)

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

Shock-induced Nanoscale Pore Collapse and Hotspot in Cyclotetramethylene Tetranitramine (HMX)

Continuum and molecular dynamics simulations of pore collapse in shocked β-tetramethylene tetranitramine (β-HMX) single crystals

Initial decomposition mechanism of solid-state HMX for steady shock waves

Anisotropic Responses and Initial Decomposition of Condensed-Phase Β-Hmx under Shock Loadings Via Molecular Dynamics Simulations in Conjunction with Multiscale Shock Technique.

A model for hot spot formation in shocked energetic materials

Phase Transitions and Mechanical Properties of Octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine in Different Crystal Phases by Molecular Dynamics Simulation

Crystal-scale Modelling of Pore Collapse in Cyclotetramethylene Tetranitramine (HMX) under Different Shock Strengths

Macro-scale sensitivity through meso-scale hotspot dynamics in porous energetic materials: Comparing the shock response of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) and 1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (HMX)

Hot Spot Formation and Chemical Reaction Initiation in Shocked HMX Crystals with Nanovoids: a Large-Scale Reactive Molecular Dynamics Study.

Intergranular Hotspots: A Molecular Dynamics Study on the Influence of Compressive and Shear Work

The Potential Energy Hotspot: Effects from Impact Velocity, Defect Geometry, and Crystallographic Orientation

Thermal-Mechanical-Chemical Responses of Polymer-Bonded Explosives Using A Mesoscopic Reactive Model under Impact Loading

Mechanisms of shock-induced initiation at micro-scale defects in energetic crystal-binder systems

Johnson–Cook yield functions for cyclotetramethylene-tetranitramine (HMX) and cyclotrimethylene-trinitramine (RDX) derived from single crystal plasticity models

Multi-scale modeling of shock initiation of a pressed energetic material III: Effect of Arrhenius chemical kinetic rates on macro-scale shock sensitivity

Modeling meso-scale energy localization in shocked HMX, Part I: machine- learned surrogate model for effect of loading and void size

Shock Initiated Thermal and Chemical Responses of HMX Crystal from ReaxFF Molecular Dynamics Simulation.

Thermomechanical analysis of shock-induced hot spot formation in condensed energetic materials

Explosive chemistry: Simulating the chemistry of energetic materials at extreme conditions