Multiscale Modeling of Shock Wave Localization in Porous Energetic Material

Mitchell A. Wood,David E. Kittell,Cole D. Yarrington,Aidan P. Thompson

DOI: https://doi.org/10.1103/PhysRevB.97.014109

2017-11-08

Abstract:Shock wave interactions with defects, such as pores, are known to play a key role in the chemical initiation of energetic materials. The shock response of hexanitrostilbene is studied through a combination of large scale reactive molecular dynamics and mesoscale hydrodynamic simulations. In order to extend our simulation capability at the mesoscale to include weak shock conditions (< 6 GPa), atomistic simulations of pore collapse are used to define a strain rate dependent strength model. Comparing these simulation methods allows us to impose physically-reasonable constraints on the mesoscale model parameters. In doing so, we have been able to study shock waves interacting with pores as a function of this viscoplastic material response. We find that the pore collapse behavior of weak shocks is characteristically different to that of strong shocks.

Materials Science,Computational Physics

What problem does this paper attempt to address?

### Problems the Paper Attempts to Solve

The main purpose of this paper is to study the behavior of shock wave localization in porous explosive materials through multi-scale modeling. Specifically:

1. **Interaction between Shock Waves and Defects**:
- Investigate how shock waves interact with defects (such as pores) in porous explosive materials.
- This interaction is crucial for the chemical initiation process.

2. **Improvement of Multi-Scale Simulation Methods**:
- Combine large-scale reactive molecular dynamics (MD) and mesoscale hydrodynamic simulations to extend simulation capabilities.
- Define a strain rate-dependent strength model by simulating pore collapse behavior at the atomic scale.

3. **Establishing Physically Reasonable Model Parameters**:
- Compare results from different simulation methods to impose physically reasonable constraints to optimize mesoscale model parameters.
- Study the differences in pore collapse behavior under weak shock conditions (<6 GPa) and strong shock conditions.

4. **Understanding Thermochemical Initiation Mechanisms**:
- Explore the mechanisms of hot spot formation during pore collapse and their impact on the sensitivity of explosive materials.
- Capture different hot spot formation mechanisms through large-scale MD simulations and use them to train mesoscale plasticity models.

In summary, this paper aims to improve mesoscale simulation methods by combining experimental data and atomic-scale simulation results, thereby better understanding and predicting the behavior of porous explosive materials under shock wave action.