Large eddy simulation on the vented gasoline-air mixture explosions in a semi-confined pipe with continuous circular hollow obstacles
-
摘要: 采用WALE模型和Zimont预混火焰模型对内置圆孔障碍物油气泄压爆炸火焰特性进行了大涡模拟,并将大涡模拟计算结果和RNG k-ε湍流模型计算结果以及实验结果进行对比分析,验证了大涡模拟的精确性。结果表明:(1)大涡模拟在预测油气爆炸超压、火焰传播速度以及火焰形态变化等方面比RNG k-ε湍流模型精确度更高,且能表现出更多流场的精细化结构;(2)障碍物诱导管道内形成湍流度较高的流场区域,导致火焰产生褶皱弯曲变形,增大火焰面积,加速火焰传播;(3)爆炸超压、火焰传播速度和火焰面积内在联系密切,具有显著的耦合性,且随时间的变化趋势存在高度的一致性。Abstract: WALE LES model coupled with Zimont premixed combustion model were applied to study the characteristics of the flame of gasoline-air mixture explosions. Results by LES simulation, RNG k-ε model and experiments were compared with each other, and the conclusions show that:(1) LES is more accurate than RNG k-ε turbulence model in predicting gasoline-air mixture explosion overpressure, flame propagation speed and the variation of flame shape, and it can show more refined structure of flow field; (2) The obstacles induce a turbulent flow field in the pipeline, which can cause the flame to produce wrinkle and bend deformation, increase the flame surface area and accelerate the flame propagation; (3) The explosion overpressure, the flame propagation speed and the flame surface area are closely related. They have significant coupling and show nearly the same changing tendency.
-
在天体物理、地球物理、固体物理、核武器以及惯性约束聚变等研究中, 常常需要精确了解极端压力(太帕)条件下物质的动态响应特性和状态方程, 而常规的动态加载技术, 如化爆、轻气炮等只能够将飞片发射至8 km/s以下的速度, 在重材料中实现的压力也不过数百万大气压, 远不能满足上述研究的要求。近年来, 为获得更高的压力, 人们发展了一些超高压加载技术, 如磁压缩、激光驱动、多级爆轰等。其中, 利用成熟的二级轻气炮加载飞片撞击炸药产生强爆轰, 从而驱动重金属飞片达到更高速度的方法越来越受到重视。由于实验所需设备简单、费用低, 而且在测试精度、数据可信度等方面具有较明显的优势, 此项技术成为开展极端高压物理研究的重要手段之一。
强爆轰是指爆轰产物中的压力在p-v平面上的位置高于炸药正常爆轰时CJ点的状态。高速飞片冲击炸药起爆等方法可以产生强爆轰状态。利用强爆轰驱动飞片获得高压力和高速度飞片的实验研究, 逐渐成为超高速发射研究的热点问题。P.K.Tang[1-3]利用强爆轰方法研究了HMX和TATB炸药, 得到了用标准的JWL状态方程计算塑形黏结炸药HMX和TATB的强爆轰状态与实验差别较大的结论。赵峰等[4]、文尚刚等[5-6]利用轻气炮(或炸药)发射(或驱动)一级飞片冲击目标炸药驱动二级飞片实验, 研究了强爆轰驱动飞片情形, 实验结果表明二级飞片的速度可达到6.0~10.0 km/s。
随着强爆轰实验和理论研究的深入发展, 以及数值模拟技术在爆轰领域的广泛应用, 要求对强爆轰产物状态的描述更精确。目前, 常用爆轰产物的状态方程有JWL、BKW、HOM等。但J. H. Kineke等[7]、L. Green等[8]的研究表明, 用于计算CJ爆轰产物的状态方程(JWL)不适用描述强爆轰的行为特征, 而这种差别对高能炸药尤其明显。Л.B.AлЬтшулер[9]提出了爆轰产物高压状态方程, 并计算了PBX-9404等4种固体炸药的强D-u曲线和p-V曲线, 与实验结果符合较好。邓全农等[10]建立了DG01A炸药的爆轰产物方程, 利用强爆轰状态方程计算的DG01A、PBX-9404炸药强爆轰状态与实验符合较好。M.van Thiel等[11]、L.Green等[12]根据强爆轰实验确定的爆轰产物Hugoniot曲线, 推导出实验等熵线, 将JWL状态方程的Grüneisen系数作为相对比容的函数。标准JWL状态方程能较精确地描述爆轰产物的膨胀驱动做功过程, P.K.Tang在标准JWL的基础上, 给出了适用于高压状态的爆轰产物的等熵膨胀线方程(JWL), 使它能够用来模拟炸药的强爆轰的流体力学过程。用该方程计算的PBX-9501和PBX-9502炸药的强爆轰Hugoniot参数和强爆轰声速与实验结果符合很好。潘昊等[13]采用Hybrid反应率结合JWLT状态方程, 研究了LX17、超细TATB等钝感炸药的冲击起爆过程, 并计算了爆轰波对碰现象, 获得了较好的结果。姚阳等[14]在DEFFL二维流体弹塑性流体力学程序中使用JWLT状态方程, 研究了飞片高速撞击PBX-9501后冲击起爆驱动二级飞片的实验模型, 并将计算结果和实验结果进行比较。
已有关于强爆轰的研究, 主要针对炸药爆轰产物的动力学行为, 对强爆轰驱动飞片的超高速发射研究相对较少, 用流体力学程序对强爆轰驱动飞片的数值模拟计算结果, 没有达到实验结果和计算结果比较一致的程度。
本文中, 利用自编一维平面弹塑性流体动力学程序, 使用JWLT状态方程描述强爆轰产物的流体力学行为, 对强爆轰驱动飞片实验进行了数值模拟, 得到二级飞片的自由面粒子速度曲线, 与实验结果比较吻合, 拟为解读实验结果和设计新的强爆轰实验提供计算依据和参考。
1. 计算方法
1.1 一维弹塑性流体力学Lagrange方程组
拉格朗日一维轴对称的运动方程、连续方程和能量方程分别为:
u=−r∂σr/∂m−v(σr−σθ)/r (1) v=r∂r/∂m (2) ˙E=−(q+p)˙v+v[˙ε1s1+˙ε2s2] (3) 式中:
为速度, r为半径, m为质量, σr径向应力, σθ为环向应力, v=1/ρ为比容, q为人为黏性力, p为状态方程计算的压力, εi(i=1, 2)为应变张量ε的分量, ρ为密度, si(i=1, 2)为应力张量s的分量。
1.2 本构关系
本构关系的法向应力、切向(环向)应力、应力偏量和应变率分别为(力的方向, 拉正压负):
σr=p+q+s1 (4) σθ=p+q+s2 (5) ˙si=2μ[εi+˙v/(2v)]i=1,2 (6) ˙ε1=−∂u/∂r,˙ε2=−u/r,˙ε3=0 (7) 式中:p为由状态方程计算得到的压力, q为人为黏性, s为应力偏量, ε为应变偏量, Y0为材料的初始屈服强度, μ为材料的剪切模量。
1.3 状态方程
对于惰性材料, 使用Grüneisen形式的状态方程。在压缩状态, 状态方程为:
p=γ(E−EH)/v+c20(v−v0)/[v0−λ(v0−v)]2 (8) 式中:γ为Grüneisen系数,E和EH分别为比内能和Hugoniot状态下的比内能,v和v0分别为比容和初始时刻的比容,c0和λ都为Grüneisen系数。可根据冲击波关系式D=c0+λup+λ'up2,由实验结果拟合得到c0和λ。
在膨胀状态, 状态方程为:
p=γ[E−Cv(v/v0−1)/(3β)]/v (9) 式中:β为线膨胀系数。
对于爆轰产物, 使用JWLT状态方程描述行为特征。P.K.Tang在对PBX-9501和PBX-9502炸药的爆轰产物状态方程进行研究时, 首先假设爆轰产物的标准JWL状态方程在低压部分是可以接受的, 而对强爆轰区中的状态进行描述时有所欠缺。为此, 他建议增加Hugoniot曲线在CJ点以上的斜率, 更好地匹配强爆轰Hugoniot实验数据, 同时也增加强爆轰产物的声速。思路类似于可压缩理想气体在高压区的概念, 只对JWL状态方程中的高压指数项做修正。
经P.K.Tang修改后的JWLT状态方程为:
p=[1+Fp(ˉv)]Ae−R1ˉv+Be−R2ˉv+Cˉv−(1+ω) (10) e=[1+Fe(v)]Ae−R1ˉv/R1+Be−R2ˉv/R2+Cˉv−ω/ω (11) 修正项为当前比容与CJ比容之差的简单多项式:
Fp(ˉv)=A0(ˉvCJ−ˉv)2+B0(ˉvCJ−ˉv)3 (12) Fe(ˉv)=(A0−3B0/R1){2[1−e−R1(ˉvCJ−ˉv)]/R21−2(ˉvCJ−ˉv)/R1+(ˉvCJ−ˉv)2}+B0(ˉvCJ−ˉv)3 (13) 修正后的等熵线压力和内能的表达式中, 只有两个新的参数A0和B0。
2. 计算模型和参数
二级轻气炮驱动一级飞片, 一级飞片经过加速后撞击直径约16 mm、厚度约4 mm的JO-9159炸药, 使其发生强爆轰, 并驱动相应的次级钢飞片(直径约16 mm、厚度约0.5 mm)达到更高的速度。为了防止强爆轰波过高压力对二级飞片的冲击破坏, 在二级飞片前放置一层厚度约0.5 mm的有机玻璃缓冲层。
共进行了3发实验, 参数见表 1。表中, L1、L2、L3和L4分别为钢飞片、JO-9159炸药、有机玻璃缓冲层和次级钢飞片的厚度。
表 1 实验模型参数Table 1. The parameter of experiment modelNo. L1/mm v0/(km·s-1) L2/mm L3/mm L4/mm 1 2.52 4.97 3.99 0.51 0.47 2 2.52 6.09 3.98 0.50 0.48 3 2.53 5.89 4.00 0.50 0.36 计算中, 使用弹塑性流体本构模型和Grüneisen状态方程描述惰性材料的力学行为, 使用JWLT状态方程描述炸药爆轰产物的行为。
材料的强度参数和状态方程参数见表 2, JO-9159炸药的参数分别为:ρ=1.860 g/cm3, D=8.862 km/s, p=37.0 GPa, Q=5.267 kJ/g, A=934.8 GPa, B=12.72 GPa, R1=4.6, R2=1.1, ω=0.37, A0=2.52 GPa, B0=36.48 GPa。
表 2 材料参数Table 2. The parameter of material材料 强度参数 状态方程参数 ρ/(g·cm-3) G/GPa Y/GPa c0/(km·s-1) λ Γ 钢 7.85 83.1 0.60 3.449 2.087 1.69 有机玻璃 1.186 2.32 0.15 2.566 1.624 0.85 钢飞片以一定的初速度撞击JO-9159炸药, 由于飞片的高速运动, 使炸药爆轰产物中压力高于其CJ压力, 形成强爆轰。实验中, 获得了二级飞片的自由面速度曲线。
3. 计算结果与实验结果的比较
计算得到实验3炸药中不同时刻的压力和密度分布, 如图 1~2所示。由图 1可知, 在0.2
s时刻, 部分炸药已经爆轰, 且压力高达84 GPa, 远远超过JO-9159炸药的CJ压力37 GPa, 随着时间的推进, 压力最大值154 GPa出现在0.6
s, 之后, 随着飞片PMMA和Fe的运动, 炸药中压力随之下降, 在1.0
s时刻, 炸药爆轰产物中压力仅为17GPa。计算结果表明, 在炸药中形成了很强的强爆轰, 驱动重金属飞片达到更高的速度是可能的。
计算得出的3个实验模型的二级飞片的自由面速度与实验结果的比较, 如图 3~5所示。由图看出, 除自由面粒子速度在第一次起跳后, 计算结果稍高于实验结果外, 计算得到的自由面速度曲线和实验测量结果几乎重合。这表明, 用JWLT状态方程描述炸药的强爆轰状态是可行的; 利用强爆轰驱动重金属飞片达到更高的速度是可行的。
4. 讨论
利用JWLT状态方程描述炸药的超压状态方程, 使用自编一维流体动力学程序, 计算了金属飞片以一定速度撞击炸药产生强爆轰驱动组合飞片实验模型, 得到了撞击过程中不同时刻炸药中的压力和密度分布和飞片自由面粒子速度曲线。
计算结果和实验测量结果比较表明, 计算结果和实验结果符合较好, 表明了自编程序使用JWLT状态方程用于计算炸药的强爆轰问题的正确性, 可以用该程序用于强爆轰问题的计算和数值模拟, 设计新实验。
然而计算中发现, 有很多的因素在影响爆轰产物的流体力学行为。如本文中使用的JWLT状态方程参数中由实验数据拟合得到的A0和B0, 是否可以使用其他的拟合方法获得, 就是目前正在研究的问题。对于强爆轰, 目前还没有公认的、较好的解决方法, 还需要进行更多、更细致、更深入的研究。
-
-
[1] 杜扬, 李国庆, 吴松林, 等.T型分支管道对油气爆炸强度的影响[J].爆炸与冲击, 2015, 35(5):729-734. http://www.bzycj.cn/CN/abstract/abstract9524.shtmlDU Yang, LI Guoqing, WU Songlin, et al. Explosion intensity of gasoline-air mixture in the pipeline containing a T-shaped branch pipe[J]. Explosion and Shock Waves, 2015, 35(5):729-734. http://www.bzycj.cn/CN/abstract/abstract9524.shtml [2] LI G, DU Y, QI S, LI Y, et al. Explosions of gasoline-air mixtures in a closed pipe containing a T-shaped branch structure[J]. Journal of Loss Prevention in the Process Industries, 2016, 43:529-536. doi: 10.1016/j.jlp.2016.07.022 [3] JOHANSEN C, CICCARELLI G. Visualization of the unburned gas flow field ahead of an accelerating flame in an obstructed square channel[J]. Combustion and Flame, 2009, 156(2):405-416. doi: 10.1016/j.combustflame.2008.07.010 [4] 陈鹏, 李艳超, 黄福军, 等.方孔障碍物对瓦斯火焰传播影响的实验与大涡模拟[J].爆炸与冲击, 2017, 37(1):21-26. http://www.bzycj.cn/CN/abstract/abstract9681.shtmlCHEN Peng, LI Yanchao, HUANG Fujun, et al. LES approach to premixed methane/air flame propagating in the closed duct with a square-hole obstacl[J]. Explosion and Shock Waves, 2017, 37(1):21-26. http://www.bzycj.cn/CN/abstract/abstract9681.shtml [5] 王公忠, 张建华, 李登科, 等.障碍物对预混火焰特性影响的大涡数值模拟[J].爆炸与冲击, 2017, 37(1):68-76. http://www.bzycj.cn/CN/abstract/abstract9687.shtmlWANG Gongzhong, ZHANG Jianhua, LI Dengke, et al. Large eddy simulation of impacted obstacles' effects on premixed flame's characteristics[J]. Explosion and Shock Waves, 2017, 37(1):68-76. http://www.bzycj.cn/CN/abstract/abstract9687.shtml [6] WEN X, YU M, JI W, et al. Methane-air explosion characteristics with different obstacle configurations[J]. International Journal of Mining Science and Technology, 2015, 25(2):213-218. doi: 10.1016/j.ijmst.2015.02.008 [7] WEN X, YU M, LIU Z, et al. Large eddy simulation of methane-air deflagration in an obstructed chamber using different combustion models[J]. Journal of Loss Prevention in the Process Industries, 2012, 25(4):730-738. doi: 10.1016/j.jlp.2012.04.008 [8] NA'INNA A M, PHYLAKTOU H N, ANDREWS G E. The acceleration of flames in tube explosions with two obstacles as a function of the obstacle separation distance[J]. Journal of Loss Prevention in the Process Industries, 2013, 26(6):1597-1603. doi: 10.1016/j.jlp.2013.08.003 [9] BLANCHARD R, ARNDT D, GRÄTZ R, et al. Explosions in closed pipes containing baffles and 90 degree bends[J]. Journal of Loss Prevention in the Process Industries, 2010, 23(2):253-259. doi: 10.1016/j.jlp.2009.09.004 [10] HISKEN H, ENSTAD G A, MIDDHA P, et al. Investigation of concentration effects on the flame acceleration in vented channels[J]. Journal of Loss Prevention in the Process Industries, 2015, 36:447-459. doi: 10.1016/j.jlp.2015.04.005 [11] 杜扬, 李国庆, 王世茂, 等.障碍物数量对油气泄压爆炸特性的影响[J].化工学报, 2017, 68(7):2946-2955. http://d.old.wanfangdata.com.cn/Periodical/hgxb201707040DU Yang, LI Guoqing, WANG Shimao, et al. Effects of obstacle number on the characteristics of vented gasoline-air mixture explosions[J]. CIESC Journal, 2017, 68(7):2946-2955. http://d.old.wanfangdata.com.cn/Periodical/hgxb201707040 [12] 李国庆, 杜扬, 齐圣, 等.障碍物对油气-空气混合气体泄压爆炸火焰传播特性影响[J].中国安全生产科学技术, 2017, 13(1):163-168. http://d.old.wanfangdata.com.cn/Periodical/zgzyaqwsgltxrz201701042LI Guoqing, DU Yang, QI Sheng, et al. Effects of obstacles on the flame propagation characteristics of vented gasoline-air mixtures explosions[J]. Journal of Safety Science and Technology, 2017, 13(1):163-168. http://d.old.wanfangdata.com.cn/Periodical/zgzyaqwsgltxrz201701042 [13] LI G, DU Y, LIANG J, et al. Characteristics of gasoline-air mixture explosions with different obstacle configurations[J]. Journal of the Energy Institute, 2018, 91:194-202. doi: 10.1016/j.joei.2017.01.001 [14] 王世茂, 杜扬, 张少波, 等.顶部开口条件下油罐油气爆炸数值模拟[J].后勤工程学院学报, 2015(4):51-56. doi: 10.3969/j.issn.1672-7843.2015.04.009WANG Shimao, DU Yang, ZHANG Shaobo, et al. Numerical simulation on fuel air mixture explosion in the oil tank with an open top[J]. Journal of Logistics Engineering University, 2015(4):51-56. doi: 10.3969/j.issn.1672-7843.2015.04.009 [15] 李国庆, 杜扬, 白洁, 等.T型分支管道内油气爆炸火焰传播特性数值模拟研究[J].中国安全生产科学技术, 2016, 12(9):120-127. http://d.old.wanfangdata.com.cn/Periodical/zgzyaqwsgltxrz201609022LI Guoqing, DU Yang, BAI Jie, et al. Numerical Simulation on flame propagation features of gasoline-vaper explosion in T-shaped branch pipe[J]. Journal of Safety Science and Technology, 2016, 12(9):120-127. http://d.old.wanfangdata.com.cn/Periodical/zgzyaqwsgltxrz201609022 [16] 蒋新生, 杜扬, 唐晓寅, 等.油料洞库油气爆炸抑制数值模拟[J].后勤工程学院学报, 2008, 24(4):13-18. doi: 10.3969/j.issn.1672-7843.2008.04.004JIANG Xinsheng, DU Yang, TANG Xiaoyin, et al. Numerical simulation on fuel-air mixture explosion suppression in underground oil depot[J]. Journal of Logistics Engineering University, 2008, 24(4):13-18. doi: 10.3969/j.issn.1672-7843.2008.04.004 [17] NICOUD F, DUCROS F. Subgrid-scale stress modelling based on the square of the velocity gradient tensor[J]. Flow, Turbulence and Combustion, 1999, 62(3):183-200. doi: 10.1023/A:1009995426001 [18] 何标, 蒋新生, 孙国骏, 等.基于大涡模拟的气体羽流分层特性数值模拟[J].后勤工程学院学报, 2015(1):38-44. doi: 10.3969/j.issn.1672-7843.2015.01.008HE Biao, JIANG Xinsheng, SUN Guojun, et al. Numerical simulation of gas plume stratification based on large eddy simulation[J]. Journal of Logistics Engineering University, 2015(1):38-44. doi: 10.3969/j.issn.1672-7843.2015.01.008 [19] ZIMONT V, BATTAGLIA V. Joint RANS/LES approach to premixed flames modelling in the context of the TFC combustion model[J]. Flow, Turbulence and Combustion, 2006, 77(1):305-331. https://www.researchgate.net/publication/286970583_Joint_RANSLES_approach_to_premixed_flames_modelling_in_the_context_of_the_TFC_combustion_model [20] ZIMONT V L. Gas premixed combustion at high turbulence. Turbulent flame closure combustion model[J]. Experimental Thermal & Fluid Science, 2000, 21(1/2/3):179-186. http://d.old.wanfangdata.com.cn/NSTLQK/NSTL_QKJJ027474251/ [21] MANNAA O, MANSOUR M S, ROBERTS W L, et al. Laminar burning velocities at elevated pressures for gasoline and gasoline surrogates associated with RON[J]. Combustion & Flame, 2015, 162(6):2311-2321. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=JJ0235325261 [22] PATEL S N D H, JARVIS S, IBRAHIM S S, et al. An experimental and numerical investigation of premixed flame deflagration in a semiconfined explosion chamber[J]. Proceedings of the Combustion Institute, 2002, 29(2):1849-1854. doi: 10.1016/S1540-7489(02)80224-3 [23] XU C, CONG L, YU Z, et al. Numerical simulation of premixed methane-air deflagration in a semi-confined obstructed chamber[J]. Journal of Loss Prevention in the Process Industries, 2015, 34:218-224. doi: 10.1016/j.jlp.2015.02.007 [24] SARLI V D, BENEDETTO A D, RUSSO G. Using large eddy simulation for understanding vented gas explosions in the presence of obstacles[J]. Journal of Hazardous Materials, 2009, 169(1/2/3):435-442. https://www.researchgate.net/publication/24396517_Using_Large_Eddy_Simulation_for_understanding_vented_gas_explosions_in_the_presence_of_obstacles 期刊类型引用(3)
1. 毛强,蒋志强. 高精度平面构件平面度及壁厚非接触测量装置研制. 机械. 2023(08): 69-74 . 百度学术
2. ZHAO YuZhe,MA TianBao,LI Jian,NING JianGuo. Study on dynamic response of multi-body structure under explosive driving. Science China(Technological Sciences). 2016(09): 1360-1369 . 必应学术
3. 李党娟,刘唐唐,孙浩,刘群华. 一种高精度重锤落体测速装置. 国外电子测量技术. 2015(10): 84-87 . 百度学术
其他类型引用(2)
-