弱点火条件下炸药燃烧转爆轰的数值模拟

计算物理 ›› 2016, Vol. 33 ›› Issue (1): 15-22.

弱点火条件下炸药燃烧转爆轰的数值模拟

董贺飞, 洪滔, 张晓立

北京应用物理与计算数学研究所, 北京 100094

收稿日期:2014-11-30 修回日期:2015-04-20 出版日期:2016-01-25 发布日期:2016-01-25
作者简介:董贺飞(1981-),男,博士,副研究员,主要从事爆轰理论与数值模拟研究,E-mail:dong_hefei@iapcm.ac.cn
基金资助:
中物院科学技术发展基金（2013B0101014，2014A09015，2015B0101021）、中物院安全弹药研发中心（RMC-2014-ZZ-PG-B01）及国防基础科研计划（B1520132012）资助项目

Numerical Simulation of Deflagration to Detonation Transition in Explosives Under Weak Ignition

DONG Hefei, HONG Tao, ZHANG Xiaoli

Institute of Applied Physics and Computational Mathematics, Beijing 100094, China

Received:2014-11-30 Revised:2015-04-20 Online:2016-01-25 Published:2016-01-25

摘要/Abstract

摘要： 考察颗粒炸药从传导燃烧到对流燃烧再到爆轰的过程.对装填密度为85%的HMX颗粒炸药的燃烧转爆轰过程进行数值模拟，分析传导燃烧、对流燃烧和爆轰的发展过程.点火早期燃烧速度很低，火焰面在8.16 ms之内只前进了不到0.2 mm；形成对流燃烧之后燃烧速度快速增加，只用了0.1 ms就形成了速度为8 165 m·s^-1的稳定爆轰.当炸药颗粒直径或点火压力减小时，形成稳定爆轰所需的时间增加.

关键词: 燃烧转爆轰, 传导燃烧, 对流燃烧, 数值模拟

Abstract: By introducing conductive burning process into classical deflagration-to-detonation transition (DDT) model, transition process from low speed conductive burning to convective burning to detonation was proposed. Transition process in HMX granular bed with 85% loading density was simulated. Development of conductive burning, convective burning and detonation was analyzed. In early stage combustion propagation rate is very slow. It propagates no more than 0. 2mm within 8. 16ms. After onset of convective burning, it tooks 20 ms to form a steady detonation with a velocity of 8 165 m·s^-1. Time to form detonation increases with decrease of particle diameter and ignition pressure.

Key words: DDT, conductive burning, convective burning, numerical simulation

中图分类号:

董贺飞, 洪滔, 张晓立. 弱点火条件下炸药燃烧转爆轰的数值模拟[J]. 计算物理, 2016, 33(1): 15-22.

DONG Hefei, HONG Tao, ZHANG Xiaoli. Numerical Simulation of Deflagration to Detonation Transition in Explosives Under Weak Ignition[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2016, 33(1): 15-22.

参考文献

[1] BERNECKER P R, SANDUSKY H W, CLAIRMONT A R. Deflagration-to-detonation transition studies of a double-based propellant[C]. 8th International Detonation Symposium, New Mexico, USA, 1985: 658-668.
[2] BAER M R, NUNZIATO J W. A two-phase mixture theory for the deflagration-to-detonation transition (DDT) in reactive granular materials[J]. Int J Multiphase Flow, 1986, 12: 861-889.
[3] POWERS J M, STEWART D S, KRIER H. Theory of two-phase detonation-Part 1: Modeling[J]. Combust Flame, 1990, 80: 264-279.
[4] GONTHIER K A, POWERS J M. A high-resolution numerical method for a two-phase model of deflagration-to-detonation transition[J], J Comput Phys, 2000, 163: 376-433.
[5] YANG T, JIN Z M, ZHANG X B. A mathematical model of deflagration to detonation transition in highly loaded propellant beds[J]. Journal of Ballistics,1992,(2): 1-9.
[6] YANG T, JIN Z M. Numerical simulation of deflagration-to-detonation transition in highly packed granular propellant beds[J]. Journal of Propulsion Technology,1994,(6): 76-82.
[7] ZHANG X B, YUAN Y X. Modeling of DDT in high loaded propellant beds[J]. Journal of Ballistics,1996, 8: 16-25.
[8] QIN G C, HOU X, CHERN L Q, et al. Numerical simulation of deflagration-to-detonation transition in high energy propellant[J]. Journal of Solid Rocket Technology,2006, 29(3): 186-189.
[9] JIA X R, LI D X, SUN J S, et al. An analysis of the deflagrantion-to-detonation transition (DDT) in NEPE propellants[J]. Acta Armamentarii, 1997, 18: 46-51.
[10] DUAN B F. Study on the deflagration to detonation transition in insensitive industry explosives[M]. Beijing: China Water & Power Press, 2009.
[11] DONG H F, ZHAO Y H, HONG T. Numerical simulation of the deflagration-to-detonation transition behavior of explosive HMX[J]. Chinese Journal of High Pressure Physics, 2012, 26(6): 601-607.
[12] JACKSON S I. Deflagration phenomena in energetic materials: an overview[M]//ASAY B W. Shock Wave Science and Technology Reference Library, Vol. 5: Non-shock Initiation of Explosives, Springer, 2010.
[13] MARGOLIN A D, CHUIKO S V. Combustion instability of a porous charge with spontaneous penetration of the combustion products into the pores[J]. Combustion, Explosion, and Shock Waves, 1966, 2(3): 72-75.
[14] SON S F, BERGHOUT H L, BOLME C A, et al. Burn rate measurements of HMX, TATB, DHT, DAAF, and BTATz[J]. Proceedings of the Combustion Institute, 2000, 28: 919-924.
[15] GLASCOE E A, MAIENSCHEIN J L, LORENZ K T, et al. Deflagration measurements of three insensitive high explosives: LLM105, TATB, and DAAF[C]. 14th International Detonation Symposium, Idaho, USA, 2010.
[16] WEN S G, WANG S Q, HUANG W B, et al. An experimental study on deflagration-to-detonation transition in high-density composition[J]. Explosion and Shock Wave, 2007, 27(6): 567-571.
[17] WANG J, WEN S G. Experimental study on deflagration-to-detonation transition in two pressed high-density explosives[J]. Chinese Journal of High Pressure Physics, 2009, 23: 441-446.
[18] CHANG S C. The method of space-time conservation element and solution element-a new approach for solving the NavierStokes and Euler equations[J]. J Comput Phys, 1995, 119: 295-324.
[19] ZHANG D L, WANG J T, WANG G. High-order CE/SE method and applications[J]. Chinese Journal of Computational Physics, 2009, 26(2): 211-310.
[20] DONG H F, HONG T, ZHANG X L. Numerical simulation of explosive dust detonation with CE/SE method[J]. Chinese Journal of Computational Physics, 2012, 29(4): 495-502.
[21] GONTHIER K A, POWERS J M. A numerical investigation of transient detonation in granulated material[J]. Shock Waves, 1996, 6: 183-195.
[22] NARIN B. One-and two-dimensional numerical simulation of deflagration to detonation transition in solid energetic materials[D]. PhD Thesis, Middle East Technical University, 2010.