Conditions for shock wave induced flame instability and detonation
-
摘要: 采用九阶WENO和十阶中心差分格式数值求解激波与火焰作用过程,考察了激波强度、火焰尺寸对激波与球形火焰作用过程的影响。结果表明,增大激波强度或火焰尺寸均可在流场中引发爆轰,但激波强度的影响更大,并且其引发的爆轰可使火焰迅速膨胀,放热率提高,从而影响燃烧特性;此外,爆轰波传播过程中会迅速消耗可燃预混气,合并原有的反射激波,并在流场中形成局部高压区,极大地改变流场结构。Abstract: A computational study of the interaction between shock waves and a spherical flame was carried out using the ninth-order WENO and the tenth-order central difference schemes, and the influence of shock intensity and flame size on the interaction process was investigated. It can be found from the results of our study that the increase of the shock intensity and the flame size can both induce detonation in the flow field, but the influence of the shock intensity is relatively stronger. Further, the detonation induced by shock wave can lead to quick flame expansion and increase its heat release rate, thereby affecting the combustion characteristics. Besides, the detonation wave will quickly burn out the combustible gas, merge the previously existing reflected shock waves in the propagation process, and form local high pressure zones, which can significantly alter the flow field structure.
-
Key words:
- shock wave /
- flame /
- detonation /
- flow field structure
-
温压炸药(thermobaric explosive,TBX)是一类能够充分利用压力效应和温度效应对目标造成毁伤的炸药。近几十年,温压炸药成为炸药研究的热点之一。A.Hahma等[1]通过测量冲击波超压研究, 比较了不同金属燃料对温压炸药TNT当量的影响;Zhang Fan等[2]利用大型爆炸罐, 研究了不同气氛条件下TNT基含铝炸药爆轰后的等静压、燃烧温度和爆炸火球状态;李秀丽等[3]采用红外热成像仪研究温压炸药的爆炸温度;阚金玲等[4-5]用红外热成像仪对温压炸药和普通炸药的火球特征参数进行了测量,发现温压炸药的能量远大于普通炸药;王晓峰等[6]根据量热法原理建立了在不同气氛条件下温压炸药爆炸能量的测量方法,用于定量评价温压炸药的爆炸总能量、爆轰能和后燃烧能。综上所述,现阶段温压炸药的研究重点是某种特定环境下能量释放规律的研究,而温压炸药是一种可以利用环境中部分氧来实现炸药能量的释放的含铝炸药,不同环境对炸药能量的释放具有不同的影响。本文中,拟采用自行设计的密闭爆炸装置对温压炸药在真空状态和空气状态下爆轰后的爆炸压力和爆炸温度进行实验测试,结合爆炸后气体产物的测试结果,分析环境状态对温压炸药爆炸性能的影响。
1. 实验
1.1 实验样品
实验原材料:奥克托今基温压炸药,铝粉质量分数为30%,理论密度为1.96 g/cm3。样品制备:将温压炸药采用模压方式压制成带8#雷管孔的药柱,药柱直径为25 mm,药柱质量为(25.000±0.050) g。
1.2 实验装置
实验装置如图 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.3 实验方法
真空环境下的实验步骤:
首先, 将点火装置短路,把实验用温压炸药样品悬挂在距离上端盖20 cm处,再将起爆雷管接到点火装置上;
然后, 将实验装置上端盖密封,用真空泵抽空爆炸罐内的空气,再向爆炸罐内缓慢充入氮气,如此循环3次,将爆炸罐内的氧气完全抽走,使爆炸罐内剩余气体的压力约为3 kPa,起爆实验样品,压力传感器和温度传感器记录50 s内的电信号数据;
最后,通过通气装置,用气体采样袋采集反应后的气体样品,利用Clarus500气相色谱仪对爆轰后的N2、CO2、CO、CH4等主要气体产物进行定量分析。
空气环境下的实验步骤:
首先, 将点火装置短路,把实验用温压炸药样品悬挂在距离上端盖20 cm处,再将起爆雷管接到点火装置上;
然后, 将实验装置上端盖密封,起爆实验样品,压力传感器和温度传感器记录50 s内的电信号数据;
最后, 通过通气装置,用气体采样袋采集反应后的气体样品,利用Clarus500气相色谱仪对爆轰后的N2、CO2、CO、CH4等主要气体产物进行定量分析。
2. 结果与分析
2.1 温压炸药爆炸压力
实验得到的温压炸药真空和空气环境下爆炸压力的电压U信号与时间t的关系曲线如图 2所示。将图 2数据处理后可得到:真空环境下的入射峰值压力和平衡压力分别为4.44和0.25 MPa, 空气环境下的入射峰值压力和平衡压力分别为8.77和0.39 MPa。由此可知,空气环境下温压炸药的爆炸入射峰值压力比真空环境下的高97%。其原因是:在真空环境下,温压炸药爆炸能量没有传播的载体,只能依靠自身反应生成的气体向外膨胀来传播,因此炸药爆炸压力衰减迅速;在空气环境下,温压炸药爆炸能量可以通过空气向外传播,因此,其爆炸压力衰减相对真空环境衰减缓慢。空气环境下的平衡压力比真空环境下的高56%,这是由于在空气环境下,空气与温压炸药爆轰后的气体产物的摩尔量总和高于真空环境下温压炸药爆轰后的气体产物的摩尔量。
2.2 温压炸药爆炸场温度
实验得到的温压炸药真空和空气环境下爆炸温度的电压信号与时间的关系曲线如图 3所示。将图 3数据处理后可得入射峰值温度和平衡温度。真空环境下,温压炸药爆炸产物入射峰值温度为943 ℃,平衡温度为283 ℃;空气环境下,温压炸药爆炸产物入射峰值温度高于1 371 ℃,平衡温度为320 ℃。温压炸药在空气环境下的爆炸场峰值温度比真空环境下的高45%以上,温压炸药在空气环境下的平衡温度比真空环境下的高13%。其原因是:真空环境下,温压炸药中的单质炸药首先爆炸生成高温高压的气体产物,高温环境下其气体产物与铝粉发生氧化还原反应并放出热量,爆炸产物温度迅速上升,与此同时高温高压的气体产物不断对外膨胀做功输出能量,使爆炸产物温度下降,由于温度的下降以及氧化反应导致的产物中氧含量的降低制约了铝粉的反应放热,致使真空环境下爆炸温度随着爆炸产物的膨胀而迅速降低,直至达到温度平衡;而在空气环境下,温压炸药中的单质炸药首先爆炸生成高温高压的气体产物,高温环境下其气体产物与铝粉发生氧化还原反应并放出热量,爆炸产物温度迅速上升,与此同时高温高压的气体产物不断对外膨胀做功输出能量,在气体产物膨胀的过程中与空气充分混合提高了氧含量,可以促进铝粉的氧化反应,提高铝粉反应的完全性及放热量,因此,在爆炸产物膨胀至温度传感器时仍能保持比真空环境下较高的温度。
2.3 气体产物测试分析
不同环境条件下,温压炸药爆炸后收集的气体产物气相色谱分析结果见表 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 将炸药爆炸后的气体视为理想气体,按照实验测得的平衡状态下的压力和温度,通过理想气体状态方程计算气体产物的物质的量,以及结合表 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 3. 结论
通过测量25 g温压炸药真空和空气条件下在容积为5.8 L的密闭爆炸罐内爆轰后的爆炸压力和爆炸温度以及气态产物分析,得到以下结论:(1)温压炸药在空气环境下爆轰后的平衡压力和平衡温度明显高于真空环境下的平衡压力和平衡温度; (2)空气中的氧气参与了温压炸药第3阶段铝粉有氧燃烧反应,证明温压炸药在空气中爆轰存在明显的后燃效应。
-
表 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 -
[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.010Gu 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/tjjs201203011Zhu 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.016Zhu 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-07Zhu 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/Y2235904Wang 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. -