• ISSN 1001-1455  CN 51-1148/O3
  • EI Compendex、CA收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

含双侧分支受限空间油气爆炸火焰行为与超压特性大涡模拟

刘冲 杜扬 梁建军 张培理 孟红

引用本文:
Citation:

含双侧分支受限空间油气爆炸火焰行为与超压特性大涡模拟

    作者简介: 刘 冲(1989- ),男,博士研究生,generalsir2013@163.com;
    通讯作者: 张培理, zpl612323@163.com
  • 中图分类号: O354;X932

Large eddy simulation of gasoline/air mixture explosion in a semi-confined space with bilateral branches

    Corresponding author: ZHANG Peili, zpl612323@163.com ;
  • CLC number: O354;X932

  • 摘要: 为研究含分支结构狭长受限空间油气爆炸特性规律,基于大涡模拟WALE模型和Zimont预混火焰模型,对横截面为100 mm×100 mm的含双侧分支管道受限空间油气泄压爆炸特性进行了数值模拟。通过对火焰形态、火焰传播速度和动态超压3个物理量的对比,验证了所建立模型对于含分支结构受限空间油气爆炸计算的适用性。基于数值模拟结果,对爆炸过程中的流场结构、火焰形态和超压变化规律进行了分析,指出了“浪花状”火焰的形成原因。结果表明:(1)火焰传播进入分支管道前,在主管道和分支管道交界处会产生旋转方向相反的对称涡旋结构,并随着火焰传播不断向分支管道内部发展;(2)当火焰传播进入分支管道后,分支管道内部前期已建立流场决定了火焰的形态,火焰锋面在涡旋结构作用下呈“浪花状”,此后火焰和流场相互影响,流场向湍流转捩,火焰锋面褶皱变形;(3)爆炸超压升压过程可划分为4个阶段,受到火焰锋面面积和分支管道泄压共同作用,表明爆炸流场、火焰行为和动态超压呈现出显著耦合性。
  • 图 1  实验系统

    Figure 1.  Experimental system

    图 2  实验管道

    Figure 2.  Experimental pipe

    图 3  网格划分

    Figure 3.  Computational grids

    图 4  实验和数值模拟的火焰形态

    Figure 4.  Experimental and simulated flame structures

    图 5  实验和数值模拟的火焰位置

    Figure 5.  Experimental and simulated flame front locations

    图 6  实验和数值模拟的火焰速度

    Figure 6.  Experimental and simulated flame speeds

    图 7  P1处实验和数值模拟的爆炸超压曲线

    Figure 7.  Experimental and simulated overpressure curvesat point P1

    图 8  P1和P2处数值模拟的爆炸超压曲线

    Figure 8.  Simulated overpressure curves at points P1 and P2

    图 9  火焰传播形态和流场速度矢量

    Figure 9.  Flame propagations and velocity vectors of flow field

    图 10  数值模拟的超压曲线和火焰传播形态

    Figure 10.  Simulated overpressure curve and flame propagation structure

    表 1  实验和数值模拟的典型超压峰值

    Table 1.  Experimental and simulated typical overpressure peaks

    方法 p1/kPa ε1/% p2/kPa ε2/% pmax/kPa εmax/%
    实验 12 0 25 0 58 0
    数值模拟 27 8.0 67 15.5
    下载: 导出CSV
  • [1] SUN S C, LIU G, LIU J X, et al. Effect of porosity and element thickness on flame quenching for in-line crimped-ribbon flame arresters [J]. Journal of Loss Prevention in the Process Industries, 2017, 50: 221–228. DOI: 10.1016/j.jlp.2017.09.017.
    [2] 王波, 杜扬, 齐圣, 等. 油气爆炸在细长密闭管道内的振荡传播特性 [J]. 振动与冲击, 2017, 36(17): 97–103, 126. DOI: 10.13465/j.cnki.jvs.2017.17.016.
    WANG B, DU Y, QI S, et al. Oscillation propagation characteristics of gasoline-air mixture explosion in elongated closed tubes [J]. Journal of Vibration and Shock, 2017, 36(17): 97–103, 126. DOI: 10.13465/j.cnki.jvs.2017.17.016.
    [3] 杜扬, 王世茂, 齐圣, 等. 油气在顶部含弱约束结构受限空间内的爆炸特性 [J]. 爆炸与冲击, 2017, 37(1): 53–60. DOI: 10.11883/1001-1455(2017)01-0053-08.
    DU Y, WANG S M, QI S, et al. Explosion of gasoline/air mixture in confined space with weakly constrained structure at the top [J]. Explosion and Shock Waves, 2017, 37(1): 53–60. DOI: 10.11883/1001-1455(2017)01-0053-08.
    [4] QI S, DU Y, WANG S M, et al. The effect of vent size and concentration in vented gasoline-air explosions [J]. Journal of Loss Prevention in the Process Industries, 2016, 44: 88–94. DOI: 10.1016/j.jlp.2016.08.005.
    [5] YU M G, WAN S J, ZHENG K, et al. Influence on the methane/air explosion characteristics of the side venting position in a pipeline [J]. Process Safety and Environmental Protection, 2017, 111: 292–299. DOI: 10.1016/j.psep.2017.07.017.
    [6] 李国庆, 杜扬, 齐圣, 等. 连续圆孔障碍物对油气泄压爆炸火焰特性影响大涡模拟 [J]. 爆炸与冲击, 2018, 38(6): 1286–1394. DOI: 10.11883/bzycj-2017-0215.
    LI G Q, DU Y, QI S, et al. Large eddy simulation on the vented gasoline-air mixture explosions in a semi-confined pipe with continuous circular hollow obstacles [J]. Explosion and Shock Waves, 2018, 38(6): 1286–1394. DOI: 10.11883/bzycj-2017-0215.
    [7] LI G Q, DU Y, WANG S M, et al. Large eddy simulation and experimental study on vented gasoline-air mixture explosions in a semi-confined obstructed pipe [J]. Journal of Hazardous Materials, 2017, 339: 131–142. DOI: 10.1016/j.jhazmat.2017.06.018.
    [8] WEN X P, YU M G, JI W T, 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.
    [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] ZHU C J, LIN B Q, JIANG B Y. Flame acceleration of premixed methane/air explosion in parallel pipes [J]. Journal of Loss Prevention in the Process Industries, 2012, 25(2): 383–390. DOI: 10.1016/j.jlp.2011.10.004.
    [11] XIAO H H, HE X C, WANG Q S, et al. Experimental and numerical study of premixed flame propagation in a closed duct with a 90° curved section [J]. International Journal of Heat and Mass Transfer, 2013, 66: 818–822. DOI: 10.1016/j.ijheatmasstransfer.2013.07.091.
    [12] ZHU C J, LIN B Q, YE Q, et al. Effect of roadway turnings on gas explosion propagation characteristics in coal mines [J]. Mining Science and Technology (China), 2011, 21(3): 365–369. DOI: 10.1016/j.mstc.2011.05.006.
    [13] NIU Y H, SHI B M, JIANG B Y. Experimental study of overpressure evolution laws and flame propagation characteristics after methane explosion in transversal pipe networks [J]. Applied Thermal Engineering, 2019, 154: 18–23. DOI: 10.1016/j.applthermaleng.2019.03.059.
    [14] ZHANG P L, DU Y, ZHOU Y, et al. Explosions of gasoline-air mixture in the tunnels containing branch configuration [J]. Journal of Loss Prevention in the Process Industries, 2013, 26(6): 1279–1284. DOI: 10.1016/j.jlp.2013.07.003.
    [15] LI G Q, DU Y, QI S, 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.
    [16] 杜扬, 李国庆, 吴松林, 等. T型分支管道对油气爆炸强度的影响 [J]. 爆炸与冲击, 2015, 35(5): 729–734. DOI: 10.11883/1001-1455(2015)05-0729-06.
    DU Y, LI G Q, WU S L, 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. DOI: 10.11883/1001-1455(2015)05-0729-06.
    [17] 杜扬, 李蒙, 李国庆, 等. 含双侧分支结构受限空间油气泄压爆炸超压特性与火焰行为 [J]. 化工进展, 2018, 37(7): 2557–2564. DOI: 10.16085/j.issn.1000-6613.2017-2522.
    DU Y, LI M, LI G Q, et al. Effects of bilateral branches structure on characteristics of gasoline-air mixtures explosion overpressure and flame behavior in a semi-confined space [J]. Chemical Industry and Engineering Progress, 2018, 37(7): 2557–2564. DOI: 10.16085/j.issn.1000-6613.2017-2522.
    [18] 刘冲, 杜扬, 李国庆, 等. 狭长密闭空间内油气爆炸火焰特性大涡模拟 [J]. 化工学报, 2018, 69(12): 5348–5358. DOI: 10.11949/j.issn.0438-1157.20180614.
    LIU C, DU Y, LI G Q, et al. Large eddy simulation of gasoline-air mixture explosion in closed narrow-long space [J]. CIESC Journal, 2018, 69(12): 5348–5358. DOI: 10.11949/j.issn.0438-1157.20180614.
    [19] 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.
    [20] ZIMONT V L, BATTAGLIA V. Joint RANS/LES approach to premixed flamemodelling in the context of the TFC combustion model [J]. Flow, Turbulence and Combustion, 2006, 77(1−4): 305–331. DOI: 10.1007/s10494-006-9048-0.
    [21] 温小萍. 瓦斯湍流爆燃火焰特性与多孔介质淬熄抑爆机理研究[D]. 大连: 大连理工大学, 2014: 63−65.
    [22] ZHENG K, YU M G, LIANG Y P, et al. Large eddy simulation of premixed hydrogen/methane/air flame propagation in a closed duct [J]. International Journal of Hydrogen Energy, 2018, 43(7): 3871–3884. DOI: 10.1016/j.ijhydene.2018.01.045.
  • [1] 张培理杜扬 . 油气爆炸的氮气非预混抑制实验. 爆炸与冲击, 2016, 36(3): 347-352. doi: 10.11883/1001-1455(2016)03-0347-06
    [2] 杜扬王世茂齐圣王波李阳超李国庆 . 油气在顶部含弱约束结构受限空间内的爆炸特性. 爆炸与冲击, 2017, 37(1): 53-60. doi: 10.11883/1001-1455(2017)01-0053-08
    [3] 张培理杜扬 . 油气爆炸过程火焰燃烧模式的实验估计. 爆炸与冲击, 2016, 36(5): 688-694. doi: 10.11883/1001-1455(2016)05-0688-07
    [4] 李国庆杜扬齐圣王世茂李蒙李润 . 连续圆孔障碍物对油气泄压爆炸火焰特性影响大涡模拟. 爆炸与冲击, 2018, 38(6): 1286-1394. doi: 10.11883/bzycj-2017-0215
    [5] 杨艺何学秋刘建章王从银 . 瓦斯爆燃火焰内部流场分形特性研究. 爆炸与冲击, 2004, 24(1): 30-36.
    [6] 李艳超毕明树高伟 . 耦合火焰不稳定的爆炸超压预测. 爆炸与冲击, 2020, 40(1): 012101-1-012101-6. doi: 10.11883/bzycj-2019-0004
    [7] 杜扬李国庆李阳超齐圣王世茂王波 . T型分支管道对油气爆炸压力的影响. 爆炸与冲击, 2017, 37(2): 323-331. doi: 10.11883/1001-1455(2017)02-0323-09
    [8] 杜扬李国庆吴松林张培理周毅齐圣王世茂 . T型分支管道对油气爆炸强度的影响. 爆炸与冲击, 2015, 35(5): 729-734. doi: 10.11883/1001-1455(2015)05-0729-06
    [9] 张陶惠君明解立峰郭学永於津 . FAE爆炸场超压与威力的实验研究. 爆炸与冲击, 2004, 24(2): 176-181.
    [10] 喻健良闫兴清 . 硅酸铝棉对火焰速度和爆炸超压的抑制作用. 爆炸与冲击, 2013, 33(4): 363-368. doi: 10.11883/1001-1455(2013)04-0363-06
    [11] 李润之司荣军 . 低温环境下甲烷爆炸流场特性模拟. 爆炸与冲击, 2015, 35(6): 901-906. doi: 10.11883/1001-1455(2015)06-0901-06
    [12] 高伟阿部俊太郎荣建忠土桥律 . 气流特征对水平长管内石松子粉尘爆炸火焰结构的影响. 爆炸与冲击, 2015, 35(3): 372-379. doi: 10.11883/1001-1455-(2015)03-0372-08
    [13] 陆晓霞李磊赵守田严晓芳王金生韩超 . 壳体约束对液体爆炸抛撒流场特性的影响. 爆炸与冲击, 2016, 36(6): 803-810. doi: 10.11883/1001-1455(2016)06-0803-08
    [14] 姚哲芳任辉启沈兆武 . 柱状装药爆炸条件下厚壁圆筒爆室内 流场及结构动力响应分析. 爆炸与冲击, 2012, 32(5): 449-456. doi: 10.11883/1001-1455(2012)05-0449-08
    [15] 浦锡锋王仲琦白春华周刚 . 用于爆炸流场与结构间相互作用分析的Euler-Lagrange耦合模拟技术. 爆炸与冲击, 2011, 31(1): 6-10. doi: 10.11883/1001-1455(2011)01-0006-05
    [16] 齐圣杜扬梁建军张培理 . 受限空间油气爆燃火焰形态. 爆炸与冲击, 2016, 36(6): 832-838. doi: 10.11883/1001-1455(2016)06-0832-07
    [17] 杜扬王世茂袁广强齐圣王波李国庆李阳超 . 含弱约束端面短管道油气爆炸特性实验研究. 爆炸与冲击, 2018, 38(2): 465-472. doi: 10.11883/bzycj-2015-0242
    [18] 陈二云赵改平杨爱玲 . 环形激波聚焦流场特性的数值研究. 爆炸与冲击, 2012, 32(3): 291-296. doi: 10.11883/1001-1455(2012)03-0291-06
    [19] 饶国宁周莉宋述忠解立峰李斌彭金华 . 云爆药剂爆炸超压测试及威力评价. 爆炸与冲击, 2018, 38(3): 579-585. doi: 10.11883/bzycj-2016-0245
    [20] 秦毅陈小伟黄维 . 密闭空间可燃气体爆炸超压预测. 爆炸与冲击, 2020, 40(3): 032202-1-032202-12. doi: 10.11883/bzycj-2019-0175
  • 加载中
图(10)表(1)
计量
  • 文章访问数:  1468
  • HTML全文浏览量:  841
  • PDF下载量:  6
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-10-23
  • 录用日期:  2020-03-25
  • 网络出版日期:  2020-06-02
  • 刊出日期:  2020-06-01

含双侧分支受限空间油气爆炸火焰行为与超压特性大涡模拟

    作者简介:刘 冲(1989- ),男,博士研究生,generalsir2013@163.com
    通讯作者: 张培理, zpl612323@163.com
  • 1. 陆军勤务学院油料系,重庆 401311
  • 2. 陆军勤务学院勤务指挥系,重庆 401311

摘要: 为研究含分支结构狭长受限空间油气爆炸特性规律,基于大涡模拟WALE模型和Zimont预混火焰模型,对横截面为100 mm×100 mm的含双侧分支管道受限空间油气泄压爆炸特性进行了数值模拟。通过对火焰形态、火焰传播速度和动态超压3个物理量的对比,验证了所建立模型对于含分支结构受限空间油气爆炸计算的适用性。基于数值模拟结果,对爆炸过程中的流场结构、火焰形态和超压变化规律进行了分析,指出了“浪花状”火焰的形成原因。结果表明:(1)火焰传播进入分支管道前,在主管道和分支管道交界处会产生旋转方向相反的对称涡旋结构,并随着火焰传播不断向分支管道内部发展;(2)当火焰传播进入分支管道后,分支管道内部前期已建立流场决定了火焰的形态,火焰锋面在涡旋结构作用下呈“浪花状”,此后火焰和流场相互影响,流场向湍流转捩,火焰锋面褶皱变形;(3)爆炸超压升压过程可划分为4个阶段,受到火焰锋面面积和分支管道泄压共同作用,表明爆炸流场、火焰行为和动态超压呈现出显著耦合性。

English Abstract

  • 汽油作为一种工业制成品,广泛应用于交通运输和工业清洗行业。但是,作为一种易燃易爆易挥发的危险品,在储存、运输及加注的全过程中,汽油时刻有着安全威胁。为有效应对安全威胁,降低人类的生命和财产损失,非常有必要针对汽油蒸气/空气混合气的爆炸特征开展深入研究。此前,已有大量针对不同可燃气体爆炸特征的研究,发现受限空间结构和形状特征对可燃气爆炸超压和火焰的变化规律有较大影响。主要影响因素包括为狭长受限空间长径比[1-2]、半受限空间开口尺寸及位置[3-5]、受限空间障碍物形状及尺寸[6-8]和所含弯管及分支结构特征[9-11]等。特别地,弯管及分支结构是储油洞库坑道及输油管路中常见的结构形式,研究表明弯管及分支结构的存在对爆炸超压和火焰传播速度有较强的激励作用,并在一定程度上可与挡板式障碍物类比[9]。Zhu等[12]发现含弯管结构受限空间可燃气爆炸规律主要受弯管结构数量影响,而受弯管方向影响较小。Niu等[13]对甲烷在输气管网中的爆炸特征进行了研究,发现爆炸超压峰值及火焰传播速度受连接平行管网的分支长度影响,两分支长度不同时爆炸更为剧烈复杂。

    尽管弯管和分支结构在形状特征上有一定相似,但由于分支结构具有更复杂的通道和多个流动终点[14],因此确有必要对含多分支结构受限空间油气爆炸特征进行研究。Zhang等[14]对5种不同分支结构受限空间爆炸进行研究,发现油气爆燃指数的增加与分支数量呈函数关系。Li等[15]、杜扬等[16]通过实验研究,发现分支结构可以增大汽油蒸气/空气混合气爆炸超压峰值、升压速率及火焰传播速度。杜扬等[17]对含双侧分支受限空间油气爆炸进行研究,发现火焰在含有双侧分支结构的管道内会呈现“浪花状”火焰。目前,人们已经对含分支结构受限空间油气爆炸进行了一定的研究,但研究多集中于实验分析和现象描述,对爆炸流场和火焰的耦合关系仍需深入研究。

    本文中,在已有实验研究[17]的基础上,通过大涡模拟方法对含双侧分支受限空间油气爆炸过程进行模拟,并将火焰形态、传播速度和爆炸动态超压等模拟结果与实验进行验证,分析含双侧分支受限空间油气爆炸流场结特征、火焰传播行为和动态超压特性,指出“浪花状”火焰形成的内在原因,揭示流场、火焰和超压的耦合机理。

    • 大涡模拟通过对小尺度湍流脉动进行过滤,求解计算量比直接模拟方法小,但空间分辨率比雷诺平均方法高,适合于捕捉油气爆炸的湍流特征[18]

      大涡模拟对小尺度湍流脉动进行Favre过滤后,利用N-S方程对大尺度湍流脉动进行直接求解。经Favre过滤后,大涡模拟控制方程的连续性方程、动量方程、能量方程、反应进程变量及其守恒方程为:

      $ \frac{{\partial \rho }}{{\partial t}} + \frac{{\partial (\rho {{\tilde u}_j})}}{{\partial {x_j}}} = 0 $

      $ \frac{{\partial (\rho {{\tilde u}_j})}}{{\partial t}} + \frac{{\partial (\rho {{\tilde u}_i}{{\tilde u}_j})}}{{\partial {x_j}}} = \frac{{\partial {\sigma _{ij}}}}{{\partial {x_i}}} - \frac{{\partial \bar p}}{{\partial {x_i}}} - \frac{{\partial {\tau _{ij}}}}{{\partial {x_j}}} $

      $ \frac{{\partial (\rho {{\tilde h}_{\rm{s}}})}}{{\partial t}} + \frac{{\partial (\rho {{\tilde u}_i}{{\tilde h}_{\rm{s}}})}}{{\partial {x_i}}} - \frac{{\partial \bar p}}{{\partial t}} - {\bar u_j}\frac{{\partial \bar p}}{{\partial {x_i}}} - \frac{\partial }{{\partial {x_i}}}\left( {\lambda \frac{{\partial \tilde T}}{{\partial {x_i}}}} \right) = - \frac{\partial }{{\partial {x_j}}}\left[ {\rho \left( {\widetilde {{u_i}{h_{\rm{s}}}} - {{\tilde u}_i}{{\tilde h}_{\rm{s}}}} \right)} \right] $

      $ c = 1 - \frac{{{Y_{\rm{f}}}}}{{Y_{\rm{\!f}}^{\rm{0}}}} $

      $ \frac{{\partial (\rho \tilde c)}}{{\partial t}} + \nabla \cdot \left( {\rho {{\tilde u}}\tilde c} \right) + \nabla \cdot \left[ {\rho \left( {{{\widetilde {{ u}c}}} - {{\tilde u}}\tilde c} \right)} \right]{\rm{ = }}\overline {\nabla \cdot \left( {\rho D\nabla \tilde c} \right)} + \overline {{{\dot \omega }_{\rm{c}}}} $

      式中:ρ为密度,t为时间,uiuj为速度分量,xixj为方向分量,τij为黏性应力张量分量,p为压力,hs为显焓,λ为热导率,T为温度,D为扩散系数,${\dot \omega _{\rm{c}}}$为归一化化学反应速率。Yf为局部未燃燃料质量分数,Yf0为初始未燃混合气中燃料质量分数;c为反应进程变量,c=0代表化学反应未进行,c=1代表化学反应完全;符号~表示Favre过滤量,—表示物理空间过滤量。

      为使大涡模拟控制方程封闭,需建立亚网格模型。壁面自适应局部涡黏模型(WALE)能够较好地捕捉流场从层流到湍流的转变,且计算量较小[19]。因此,采用WALE模型作为计算的亚网格模型。涡黏模型方程为:

      $ {\tau _{ij}} - \frac{1}{3}{\tau _{kk}}{\delta _{ij}} = - 2{\mu _{{\rm{S}}}}{\bar S_{ij}} $

      式中:δij为单位张量分量;μS为亚网格黏性系数,为:

      $ {\mu _{{\rm{S}}}} = \rho L_{\rm{S}}^2\frac{{{{\left( {{S_{ij}}{S_{ij}}} \right)}^{3/2}}}}{{{{\left( {{{\bar S}_{ij}}{{\bar S}_{ij}}} \right)}^{3/2}} + {{\left( {{S_{ij}}{S_{ij}}} \right)}^{5/4}}}} $

      式中:LS为亚网格尺度混合长度,${S_{ij}}$为旋率张量分量。

    • 湍流燃烧过程中,火焰锋面厚度较小,为取得精确的结果,必须采用更精细的网格。但这样会使网格量增大,运算成本陡增。因此,采用Zimont模型[20],对火焰锋面进行增厚处理,确保在有限运算量下精确模拟火焰锋面。在Zimont模型中,湍流火焰速度公式为:

      $ {u_{\rm{t}}} = {A^ * }{\left( {{{u'_\Delta} }} \right)^{3/4}}S_{\rm{L}}^{1/2}{\chi ^{ - 1/4}}L_\Delta ^{1/4} $

      式中:A*为模型常数,计算中取值0.5[20]${u'_\Delta }$为亚网格脉动速度;SL为层流火焰速度,计算中取0.43 m/s[6]$\chi $为未燃混合物的摩尔传热系数;$L_\Delta ^{}$为亚网格湍流长度尺度。

    • 以文献[6]中的实验数据为基础开展数值模拟。其实验系统如图1所示,由实验管道、点火系统、压力采集系统、高速摄影系统、油气产生系统、油气浓度测试系统、同步触发系统和数据采集系统等部分组成。图2为实验管道尺寸,数值模拟几何模型与实验管道尺寸一致。主管道及两侧分支管道内部横截面尺寸均为100 mm×100 mm,主管道左端面密闭,中心处为点火点,右端面开口率为50%。实验前,用聚乙烯薄膜密封右端开口,但爆炸破膜压力很小,对爆炸发展演变影响较小,因此在数值计算过程中未考虑薄膜的影响[8]。为保证计算结果的准确性,数值计算中在管道右端开口处设置了一个700 mm×500 mm×800 mm的外部空间模拟实验外场,以抑制爆炸回流对管道内部的影响,如图2所示。沿爆炸传播方向,在分支前、后各275 mm的主管道处分别设置两个压阻型传感器P1、P2,对应数值计算相同位置两个压力监测点。实验时,管道充装汽油体积分数1.73%的汽油-空气混合气,环境压力约为1.013×105 Pa,环境温度约为300 K。高速摄影仪记录爆炸过程火焰发展过程,帧率为500 Hz。

      图  1  实验系统

      Figure 1.  Experimental system

      图  2  实验管道

      Figure 2.  Experimental pipe

    • 大涡模拟计算过程中,对较大尺寸涡团直接求解、对较小尺寸涡团过滤的特性决定了利用大涡模拟进行数值计算没有“网格无关性”,越小的网格尺寸和时间步长会获得越精确的计算结果,这与雷诺平均方法不同。但随着网格尺寸趋近更小,计算量将趋近于直接模拟方法,因此在网格划分上要尽可能平衡计算精度和计算成本[21]。经过调试验证后,本文中网格划分如图3所示,图3(a)为计算域整体网格划分,图3(b)为分支处局部加密网格。因计算区域规则,内场和外场网格均为六面体结构网格。在内场分支结构区域加密网格,最小为2 mm;外场沿爆炸传播方向网格逐渐增大,最大为10 mm:计算域网格总数约为490万个。

      图  3  网格划分

      Figure 3.  Computational grids

    • 为与实验条件相对应,对含双侧分支结构受限空间充装汽油体积分数1.73%的汽油-空气混合气的爆炸过程进行模拟。主管道和分支管道壁面设置为绝热无滑移壁面且温度为300 K,外场空间壁面设置为压力出口且压力为1.013×105 Pa。环境温度为300 K,流场初始速度和反应进程变量为零。在实验管道左端面中心处设置一个半径5 mm的半球形区域,将其内部反应进程变量设为1,用以模拟实验点火过程。

    • 计算过程中,采用二阶迎风格式对对流项进行离散,采用二阶中心差分格式对扩散项进行离散,采用SIMPLE算法求解压力和速度耦合方程组。时间步长为10−4 s,每时间步迭代20次,保证能量残差小于10−6,动量残差小于10−5,质量残差小于10−6,反应进程变量残差小于10−4。所有计算在Think Station P700工作站完成,耗时约100 h。

    • 油气爆炸关键特征通常包括火焰传播行为和超压变化规律,为验证数值模拟结果的准确性,选取火焰传播形态、火焰锋面传播速度和超压曲线三方面进行数值模拟,并和文献[17]实验数据对比。

    • 图4为火焰形态的实验和数值模拟结果。在初始阶段(32 ms前),火焰受到管道侧壁的约束沿轴向拉伸,以典型的“指形”火焰传播,且此时火焰锋面未受到分支结构的扰动作用,火焰锋面包络面清晰光滑,为层流状。当火焰传播至分支管道附近时,因壁面约束消失和分支管道尖角的障碍物作用,致使流场波动剧烈,火焰从层流状向湍流状转捩,火焰锋面褶皱弯曲。火焰传播至分支管道(54~64 ms)时,火焰锋面形成了三路传播路径,即火焰继续沿主管道向末端传播,同时在两侧分支管道中对称传播。在分支管道中,火焰锋面因受到右侧壁面的约束,在传播中逐渐向左侧壁面弯曲,并形成了实验中对称的“浪花状”火焰。并且,随着反应的不断进行,狭长火焰锋面两侧大量未燃气被同时引燃,火焰锋面严重变形并形成明显褶皱。

      图  4  实验和数值模拟的火焰形态

      Figure 4.  Experimental and simulated flame structures

      图4可以看出,火焰在主管道中传播时,数值模拟能够直观准确地反映实验中的火焰形态变化。当火焰传播进入分支管道时,数值模拟能够较准确地反映实验中“浪花状”火焰的轮廓和锋面位置,但在细节上和实验仍存在一定差别。究其原因有两点:一方面,高速摄影仪在拍摄较大尺寸物体时像素数会呈指数降低,实验照片对火焰的细节捕捉不够;另一方面,数值模拟采用了Zimont模型计算火焰传播锋面,而该模型为了更好地捕捉到火焰锋面位置,对火焰锋面进行了增厚处理,导致在60和62 ms时对分支管道中火焰锋面末端弯曲细节表达不精准。鉴于本文的研究目的,是为了分析火焰与流场耦合规律,因此可以忽略。

      在火焰形态对比基础上,对火焰位置和传播速度进行了对比,如图56所示。由图5可以发现,数值模拟的火焰位置能够较好地反应火焰传播位置随时间的变化,最大误差为14.436%,发生在42 ms处。这是因为,此时火焰锋面传播至分支管道附近,爆炸流场复杂,火焰锋面结构出现严重溃散,扰动随机性较大,导致误差增大。由图6可以发现,数值模拟的火焰速度与实验的变化规律整体一致,在爆炸初期(28 ms前),因参与反应的油气较少,反应速率较小,火焰传播速度缓慢升高;随着化学反应持续进行(实验28~40 ms,数值模拟28~45 ms),火焰锋面向末端传播过程中沿轴向不断拉伸,火焰表面积增大,参与反应的单位时间油气增多,化学反应速率增大,进而火焰传播速度迅速增大;当火焰锋面传播到分支管道附近时,初期分支管道的泄压作用较大,火焰传播速度增长趋势出现停滞,其后随着分支管道内油气被引燃,火焰锋面表面积呈指数增加,火焰传播速度继续增大。数值模拟误差最大为24%,发生在45 ms处。

      图  5  实验和数值模拟的火焰位置

      Figure 5.  Experimental and simulated flame front locations

      图  6  实验和数值模拟的火焰速度

      Figure 6.  Experimental and simulated flame speeds

    • 在压力传感器P1处实验和数值模拟的动态超压曲线,如图7所示。表1为实验和数值模拟的3个典型超压峰值的对比。可见,两者超压曲线的规律基本一致,分别在65和68 ms达到爆炸超压峰值,峰值误差为8.6%。但在数值模拟中,未考虑薄膜破裂时的破膜压力对管道内压力场的影响,因此在爆炸超压上升阶段,实验存在3个峰值,而数值模拟只存在2个峰值。

      方法 p1/kPa ε1/% p2/kPa ε2/% pmax/kPa εmax/%
      实验 12 0 25 0 58 0
      数值模拟 27 8.0 67 15.5

      表 1  实验和数值模拟的典型超压峰值

      Table 1.  Experimental and simulated typical overpressure peaks

      图  7  P1处实验和数值模拟的爆炸超压曲线

      Figure 7.  Experimental and simulated overpressure curvesat point P1

      在不考虑破膜压力的情况下,爆炸超压升压过程可以划分为4个阶段,以数值模拟得到的曲线为例:第Ⅰ阶段为缓慢升压阶段(0~20 ms),该阶段从模拟点火开始至火焰拉伸呈“指形”终止,在该阶段爆炸集中于点火源附近,参与反应未燃气量少,反应释能小,超压随之上升缓慢;第Ⅱ阶段为超压上升阶段(20~50 ms),随着爆炸进行,参与反应的未燃气逐渐增多,单位时间释能迅速增大,超压开始明显上升;第Ⅲ阶段为超压平衡阶段(50~56 ms),主管道泄流效应和分支管道二次泄压作用弱化了内场超压上升,超压上升速率减慢;第Ⅳ阶段为快速升压阶段(56~65 ms),分支管道增大了燃烧面积,反应速率迅速上升,超压再一次快速上升,直至可燃气消耗殆尽转入下降阶段。

      在小尺度狭长空间发生油气泄爆时,通常管道内压力波动以静压为主,动压对超压的影响非常小,因此空间内部各点的超压基本一致。而对于含双侧分支的受限空间,当压力波传播经过分支附近时,内部空间横截面的变化必然影响压力波的传播,因此在分支管道前后测点的超压一定会有所差异。为此,对比了分支管道前后测点P1和P2处的超压曲线,如图8所示。由图8可以看出,P1和P2处的超压曲线演变规律基本一致,仅在局部有数值大小等方面的差异,具体有3点:(1)在缓慢升压阶段和超压上升阶段,P2处超压上升滞后于P1;(2)P1处最大超压峰值大于P2,形成时刻略滞后于P2;(3)P2处的超压振荡比P1剧烈。P1和P2处的超压特性差异说明,分支管道导致主管道横截面突扩,在分支管道附近诱导稀疏波导致压力波衰减,分支管道对主管道起到了泄压作用,并且P2位置更靠近泄爆口,受泄流效应的影响也更大。

      图  8  P1和P2处数值模拟的爆炸超压曲线

      Figure 8.  Simulated overpressure curves at points P1 and P2

    • 由文献[22],在爆炸发展过程中,流场变化会影响火焰传播行为,火焰传播行为的改变亦会对爆炸流场产生影响,火焰传播和湍流流场间存在一定的正反馈耦合机制。此外,火焰结构在一定程度上表征反应速率的快慢,进而影响爆炸超压的波动。

      图9为火焰传播形态和流场速度矢量图,白色箭头代表速度矢量,颜色越白速度数值越大,速度矢量所在平面为管道中轴面(y=0.05 m)。由于重点关注分支管道附近的火焰和流场形态,因此火焰传播至分支管道前(20、32 ms)的图片包含点火端至分支管道;而为了清晰显示火焰在分支管道中传播形成“浪花状”的过程,这个阶段(48、60、64 ms)的图片仅包含分支管道区域。

      图  9  火焰传播形态和流场速度矢量

      Figure 9.  Flame propagations and velocity vectors of flow field

      在20 ms时,火焰以层流状态传播,火焰锋面前端附近速度矢量沿半球面指向未燃区,主管道未燃区气体在火焰膨胀作用下沿水平方向向管道末端流动。当流动经过分支管道时,两侧壁面对未燃气的约束消失,未燃气对称地向两侧分支管道膨胀,在主管道与分支管道交界处建立了层流流场,此时分支管道内部流场速度趋近于零。随着反应进行,在32 ms时火焰以“指形”传播,火焰面积增大,主管道轴向速度场迅速增强,由于受到分支管道尖角的扰动流场在分支管道口附近向湍流转捩,且因分支管道右侧壁面的约束,在主管道与分支管道交界处形成了对称的涡旋结构,上侧旋转方向为逆时针,下侧为顺时针,此外分支管道内部流场速度逐渐增大。当反应进行到48 ms、火焰传播至分支管道口时,火焰锋面在复杂湍流流场作用下显著变形,而早期建立的对称涡旋结构在火焰冲击作用下向分支管道内部发展且旋转方向维持不变,但涡旋结构沿分支管道轴向拉伸成为较大的椭圆形,此时火焰面积较小且受到较强的泄压作用,导致火焰传播速度上升趋势减缓。此后,火焰传播进入分支管道,在漩涡结构作用下火焰锋面沿分支管道右侧壁面拉伸,而后向左侧壁面弯曲形成60 ms时的“浪花状”火焰。与此同时,火焰锋面前侧未燃气不断被引燃,火焰膨胀促进分支管道内部速度矢量指向主管道,分支管道内未燃气在此流场作用下进入主管道,导致主管道与分支管道交界处形成新的燃烧面,如64 ms时。

      为了便于分析火焰形态和动态超压的耦合关系,将P1处数值模拟的动态超压曲线和典型时刻火焰形态结构一起呈现,如图10所示。根据前述分析,爆炸超压变化过程可划分为4个升压阶段和1个振荡下降阶段。第Ⅰ阶段,在点火初期,只有点火头周围少量油气参与化学反应,而后火焰以半球形传播,参与反应的可燃气有限,爆炸缓慢发展,超压上升趋势不明显。第Ⅱ阶段,爆炸进入快速升压阶段,初期火焰以“指形”结构传播(32 ms),火焰面积随着锋面拉伸不断增大,参与反应的油气增多,超压逐渐上升,当火焰传播至分支管道附近时,因流场复杂导致的火焰锋面褶皱严重进一步导致反应更加剧烈,超压上升趋势更明显。第Ⅲ阶段,火焰传播进入分支管道形成“伪受限空间”燃烧(54、56 ms),初始火焰传播速度较慢,火焰面积变化不明显,参与反应的油气量没有显著增加,而火焰锋面靠近末端泄爆口,主管道泄压作用增强,导致这个阶段升压速率下降。第Ⅳ阶段,当火焰在分支管道中涡旋结构的作用下充分发展时,火焰锋面不规则溃散严重,同时参与反应的油气量也大幅增多,爆炸剧烈进行,升压速率较前3阶段都大。振荡下降阶段,火焰传播到64 ms时,由于主管道末端的泄压效应,超压迅速下降,而后在压力波来回传播和管道内残余油气持续爆炸的两方面作用下形成超压振荡。

      图  10  数值模拟的超压曲线和火焰传播形态

      Figure 10.  Simulated overpressure curve and flame propagation structure

    • 利用大涡模拟WALE模型,模拟了双侧分支结构受限空间油气泄压爆炸的超压特性和火焰行为,与实验数据对比了火焰传播形态、火焰传播速度和动态超压曲线,验证了所选取模型对含多分支结构受限空间泄压爆炸计算的有效性。然后,运用该模型对双侧分支结构受限空间油气泄压爆炸的空间流场结构、火焰行为和超压特性进行了研究。得出以下结论:

      (1)火焰传播进入分支管道前,火焰行为会影响分支管道内部流场,在主管道和分支管道交界处会产生方向相反的对称涡旋结构,随着火焰进一步传播,涡旋结构会向分支结构内部发展,建立旋转流场;

      (2)当火焰传播进入分支管道后,分支管道内部前期已建立流场决定了火焰的形态,火焰锋面在涡旋结构作用下形成“浪花状”火焰,此后火焰和流场相互影响,流场向湍流转捩,火焰锋面溃散;

      (3)爆炸超压升压过程可划分为4个阶段,受到火焰锋面面积、主管道泄流效应和分支管道二次泄压共同作用,爆炸流场、火焰行为和动态超压具有显著耦合性。

参考文献 (22)

目录

    /

    返回文章
    返回