LES approach to premixed methane/air flame propagating in the closed duct with a square-hole obstacle
-
摘要: 为揭示置障管道内甲烷/空气预混火焰传播特性,运用高速摄影技术对甲烷/空气预混火焰的形状变化和火焰前锋的速度特性进行实验,并利用大涡模拟对管道内的流场结构进行数值分析。结果表明:置障管道内依次出现了球形火焰、指尖形火焰及“蘑菇”状火焰,且“蘑菇”状火焰出现之后,火焰开始反向传播;“蘑菇”状火焰是双涡旋结构与火焰前锋面相互作用的结果,而火焰的反向传播是由流场中出现逆流结构引起的;障碍物对火焰前锋有明显的加速作用;大涡模拟成功再现了实验中观察到的火焰形状、火焰前锋速度及流场结构,说明大涡模拟适用于置障管道内预混火焰传播特性的研究。Abstract: Aiming at revealing the characteristics of premixed methane/air flame propagating in an obstructed duct.A 4 mm thick obstacle with a square hole of 50 mm×50 mm was equipped at 210 mm from the ignition source. In the experiment, the high-speed video photography was used to study the flame shape changes and flame front speed. In the numerical simulation, the large eddy simulation (LES) was applied to investigating the flow structure. The results demonstrate that the flame-tip successively takes on a spherical, finger and mushroom-like shape, and the flame begins to propagate in the backward direction after the mushroom-like flame appears. The mushroom-like flame can be explained by the interaction of the flame with two vortexes, and the reverse flow emerged in the flow field leads to the backward motion of the premixed flame. The flame speed is accelerated significantly due to the obstacle and the flame tip speed reaches the maximum value of 17 m/s when the flame passes through the square hole of the obstacle. The flame shape changes and the flow structure observed in the experiments can be well reproduced in the numerical simulations using the LES model. It is indicated that the LES model can be used to predict the premixed flame propagating in an obstructed duct.
-
瓦斯爆炸是发生在煤矿井下最严重的灾害之一。在实际矿井巷道中不可避免地存在各种障碍物,如通风设备、支护设施及各种管道等,当火焰与障碍物相互作用时,火焰的形态、火焰传播速度、流场结构及爆炸压力等都会发生显著变化[1-4]。因此,研究火焰与障碍物的相互作用机制及由此引发的流场结构变化对于有效预防瓦斯爆炸事故的发生及降低灾害损失有着重大意义。
学者们对置障条件下的瓦斯爆炸机理、传播特性以及火焰和流场相互作用进行了深入研究[5-6]。A.R.Masri等[7]进行了障碍物孔型对火焰传播速度影响实验,指出与三角形和圆形障碍物相比,方形障碍物对火焰的加速作用更明显。R.Hall等[8]分析了障碍物数量对湍流预混火焰传播特性的影响,发现爆炸超压随障碍物数量的增加而增大,但存在上限,超过上限后,爆炸超压随障碍物数量的增加而减小。P.Chen等[9]通过改变障碍物阻塞比证明了火焰传播速度和压力增长速度随阻塞比的增大而增加,并指出火焰穿越障碍物后速度的改变是火焰和湍流相互作用的结果。障碍物能显著增强管道内流场的湍流强度,提高火焰燃烧效率[10]。C.T.Johansen等[11]利用新型的纹影技术捕捉到了障碍物后的涡旋结构,并指出火焰传播特性变化的根本原因在于火焰与涡旋的相互作用。V.D.Sarli等[12]利用粒子示踪技术得到了障碍物后的涡旋结构,实验证明了火焰与涡旋的相互作用是火焰加速及火焰形态变化的主要原因。实验手段很难得到气体爆燃过程中火焰和流场的精细结构,大涡模拟(large eddy simulation, LES)可以揭示火焰传播过程中的流动现象,有助于理解层流火焰向湍流火焰转变、燃烧不稳定性和火焰形态变化等[13]。本文中通过实验和大涡模拟对含障碍物管道内甲烷/空气预混火焰传播特性进行研究,研究置障管道内火焰传播的动力学特性,通过对比实验和模拟结果,验证大涡模拟的有效性。
1. 实验系统
图 1所示为甲烷/空气预混火焰传播实验系统示意图。该系统由定容燃烧管道、高速摄影系统、混气系统、高压点火系统和同步控制系统等组成。实验管道为一方形直管,长度为500 mm,横断面尺寸为110 mm×80 mm。障碍物(厚度4 mm)与点火点相距200 mm,方孔截面积为50 mm×50 mm,障碍物的材质为阻燃有机玻璃板,实验前后,方孔玻璃板无结构损坏和形变,障碍物的上下端与固定在管道内的凹槽相连接。为提供高速摄影所需光路,管道前后两侧由石英玻璃板制成。实验中的点火器和高速摄像系统由同步控制器进行同步控制。高速摄像机的拍摄速度设置为960 s-1。实验气体的甲烷的体积分数为7%,管道的初始温度和压力分别为298 K和0.1 MPa。可燃气被高能点火器点燃之前,静置60 s,点火电极安装在管道左端的法兰上,高能点火器点火能可调范围是3~20 J。
2. 数学模型
2.1 控制方程
随着计算机性能的提高,计算流体力学在燃烧科学领域应用日益广泛[14-15]。本文中大涡模拟所采用的控制方程为滤波后的三维瞬态质量守恒方程、动量守恒方程、能量守恒方程和反应进程变量方程[3, 16]。大涡模拟放弃对全尺度范围上涡结构瞬时运动的模拟,只将比网格尺度大的湍流运动通过瞬时N-S方程直接计算,而小尺度涡对大涡运动的影响则通过亚格子模型(sub-gridscale model)进行建模处理。其中,反应进程变量方程控制着甲烷/空气反应物的化学反应状态,通过控制反应进程变量捕捉火焰燃烧面,大涡模型的控制方程为[17]:
∂ρ∂t+∂∂xi(ρ˜ui)=0 (1) ∂∂t(ρ˜ui)+∂∂xj(ρ˜ui˜uj)=∂σij∂xj−∂ˉp∂xi−∂τij∂xj (2) ∂ρ˜hs∂t+∂ρ˜ui˜hs∂xi−∂ˉp∂t−ˉuj∂ˉp∂xi−∂∂xi(λ∂˜T∂xi)=−∂∂xj[ρ(~uihs−˜ui˜hs)] (3) 式中:LES滤波的参量标注为横线上标,质量权重滤波的参量标注为波浪线上标;ρ为密度,p为压力,ui、uj为速度分量,t为时间,应力张量σij由分子黏度μ决定,可表示为σij=μ(∂˜ui∂xj+∂˜uj∂xi)−23μ∂˜uk∂xkδij,τij为亚网格尺度应力,将其定义为τij=ρ~uiuj−ρ˜ui˜uj; hs为显焓,λ为热导率,亚网格热焓通量可通过梯度假设近似为:ρ(~uihs−˜ui˜hs)=−μSGScpPrSGS∂˜T∂xj,其中:μSGS为亚网格黏度,PrSGS为亚网格的普朗特数,cp为定压比热,T为温度。
2.2 初始条件和边界条件
点火前,管道内甲烷/空气混合气保持静止,初温和初压分别为298 K和0.1 MPa。设反应过程变量c为量纲一化的产物质量分数,相应地,c=0对应为未燃反应物,c=1对应为已燃产物。将初始化学反应过程变量在整个流体区域设置为c=0,并通过将c设置为0.1来实现点火,点火半径为10 mm。管道边界条件为无滑移、绝热边界条件。网格为结构网格,单元特征尺寸为2 mm×2 mm×2 mm。应用ANSYS Fluent软件作为计算平台,压力速度耦合采用SIMPLE算法,对流项采用二阶逆风格式,扩散项采用中心差分格式。
3. 结果与讨论
3.1 火焰前锋发展的高速摄影图片
图 2所示为置障条件下甲烷/空气预混火焰传播的高速摄影图像。t=10 ms时,由于未受到障碍物和管道侧壁影响,火焰阵面以球形向前传播;之后,管道侧壁开始阻碍火焰的自由发展,火焰形状由球形转变成指尖形,且火焰裙边在t=32 ms时接触管道侧壁,火焰燃烧表面积开始减小;t=42 ms时,障碍物开始作用于火焰前锋,平滑的弧面火焰前沿开始变尖;t=44 ms时,火焰恰好穿过障碍物的方孔;t=48 ms时,穿过方孔的火焰再次发生变形,并在t=50 ms时形成“蘑菇”状的火焰前端;50 ms之后“蘑菇”状火焰开始破碎,值得注意的是,“蘑菇”状靠近上下壁面的部分开始反向传播,并在t=53 ms时反向穿过方孔结构。从火焰形态特征而言,置障条件下管道内的预混火焰经历了层流向湍流燃烧的转变。
3.2 火焰前锋发展的大涡模拟
图 3所示为置障条件下甲烷/空气预混火焰传播的大涡模拟,其中燃烧面的反应进程变量c=0.1。
从图 3中可以发现,基于大涡模型的数值模拟成功再现了实验观察到的火焰形态及传播现象,如球形火焰、指尖形火焰、“蘑菇”状火焰及其反向流动等。这也证明了大涡模型对于模拟置障条件下甲烷/空气预混火焰传播的可行性。
3.3 流场结构
图 4所示为置障条件下甲烷/空气预混火焰传播的流场结构。数值模拟中红色代表已燃气体,蓝色代表未燃气体,白色箭头线代表速度矢量。
在t=42 ms之前,管道内的气体流动是层流状态,流线分布均匀无交叉;t=49 ms时,未燃区域出现双涡旋结构,且涡旋范围不断扩大,涡旋中心逐渐向障碍物移动;t=50 ms时,已燃区的流动方向由原来的正向流动反转成逆向;t=51 ms时,未燃区域方孔障碍物附近也出现逆流。可以推断,“蘑菇”状火焰的形成与双涡旋结构密切相关,而燃烧后期的火焰反向传播是由流场的逆流结构造成的。
3.4 火焰前锋的动力学变化
图 5所示为火焰前锋位置随时间的变化关系的数值模拟和实验结果。可见,t=42 ms时刻是火焰前锋位置随着时间变化的转折点,在42 ms之前,变化速率较为缓慢,在42 ms之后,变化速率突然增大。图 6所示为火焰前锋速度随时间的变化特性。从图 6中可以看出,在32 ms之前,火焰前锋速度随着燃烧表面积的增大而缓慢上升;当t=32 ms时,火焰裙边触壁,燃烧表面积突然减小,进而引起火焰传播速度的小幅度减小;当t=44 ms时,火焰前锋速度达到最大值17 m/s,此时火焰前锋正在穿越障碍物的方孔。
4. 结论
本文中通过实验和数值模拟研究了方孔障碍物对管道内甲烷/空气预混火焰传播特性的影响。采用高速摄影技术揭示了火焰过程中的形状变化和火焰传播速度特性,并应用大涡模拟对火焰传播行为做了进一步分析。得到结论如下:
(1) 由于未受到方孔障碍物的影响,前期火焰经历了球形火焰和指尖形火焰两个阶段,指形火焰触壁之前,火焰燃烧面积快速增大;触壁之后,火焰燃烧面积迅速减小。临近方孔时,平滑的弧面火焰前沿开始变尖,并且火焰前沿在穿过方孔障碍物后形成“蘑菇”状火焰;火焰前沿接近管道尾端时,“蘑菇”状火焰的靠近上下壁面的部分开始反向传播,并在t=53 ms时反向穿过方孔结构。
(2) 火焰前沿形状的变化与管道内的流场结构密切相关。通过分析流场结构可以发现,“蘑菇”状火焰的形成与双涡旋结构有直接联系,而燃烧后期的火焰反向传播是由流场中出现的逆流造成的。
(3) 方孔障碍物对预混火焰的加速作用显著,火焰未穿过障碍物之前,火焰裙边触壁时刻,火焰传播速度达到5.72 m/s;当火焰穿越方孔时,火焰传播速度再次达到最大值17 m/s。
(4) 大涡模拟成功再现了实验中观察到的火焰形状、火焰前锋速度及流场结构。说明大涡模拟对于模拟置障管道内甲烷/空气预混火焰传播的可行性。
-
-
[1] Dorofeev S B. Flame acceleration and explosion safety applications[J]. Proceedings of the Combustion Institute, 2011, 33(2):2161-2175. doi: 10.1016/j.proci.2010.09.008 [2] Alharbi A, Masri A R, Ibrahim S S. Turbulent premixed flames of CNG, LPG, and H2 propagating past repeated obstacles[J]. Experimental Thermal and Fluid Science, 2014, 56(7):2-8. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=cab80255e3dd2ac1e24881d2a44ea68e [3] Johansen C T, Ciccarelli G. Modeling the initial flame acceleration in an obstructed channel using large eddy simulation[J]. Journal of Loss Prevention in the Process Industries, 2013, 26(4):571-585. doi: 10.1016/j.jlp.2012.12.005 [4] 陈志华, 叶经方, 范宝春, 等.方形管内楔形障碍物对火焰结构与传播的影响[J]. 爆炸与冲击, 2006, 26(3):208-213. doi: 10.3321/j.issn:1001-1455.2006.03.003Chen Zhihua, Ye Jingfang, Fan Baochun, et al. Effects of a wedge obstacle on flame propagation and its structure[J]. Explosion and Shock Waves, 2006, 26(3):208-213. doi: 10.3321/j.issn:1001-1455.2006.03.003 [5] Ciccarelli G, Dorofeev S. Flame acceleration and transition to detonation in ducts[J]. Progress in Energy & Combustion Science, 2008, 34(4):499-550. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=718edb9b92782b6bc9b5758162d9bbdb [6] Kundu S, Zanganeh J, Moghtaderi B. A review on understanding explosions from methane-air mixture[J]. Journal of Loss Prevention in the Process Industries, 2016(40):507-523. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=bd0c31f6cbd2e040995370fc181d2079 [7] Masri A R, Ibrahim S S, Nehzat N, et al. Experimental study of premixed flame propagation over various solid obstructions[J]. Experimental Thermal and Fluid Science, 2000, 21(1/2/3):109-116. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=ff5b56c451e0a34d392f7a4e791e637b [8] Hall R, Masri A R, Yaroshchyk P, et al. Effects of position and frequency of obstacles on turbulent premixed propagating flames[J]. Combustion and Flame, 2009, 156(2):439-446. doi: 10.1016/j.combustflame.2008.08.002 [9] Chen P, Li Y C, Huang F J, et al. Experimental and LES investigation of premixed methane/air flame propagating in a chamber for three obstacle BR configurations[J]. Journal of Loss Prevention in the Process Industries, 2016, 41(5):48-54. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=a21f326d49cd196d10940be6b2b2ce86 [10] Wen X P, Yu M G, Liu Z C, 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(25):730-738. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=01e1b3eb9b77325460e27f3d572d02aa [11] Johansen C T, 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 [12] Sarli V D, Benedetto A D, Russo G, et al. Large eddy simulation and PIV measurements of unsteady flames accelerated by obstacles[J]. Flow Turbulence and Combustion, 2009, 83(2):227-250. doi: 10.1007/s10494-008-9198-3 [13] Sarli V D, Benedetto A D, Russo G. Sub-grid scale combustion models for large eddy simulation of unsteady premixed flame propagation around obstacles[J]. Journal of Hazardous Materials, 2010, 180(1/2/3):71-78. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=8e2c1f3fb6857178dce5a963fa5e1fac [14] Ibrahim S S, Gubba S R, Malalasekera W, et al. Calculations of explosion deflagration flames using a dynamic flame surface density model[J]. Combustion, Explosion, and Shock Waves, 2012, 48(4):393-405. doi: 10.1134/S0010508212040041 [15] 马秋菊, 张奇, 庞磊.巷道壁面与瓦斯爆炸相互作用的数值模拟[J]. 爆炸与冲击, 2014, 34(1):23-27. doi: 10.3969/j.issn.1001-1455.2014.01.005Ma Qiuju, Zhang Qi, Pang Lei. Numerical simulation on interaction between laneway surface and methane explosion[J]. Explosion and Shock Waves, 2014, 34(1):23-27. doi: 10.3969/j.issn.1001-1455.2014.01.005 [16] Gubba S R, Ibrahim S S, Malalasekera W, et al. Measurements and LES calculations of turbulent premixed flame propagation past repeated obstacles[J]. Combustion and Flame, 2011, 158(12):2465-2481. doi: 10.1016/j.combustflame.2011.05.008 [17] Zimont V L, Battaglia V. Joint RANS/LES approach to premixed flame modelling in the context of the TFC combustion model[J]. Flow Turbulence and Combustion, 2006, 77(1):305-331. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dc0a05d36d25fa51a6a6e96a1f6b0ad9 期刊类型引用(22)
1. 张浩,郝斌. 密闭约束空间不同预混可燃气体火焰动力学对比数值模拟研究. 化工安全与环境. 2024(06): 41-45 . 百度学术
2. 刘佳佳,张扬,张翔,聂子硕. Y型通风采煤工作面瓦斯爆炸传播规律模拟研究. 爆炸与冲击. 2023(08): 182-194 . 本站查看
3. 王皓楠,戚承志,陈昊祥,班力壬,段秋宇,马啸宇,封焱杰,罗伊. 障碍物影响下可燃气体爆炸火焰传播规律的研究进展. 市政技术. 2022(01): 139-147+190 . 百度学术
4. 马云龙. 圆孔障碍物对压力重叠影响的实验研究. 煤矿安全. 2022(07): 70-73 . 百度学术
5. 徐景德,田思雨,张延炜,张亮. 瓦斯爆炸主动式阻隔爆技术研究进展. 华北科技学院学报. 2022(04): 71-77 . 百度学术
6. 王磊. 球形障碍物对瓦斯爆燃火焰影响试验研究. 中国安全科学学报. 2021(03): 54-59 . 百度学术
7. 胡洋,吴秋遐,张延炜,秦汉圣,庞磊,梁金虎,杨雨欣,陈明虎. 激光纹影技术在煤矿瓦斯爆燃研究中的应用. 煤矿安全. 2021(05): 66-71 . 百度学术
8. 高伟,周永浩,李智勇,李艳超,姜海鹏,黄磊,张锎,毕明树. 障碍物对开敞空间贫燃天然气爆炸特性的影响研究. 安全与环境学报. 2021(04): 1539-1544 . 百度学术
9. 余明高,阳旭峰,郑凯,栾鹏鹏. 我国煤矿瓦斯爆炸抑爆减灾技术的研究进展及发展趋势. 煤炭学报. 2020(01): 168-188 . 百度学术
10. 李国庆,杜扬,白洁,武军,李孟源,吴晓澍,朱亮. 平板障碍物通道形状对油气爆炸传播特性影响. 化工学报. 2020(04): 1912-1921 . 百度学术
11. 景国勋,张胜旗,段新伟,郭绍帅,刘闯. 竖直管道内煤尘浓度对瓦斯爆炸特性影响研究. 中国安全科学学报. 2020(03): 15-20 . 百度学术
12. 程方明,常助川,史合,高彤彤,罗振敏,葛天姣. 多孔障碍物对预混火焰传播的影响. 爆炸与冲击. 2020(08): 33-41 . 本站查看
13. 贾海林,翟汝鹏,李第辉,项海军,杨永钦. 三种盐类超细水雾抑制管道内甲烷-空气预混气爆炸的差异性. 爆炸与冲击. 2020(08): 42-53 . 本站查看
14. 李国庆,张笈玮,武军,王世茂,齐晓光,张培理,齐圣. 方管内汽油-空气混合气体密闭爆炸和泄爆特性研究. 爆炸与冲击. 2020(10): 4-17 . 本站查看
15. 程方明,常助川,高彤彤,纪鹏飞,罗振敏. 不同位置多孔障碍物对预混火焰传播的影响. 中国安全科学学报. 2020(11): 114-120 . 百度学术
16. 胡洋,尹尚先,ARNTZENJ.Bj?rn,朱建芳,李雪冰. 矿井瓦斯/空气预混气体爆燃的激光纹影测试系统设计. 光学精密工程. 2019(05): 1045-1051 . 百度学术
17. 程科,叶文华,冯阳,肖磊明,陈阳明,王一茹. 破碎机锤轴转动对爆炸过程影响的数值模拟. 工程爆破. 2019(04): 20-25 . 百度学术
18. 胡洋,尹尚先,Bj?rn J.ARNTZEN,朱建芳,李雪冰,Ragnhild Dybdal OIE,秦汉圣. 瓦斯/空气预混气体爆燃流场测试系统多参数耦合同步控制实验方法. 爆炸与冲击. 2019(09): 121-127 . 本站查看
19. 景国勋,吴昱楼,郭绍帅,邵泓源,刘闯,张胜旗. 障碍物对瓦斯煤尘爆炸火焰传播规律的影响. 中国安全生产科学技术. 2019(09): 99-104 . 百度学术
20. 刘希亮,李烨,王新宇,郭佳奇. 不同截面管廊在燃气爆炸作用下的动力响应分析. 应用力学学报. 2019(05): 1111-1115+1260 . 百度学术
21. 周宁,王文秀,张国文,宗永迪,赵会军,袁雄军. 障碍物对丙烷-空气爆炸火焰加速的影响. 爆炸与冲击. 2018(05): 1106-1114 . 本站查看
22. 李国庆,杜扬,齐圣,王世茂,李蒙,李润. 连续圆孔障碍物对油气泄压爆炸火焰特性影响大涡模拟. 爆炸与冲击. 2018(06): 1286-1294 . 本站查看
其他类型引用(9)
-