MOLECULAR DYNAMICS SIMULATIONS OF HIGH VELOCITY SHOCK COMPRESSED TNT

Liu Hai; Li Qikaiy; He Yuanhang

doi:10.6052/0459-1879-14-141

Volume 47 Issue 1

Turn off MathJax

Article Contents

Article Navigation > Chinese Journal of Theoretical and Applied Mechanics > 2015 > 47(1): 174-179

Liu Hai, Li Qikaiy, He Yuanhang. MOLECULAR DYNAMICS SIMULATIONS OF HIGH VELOCITY SHOCK COMPRESSED TNT[J]. Chinese Journal of Theoretical and Applied Mechanics, 2015, 47(1): 174-179. doi: 10.6052/0459-1879-14-141

Citation:

Liu Hai, Li Qikaiy, He Yuanhang. MOLECULAR DYNAMICS SIMULATIONS OF HIGH VELOCITY SHOCK COMPRESSED TNT[J]. Chinese Journal of Theoretical and Applied Mechanics, 2015, 47(1): 174-179. doi: 10.6052/0459-1879-14-141

Citation:

PDF( 10289 KB)

MOLECULAR DYNAMICS SIMULATIONS OF HIGH VELOCITY SHOCK COMPRESSED TNT

doi: 10.6052/0459-1879-14-141

1.
State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology, Beijing 100081, China
2. School of Materials Science and Engineering, Tsinghua University, Beijing 100086, China

Received Date: 2014-08-13
Rev Recd Date: 2014-09-24
Publish Date: 2015-01-18

Abstract

Abstract

We simulate the shock compression behavior of TNT with ReaxFF-MD. When shock compression is complete, all of the TNT molecules are decomposed, and when volume compression is up to the 40% of original volume, pressure of the system reaches a peak. Close behind is rarefaction wave reverse stretching the compressed energetic materials and leading to a large number of atoms or molecules group splash to the downstream, pressure begin to unload at the same time. Density and particle wave velocity profile show a greater density in the compressed region, and the particles in a stationary state, but sharp velocity gradient in the region of compression wave. In the earlier chemical characteristics, TNT molecules shed the H, O atoms under the effect of shock compression, and then the residues aggregate to the larger clusters, and this phase associated with translational-vibrational relaxation processes. The rotational mode is subsequently transferred into vibrational modes with a time scale of 0.5 ps. Fragment analysis shows that a large number of C—H, O＝N bonds rupture to form the OH, H₂, H₂O, N₂ groups and parts of H, O atoms are free in the system. The molar mass of the carbon-containing clusters under the joint actions of compressional wave ahead and rear compression is accumulating gradually from the analysis. Atomic ratio in the carbon-containing clusters tends to balance (O/C=0.680, H/C=0.410, N/C=0.284), but less than the ratio in the initial structure.
- shock compression,
- TNT,
- reactive force field,
- molecular dynamics

FullText(HTML)

References(8)

References

Van Duin ACT, Dasgupta S, Lorant F, et al. ReaxFF: a reactive force field for hydrocarbons. J Phys Chem A, 2001, 105(41): 9396-9409

Strachan A, Kober EM, van Duin ACT, et al. Thermal decomposition of RDX from reactive molecular dynamics. J Chem Phys , 2005, 122: 054502

Agrawalla S, van Duin ACT. Development and application of a ReaxFF reactive force field for hydrogen combustion. J Phys Chem A, 2011, 115(6): 960-972

Nomura K, Kalia RK, Nakano A, et al. Dynamic transition in the structure of an energetic crystal during chemical reactions at shock front prior to detonation. Phys Rev Lett , 2007, 99: 148303

Liu LC, Liu Y, Zybin SV, et al. ReaxFF-lg: correction of the ReaxFF reactive force field for london dispersion, with applications to the equations of state for energetic materials. J Phys Chem A, 2011, 115(40): 11016-11022

Wen YS, Xue XG, Zhou XQ, et al. Twin induced sensitivity enhancement of hmx versus shock: a molecular reactive force field simulation. J Phys Chem C , 2013, 117: 24368-24374

Plimpton S. Fast parallel algorithms for short-range molecular dynamics. J Comp Phys 1995, 117(1): 1-19 http://lammps.sandia.gov.

Dremin AN. Shock discontinuity zone effect: the main factor in the explosive decomposition detonation process. Philosophical Transactions of the Royal Society of London. Series A: Physical and Engineering Sciences, 1992, 339(1654): 355-364

Relative Articles

Supplements(0)

Cited By

Proportional views