• ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

激波诱导火焰失稳与爆轰的条件研究

朱跃进 于蕾 张彭岗 潘振华 潘剑锋 董刚

黄亚峰, 田轩, 冯博, 王晓峰. 温压炸药爆炸性能实验研究[J]. 爆炸与冲击, 2016, 36(4): 573-576. doi: 10.11883/1001-1455(2016)04-0573-04
引用本文: 朱跃进, 于蕾, 张彭岗, 潘振华, 潘剑锋, 董刚. 激波诱导火焰失稳与爆轰的条件研究[J]. 爆炸与冲击, 2017, 37(4): 741-747. doi: 10.11883/1001-1455(2017)04-0741-07
Huang Yafeng, Tian Xuan, Feng Bo, Wang Xiaofeng. Experimental study on explosion performance of thermobaric explosive[J]. Explosion And Shock Waves, 2016, 36(4): 573-576. doi: 10.11883/1001-1455(2016)04-0573-04
Citation: Zhu Yuejin, Yu Lei, Zhang Penggang, Pan Zhenhua, Pan Jianfeng, Dong Gang. Conditions for shock wave induced flame instability and detonation[J]. Explosion And Shock Waves, 2017, 37(4): 741-747. doi: 10.11883/1001-1455(2017)04-0741-07

激波诱导火焰失稳与爆轰的条件研究

doi: 10.11883/1001-1455(2017)04-0741-07
基金项目: 

国家自然科学基金项目 11402102

国家自然科学基金项目 11372140

江苏省自然科学基金青年项目 BK20140524

江苏省博士后基金项目 1402013B

江苏大学高级专业人才科研启动基金项目 14JDG031

详细信息
    作者简介:

    朱跃进(1986-),男,博士,讲师,zyjwind@163.com

  • 中图分类号: O381

Conditions for shock wave induced flame instability and detonation

  • 摘要: 采用九阶WENO和十阶中心差分格式数值求解激波与火焰作用过程,考察了激波强度、火焰尺寸对激波与球形火焰作用过程的影响。结果表明,增大激波强度或火焰尺寸均可在流场中引发爆轰,但激波强度的影响更大,并且其引发的爆轰可使火焰迅速膨胀,放热率提高,从而影响燃烧特性;此外,爆轰波传播过程中会迅速消耗可燃预混气,合并原有的反射激波,并在流场中形成局部高压区,极大地改变流场结构。
  • 温压炸药(thermobaric explosive,TBX)是一类能够充分利用压力效应和温度效应对目标造成毁伤的炸药。近几十年,温压炸药成为炸药研究的热点之一。A.Hahma等[1]通过测量冲击波超压研究, 比较了不同金属燃料对温压炸药TNT当量的影响;Zhang Fan等[2]利用大型爆炸罐, 研究了不同气氛条件下TNT基含铝炸药爆轰后的等静压、燃烧温度和爆炸火球状态;李秀丽等[3]采用红外热成像仪研究温压炸药的爆炸温度;阚金玲等[4-5]用红外热成像仪对温压炸药和普通炸药的火球特征参数进行了测量,发现温压炸药的能量远大于普通炸药;王晓峰等[6]根据量热法原理建立了在不同气氛条件下温压炸药爆炸能量的测量方法,用于定量评价温压炸药的爆炸总能量、爆轰能和后燃烧能。综上所述,现阶段温压炸药的研究重点是某种特定环境下能量释放规律的研究,而温压炸药是一种可以利用环境中部分氧来实现炸药能量的释放的含铝炸药,不同环境对炸药能量的释放具有不同的影响。本文中,拟采用自行设计的密闭爆炸装置对温压炸药在真空状态和空气状态下爆轰后的爆炸压力和爆炸温度进行实验测试,结合爆炸后气体产物的测试结果,分析环境状态对温压炸药爆炸性能的影响。

    实验原材料:奥克托今基温压炸药,铝粉质量分数为30%,理论密度为1.96 g/cm3。样品制备:将温压炸药采用模压方式压制成带8#雷管孔的药柱,药柱直径为25 mm,药柱质量为(25.000±0.050) g。

    实验装置如图 1所示,密闭爆炸装置为一钢结构的圆柱型弹体,其高为400 mm,外径为270 mm,内径为188 mm,内容积为5.8 L。本实验装置的温度传感器采用具有自恢复能力的快速反应钨铼热电偶,布置在距离端盖中心40 mm处,下端距离上端盖底部180 mm,响应时间达10-5 s;最大可耐压力达135 MPa;测温系统频带宽度为200 kHz;放大倍数为100,温度范围为-240~1 200 ℃, 精度小于1%。本实验装置的压力传感器采用超高温硅压阻传感器,布置在距离端盖中心40 mm处,其压力范围为0~140 MPa,精度小于1%。

    图  1  实验装置示意图
    Figure  1.  Schematics of experimental device

    真空环境下的实验步骤:

    首先, 将点火装置短路,把实验用温压炸药样品悬挂在距离上端盖20 cm处,再将起爆雷管接到点火装置上;

    然后, 将实验装置上端盖密封,用真空泵抽空爆炸罐内的空气,再向爆炸罐内缓慢充入氮气,如此循环3次,将爆炸罐内的氧气完全抽走,使爆炸罐内剩余气体的压力约为3 kPa,起爆实验样品,压力传感器和温度传感器记录50 s内的电信号数据;

    最后,通过通气装置,用气体采样袋采集反应后的气体样品,利用Clarus500气相色谱仪对爆轰后的N2、CO2、CO、CH4等主要气体产物进行定量分析。

    空气环境下的实验步骤:

    首先, 将点火装置短路,把实验用温压炸药样品悬挂在距离上端盖20 cm处,再将起爆雷管接到点火装置上;

    然后, 将实验装置上端盖密封,起爆实验样品,压力传感器和温度传感器记录50 s内的电信号数据;

    最后, 通过通气装置,用气体采样袋采集反应后的气体样品,利用Clarus500气相色谱仪对爆轰后的N2、CO2、CO、CH4等主要气体产物进行定量分析。

    实验得到的温压炸药真空和空气环境下爆炸压力的电压U信号与时间t的关系曲线如图 2所示。将图 2数据处理后可得到:真空环境下的入射峰值压力和平衡压力分别为4.44和0.25 MPa, 空气环境下的入射峰值压力和平衡压力分别为8.77和0.39 MPa。由此可知,空气环境下温压炸药的爆炸入射峰值压力比真空环境下的高97%。其原因是:在真空环境下,温压炸药爆炸能量没有传播的载体,只能依靠自身反应生成的气体向外膨胀来传播,因此炸药爆炸压力衰减迅速;在空气环境下,温压炸药爆炸能量可以通过空气向外传播,因此,其爆炸压力衰减相对真空环境衰减缓慢。空气环境下的平衡压力比真空环境下的高56%,这是由于在空气环境下,空气与温压炸药爆轰后的气体产物的摩尔量总和高于真空环境下温压炸药爆轰后的气体产物的摩尔量。

    图  2  温压炸药爆炸压力随时间的变化
    Figure  2.  Explosion pressure of thermobaric explosive varying with time

    实验得到的温压炸药真空和空气环境下爆炸温度的电压信号与时间的关系曲线如图 3所示。将图 3数据处理后可得入射峰值温度和平衡温度。真空环境下,温压炸药爆炸产物入射峰值温度为943 ℃,平衡温度为283 ℃;空气环境下,温压炸药爆炸产物入射峰值温度高于1 371 ℃,平衡温度为320 ℃。温压炸药在空气环境下的爆炸场峰值温度比真空环境下的高45%以上,温压炸药在空气环境下的平衡温度比真空环境下的高13%。其原因是:真空环境下,温压炸药中的单质炸药首先爆炸生成高温高压的气体产物,高温环境下其气体产物与铝粉发生氧化还原反应并放出热量,爆炸产物温度迅速上升,与此同时高温高压的气体产物不断对外膨胀做功输出能量,使爆炸产物温度下降,由于温度的下降以及氧化反应导致的产物中氧含量的降低制约了铝粉的反应放热,致使真空环境下爆炸温度随着爆炸产物的膨胀而迅速降低,直至达到温度平衡;而在空气环境下,温压炸药中的单质炸药首先爆炸生成高温高压的气体产物,高温环境下其气体产物与铝粉发生氧化还原反应并放出热量,爆炸产物温度迅速上升,与此同时高温高压的气体产物不断对外膨胀做功输出能量,在气体产物膨胀的过程中与空气充分混合提高了氧含量,可以促进铝粉的氧化反应,提高铝粉反应的完全性及放热量,因此,在爆炸产物膨胀至温度传感器时仍能保持比真空环境下较高的温度。

    图  3  温压炸药爆炸温度与时间关系
    Figure  3.  Explosion temperature of thermobaric explosive varying with time

    不同环境条件下,温压炸药爆炸后收集的气体产物气相色谱分析结果见表 1,表中各气体的含量是其摩尔量的百分比含量。由表 1中数据可知,空气环境下温压炸药的气体产物中CO2的含量明显低于真空环境下。按照温压炸药配方组成计算,25.00 g温压炸药,氧元素的物质的量为0.425 7 mol。温压炸药在空气环境下爆炸时,系统的氧元素的物质的量为0.519 7 mol。

    表  1  不同氛围下气体产物的摩尔分数
    Table  1.  Mole fraction of gas product under different conditions
    氛围 x(CH4)/% x(CO2)/% x(N2)/% x(CO)/%
    真空 0.625 0 3.617 5 42.764 1 25.803 3
    空气 1.551 5 0.029 2 33.556 9 26.283 2
    下载: 导出CSV 
    | 显示表格

    将炸药爆炸后的气体视为理想气体,按照实验测得的平衡状态下的压力和温度,通过理想气体状态方程计算气体产物的物质的量,以及结合表 1中数据计算得表 2。由表 2中数据可知,虽然空气环境下剩余氧元素的物质的量nrem(O)比真空环境下的高0.017 3 mol,由于空气环境下氧元素的总摩尔量nal(O)比真空环境下的高0.094 0 mol,因此,空气环境下实际参加氧化反应的氧元素的物质的量nrea(O)比真空环境下的多0.076 7 mol,正是由于这部分氧气参加氧化反应放出的热量使空气环境下的平衡温度比真空环境下的高37 ℃。

    表  2  不同氛围下气体产物的摩尔量
    Table  2.  Mole of gas product under different conditions
    氛围 nal(O)/mol n(CH4)/mol n(CO2)/mol n(N2)/mol n(CO)/mol nrem(O)/mol nrea(O)/mol
    真空 0.425 7 0.002 0 0.011 3 0.134 2 0.080 9 0.103 5 0.322 2
    空气 0.519 7 0.007 1 0.000 1 0.154 0 0.120 6 0.120 8 0.398 9
    下载: 导出CSV 
    | 显示表格

    通过测量25 g温压炸药真空和空气条件下在容积为5.8 L的密闭爆炸罐内爆轰后的爆炸压力和爆炸温度以及气态产物分析,得到以下结论:(1)温压炸药在空气环境下爆轰后的平衡压力和平衡温度明显高于真空环境下的平衡压力和平衡温度; (2)空气中的氧气参与了温压炸药第3阶段铝粉有氧燃烧反应,证明温压炸药在空气中爆轰存在明显的后燃效应。

  • 图  1  计算区域示意图

    Figure  1.  Sketch of computational domain

    图  2  实验纹影[4]与计算密度的比较

    Figure  2.  Comparison between experimental schlieren images[4] and computational density images at selected times

    图  3  计算密度(Case 2)

    Figure  3.  Computational density images (Case 2)

    图  4  计算密度与质量分数(Case 3)

    Figure  4.  Images of computational density and mass fraction (Case 3)

    图  5  计算密度与质量分数(Case 5)

    Figure  5.  Images of computational density and mass fraction (Case 5)

    图  6  火焰有效面积随时间的变化

    Figure  6.  Time histories of flame effective area

    图  7  火焰平均反应放热率随时间的变化

    Figure  7.  Time histories of average reaction heat release rate

    表  1  不同激波马赫数和火焰尺寸的4组算例

    Table  1.   Four cases with different shock Mach numbers and flame sizes

    算例 Ma R0/m
    Case 1 1.7 0.019
    Case 2 2.1 0.019
    Case 3 2.5 0.019
    Case 4 1.7 0.024
    Case 5 2.1 0.024
    Case 6 2.5 0.024
    下载: 导出CSV
  • [1] Marble F E, Hendricks G J, Zukoski E E. Progress toward shock enhancement of supersonic combustion processes[C]//23rd Joint Propulsion Conference. San Diego, CA, 1987.
    [2] Oran E S, Gamezo V N. Origins of the deflagration-to-detonation transition in gas-phase combustion[J]. Combustion and Flame, 2007, 148(1/2):4-47. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=4cf0d5c02c507b4281e25a3cd95da815
    [3] Markstein G H. A shock-tube study of flame front-pressure wave interaction[J]. Symposium (International) on Combustion, 1957, 6(1):387-398. doi: 10.1016/S0082-0784(57)80054-X
    [4] Thomas G O, Bambrey R, Brown C. Experimental observations of flame acceleration and transition to detonation following shock-flame interaction[J]. Combustion Theory and Modelling, 2001, 5(4):573-594. doi: 10.1088/1364-7830/5/4/304
    [5] Ju Y, Shimano A, Inoue O. Vorticity generation and flame distortion induced by shock flame interaction[J]. Symposium (International) on Combustion, 1998, 27(1):735-741. doi: 10.1016/S0082-0784(98)80467-0
    [6] Khokhlov A M, Oran E S, Thomas G O. Numerical simulation of deflagration-to-detonation transition: The role of shock-flame interactions in turbulent flames[J]. Combustion and Flame, 1999, 117(1/2):323-339.
    [7] Khokhlov A M, Oran E S. Numerical simulation of detonation initiation in a flame brush: The role of hot spots[J]. Combustion and Flame, 1999, 119(4):400-416. doi: 10.1016/S0010-2180(99)00058-9
    [8] Gamezo V N, Oran E S, Khokhlov A M. Three-dimensional reactive shock bifurcations[J]. Proceedings of the Combustion Institute, 2005, 30(2):1841-1847. doi: 10.1016/j.proci.2004.08.259
    [9] Teng H H, Jiang Z L, Hu Z M. Detonation initiation developing from the Richtmyer-Meshkov instability[J]. Acta Mechanica Sinica, 2007, 23(4):343-349. doi: 10.1007/s10409-007-0085-2
    [10] 谷壮志, 王超, 施红辉, 等.激波诱导火焰变形的数值模拟[J].浙江理工大学学报, 2011, 28(4):529-533. doi: 10.3969/j.issn.1673-3851.2011.04.010

    Gu Zhuangzhi, Wang Chao, Shi Honghui, et al. The numerical simulation of the flame deformation induced by shock wave[J]. Journal of Zhejiang Sci-Tech University, 2011, 28(4):529-533. doi: 10.3969/j.issn.1673-3851.2011.04.010
    [11] Dong G, Fan B, Ye J. Numerical investigation of ethylene flame bubble instability induced by shock waves[J]. Shock Waves, 2008, 17(6):409-419. doi: 10.1007/s00193-008-0124-3
    [12] 朱跃进, 董刚, 范宝春.受限空间内激波与火焰作用的三维计算[J].推进技术, 2012, 33(3):405-411. http://d.old.wanfangdata.com.cn/Periodical/tjjs201203011

    Zhu Yuejin, Dong Gang, Fan Baochun. Three-dimensional computation of the interactions between shock waves and flame in a confined space[J]. Journal of Propulsion Technology, 2012, 33(3):405-411. http://d.old.wanfangdata.com.cn/Periodical/tjjs201203011
    [13] 朱跃进, 董刚, 刘怡昕, 等.激波诱导火焰变形、混合和燃烧的数值研究[J].爆炸与冲击, 2013, 33(4):430-437. doi: 10.3969/j.issn.1001-1455.2013.04.016

    Zhu Yuejin, Dong Gang, Liu Yixin, et al. A numerical study on shock induced distortion, mixing and combustion of flame[J]. Explosion and Shock Waves, 2013, 33(4):430-437. doi: 10.3969/j.issn.1001-1455.2013.04.016
    [14] Zhu Y J, Dong G, Liu Y X. Three-dimensional numerical simulations of spherical flame evolutions in shock and reshock accelerated flows[J]. Combustion Science and Technology, 2013, 185(10):1415-1440. doi: 10.1080/00102202.2013.798656
    [15] 朱跃进, 董刚.激波冲击火焰的涡量特性研究[J].爆炸与冲击, 2015, 35(6):839-845. doi: 10.11883/1001-1455(2015)06-0839-07

    Zhu Yuejin, Dong Gang. A study of vorticity characteristics of shock-flame interaction[J]. Explosion and Shock Waves, 2015, 35(6):839-845. doi: 10.11883/1001-1455(2015)06-0839-07
    [16] Jiang G S, Shu C W. Efficient implementation of weighted ENO schemes[J]. Journal of Computational Physics, 1996, 126:202-228. doi: 10.1006/jcph.1996.0130
    [17] 王革, 关奔.激波作用下R22气泡射流现象研究[J].力学学报, 2013, 45(5):707-715. http://d.old.wanfangdata.com.cn/Thesis/Y2235904

    Wang Ge, Guan Ben. A study on jet phenomenon of R22 gas cylinder under the impact of shock wave[J]. Chinese Journal of Theoretical and Applied Mechanics, 2013, 45(5):707-715. http://d.old.wanfangdata.com.cn/Thesis/Y2235904
    [18] Gordon S, Mcbride B J. Computer program for calculation of complex chemical equilibrium compositions and application Ⅰ: analysis[Z]. NASA Reference Publication 1311, 1994.
  • 加载中
图(7) / 表(1)
计量
  • 文章访问数:  4626
  • HTML全文浏览量:  1251
  • PDF下载量:  382
  • 被引次数: 0
出版历程
  • 收稿日期:  2015-12-23
  • 修回日期:  2016-05-03
  • 刊出日期:  2017-07-25

目录

/

返回文章
返回