Numerical simulation of methane-air explosion in a connected device with volume fraction gradient
-
摘要: 为了研究具有体积分数梯度的连通装置内甲烷-空气爆炸特性，以60 L圆柱体容器和20 L圆柱体容器通过3 m长，截面为0.035 m×0.035 m的方形管道而连接形成的容器管道连通装置作为研究对象，利用Fluidyn软件对均一体积分数的连通装置以及具有体积分数梯度的连通装置中甲烷-空气爆炸的特性进行了数值模拟。结果表明：连通装置中甲烷的均一体积分数为6.517%～8.067%时，并由大容器中心点火工况时，最大爆炸压力、最大爆炸压力上升速率、最高温度和最大速度，以及这些爆炸参数达到最大值时的时刻值随体积分数的变化约呈线性关系；连通装置大容器甲烷体积分数6.0%体积分数梯度为2.0%～8.0%且大容器中心点火时，最大爆炸压力、最大爆炸压力上升速率、最高温度和最大速度随体积分数梯度总体上呈现先增大后减小趋势；大容器中心点火时，最大爆炸压力位于小容器，最大压力上升速率位于管道1或管道2，最大速度位于管道3，速度值可达400～600m/s。本研究可为连通装置内可燃气体爆炸事故防控提供理论指导。Abstract: A connected vessel is a common typical chemical plant, and its explosion hazard is much higher than that of an independent vessel. In an actual explosion accident, the combustible gas volume fraction in the connected device presents a non-uniform state, and there is a volume fraction gradient. A connected device was chosen as the research object. The device was formed by connecting two cylindrical vessels with the volumes of 60 litres and 20 litres, respectively, through a square pipe as long as 3 meters, with a cross section of 35 mm×35 mm. To explore the methane-air explosion characteristics in the connected device with combustible gas volume fraction gradient, the Fluidyn software was applied to numerically simulate the methane-air explosions in the connected devices with uniform and non-uniform combustible gas volume fractions, respectively. The results show as follows. When the volume fraction of the methane in the connected device is uniform and ranges from 6.517% to 8.067% and the ignition is located in the center of the large vessel, the maximum explosion pressure, the maximum explosion pressure rise rate, the maximum temperature and the maximum velocity as well as their arrival times change linearily with the volume fraction of the methane. When the volume fraction of the methane in the large vessel of the connected device is 6.0%, the volume fraction gradient of the methane is 2.0% to 8.0%, and the ignition is located in the center of the large vessel, the maximum values of the parameters, including explosion pressure, explosion pressure rise rate, flame temperature and velocity, increase firstly and then decrease with increasing volume fraction gradient. When the ignition is located in the center of the large vessel, the maximum explosion pressure is in the small vessel, the maximum pressure rising rate is in the pipe connected to the large vessel, and the maximum flame velocity is in the pipe connected to the small vessel, and the flame velocity can reach 400-600 m/s. The research results can provide a theoretical guidance for preventing and controling combustible gas explosion accident in connected devices.
-
Key words:
- connected device /
- volume fraction gradient /
- methane explosion /
- pressure rise rate
-
随着水下精确制导武器的发展，舰船遭受近场水下爆炸攻击的概率越来越大[1-2]。近场水下爆炸能对舰船等海上作战平台造成毁灭性的打击[3-5]，因此，研究近场水下爆炸对现代海洋战争有着重要的战略意义。近场水下爆炸载荷主要可分为冲击波阶段及气泡射流阶段[6-9]，20世纪80年代以前，水下爆炸相关研究主要集中于冲击波对结构造成的破坏。从20世纪80年代开始，人们逐渐意识到气泡脉动对结构的损伤可能不亚于冲击波，进行了大量的水下爆炸实验，得到了广泛使用的水下爆炸经验公式[10]。一般认为，气泡脉动带来的滞后流及气泡脉动压力能造成舰船的总体毁伤[11]，气泡收缩的最后阶段，气泡坍塌继而形成的高速水射流[12]则能严重毁伤舰船局部。近年来，水下爆炸中的水射流现象及其载荷特性已成为关注的热点。
在开展水下爆炸气泡与边界相互作用的研究时，因炸药爆炸具有一定的危险性且实验成本较高，多采用电火花[13-16]或激光[17-19]生成气泡的方法。但电火花生成的气泡内部为可冷凝的水蒸气，而炸药水下爆炸生成的气泡中心为高温高压气体、外部为水蒸气，二者存在物理本质上的差异，所以有必要开展炸药水下爆炸实验。Klaseboer等[20]通过在水池中进行水下爆炸系列实验，拍摄了完整的气泡坍塌形成水射流的过程，并基于边界元方法模拟了爆炸气泡与水下弹性结构相互作用的过程。牟金磊等[21]在水箱内开展了水下爆炸系列实验，测得气泡射流阶段的压力时程曲线呈马鞍形，并对Blake判断准则进行了修正。Gan等[22]在水箱内开展了简化的细长浮式结构在水下爆炸作用下的实验和数值模拟研究，并采用CEL方法对第1次气泡脉动期间气泡与水面结构作用的全过程开展了数值模拟。Cui等[23]在水箱内开展了小当量炸药在不同边界条件下的水下爆炸气泡特性实验研究，得到了不同边界条件下气泡的形态及压力特性。然而，目前的研究主要集中在近场水下爆炸气泡水射流在不同边界条件下的演化过程和压力特性，对水射流速度和时间特性的研究不够深入。
舰船遭遇近场水下爆炸攻击时，爆心的位置多为舰底或舷侧。本文中，为了研究底部爆炸时气泡水射流阶段的载荷特性，首先开展TNT当量2.5 g的炸药在固支方板底部不同爆距下起爆的水下爆炸系列实验，通过高速摄像拍摄不同爆距下气泡与钢板相互作用形成水射流的图像。实验现象表明，固支方板底部水下爆炸形成的射流可分为吸附式射流和非吸附式射流。接着，基于ABAQUS软件，采用CEL方法开展了数值模拟研究，先验证数值模拟的准确性，再通过系列数值模拟探究吸附式射流向非吸附式射流转化的临界点，同时给出临界点附近爆距下的气泡演化过程。最后，以速度云图形式给出不同爆距下水射流形成的过程，并讨论爆距对水射流开始形成、达到最大速度和击中固支方板时速度特性和时间特性的影响。可为利用气泡射流现象对舰船等典型海上目标实施毁伤打击提供参考。
1. 固支方板底部水下爆炸实验
1.1 实验方案
固支方板底部水下爆炸实验在江苏永丰机械有限责任公司的2 m×2 m×2.2 m水箱中进行。图1为实验系统，主要由敞口水箱、高速摄像机、水下压力测试系统组成。水下压力测试系统由NI公司的PXle-4480数据采集卡、PXle-1082机箱、PXle-8880控制器以及PCB公司的138A06水下压力传感器组成，测试系统的采样率为2×106 s−1。高速摄像机帧频为3 200 s−1。炸药为TNT，药量为2.5 g。图2为固支方板实物，固支钢板尺寸为70 cm×70 cm×1 cm，材料为Q235钢，钢板的四周焊上挡水板，挡水板尺寸为70 cm×10 cm×1 cm，钢板的四角焊上支架，支架材料为直径30 mm的钢棒，底座为200 mm×200 mm×8 mm的钢板，可使钢板在冲击波作用下保持稳定。在距炸药70.0 cm、距水面20.0 cm处布置水下压力传感器，并在传感器正下方悬挂配重为1.5 kg的秤砣。炸药由钢丝固定在钢板中心的正下方，水面与水箱底部的距离为1.2 m，没过钢板。
为探究炸药到靶板的距离对气泡射流现象的影响，根据汪斌等[24]的实验，本文中设计炸药到钢板的距离d分别为15.0、17.5、20.0和28.0 cm，定义距离参数γ，并有经验公式[25]：
γ = d/Rm (1) Rm=3.5(Wh+10)1/133 (2) 式中：Rm为气泡最大半径的理论值（m），W为炸药质量（kg），h为炸药深度（m）。
本文中对应的距离参数γ分别为0.684、0.798、0.913和1.282。
1.2 实验结果
1.2.1 水下传感器测试结果
在距爆心70.0 cm、距水面20.0 cm处布置了水下压力传感器，图3为水下传感器记录的实验1压力曲线，其中冲击波压力峰值为7.64 MPa，气泡脉动压力峰值出现在38.39 ms，其值为1.77 MPa。对自由场中水下爆炸冲击波的传递，Cole[10]通过对大量实验数据的整理，总结了水下爆炸冲击波的经验公式，并得到了广泛认可，Zamyshlyayev等[26]在其基础上改进为：
pm1={44.1(W1/3Rc)1.56≤Rc/r≤1252.4(W1/3Rc)1.13 12＜Rc/r≤240 (3) 式中：pm1为冲击波超压（MPa），Rc为测点到爆心的距离（m），r为药包半径（m）。王高辉等[27]指出，自由边界和平板边界对冲击波传播特性的影响主要集中于入射波峰值之后，对入射波峰值的影响几乎可以忽略不计。因此，可将本文中测得的冲击波峰值压力与经验公式计算结果进行比较。表1为不同距离参数时水下压力传感器的冲击波峰值压力与2.5 g TNT 水下爆炸经验公式计算结果的比较。可以看出，实验冲击波压力峰值与经验公式计算结果的偏差δp在10% 以内，并且相同距离参数时偏差之间不超过3%。
表 1 冲击波峰值压力的实验结果和经验公式结果Table 1. Experimental results and empirical formula results of shock wave peak pressures实验 γ pm1/MPa δp/% 实验 经验公式 1 0.684 7.64 8.19 6.72 2 0.684 7.49 8.19 8.55 3 0.798 7.72 8.20 5.85 4 0.798 7.79 8.20 5.00 5 0.913 7.81 8.21 4.87 6 0.913 7.77 8.21 5.36 7 0.913 7.91 8.21 3.65 8 1.282 7.54 8.15 7.49 9 1.282 7.61 8.15 6.63 10 1.282 7.58 8.15 6.99 1.2.2 高速摄像测试结果
图4为距离参数γ=0.684时气泡演变过程的高速摄影图像。在0～3.75 ms时，气泡受到的Bjerknes力较小，其形状为球形。在3.75～19.06 ms时，因钢板阻碍了气泡的发展，气泡上半部分明显呈扁平状，下半部分仍保持球状。在气泡收缩阶段，气泡在Bjerknes力及浮力作用下向上迁移，气泡与钢板发生接触，吸附在钢板上，在25.63 ms时能明显看到吸附在钢板上的气泡在钢板与水面间形成盘状区域。接着，气泡开始进一步吸附在钢板上，在34.36 ms时，已吸附在钢板上的气泡呈漏斗形，而未吸附在钢板上的部分仍保持为球形。在38.75 ms时，气泡完全吸附在钢板上，整体呈漏斗形。由于漏斗形气泡底部曲率半径最小，最容易发生失稳，导致气泡底部发生凹陷并产生向上的射流。在44.69 ms时，气泡在钢板附近溃散形成气泡群，且在气泡群下方能看到，爆炸残留物在水流的作用下呈放射状指向钢板。
图 4 当γ=0.684时气泡演变过程的实验图像Figure 4. Experimental images of bubble evolution process when γ=0.684图5为γ=0.798时气泡与钢板作用演化过程，其现象基本与γ=0.684时的一致。
图 5 当γ=0.798时气泡演变过程的实验图像Figure 5. Experimental images of bubble evolution process when γ=0.798图6为γ=0.913时气泡与钢板作用演化过程。在0～9.38 ms时，气泡近似自由膨胀，其形状为球形。在12.19 ms时，气泡上表面在钢板Bjerknes力排斥作用下，上半部分呈扁平状。在19.06 ms时，气泡膨胀过程结束，在浮力及收缩过程中钢板Bjerknes力吸引作用下，气泡上表面与钢板接触，形成盘状区域，见28.13 ms时图像。与γ=0.684, 0.798时不同，当γ=0.913时，气泡在收缩阶段并未完全吸附在钢板上，而是与钢板间存在一层水层，见γ=0.798在34.36 ms时图像。在收缩过程中，气泡上表面位置保持不变，而下表面向上表面靠拢，在39.69 ms时能看到气泡半径收缩至最小，形成指向钢板的射流。在43.13 ms时，气泡发生溃散，形成小气泡群。
图 6 当γ=0.913时气泡演变过程的实验图像Figure 6. Experimental images of bubble evolution process when γ=0.913图7为γ=1.282时气泡与钢板作用演化过程。在0～9.38 ms时，气泡形态受Bjerknes力影响较弱，其形状为球形。在7.81～19.06 ms时，随着气泡的膨胀，其上表面与钢板间的距离减小，气泡受钢板的阻力变大，气泡上半部分难以自由膨胀，导致上半部分呈扁平状，下半部分为球形。在气泡收缩阶段(19.06～40.63 ms)，在Bjerknes力的影响下，气泡恢复球形，同时质心位置发生上移。在第1次气泡脉动结束后，在40.93 ms时能看到一束气液混合的高速射流冲向钢板，由相邻两帧高速摄像可以估算射流速度为65 m/s。同时，随着气泡的收缩，能清晰看到爆炸残留物在水流的带动下呈放射状指向钢板，见42.19 ms时图像。
图 7 当γ=1.282时气泡演变过程的实验图像Figure 7. Experimental images of bubble evolution process when γ=1.282气泡脉动周期和气泡最大半径是表征气泡脉动现象的两个主要参数，对TNT有经验公式[25]：
T=2.11W1/133(h+10)5/566 (4) 式中：T为第1次气泡脉动周期（s）。表2为气泡脉动周期和气泡最大半径的实验值及经验公式计算结果，其中气泡半径的表征方法同文献[22]。气泡脉动周期及气泡最大半径均小于经验公式计算结果。由于固支方板与气泡间的Bjerknes效应，使气泡坍塌速度加快，导致气泡脉动周期偏小。由于顶部固支方板阻碍了气泡自由膨胀，而经验公式在自由液面的前提下，所以实验测得的气泡最大半径小于经验公式计算结果。
表 2 气泡脉动周期和气泡最大半径的实验结果和经验公式结果Table 2. Experimental results and empirical formula results of bubble pulsation periods and bubble maximum radii实验 γ T/ms δT/% Rm/cm δR/% 实验 经验公式 实验 经验公式 1 0.684 37.81 41.52 8.94 21.50 21.94 2.01 2 0.684 38.44 41.52 7.42 21.50 21.94 2.01 3 0.798 39.69 41.43 4.20 21.60 21.92 1.46 4 0.798 39.37 41.43 4.97 21.40 21.92 2.37 5 0.913 39.81 41.35 3.72 21.50 21.90 1.83 6 0.913 39.75 41.35 3.87 21.60 21.90 1.37 7 0.913 39.69 41.35 4.01 21.60 21.90 1.37 8 1.282 38.75 41.08 5.67 21.00 21.85 3.89 9 1.282 39.75 41.08 3.24 21.10 21.85 3.43 10 1.282 39.56 41.08 3.70 21.00 21.85 3.89 1.2.3 吸附式射流和非吸附式射流现象
高速摄影图像结果表明，炸药在固支方板底部起爆时气泡射流的演化过程可以分为两种。（1）当爆距较小（γ=0.684, 0.798）时，炸药起爆后高温高压气体产物先发生膨胀，因钢板限制了气泡的发展，气泡在膨胀至最大体积时呈上扁下圆的不规则球形。接着气泡开始收缩，由于爆距较小，气泡在收缩过程中逐渐吸附在钢板上，且已吸附的部分呈漏斗形，而未吸附的部分保持球形，吸附过程结束后，气泡整体呈漏斗形。接着气泡开始坍塌，由于漏斗形气泡底部曲率半径最小，容易发生失稳，气泡从底部开始塌陷，形成指向钢板的射流。在这种情况下，射流作用在钢板上几乎没有穿过水层，能量损耗小，此时钢板对射流的作用效果较强。（2）爆距较大（γ=0.913, 1.282）时，炸药爆炸后气泡开始膨胀，膨胀至最大体积时，气泡呈上扁下圆的不规则球形。接着气泡开始收缩，由于爆距较大，钢板对气泡的Bjerknes引力较弱，气泡并未吸附在钢板上，而与钢板相隔一定距离继续收缩。在气泡收缩过程中，气泡底部发生坍塌，形成至下而上的射流。在这种情况下，射流在作用到钢板前需要穿过水层，因水中阻力较大，且气泡收缩过程中气泡顶部的水也以一定速度流向气泡，底部射流与顶部水流发生对流，进一步阻碍了底部射流的运动，此时射流对钢板的作用效果较弱。
2. 数值模拟
为了进一步探究固支方板底部水下爆炸气泡射流的演化过程，采用ABAQUS软件开展数值模拟。欧拉法允许材料在网格中流动，一般适用于气体及液体材料的数值模拟，拉格朗日法则常常用来描述结构的变形，而CEL方法结合了拉格朗日法和欧拉法的优点[28]。因此，本文中采用CEL方法开展固支方板底部水下爆炸的数值模拟。
2.1 流体的控制方程
流体的运动都遵循流体力学的基本控制方程[29]：连续性方程、动量方程和能量方程。
认为水和空气可压缩，可压缩流体的连续性方程为[29]：
∂ρ/∂t+∇⋅(ρV)=0 (5) 式中：ρ为流体的密度，t为时间，流体的速度矢量
V = (u,v,w) 。可压缩流体的动量方程为[29]：
∂(ρV)/∂t+∇⋅(ρVV)=−∇p+∇⋅τ+f (6) 式中：f为体积力，p为压力，
τ 为黏性应力张量。对于牛顿流体，动量方程可以展开[29]：
∂(ρV)/∂t+∇⋅(ρVV)=∇⋅(μ∇V)−∇p+∇⋅(μ(∇V)T)+∇(λ∇⋅V)+f (7) 式中：μ为动黏性系数，λ为膨胀黏性系数，通常λ＝−(2/3) μ。
能量方程为[29]：
∂(ρe)/∂t+∇⋅(ρVe)=−∇⋅˙qs−∇⋅(pV)+∇⋅(τ⋅V)+f⋅V+˙qV (8) 式中：e为质量总能量，
˙qs 为黏性耗散功率，˙qV 为黏性正应力的膨胀功率。在本文中，设水下爆炸过程为绝热过程，于是˙qs = 0,˙qV=0 。2.2 流体的状态方程
在三维情况下，可压缩流体的控制方程包含6个变量（如ρ、u、v、w、p和e），但控制方程只有5个，所以需加入描述压力p和密度ρ关系的状态方程。而一种材料存在多种不同形式的状态方程[30]，因此选取材料的合适状态方程对数值模拟至关重要，下面给出数值模拟中的状态方程。
2.2.1 水的状态方程
水的状态方程选用Mie-Grüneisen状态方程，该方程中压力是质量和比内能的函数，最常见形式为：
p−pH=Γρ(Em−EH) (9) 式中：pH为Hugoniot压力，EH为Hugoniot能量，Em为比内能，Mie-Grüneisen系数
Γ=Γ0ρ0/ρ ，ρ0为参考密度。Hugoniot能量EH和Hugoniot压力pH的关系为：
EH=ηpH2ρ0 (10) 式中：
η=1−ρ0/ρ 。Hugoniot压力通常为：pH=ρ0c20η(1−sη)2 (11) 设激波速度us和粒子速度up呈线性关系，即us=c0+sup，则Mie-Grüneisen状态方程为：
p=Γ0ρ0Em+ρ0c20η(1−sη)2(1−Γ0η2) (12) 式中：c0为介质中的声速，Г0为材料常数，s为us-up方程中的常数。
本文中，水的密度为1 000 kg/m3，水中的声速为1 450 m/s。
2.2.2 空气的状态方程
理想气体的状态方程为：
p+pa=ρRT (13) 式中：pa为环境压力， T为温度，R为气体常数。
本文中，空气的环境压力pa为101 300 Pa，空气密度ρ为1.225 kg/m3，此外，还定义空气的比热为717.4 J/(kg·K)。
2.2.3 爆轰产物的状态方程
爆轰产物的JWL状态方程为：
p=A(1−ωρR1ρ0)exp(−R1ρ0/ρ0ρρ)+B(1−ωρR2ρ0)exp(−R2ρ0/ρ0ρρ)+ωρEm (14) 式中：A、B、R1、R2和ω为材料参数，ρ0为装药密度，ρ为爆轰产物密度，Em为比内能。
本文中，材料参数分别为：ρ0=1 630 kg/m3，A=373.77 GPa，B=3.747 1 GPa，R1=4.15，R2=0.9，ω=0.35，Em=3.6 MJ/kg。
2.3 Johnson-Cook材料模型
固支方板材料为Q235钢，其材料密度为7 800 kg/m3，杨氏模量为200 GPa，泊松比为0.3。Q235钢的塑性力学性能采用Johnson-Cook材料模型，因固支方板在实验中未发生破坏，为简化计算不引入Johnson-Cook失效模型，因此Johnson-Cook材料模型为：
σeq=(A+Bεneq)(1+Cln˙ε∗eq)(1−T∗m) (15) 式中：A、B、n、C和m为模型参数，σeq为等效应力，εeq为等效塑性应变，
˙ε∗eq = ˙εeq/˙ε0 为无量纲等效塑性应变率，˙ε0 为参考应变率，T∗ = (T−Tr)/(Tm−Tr) 为无量纲温度，其中Tr为参考温度，Tm为材料的熔点，T为当前温度。Q235钢的材料参数分别为[31]：A=244.8 MPa，B=899.7 MPa，n=0.940，m=0.757，C=0.039 1。
2.4 数值模拟模型
图8为数值模拟模型示意图，采用1/4模型，水、空气、炸药采用欧拉单元建模，采取欧拉体积分数法定义欧拉域中水、空气、炸药的尺寸及相对位置。炸药设置为直径12.5 mm、高12.5 mm的圆柱。欧拉域顶部设置为无反射边界（Eulerian boundary>Outflow>Nonreflecting），欧拉域侧面和底面设置为固定边界（Symmetry/Antisymmetry/Encastre>Encastre），对称面设置为对称边界条件（Symmetry/Antisymmetry/Encastre>XSYMM /YSYMM）。欧拉域的网格采用非等距网格，爆心处网格尺寸为0.25 cm[30]，最外侧网格尺寸为4 cm，欧拉域网格总数为141万。固支方板采用壳单元建模，在固支方板角点处设置全约束边界条件，网格尺寸为4 cm。表3为数值模拟的距离参数。
表 3 数值模拟的距离参数Table 3. Distance parameters in numerical simulation数值模拟 d/cm γ 1 12.0 0.546 2 15.0 0.684 3 16.0 0.730 4 17.5 0.798 5 18.0 0.821 6 19.0 0.867 7 20.0 0.913 8 21.0 0.959 9 25.0 1.143 10 28.0 1.282 2.5 数值模拟结果分析与讨论
2.5.1 数值模拟与实验的对比
图9为γ=0.684时气泡演变过程的数值模拟结果与实验结果的对比。在数值模拟28.30～33.00 ms时和实验结果31.88～35.63 ms时，气泡吸附在钢板上的过程中，其上半部分呈漏斗形，下半部分保持球形；吸附过程结束（数值模拟34.00 ms）时，气泡整体呈漏斗形；随后 （数值模拟36.00～37.20 ms），气泡底部坍塌形成射流，最后（数值模拟41.90 ms）气泡在钢板附近溃散。
图 9 当γ=0.684时气泡演变过程的实验和数值模拟图像Figure 9. Experimental and numerical simulation images of bubble evolution process when γ=0.684图10为γ=0.798时气泡演变过程的数值模拟结果与实验结果的对比。实验和数值模拟结果均显示，气泡在收缩过程中吸附在钢板上，同时气泡吸附的过程与γ=0.684时相似。
图 10 当γ=0.798时气泡演变过程的实验和数值模拟图像Figure 10. Experimental and numerical simulation images of bubble evolution process when γ=0.798图11为γ=0.913时气泡演变过程的数值模拟结果与实验结果的对比。数值模拟和实验结果基本一致，气泡在收缩过程中有吸附在钢板上的趋势。在数值模拟36.00 ms时，钢板对气泡的Bjerknes力无法使气泡完全吸附在钢板上。在数值模拟36.80 ms时，射流发生，气泡的上表面并未贴合在钢板上。
图 11 当γ=0.913时气泡演变过程的实验和数值模拟图像Figure 11. Experimental and numerical simulation images of bubble evolution process when γ=0.913图12为γ=1.282时气泡演变过程的数值模拟结果与实验结果的对比。气泡在膨胀阶段(0～19.06 ms)由球形向上边下圆的不规则球形演化。在气泡收缩阶段，气泡上表面在钢板Bjerknes力作用下基本不发生移动，而下表面向上表面靠近，最终形成射流。在数值模拟37.00 ms时，射流贯穿气泡，气泡演化为环状气泡。
图 12 当γ=1.282时气泡演变过程的实验和数值模拟图像Figure 12. Experimental and numerical simulation images of bubble evolution process when γ=1.282按照图13的气泡上、下表面位置选取方式，图14为γ分别为0.684、0.798、0.913和1.282时气泡上、下表面位置曲线。气泡上、下表面位移曲线位移幅值的数值模拟和实验结果基本一致，但数值模拟的气泡脉动时间较短，误差在10%以内。在10～30 ms，气泡上表面位移曲线呈马鞍状，这是因为：钢板受冲击波作用后开始振动，顶部固支方板先向下运动，使气泡上表面下移，随着钢板的回弹，气泡上表面发生上移。
2.5.2 吸附式射流向非吸附式射流转化的临界点
图15～16分别为γ=0.821, 0.867时气泡演化过程的数值模拟图像。当γ=0.821时，在气泡收缩阶段（19～37 ms），气泡在钢板Bjerknes引力作用下整体吸附在钢板上，产生吸附式射流。当γ=0.867时，在气泡收缩阶段（19～37 ms），气泡受到钢板Bjerknes引力的作用，部分气泡向钢板运动，呈现吸附的趋势，但钢板对气泡的Bjerknes引力不足以使气泡完全吸附在钢板上，在气泡坍塌阶段并未发生吸附式射流，见37 ms时图像。因此，吸附式射流向非吸附式射流转化的临界点在γ=0.821～0.867。炸药起爆时，因钢板的阻挡，气泡上半部分难以自由膨胀，使气泡发生吸附时对应的爆距小于1倍最大气泡半径，依据本文的数值模拟结果，临界点在γ=0.821～0.867。
图 15 当γ=0.821时气泡演变过程的数值模拟图像Figure 15. Numerical simulation images of bubble evolution process when γ=0.821
图 16 当γ=0.867时气泡演变过程的数值模拟图像Figure 16. Numerical simulation images of bubble evolution process when γ=0.8672.5.3 爆距对射流速度特性的影响
图17为水射流顶部节点（气泡底部节点）的速度曲线，图18为不同爆距下射流速度云图演化的数值模拟图像。当γ=0.546～0.821时，射流均从气泡底部开始形成，随着爆距的增加，气泡在射流形成阶段的形状越来越细长，气泡底部的曲率越来越大，导致射流速度随着爆距增大而增加。当γ=0.867时，随着爆距的增加，钢板对气泡的作用效果减弱，在射流阶段气泡明显脱离钢板，同时射流速度减小。当γ=0.913时，能看到36.60 ms时气泡顶部形成一束射向下方的射流，但因钢板Bjerkness引力的限制，顶部向下的射流难以继续演化，同时，气泡底部在钢板Bjerkness引力和浮力的作用下形成由气泡底部射向气泡顶部的射流，在37.00 ms时能看到两股射流相互作用，形成了对射流。当γ=0.959时，看到36.60 ms时从气泡顶部形成一束裙摆状射流向下运动，同时，气泡底部发生失稳形成一束更强的射流向上运动，两股射流形成了对射流。当γ=1.143, 1.282时，现象较相近，在射流形成阶段，气泡上表面位置保持不变，气泡下表面向钢板运动，在浮力和Bjerkness引力作用下底部进一步坍塌，最终形成射流，但由于爆距较大，Bjerkness引力较弱，形成的射流速度也较低。综上所述，随着爆距的增大，射流速度先增大后减小，在吸附式射流和非吸附式射流转化的临界点附近，射流速度达到最大。
图 17 在不同距离参数下射流顶部节点的速度曲线Figure 17. Velocity curves of jet top nodes under different burst distances
18 在不同距离参数下射流速度演化过程的数值模拟图像18. Numerical simulation images of jet velocity evolution processes at different burst distances图19为不同距离参数下射流开始形成的时间，以图17中速度曲线起始阶段的拐点作为射流形成时间。当γ<0.867时，射流形成时间随着爆距增加而增加，当γ≥0.867时，射流形成时间随着爆距增加基本不变。这说明，顶部固支方板加速了气泡坍塌形成射流的过程，与张阿漫等[32]的观点一致。本文中还发现：发生吸附式射流时，爆心离固支方板越近，射流形成时间越早；发生非吸附式射流时，射流形成时间不受爆距影响。拟合可得射流形成时间Tf（ms）与距离参数γ的关系（相关系数为0.993 7）：
图 19 射流形成时间与距离参数的关系Figure 19. Relationship between jet formation time and distance parameterTf=36.26−2.77/(1+exp((γ−0.715)/0.05581)) (16) 图20为不同距离参数下射流的最大速度及达到最大速度时对应的时间。发生吸附式射流时，射流最大速度随着爆距的增加而增加，且射流达到最大速度出现的时间越来越迟，这是因射流形成时间变迟导致的；发生非吸附式射流时，射流最大速度随着爆距的增加逐渐减小，当γ≥0.959时，射流最大速度基本保持不变，而射流最大速度出现的时间越来越早。可拟合得到射流最大速度Vts（m/s）与距离参数γ的关系（相关系数为0.896 0）：
图 20 不同距离参数下最大射流速度和最大射流速度时间Figure 20. Maximum jet velocities and maximum jet velocity times at different distance parametersVts=74.31+546.3exp(−0.5((γ−0.821)/0.058 88)2)0.546≤γ≤1.282 (17) 和射流达到最大速度的时间Tts（ms）与距离参数γ的关系（相关系数为0.948 5）：
Tts=36.80−2.3/(1+(γ/0.78)16.3)0.546≤γ≤1.282 (18) 图21为不同爆距下射流击中钢板时的速度，当γ≥0.959时射流未能到达钢板，因此图中没有显示。当γ≤0.821时，射流击中钢板时的速度和射流击中钢板时的时间都随着爆距的增大而增大；当0.821＜γ≤0.913时，射流击中钢板时的速度随着爆距的增加迅速减小。这是因为，发生非吸附式射流时，钢板与气泡间存在较厚的水层，射流在穿过水层时速度迅速衰减，导致射流速度降低。当γ＞0.959时，因射流在水层中的衰减效应和对射流现象，导致射流未能到达钢板。拟合可得射流击中钢板时的速度Vhit（m/s）与距离参数γ的关系（相关系数为0.995 5）：
图 21 不同距离参数下射流击中钢板的速度和时间Figure 21. Velocities and times of jet hitting steel plate at different distance parametersVhit={−2.27+269.62exp(−(γ−0.842)2/(2×0.149592))0.546≤γ＜0.842−2.27+269.62exp(−(γ−0.842)2/(2×0.031782))0.842≤γ＜0.9130γ≥0.913 (19) 和射流击中钢板的时间Thit（ms）与距离参数γ的关系（相关系数为0.990 6）：
Thit = 36.31 + 9.167−11exp(γ/0.037 26) 0.546≤γ≤0.913 (20) 3. 结 论
开展了TNT炸药在水面固支方板底部不同爆距下的水下爆炸实验，观测了距离参数γ=0.684～1.282的气泡射流现象，通过水下压力传感器测得了冲击波和气泡脉动的压力数据，基于ABAQUS软件，采用CEL方法开展了系列数值模拟。在此基础上，得到如下结论。
（1）高速摄像图像显示：当距离参数γ=0.684, 0.798时，气泡吸附在钢板上，形成吸附式射流；当γ=0.913, 1.282时，气泡未吸附在钢板上，形成非吸附式射流。
（2）数值模拟结果表明，吸附式射流向非吸附式射流转化的临界点在距离参数为0.821～0.867。这是由于，在钢板的阻挡下，气泡上半部分难以自由膨胀，使气泡发生吸附时对应的爆距小于1倍最大气泡半径。
（3）在吸附式射流对应的爆距范围内，爆距越小，射流形成时间越早，说明顶部固支方板加速了气泡坍塌形成射流的过程，且气泡距离钢板越近，这种加速效果越明显。但在非吸附式射流对应的爆距范围，固支方板对气泡坍塌过程的影响较弱。
（4）在吸附式射流对应的爆距范围，随着爆距的增加，射流形成过程中的最大速度越大，射流达到最大速度的时间越迟，射流速度最大值可达到621 m/s。这是因为，爆距越大气泡坍塌时间越迟，气泡能量越集中，因此形成的射流速度越大。而在非吸附式射流对应的爆距范围，随着爆距的增大，射流最大速度迅速下降，当γ≥0.959时，射流最大速度基本为约100 m/s，最大速度出现的时间在小范围内变动。
（5）射流击中钢板时的速度随着爆距的增大先升后降，在临界点附近最大值可达到269 m/s，射流击中钢板的时间随着爆距的增加持续增加。由于水流的阻挡，当γ≥0.959时，射流未击中钢板。
-
图 4 实验装置和物理模型的比较
Figure 4. Comparison between the experimental apparatus and the physical model
图 5 不同工况条件下最大爆炸压力随体积分数和体积分数梯度变化曲线
Figure 5. Changes of the maximum explosion pressure with volume fraction and volume fraction gradient under different working conditions
图 6 不同工况条件下最大爆炸压力时刻值随体积分数和体积分数梯度变化曲线
Figure 6. Changes of the arrival time of the maximum explosion pressure with volume fraction and volume fraction gradientunder different working conditions
图 7 不同工况下最大压力上升速率随体积分数和体积分数梯度的变化
Figure 7. Changes of the maximum pressure rise rate with volume fraction and volume fraction gradient under different working conditions
图 8 不同工况下最大压力上升速率时刻随体积分数和体积分数梯度的变化
Figure 8. Changes of the arrival time of the maximum pressure rise rate with volume fraction and volume fraction gradient under different working conditions
图 9 不同工况下最高温度随体积分数和体积分数梯度的变化
Figure 9. Changes of the maximum temperature with volume fraction and volume fraction gradient under different working conditions
图 10 不同工况下最高温度时刻值随体积分数和体积分数梯度的变化
Figure 10. Changes of the arrival time of the maximum temperature with volume fraction and volume fraction gradient under different working conditions
表 1 模拟工况
Table 1. Simulated working conditions
工况 体积分数梯度/% 体积分数分布/% 区域1 区域2 区域3 区域4 区域5 1 0 6.517 2 2 6 6.5 7 7.5 8 3 0 7.034 4 4 6 7 8 9 10 5 0 7.551 6 6 6 7.5 9 10.5 12 7 0 7.697 8 6.5 6 7.625 9.25 10.875 12.5 9 0 7.827 10 7 6 7.75 9.5 11.25 13 11 0 7.958 12 7.5 6 7.875 9.75 11.625 13.5 13 0 8.067 14 8 6 8 10 12 14 
下载: 导出CSV
表 2 连通装置不同体积分数时最大爆炸压力相关参数
Table 2. Related parameters of the maximum explosion pressures for different volume fractions of connected devices
工况 体积分数/% 最大爆炸压力 位置 时刻/s 数值/MPa 1 6.517 区域5 0.157 0.716 3 7.034 区域5 0.147 0.758 5 7.551 区域5 0.139 0.797 7 7.696 区域5 0.137 0.810 9 7.827 区域5 0.135 0.820 11 7.958 区域5 0.133 0.829 13 8.067 区域5 0.132 0.834
下载: 导出CSV
表 3 连通装置不同体积分数梯度时最大爆炸压力相关参数
Table 3. Related parameters of the maximum explosion pressures for different volume fraction gradients of connected devices
工况 体积分数梯度/% 最大爆炸压力 位置 时刻/s 数值/MPa 2 2.0 区域5 0.154 0.754 4 4.0 区域5 0.147 0.820 6 6.0 区域5 0.151 0.814 8 6.5 区域5 0.154 0.817 10 7.0 区域5 0.156 0.821 12 7.5 区域5 0.158 0.821 14 8.0 区域5 0.160 0.791
下载: 导出CSV
表 4 连通装置不同体积分数时最大压力上升速率相关参数
Table 4. Related parameters of the maximum pressure rise rates for different volume fractions of connected devices
工况 体积分数/% 最大压力上升速率 位置 时刻/s 数值/（GPa·s−1） 1 6.517 区域2 0.150 0.194 3 7.034 区域3 0.138 0.220 5 7.551 区域2 0.132 0.247 7 7.697 区域2 0.130 0.255 9 7.830 区域2 0.129 0.260 11 7.958 区域2 0.127 0.265 13 8.067 区域3 0.123 0.267
下载: 导出CSV
表 5 连通装置不同体积分数梯度时最大压力上升速率相关参数
Table 5. Parameters related to the maximum pressure rise rates in the connected devices with different volume fraction gradients
工况 体积分数梯度/% 最大压力上升速率 位置 时刻/s 数值/（GPa·s−1） 2 2.0 区域3 0.144 0.220 4 4.0 区域2 0.141 0.263 6 6.0 区域2 0.144 0.240 8 6.5 区域2 0.146 0.231 10 7.0 区域2 0.148 0.234 12 7.5 区域2 0.151 0.250 14 8.0 区域2 0.153 0.247
下载: 导出CSV
表 6 连通装置不同均一体积分数下最高温度相关参数
Table 6. Related parameters of the maximum temperatures in connected devices with different volume fractions
工况 体积分数/% 最高温度 位置 时刻/s 数值/℃ 1 6.517 区域1 0.189 1 917.58 3 7.034 区域1 0.182 2 023.45 5 7.551 区域1 0.179 2 121.45 7 7.697 区域1 0.176 2 156.08 9 7.827 区域1 0.175 2 181.71 11 7.958 区域1 0.173 2 209.39 13 8.067 区域1 0.172 2 223.58
下载: 导出CSV
表 7 连通装置不同体积分数梯度时最高温度相关参数
Table 7. Rrelated parameters of the maximum temperatures in connected devices with different volume fraction gradients
工况 体积分数梯度/% 最高温度 位置 时刻/s 数值/℃ 2 2.0 区域5 0.208 1 946.72 4 4.0 区域5 0.201 2 250.03 6 6.0 区域5 0.205 2 248.07 8 6.5 区域5 0.207 2 247.44 10 7.0 区域5 0.208 2 291.28 12 7.5 区域5 0.217 2 302.43 14 8.0 区域5 0.242 2 251.81
下载: 导出CSV
表 8 连通装置不同体积分数时最大速度相关参数
Table 8. Related parameters of the maximum velocities in connected devices with different volume fractions
工况 体积分数/% 最大速度 位置 时刻/s 数值/（m·s−1） 1 6.517 区域4 0.141 450.01 3 7.034 区域4 0.133 486.93 5 7.551 区域4 0.125 517.60 7 7.696 区域4 0.122 527.40 9 7.827 区域4 0.121 533.92 11 7.957 区域4 0.119 540.47 13 8.067 区域4 0.118 543.77
下载: 导出CSV
表 9 连通装置存在不同体积分数梯度时最大速度相关参数
Table 9. Related parameters of the maximum velocities in connected devices with different volume fraction gradients
工况 体积分数梯度/% 最大速度 位置 时刻/s 数值/（m·s−1） 2 2.0 区域4 0.138 445.20 4 4.0 区域4 0.133 456.60 6 6.0 区域4 0.136 448.40 8 6.5 区域4 0.138 443.01 10 7.0 区域4 0.141 437.55 12 7.5 区域4 0.143 431.70 14 8.0 区域4 0.145 427.08
下载: 导出CSV
-
[1] BARTKNECHT W. Explosion course prevention protection [M]. Berlin: Springer-Verlag, 1981. [2] PHYLAKTON H, ANDREWS G E. Gas explosions in linked vessels [J]. Journal of Loss Prevention in the Process Industries, 1993, 6(1): 15–19. DOI: 10.1016/0950-4230(93)80015-E. [3] LUNN G A, HOLBROW P, ANDRES S, et al. Dust explosions in totally enelosed interconnected vessel systems [J]. Journal of Loss Prevention in the Proeess Industries, 1995, 9(1): 45–58. [4] HOLBROW P, ANDRES S, LUNN G A. Dust explosions in interconnected vented vessels [J]. Journal of Loss Prevention in the Process Industries, 1996, 9(1): 91–103. DOI: 10.1016/0950-4230(95)00055-0. [5] HOLBROW P, LUNN G A. Dust explosion protection in linked vessels: guidance for containment and venting [J]. Journal of Loss Prevention in the Process Industries, 1999, 12: 227–234. DOI: 10.1016/S0950-4230(98)00050-3. [6] 严建骏, 蒋军成, 王志荣, 等. 连通容器内预混气体爆炸过程的实验研究 [J]. 中国安全生产科学技术, 2009, 4(6): 10–14.YAN J J, JIANG J C, WANG Z R, et al. Experimental investigation into explosion of premixed gases in linked vessels [J]. Safety Production Science and Technology in China, 2009, 4(6): 10–14. [7] 尤明伟, 蒋军成, 王志荣, 等. 连通容器中不同连通管径爆炸数值模拟分析 [J]. 工业安全与环保, 2010, 36(12): 25–26.YOU M W, JIANG J C, WANG Z R, et al. Numerical simulation of gas explosion in linked vessels with different pipe diameters [J]. Industrial Safety and Environmental Protection, 2010, 36(12): 25–26. [8] 王志荣, 蒋军成, 郑杨艳. 连通容器内气体爆炸过程的数值分析 [J]. 化学工程, 2006, 34(10): 13–16. DOI: 10.3969/j.issn.1005-9954.2006.10.004.WANG Z R, JIANG J C, ZHENG Y Y. Numerical analysis of gas explosion process in linked vessels [J]. Chemical Engineering, 2006, 34(10): 13–16. DOI: 10.3969/j.issn.1005-9954.2006.10.004. [9] 尤明伟, 喻源, 蒋军成, 等. 不同管长条件下连通容器预混气体的爆炸 [J]. 燃烧科学与技术, 2012, 18(3): 256–259.YOU M W, YU Y, JIANG J C, et al. Premixed flammable gas explosion in containers connected by pipes with different lengths [J]. Journal of Combustion Science and Technology, 2012, 18(3): 256–259. [10] 王志荣, 周超, 师喜林, 等. 连通容器内预混气体泄爆过程 [J]. 化工学报, 2011, 62(1): 287–291.WANG Z R, ZHOU C, SHI X L, et al. Gas explosion venting of premixed gases in linked vessels [J]. CIESC Journal, 2011, 62(1): 287–291. [11] THOMAS G O, SUTTON P, EDWARDS D H. The behavior of detonation waves at concentration gradients [J]. Combustion Flame, 1991(84): 312–322. [12] KUZNETSOV M S, ALEKSEEV V I, DOROFEEV S B, et al, Detonation propagation, decay, and reinitiation in nonuniform gaseous mixtures [C]// Twenty-Seventh International Symposium on Combustion/The Combustion Institute, 1998:2241–2247. [13] VOLLMER, K G, ETTNER F, SATTELMAYER T. Influence of concentration gradients on flame acceleration in tubes [J]. Energetic Matter Science Technology, 2011, 72: 74–77. [14] VOLLMER K G, ETTNER F, SATTELMAYER T. Deflagration-to-detonation transition in hydrogen-air mixtures with a concentration gradient [J]. Combust Science & Technology, 2012, 184(10-11): 1903–1915. DOI: 10.1080/00102202.2012.690652. [15] KESSLER D A, GAMEZO V N, ORAN E S. Wave structures and irregular detonation cells in methane-air mixtures with concentration gradients[C]//49th AIAA Aerospace Science Meeting. Orlando,Florida, 2011. [16] ISHII K, KOJIMA M. Behavior of detonation propagation in mixtures with concentration gradients [J]. Shock Waves, 2007, 17(1): 95–102. [17] WANG C J, WEN J X. Numerical simulations of hydrogen-air detonation wave propagation in a non-uniform semi-confined flat layer [C]//The Sixteenth International Colloquium on Dust Explosions and the Eleventh Colloquium on Gas, Vapour, Liquid and Hybrid Explosions. Bergen, Norway, 2014. [18] HAN W H, WANG C, CHUNG K. Role of transversal concentration gradient in detonation propagation [J]. Journal of Fluid Mechanics, 2019, 865: 602–649. DOI: 10.1017/jfm.2019.37. [19] 陈昊驰. 柱状空间内不同浓度甲烷爆炸传播特性的实验研究[D]. 石家庄: 华北科技学院, 2018: 38−44. [20] 王信群, 王婷, 徐海顺, 等. BC粉体抑爆剂改性及抑制甲烷/空气混合物爆炸 [J]. 化工学报, 2015, 66(12): 5171–5178.WANG X Q, WANG T, XU H S, et al. Modification of commercial BC dry chemical powder suppressant and experiments on suppression of methane-air explosion [J]. CIESC Journal, 2015, 66(12): 5171–5178. [21] XIAO H H, HE X C, DUAN Q L, et al. An investigation of premixed flame propagation in a closed combustion duct with a 90° blend [J]. Applied Energy, 2014, 134: 248–256. DOI: 10.1016/j.apenergy.2014.07.071. [22] 尤明伟. 连通容器气体爆炸及泄爆动力学过程研究[D]. 南京: 南京工业大学,2011. [23] 王志荣. 受限空间气体爆炸传播及其动力学过程研究[D]. 南京: 南京工业大学, 2005. 期刊类型引用(17)
1. 周刚，孔阳，崔洋洋，钱新明，傅砺烨，张琦. 城市地下排水管道中燃气爆炸及气-液两相耦合作用规律. 爆炸与冲击. 2024(03): 90-104 . 
本站查看2. 谢洋洋，孙岩. 水下爆炸多气泡耦合形态及载荷分布特性. 兵器装备工程学报. 2024(06): 10-18 . 
百度学术3. 李庆海，于福临，郭文琦，苏超，李战全，贾风光，龚月滢. 层冰边界条件下多气泡破冰试验研究. 山东交通学院学报. 2024(03): 104-111 . 百度学术4. 方厚林，卢强，郭权势，李国亮，刘存旭，陶思昊，张德志. 水下爆炸冲击波和气泡行为自由面效应的实验研究. 爆炸与冲击. 2024(08): 149-160 . 本站查看5. 陈钰帆，黄锡堤，赖志超，秦健，孟祥尧，文彦博，黄瑞源. 板结构水下爆炸下的空化特性. 火炸药学报. 2024(08): 714-729 . 百度学术6. 许江，陈涛，王立志，郭君. 不同水底介质浅水爆炸载荷特性试验研究. 中国造船. 2024(04): 180-188 . 百度学术7. 刘靖晗，唐廷，韦灼彬，高屹，李凌锋，张星. 水下爆炸作用下高桩码头损伤特性数值模拟研究. 兵器装备工程学报. 2024(10): 53-60 . 百度学术8. 赖志超，邓硕，秦健，迟卉，孟祥尧，文彦博，黄瑞源. 不同类型炸药近场水下爆炸下固支方板动态响应研究. 工程力学. 2024(11): 179-194 . 百度学术9. 叶记飞，李华，李振瀚，宋执权，傅鹏. 15 kV爆炸开关触头开断过程研究. 强激光与粒子束. 2024(11): 140-147 . 百度学术10. 剡小军，孙浩，马林，张旭辉，吴曦，雍顺宁. 水下爆炸作用下典型舰船目标毁伤特性. 兵工学报. 2024(S2): 215-221 . 百度学术11. 高龙翔，高涵，潘文，薛乐星，冯晓军. 光电测试技术在炸药爆轰性能研究中的应用进展. 火炸药学报. 2024(12): 1055-1073 . 百度学术12. 高涵，宋振伟，孔祥韶，周沪，郑成，吴卫国. 加筋板结构参数对其抗水下爆炸能力的影响分析. 中国舰船研究. 2023(05): 180-193 . 百度学术13. 巨圆圆，熊展，张磊，赵鹏铎，杜志鹏. 水下高压气泡动态特性数值模拟. 舰船科学技术. 2023(20): 14-17 . 百度学术14. 邓硕，赖志超，秦健，孟祥尧，迟卉，黄瑞源. 复杂边界条件下近场水下爆炸对固支方板的毁伤效应. 爆炸与冲击. 2023(11): 40-57 . 本站查看15. 刘元凯，秦健，迟卉，孟祥尧，文彦博，黄瑞源. 近场水下爆炸气泡与目标尺寸匹配关系研究. 爆破. 2023(04): 142-153 . 百度学术16. 齐贺阳，田永卫. 基于ABAQUS CEL的高速水冲击载荷仿真与试验研究. 计测技术. 2022(03): 36-42 . 百度学术17. 高文博，赵振宇，任建伟，卢天健. 表面覆水对浅埋爆炸冲量传递的影响. 水下无人系统学报. 2022(03): 292-299 . 百度学术其他类型引用(12)
-
下载:
百度学术
计量
- 文章访问数: 758
- HTML全文浏览量: 323
- PDF下载量: 74
- 被引次数: 29