Detonation Characteristics of RDX-Liquid Fuel-Air Three-Phase System
-
摘要: 在内径为0.2 m、高5.4 m的大型立式激波管内,同时对液体燃料和黑索金(RDX)粉末进行抛撒,采用底部直接起爆,得到了RDX-液体燃料-空气三相体系的爆压和爆速,利用烟熏技术得到了爆轰波的胞格结构、尺寸及长宽比。结果表明:90#溶剂油-空气两相体系的平均爆压为5~6 MPa,可以将其作为燃料-空气炸药的主燃料;向该燃料中添加硝酸异丙酯(IPN)可以降低其临界起爆能,添加RDX能够提高体系的爆速和爆压。此外,通过理论分析结合实验验证得出,多相爆轰的临界起爆能和爆轰胞格尺寸之间存在与气相爆轰类似的关系。Abstract: Liquid fuel and RDX dust were dispersed at the same time in a detonation tube with an inner diameter of 200 mm and length of 5.4 m, and the cloud formed was initiated by high explosive from the bottom.Detonation pressure and velocity were obtained by the pressure measurement system, and photos and detonation cell were got by smoke trace technology.The results indicate that the 90# solvent oil-air cloud detonation pressure is 5-6 MPa, which means the 90# solvent oil can be used as a main fuel of fuel-air explosive.By adding isopropyl nitrate (IPN), its critical initiation energy can be lowered, and by adding RDX the detonation pressure and velocity can be raised.The calculation and experimental verification also found a link between the critical initiation energy and detonation cell size of multi-phase detonation which is similar with the gaseous detonation.
-
1. 引言
液体燃料抛撒成雾后能和空气形成爆炸性混合物。通常该爆炸性混合物有两种截然不同的化学反应模式:一种是速度缓慢的层流燃烧模式, 火焰-未反应燃料的相对速度约为1 m/s; 另一种是伴随着高温高压的高速(约2 km/s)化学反应模式, 称之为爆轰。与气相爆轰[1-5]相比, 气-液两相或气-液-固三相爆轰更复杂, 伴随液-固燃料的抛撒、破碎等物理过程, 其反应特征时间和反应区尺度远大于气相爆轰。相比于凝聚相炸药, 多相爆轰的容积能量密度要低很多, 但其质量能量密度却远超过一般凝聚相炸药。
近年来, 关于多相爆轰的研究持续升温。一个主要原因是出于对燃料-空气炸药(Fuel-Air Explosive, FAE)的关注, 通过采用不同种类和配比的燃料, 寻找高能燃料以提高FAE的威力。另一方面, 是出于对工业生产中安全问题的考虑。煤矿中的煤尘-瓦斯-空气爆炸, 食品制药行业中的可燃粉尘-空气爆炸, 石油化工企业中的燃料-空气爆炸等多相爆轰现象都能造成重大的灾难性事故。
国内已进行了很多关于多相爆轰的研究。刘光烈[6]设计了一台多相爆轰试验管, 并利用它测量了环氧丙烷-铝粉-空气混合物的爆轰压力和正压作用时间, 测量结果与外场试验结果基本一致。刘庆明等人[7]对多相FAE的分散、爆轰过程进行了光学测量, 并测量了FAE爆炸压力场, 分析了气-液-固多相爆轰的特征和压力波形的特点, 研究了冲击波峰值超压及比冲量随传播距离变化的规律。汤明钧[8]归纳总结了多相爆轰的特点, 并对影响多相爆轰的因素进行了分析。贵大勇等人[9]研究了典型液态燃料和固-液混合燃料配方对FAE爆轰性能的影响规律, 并对各典型液态燃料的FAE爆轰威力进行了排序。徐晓峰等人[10]以立式激波管为主要实验手段, 对4种燃料(C5H8.68、环氧丙烷、正己烷和癸烷)与空气混合物的云雾爆轰性能进行了较为系统的研究。姚干兵等人[11]采用烟迹技术, 在立式激波管中测定了环氧丙烷、90号汽油、硝酸异丙酯、庚烷、癸烷、戊二烯等几种燃料气-液两相云雾爆轰的胞格尺寸, 认为液滴的碎解、汽化过程以及燃烧区前导是控制气-液两相云雾爆轰的主要因素。陈默等人[12]在长为32.4 m、内径为0.199 m的大型长直水平管道中, 对环氧丙烷-铝粉-空气三相流云雾的爆燃转爆轰(Deflagration to Detonation, DDT)过程进行了实验研究, 实现了多相体系爆燃向爆轰转变的过程, 并测得其胞格尺寸。
本研究以立式激波管为载体, 以RDX-混合燃料-空气三相体系为对象, 借助压力测试系统、烟熏技术, 进行爆轰波压力、速度和胞格尺寸的测量, 目的是获得气-液-固三相爆轰的相关参数, 以期了解多相爆轰波的结构及其传播规律。
2. 实验装置及测试系统
2.1 实验设备
整套立式激波管设备由激波管管体、喷粉喷雾系统和点火系统组成。其结构如图 1所示。激波管内径为200 mm, 高5.4 m, 管体分为3个部分:主体部分, 长4 m; 底端为起爆端, 长0.7 m; 顶端为观察窗, 长0.7 m。不使用连着观察窗的4个喷头, 以防止实验中观察窗遭到破坏。故仅利用4.7 m的管体, 有效容积为147.6 L。喷粉喷雾系统由空气压缩机、储气罐、U型管储液储粉装置、喷头等组成。点火系统由延时点火器、起爆线、雷管基座等组成。
实验开始前, 先将烟熏板置于激波管内指定位置(见图 1), 安装好起爆装置; 然后向U型管内注入燃料和粉尘, 开启空气压缩机, 向储气罐内注入一定量的压缩空气; 最后开启压力采集系统, 设置点火延时器。准备就绪后, 打开电磁阀使压缩空气冲出, 夹带着液体和粉尘的空气由喷头抛撒进入激波管, 至云团达到最佳状态点火, 记录压力曲线并取下烟熏板保存。点火延迟时间是指从开启电磁阀到点火的这段时间。根据多次重复实验, 液体-空气云团在电磁阀开启1 s后发展到最大, 故将点火延迟时间取为1 s。充入储气罐中的气压值经过多次尝试, 取0.4 MPa时, 能在500 ms内将U型管内的液体、粉尘喷撒干净[10]。为确保将U型管内粉尘吹干净, 设定电磁阀开启持续时间为1 s。电磁阀关闭后即刻开始点火。即吹粉时间为1 s, 点火延迟时间也为1 s。
本研究采用了两种液体燃料:一种是90#溶剂油, 另一种是90#溶剂油和硝酸异丙酯(IPN)按质量比4:1制成的混合物。液体燃料和RDX粉尘交叉注入26个U型管中, 即:若液体燃料注入某一个U型管, 则RDX粉尘注入其对面和临近的U型管中。除特别说明外, 起爆物均为一发雷管加3 g塑性炸药, 起爆能量为23.52 kJ。
2.2 测试系统
测试系统由压力测试系统、胞格记录装置等组成。其中, 压力测试系统由传感器、电荷放大器、数据采集卡、微机等组成。采用PCB压电式石英传感器, 编号为S1~S8, 分布如图 1所示。每两个传感器间距为0.5 m, 传感器S2距离起爆端1.4 m。
胞格记录装置主要由烟熏板、固定支架组成。其中, 烟熏板为半圆形铝板, 使用前将板清洗干净并晾干, 表面涂一层很薄的机油, 用煤油灯烟熏黑其内表面, 出现一层薄而均匀的烟灰, 用0.4 MPa压缩空气吹之, 烟灰不脱落者为合格。烟熏板用支架固定在法兰上, 板中心距离上法兰1.5 m。
3. 实验结果与讨论
表 1给出了各传感器的具体位置和参数。经过多次测试发现, 距离起爆点较近处压力大, 较远处压力稍小。故将量程大的传感器安置在激波管底端附近, 量程小的传感器安置在管体顶端附近。图 2为典型的爆轰波压力-时间曲线, 传感器同时记录了各测试点处的爆压和爆轰波到达该点的时刻, 以两个传感器之间的距离除以爆轰波到达的时间差, 即可推算出爆轰波在这两个传感器间的平均速度。
表 1 传感器位置和量程Table 1. Position and range of pressure sensorsNo. Distance from the
detonation point/
(m)Maximum
pressure/
(MPa)S2 1.4 15.330 S3 1.9 15.270 S4 2.4 6.784 S5 2.9 6.757 S6 3.4 6.220 S7 3.9 6.000 实验所用RDX为白品造粒, 经BT-9300S粒度分析仪测得其平均粒径(D50)为85 μm。90#溶剂油和IPN的理化性质如表 2所示, 二者按质量比4:1制成的混合物密度为0.736 g/cm3。如前所述, 液体燃料和RDX粉尘是交叉注入26个U型管的, 即13个装液体, 13个装粉尘。经多次重复实验得到, 当单个U型管中加入4 mL、即总共加入52 mL燃料时, 刚好发生爆轰, 故纯液体燃料的点火浓度下限为260 g/m3。在此基础上, 改变RDX粉尘的加入量, 每次向单个U型管中分别加入1、2、3、4、5 g, 对应粉尘浓度依次为88、176、264、352、440 g/m3。在激波管顶部加观察窗, 记录云雾形成、爆轰过程, 得到如图 3所示的时序照片。
表 2 90#溶剂油和IPN的理化性质Table 2. Physicochemical properties of 90# solvent oil and IPNLiquid
fuelDensity/
(g/cm3)Composition Boiling
point/
(℃)Flash
point/
(℃)Ignition
temperature/
(℃)Upper
explosive
limit/(%)Lower
explosive
limit/(%)90# solvents oil 0.65 C5H12, C6H14 35-60 13 228 5.9 1.1 IPN 1.04 C3H7NO3 98-102 11 - 100 2 3.1 爆轰压力
不同配比的RDX-90#溶剂油及RDX-混合物的爆轰压力(简称爆压)随距起爆点距离的变化趋势如图 4、图 5所示, 图中:+n g RDX(n=0, 2, 5)意为每次向单个U型管中加入n g RDX, 下同。
多相爆轰为不稳定爆轰, 爆轰波的压力、速度并非一成不变。由图 5、图 6可知, 同一组实验中, 随着爆轰波的传播, 爆轰压力的波动较大, 整体呈下降趋势。分析其原因, 可能是因为液滴、粉尘喷出之后, 较大颗粒的液滴、粉尘部分沉降, 出现激波管底部粉尘浓度较大、顶部粉尘浓度较小的情况, 所以压力呈现递减趋势。但不论是纯溶剂油还是混合溶剂体系, 其气云爆轰的最低压力均在5 MPa左右, 比一般碳氢液体燃料-空气的气云爆轰压力高出很多。这主要是因为90#溶剂油的挥发性较强, 且含碳量较高, 能量密度较大, 与空气形成爆炸性混合物后, 反应较为充分。
由图 4、图 5得知, 随着RDX粉尘浓度的增加, RDX-90#溶剂油和RDX-混合溶剂三相体系的平均爆压均呈整体上升趋势, 相比而言, 前者爆压的上升趋势不如后者明显。
3.2 爆轰速度
不同配比的RDX-90#溶剂油、RDX-混合溶剂的爆轰速度(简称爆速)随距起爆点距离的变化趋势如图 6、图 7所示。
距离起爆点较近(1.5~2.5 m)处, 爆速起伏上升; 距离大于2.5 m以后, 爆速趋于稳定, 并达到最大值。从图 6、图 7中可以得知, 随着RDX粉尘浓度的增加, RDX-90#溶剂油和RDX-混合溶剂的爆速均呈现整体上升的趋势。
图 8为RDX-燃料-空气三相体系在距离起爆点3.15 m处的平均爆速随RDX质量分数的变化趋势。无论何种配比, RDX-90#溶剂油体系的平均爆速均不小于RDX-混合溶剂。这一点也不难理解, 由表 2可知, 90#溶剂油的密度较小, 沸点较低, 饱和蒸汽压较高, 容易挥发汽化, 在激波阵面到达以前, 体系中就有一定的燃料蒸汽存在, 故而其破碎雾化的时间较短, 爆轰压力和爆轰速度也较大。
3.3 RDX参与了爆轰反应的验证
由于90#溶剂油和IPN混合物原本的爆速、爆压就比较高, 故对RDX是否参与了爆轰有所怀疑, 设计如下实验进行验证。其它条件(如:加入液体和粉尘的量、电磁阀开启持续时间、起爆能量等)不变, 将液体燃料换成环氧丙烷(Propylene Oxide, PO), 且注入U型管后不进行静置, 不待其发生汽化, 直接进行抛撒起爆, 进行3次重复实验。传感器S3~S7的参数和位置见表 1。
单纯注入PO时, 压力-时间曲线如图 9所示。由图 9可以看出, 波形出现双峰结构, 说明云团内发生了二次反应, 且二次反应比第一次猛烈。第一和第二个峰之间正是液滴发生破碎雾化的过程, 特征时间Δt=0.2 ms。从爆压和爆速(1 km/s以下)的数据来看, PO-空气体系并没有发生爆轰。
加入RDX后, 压力-时间曲线如图 10所示。由图 10可以看出, 爆压达到5~6 MPa, 爆速达到1.5 km/s以上。加入RDX粉尘之后, PO由爆燃达到了爆轰, 压力和速度都有很大突跃, 可见RDX确实参与了爆轰, 而不仅仅是燃烧。
3.4 多相爆轰胞格的测量与爆压分析
3.4.1 胞格尺寸和长宽比
所谓爆轰胞格就是爆轰波在管道内传播过程中, 横波、入射波、马赫杆三波交合点的运动轨迹。爆轰的胞格尺寸是判断可燃混合气体的最基本特征参数, 爆轰参数诸如临界管径、爆轰极限以及直接气爆能量等, 都可以通过胞格尺寸进行预测[13]。如图 11所示, 爆轰波从左至右传播, 其前导冲击波是由马赫杆和入射冲击波交替组合的波阵面。自1959年Denisov和Troshin发现爆轰波具有胞格结构以来, 实验中都使用烟迹技术观测胞格结构。本研究所用的烟迹载体为半圆形铝制板, 利用支架将其固定在激波管的内壁上, 以最大限度减少对爆轰波传播过程的干扰。
图 12是用上述铝制板得到的RDX-燃料-空气三相爆轰的典型胞格, 其中L为胞格长度, λ为胞格宽度。爆轰波从左至右传播, 图 12中用红色线标记了几个典型的胞格。与气相爆轰的胞格(见图 11)相比, 三相爆轰胞格的长宽比(L/λ)要小得多, 有些胞格长、宽相等甚至长比宽还要小。
多相爆轰为非均相爆轰, 胞格尺寸、形状不统一, 波动性较大(见图 12)。经过重复实验, 取胞格尺寸(此处特指胞格的宽度)平均值, 并利用数据处理软件计算其误差线, 如图 13所示。
图 13分别为RDX-90#溶剂油-空气体系及RDX-混合溶剂-空气体系的爆轰胞格尺寸随RDX质量分数的变化趋势。两种体系相比, 单纯90#溶剂油体系的平均胞格尺寸比加入IPN的混合溶剂体系要大, 也即前者的临界起爆能量较高。这一点很好解释, 因为IPN中带有-O-NO2基团, 相当于“敏化剂”, 降低了体系的直接起爆能量。另外, 两种体系的爆轰胞格尺寸都在RDX质量分数为55%左右时达到最大, 也即当RDX粉尘质量与液体燃料质量之比为55:45时, 体系最不敏感, 需要较大的能量才能起爆。在爆轰过程中, 液滴、粉尘的破碎相变过程需要吸收较多能量。RDX粉尘的比表面积较大, 增加的RDX粉尘自身蕴含的能量不足以抵消加热它们所需的能量, 所以随着粉尘量的增加, 起爆所需要的能量也增加。但是, RDX的能量密度远大于液体燃料, 当RDX质量分数达到一定程度(约55%)时, 体系增加的能量开始大于需要多消耗的能量, 故此后, 起爆所需的能量就迅速减少。如图 13所示, 当RDX质量分数大于55%时, 体系的爆轰胞格尺寸迅速减小。
在测量胞格宽度的同时, 也测定了胞格的长度, 并计算其长宽比, 用数据处理软件计算得到误差线, 如图 14所示。Strehlow等人发现, 低爆轰压力下, 气相爆轰波的胞格长宽比近似满足L/λ≈1.61。不过正如图 12和图 14所示, 多相爆轰波的胞格长宽比较小, 仅为1.2左右。胞格长宽比反映的是爆轰波和横波平均传播速度的比值, 其中:爆轰波的平均传播速度为定值, 由总放热量决定; 而横波的平均传播速度与波后声速相关[16]。如图 8所示, 90#溶剂油体系的爆轰波平均传播速度比混合溶剂体系大; 又由表 3知, 二者介质密度接近, 故其波后声速也接近, 因此前者的胞格长宽比大于后者。
在二次引爆型FAE的云雾起爆阶段, 若发生蹿火将大大降低FAE的作用威力。如果在液固云爆药剂中采用前述配方(RDX的质量分数为55%), 则FAE云团的直接气爆能量增加, 能够降低发生蹿火的风险。不过其效果还需要进一步的实验验证。
3.4.2 胞格尺寸与临界起爆能
根据Lee提出的表面积能量模型[13-15], 即球形爆轰波表面所含的能量等于平面爆轰在临界管径内的能量, 此时所需满足的临界条件为
4πR∗2=πd2c4 (1) 式中:R*为爆轰内核半径, dc为临界管径。
此外, 由强爆轰波理论可知, 爆轰波能量为
Ec=4πIγp0Ma2R3Sh (2) 式中:Ec为临界起爆能量, γ为气体绝热指数, I为定常数(当γ=1.4时, I=0.423), p0为压强, Ma为球形爆轰波的马赫数, RSh是球形爆轰波的半径。
依据Zeldovich准则, 当爆炸波衰减到Ma=MaCJ(MaCJ为CJ爆轰波的马赫数)时, 有RSh=R*。对于未稀释的混合气体, 临界管径与爆轰胞格尺寸之间存在如下经验公式
dc=13λ (13) 因此, 联合(1)式~(3)式, 可得气相爆轰参数中, 临界起爆能量与胞格尺寸之间的关系为
Ec=4πIγp0Ma2CJ(13λ4)3=2197πρ0v2CJIλ316 (4) 式中:ρ0为混合物密度; vCJ为CJ爆速, 近似取实测爆速。根据(4)式, 即可由气相爆轰胞格尺寸预测其直接起爆的临界起爆能量。
多相爆轰中, 无论液体还是固体粉尘, 都要先经历物理破碎、相变过程, 之后的阶段就和气相爆轰一样, 所以推测多相爆轰也能采用气相爆轰的方法预测其临界起爆能量。为了验证(4)式是否能应用于多相爆轰中, 测试了RDX-90#溶剂油-空气和RDX-混合溶剂-空气体系的临界起爆能量, 起爆物分别为一发雷管加1.5 g塑性炸药和一发雷管。测试结果与理论值相比较, 列于表 3。
表 3 临界起爆能实验值与理论值的比较Table 3. Comparison between experimental and theoretical value of critical initiation energyConstituents ρ0/
(kg/m3)vCJ/
(m/s)λ/
(mm)Theoretical
Ec/(kJ)Experimental
Ec/(kJ)RDX-90# solvent oil-air 1.633 2 1 794 24.6 14.279 14.735 RDX-mixed fuel-air 1.659 5 1 808 19.1 6.893 5.945 由此可见, 多相爆轰参数——临界起爆能和胞格尺寸之间确实存在类似(4)式的关系, 但还需要进一步的实验和理论验证。
4. 结论
(1) 立式激波管可以对液、固两相燃料进行同时抛撒, 实现气-液-固三相爆轰, 并能对爆压、爆速、最小起爆能、胞格尺寸等爆轰参数进行测量。
(2) 90#溶剂油-空气体系的平均爆压达到5~6 MPa, 可以将其作为燃料-空气炸药的主燃料。
(3) 添加RDX有助于提高90#溶剂油-IPN混合燃料体系的爆速和爆压。
(4) 加入IPN后, 混合溶剂的直接起爆能量降低, 较易实现爆轰, 故IPN可以作为一种敏化剂加入FAE燃料之中。此外, 加入IPN的混合溶剂爆轰胞格也较小, 可见临界起爆能量和爆轰胞格之间存在与气相爆轰类似的关系, 但还需要进一步的实验和理论验证。
-
表 1 传感器位置和量程
Table 1. Position and range of pressure sensors
No. Distance from the
detonation point/
(m)Maximum
pressure/
(MPa)S2 1.4 15.330 S3 1.9 15.270 S4 2.4 6.784 S5 2.9 6.757 S6 3.4 6.220 S7 3.9 6.000 表 2 90#溶剂油和IPN的理化性质
Table 2. Physicochemical properties of 90# solvent oil and IPN
Liquid
fuelDensity/
(g/cm3)Composition Boiling
point/
(℃)Flash
point/
(℃)Ignition
temperature/
(℃)Upper
explosive
limit/(%)Lower
explosive
limit/(%)90# solvents oil 0.65 C5H12, C6H14 35-60 13 228 5.9 1.1 IPN 1.04 C3H7NO3 98-102 11 - 100 2 表 3 临界起爆能实验值与理论值的比较
Table 3. Comparison between experimental and theoretical value of critical initiation energy
Constituents ρ0/
(kg/m3)vCJ/
(m/s)λ/
(mm)Theoretical
Ec/(kJ)Experimental
Ec/(kJ)RDX-90# solvent oil-air 1.633 2 1 794 24.6 14.279 14.735 RDX-mixed fuel-air 1.659 5 1 808 19.1 6.893 5.945 -
[1] 姜宗林, 滕宏辉, 刘云峰.气相爆轰物理的若干研究进展[J].力学进展, 2012, 42(2): 129-140. http://d.wanfangdata.com.cn/Periodical/lxjz201202002Jiang Z L, Teng H H, Liu Y F. Some research progress on gaseous detonation physics[J]. Advances in Mechanics, 2012, 42(2): 129-140. (in Chinese) http://d.wanfangdata.com.cn/Periodical/lxjz201202002 [2] Lee J H S. Dynamic parameters of gaseous detonations[J]. Annu Rev Fluid Mech, 1984, 16: 311-336. doi: 10.1146/annurev.fl.16.010184.001523 [3] Eaton R, Zhang B, Bergthorson J M, et al. Measurement and chemical kinetic model predictions of detonation cell size in methanol-oxygen mixtures[J]. Shock Waves, 2012, 22: 173-178. doi: 10.1007/s00193-012-0359-x [4] 王昌建, 徐胜利, 郭长铭.气相爆轰波在半圆形弯管中传播现象的实验研究[J].爆炸与冲击, 2003, 23(5): 448-453. http://www.cqvip.com/Main/Detail.aspx?id=8384401Wang C J, Xu S L, Guo C M. Experimental investigation on gaseous detonation propagation in a semi-circle bend tube[J]. Explosion and Shock Waves, 2003, 23(5): 448-453. (in Chinese) http://www.cqvip.com/Main/Detail.aspx?id=8384401 [5] 王昌建, 郭长铭, 徐胜利, 等.气相爆轰在T形管中传播现象的实验研究[J].力学学报, 2004, 36(1): 16-23. http://www.cqvip.com/Main/Detail.aspx?id=9081118Wang C J, Guo C M, Xu S L, et al. Experimental investigation on gaseous detonation propagation through a T-shape bifurcated tube[J]. Acta Mechanica Sinica, 2004, 36(1): 16-23. (in Chinese) http://www.cqvip.com/Main/Detail.aspx?id=9081118 [6] 刘光烈.多相爆轰试验管结构设计及试验验证[J].含能材料, 1997, 5(1): 36-38.Liu G L. Design of multiphase detonation tube and its experimental verification[J]. Energetic Materials, 1997, 5(1): 35-38. (in Chinese) [7] 刘庆明, 白春华, 李建平.多相燃料空气炸药爆炸压力场研究[J].实验力学, 2008, 23(4): 360-370.Liu Q M, Bai C H, Li J P. Study on blast field characteristics of multiphase fuel air explosive[J]. Journal of Experimental Mechanics, 2008, 23(4): 360-370. (in Chinese) [8] 汤明钧.多相爆轰与生产安全[J].力学与实践, 1985, 7(5): 62-63. http://www.cnki.com.cn/Article/CJFDTotal-LXYS198505017.htmTang M J. Multiphase detonation and production safety[J]. Mechanics in Engineering, 1985, 7(5): 62-63. (in Chinese) http://www.cnki.com.cn/Article/CJFDTotal-LXYS198505017.htm [9] 贵大勇, 刘吉平, 冯顺山.几种典型燃料空气炸药威力性能研究[J].含能材料, 2002, 10(3): 121-124. http://www.cqvip.com/QK/90247X/20023/6819445.htmlGui D Y, Liu J P, Feng S S. Research of power performace of several typical fuel-air explosives[J]. Energetic Materials, 2002, 10(3): 121-124. (in Chinese) http://www.cqvip.com/QK/90247X/20023/6819445.html [10] 徐晓峰.碳氢燃料爆轰特性的研究[D].南京: 南京理工大学, 2003.Xu X F. Detonation properties of several hydrocarbon fuels[D]. Nanjing: Nanjing University of Science and Technology, 2003. (in Chinese) [11] 姚干兵, 解立峰, 刘家骢.立式激波管内云雾爆轰胞格尺寸的测定与分析[J].爆炸与冲击, 2007, 27(4): 312-318.Yao G B, Xie L F, Liu J C. Measurement and analysis of the cell size caused by fuel-air detonation in a vertical shock tube[J]. Explosion and Shock Waves, 2007, 27(4): 312-318. (in Chinese) [12] 陈默, 白春华, 刘庆明.大型水平管道中环氧丙烷-铝粉-空气混和物爆燃转爆轰的实验研究[J].高压物理学报, 2011, 25(4): 359-364. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gywlxb201104013Chen M, Bai C H, Liu Q M. Experimental research of deflagration to detonation transition of epoxypropane-aluminum-air mixtures in large-scale tube[J]. Chinese Journal of High Pressure Physics, 2011, 25(4): 359-364. (in Chinese) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gywlxb201104013 [13] 张博, 白春华.气相爆轰动力学[M].北京: 科学出版社, 2012.Zhang B, Bai C H. Gaseous Detonation Mechanics[M]. Beijing: Science Press, 2012. (in Chinese) [14] Lee J H S. Initiation of gaseous detonation[J]. Annu Rev Phys Chem, 1977, 28(1): 75-104. [15] Lee J H S, Knystautas R, Guirao C. The link between cell size, critical tube diameter, initiation energy and detonability limits[C]//Fuel-Air Explosions. Waterloo: University of Waterloo Press, 1982: 157-187. [16] 刘世杰, 林志勇, 孙明波, 等.采用不同化学反应源项处理方法的胞格爆轰数值研究[J].国防科技大学学报, 2010, 32(5): 1-6.Liu S J, Lin Z Y, Sun M B, et al. Numerical Simulation of cellular detonation using different chemical reaction source term methods[J]. Journal of National University of Defense Technology, 2010, 32(5): 1-6. (in Chinese) -