DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Sensitivity of pore collapse heating to the melting temperature and shear viscosity of HMX

Abstract

A multiscale modeling strategy is used to quantify factors governing the temperature rise in hot spots formed by pore collapse from supported and unsupported shock waves in the high explosive HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). Two physical aspects are examined in detail, namely the melting temperature and liquid shear viscosity. All-atom molecular dynamics simulations of phase coexistence are used to predict the pressure-dependent melting temperature up to 5 GPa. Equilibrium simulations and the Green–Kubo formalism are used to obtain the temperature- and pressure-dependent liquid shear viscosity. Starting from a simplified continuum-based grain-scale model of HMX, in this study we systematically increase the complexity of treatments for the solid–liquid phase transition and liquid shear viscosity in simulations of pore collapse. Using a realistic pressure-dependent melting temperature completely suppresses melting for supported shocks, which is otherwise predicted when treating it as a constant determined at atmospheric pressure. Alternatively, melt pools form around collapsed pores when the pressure (and melting temperature) are reduced during the release stage of unsupported shocks. Capturing the pressure dependence of the shear viscosity increases the peak temperature of melt pools by hundreds of Kelvin through viscous work. The complicated interplay of the solid-phase plastic work, solid–liquid phase transition, and liquid-phase viscousmore » work identified here motivate taking a systematic approach to building increasingly complex grain-scale models.« less

Authors:
ORCiD logo [1]; ORCiD logo [1]
  1. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Publication Date:
Research Org.:
Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA)
OSTI Identifier:
1769080
Alternate Identifier(s):
OSTI ID: 1776453
Report Number(s):
LLNL-JRNL-809399
Journal ID: ISSN 0167-6636; 1015388
Grant/Contract Number:  
AC52-07NA27344
Resource Type:
Accepted Manuscript
Journal Name:
Mechanics of Materials
Additional Journal Information:
Journal Volume: 152; Journal ID: ISSN 0167-6636
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
36 MATERIALS SCIENCE; hot spots; melting; multiscale modeling; molecular crystals; shock waves; momentum transport

Citation Formats

Kroonblawd, Matthew P., and Austin, Ryan A. Sensitivity of pore collapse heating to the melting temperature and shear viscosity of HMX. United States: N. p., 2020. Web. doi:10.1016/j.mechmat.2020.103644.
Kroonblawd, Matthew P., & Austin, Ryan A. Sensitivity of pore collapse heating to the melting temperature and shear viscosity of HMX. United States. https://doi.org/10.1016/j.mechmat.2020.103644
Kroonblawd, Matthew P., and Austin, Ryan A. Mon . "Sensitivity of pore collapse heating to the melting temperature and shear viscosity of HMX". United States. https://doi.org/10.1016/j.mechmat.2020.103644. https://www.osti.gov/servlets/purl/1769080.
@article{osti_1769080,
title = {Sensitivity of pore collapse heating to the melting temperature and shear viscosity of HMX},
author = {Kroonblawd, Matthew P. and Austin, Ryan A.},
abstractNote = {A multiscale modeling strategy is used to quantify factors governing the temperature rise in hot spots formed by pore collapse from supported and unsupported shock waves in the high explosive HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). Two physical aspects are examined in detail, namely the melting temperature and liquid shear viscosity. All-atom molecular dynamics simulations of phase coexistence are used to predict the pressure-dependent melting temperature up to 5 GPa. Equilibrium simulations and the Green–Kubo formalism are used to obtain the temperature- and pressure-dependent liquid shear viscosity. Starting from a simplified continuum-based grain-scale model of HMX, in this study we systematically increase the complexity of treatments for the solid–liquid phase transition and liquid shear viscosity in simulations of pore collapse. Using a realistic pressure-dependent melting temperature completely suppresses melting for supported shocks, which is otherwise predicted when treating it as a constant determined at atmospheric pressure. Alternatively, melt pools form around collapsed pores when the pressure (and melting temperature) are reduced during the release stage of unsupported shocks. Capturing the pressure dependence of the shear viscosity increases the peak temperature of melt pools by hundreds of Kelvin through viscous work. The complicated interplay of the solid-phase plastic work, solid–liquid phase transition, and liquid-phase viscous work identified here motivate taking a systematic approach to building increasingly complex grain-scale models.},
doi = {10.1016/j.mechmat.2020.103644},
journal = {Mechanics of Materials},
number = ,
volume = 152,
place = {United States},
year = {Mon Nov 02 00:00:00 EST 2020},
month = {Mon Nov 02 00:00:00 EST 2020}
}

Works referenced in this record:

Non-Schmid effect of pressure on plastic deformation in molecular crystal HMX
journal, June 2019

  • Pal, Anirban; Picu, Catalin R.
  • Journal of Applied Physics, Vol. 125, Issue 21
  • DOI: 10.1063/1.5092285

Theoretical study of the defect evolution for molecular crystal under shock loading
journal, February 2019

  • Long, Yao; Chen, Jun
  • Journal of Applied Physics, Vol. 125, Issue 6
  • DOI: 10.1063/1.5067284

Multiscale modeling of shock wave localization in porous energetic material
journal, January 2018


Ultrafast dynamic response of single-crystal β -HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine)
journal, May 2018

  • Zaug, Joseph M.; Austin, Ryan A.; Armstrong, Michael R.
  • Journal of Applied Physics, Vol. 123, Issue 20
  • DOI: 10.1063/1.5029923

Elastic–plastic wave profiles in cyclotetramethylene tetranitramine crystals
journal, July 2004

  • Dick, J. J.; Hooks, D. E.; Menikoff, R.
  • Journal of Applied Physics, Vol. 96, Issue 1
  • DOI: 10.1063/1.1757026

A comparison of methods for melting point calculation using molecular dynamics simulations
journal, April 2012

  • Zhang, Yong; Maginn, Edward J.
  • The Journal of Chemical Physics, Vol. 136, Issue 14
  • DOI: 10.1063/1.3702587

Direct numerical simulation of shear localization and decomposition reactions in shock-loaded HMX crystal
journal, May 2015

  • Austin, Ryan A.; Barton, Nathan R.; Reaugh, John E.
  • Journal of Applied Physics, Vol. 117, Issue 18
  • DOI: 10.1063/1.4918538

Tandem Molecular Dynamics and Continuum Studies of Shock‐Induced Pore Collapse in TATB
journal, November 2019

  • Zhao, Puhan; Lee, Sangyup; Sewell, Tommy
  • Propellants, Explosives, Pyrotechnics, Vol. 45, Issue 2
  • DOI: 10.1002/prep.201900382

Initiation of Detonation by the Interaction of Shocks with Density Discontinuities
journal, January 1965


Constituent properties of HMX needed for mesoscale simulations
journal, March 2002


Understanding the shock and detonation response of high explosives at the continuum and meso scales
journal, March 2018

  • Handley, C. A.; Lambourn, B. D.; Whitworth, N. J.
  • Applied Physics Reviews, Vol. 5, Issue 1
  • DOI: 10.1063/1.5005997

Fast Parallel Algorithms for Short-Range Molecular Dynamics
journal, March 1995


Visualization and analysis of atomistic simulation data with OVITO–the Open Visualization Tool
journal, December 2009


Atomic-level view of inelastic deformation in a shock loaded molecular crystal
journal, August 2007


High Explosive Ignition through Chemically Activated Nanoscale Shear Bands
journal, May 2020


A generalized crystal-cutting method for modeling arbitrarily oriented crystals in 3D periodic simulation cells with applications to crystal–crystal interfaces
journal, October 2016

  • Kroonblawd, Matthew P.; Mathew, Nithin; Jiang, Shan
  • Computer Physics Communications, Vol. 207
  • DOI: 10.1016/j.cpc.2016.07.007

Predicted melt curve and liquid-state transport properties of TATB from molecular dynamics simulations
journal, January 2018


Time-Correlation Functions and Transport Coefficients in Statistical Mechanics
journal, October 1965


Temperature-dependent shear viscosity coefficient of octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX): A molecular dynamics simulation study
journal, April 2000

  • Bedrov, Dmitry; Smith, Grant D.; Sewell, Thomas D.
  • The Journal of Chemical Physics, Vol. 112, Issue 16
  • DOI: 10.1063/1.481285

Dynamic fracture and hot-spot modeling in energetic composites
journal, February 2018

  • Grilli, Nicolò; Duarte, Camilo A.; Koslowski, Marisol
  • Journal of Applied Physics, Vol. 123, Issue 6
  • DOI: 10.1063/1.5009297

Hot spot ignition mechanisms for explosives
journal, November 1992


Thermal decomposition and phase transitions in solid nitramines
journal, January 1971


Modeling The Effects of Shock Pressure and Pore Morphology on Hot Spot Mechanisms in HMX
journal, June 2018

  • Springer, H. Keo; Bastea, Sorin; Nichols, Albert L.
  • Propellants, Explosives, Pyrotechnics, Vol. 43, Issue 8
  • DOI: 10.1002/prep.201800082

Rattle: A “velocity” version of the shake algorithm for molecular dynamics calculations
journal, October 1983


Defect evolution and pore collapse in crystalline energetic materials
journal, February 2009

  • Barton, Nathan R.; Winter, Nicholas W.; Reaugh, John E.
  • Modelling and Simulation in Materials Science and Engineering, Vol. 17, Issue 3
  • DOI: 10.1088/0965-0393/17/3/035003

Thermodynamics of HMX Polymorphs and HMX/RDX Mixtures
journal, December 2016

  • Myint, Philip C.; Nichols, Albert L.
  • Industrial & Engineering Chemistry Research, Vol. 56, Issue 1
  • DOI: 10.1021/acs.iecr.6b03697

Pressure‐dependent Elastic Coefficients of β‐HMX from Molecular Simulations
journal, February 2018

  • Mathew, Nithin; Sewell, Tommy
  • Propellants, Explosives, Pyrotechnics, Vol. 43, Issue 3
  • DOI: 10.1002/prep.201700286

Quantum Chemistry Based Force Field for Simulations of HMX
journal, May 1999

  • Smith, Grant D.; Bharadwaj, Rishikesh K.
  • The Journal of Physical Chemistry B, Vol. 103, Issue 18
  • DOI: 10.1021/jp984599p

Generalized Lindemann Melting Law
journal, August 1969


The crystal structure of α-HMX and a refinement of the structure of β-HMX
journal, July 1963