Processing math: 100%
  • ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

耦合火焰不稳定的爆炸超压预测

李艳超 毕明树 高伟

李艳超, 毕明树, 高伟. 耦合火焰不稳定的爆炸超压预测[J]. 爆炸与冲击, 2020, 40(1): 012101. doi: 10.11883/bzycj-2019-0004
引用本文: 李艳超, 毕明树, 高伟. 耦合火焰不稳定的爆炸超压预测[J]. 爆炸与冲击, 2020, 40(1): 012101. doi: 10.11883/bzycj-2019-0004
LI Yanchao, BI Mingshu, GAO Wei. Explosion pressure prediction considering the flame instabilities[J]. Explosion And Shock Waves, 2020, 40(1): 012101. doi: 10.11883/bzycj-2019-0004
Citation: LI Yanchao, BI Mingshu, GAO Wei. Explosion pressure prediction considering the flame instabilities[J]. Explosion And Shock Waves, 2020, 40(1): 012101. doi: 10.11883/bzycj-2019-0004

耦合火焰不稳定的爆炸超压预测

doi: 10.11883/bzycj-2019-0004
基金项目: 国家自然科学基金(51674059)
详细信息
    作者简介:

    李艳超(1989- ),男,博士研究生,lyc092451@mail.dlut.edu.cn

    通讯作者:

    高 伟(1984- ),男,博士,教授,gaoweidlut@dlut.edu.cn

  • 中图分类号: O389;TK91

Explosion pressure prediction considering the flame instabilities

  • 摘要: 基于火焰不稳定和爆炸超压的耦合机制,通过向光滑火焰模型中引入褶皱因子,建立了褶皱火焰模型和湍流火焰模型,对密闭燃烧室内爆炸超压进行理论预测,且对比了绝热压缩和等温压缩对爆炸超压预测的影响规律。结果表明:在增强的流体动力学不稳定作用下,膨胀火焰失稳加剧,且在定容燃烧阶段形成胞状火焰;光滑火焰模型忽略了火焰不稳定,爆炸超压理论预测值比实验值偏低,且等温压缩下超压预测值低于绝热压缩下的预测值;湍流火焰模型高估了火焰褶皱程度,超压预测值远高于实验值;褶皱火焰模型可成功预测丙烷/空气爆炸压力和燃烧室体积V=25.6 m3的甲烷/空气爆炸压力;对于甲烷/空气爆炸,燃烧室体积V≤1.25 m3时,实验压力值介于褶皱火焰模型和绝热光滑火焰模型预测值之间。
  • 为缓解化石类能源引起的环境问题,碳氢类气体燃料(如天然气和液化石油气)正广泛应用于内燃机和涡轮机等供能系统。但在制备、存储、输运及使用碳氢类气体燃料过程中,不可避免地会发生燃料泄漏,进而引发燃烧和爆炸事故。由于火焰不稳定的存在,膨胀火焰表面会出现胞格结构,进而引起火焰加速并导致爆炸超压的增强[1-3]。为降低爆炸灾害后果,建立考虑火焰加速效应的爆炸超压预测模型势在必行。

    为此,国内外学者对爆炸超压预测进行了深入的研究。基于三点假设:(1) 爆炸过程中火焰表面始终光滑、(2) 火焰前端未燃气体被绝热压缩、(3) 分压升高正比于已燃气体质量分数,Dahoe等[4]建立了预测球形容器爆炸超压的光滑火焰模型:

    dpdt=3(padia,maxp0)R(pp0)1/γu[1(p0p)1/γupadia,maxppadia,maxp0]2/3SL (1)

    式中:p为预测爆炸压力,padia,max为绝热最大爆炸压力,p0为初始压力,R为密闭腔室半径,γbγu分别为已燃区域和未燃区域的比热比,SL为层流燃烧速度。

    假设火焰前端未燃气体被等温压缩,且已燃区域和未燃区域比热比相等(γu=γb=1),Bradley等[5]对微分方程(1)进行简化:

    dpdt=3SLRp2/3adia,maxpp0(padia,maxp0)1/3(1p0p)2/3 (2)

    Lautkaski[6]对偏微分方程(2)进行简化积分,推导出了受限空间内预测爆炸超压的立方根定律,但该模型仅适用于Δpp0的情况,具体的:

    Δp=(padia,maxp0)(padia,maxp0)2(SLtR)3 (3)

    值得注意的是,上述爆炸超压预测模型均假设爆炸过程中膨胀火焰始终光滑。Kuznetsov等[7]对比了常压阶段和升压阶段的富燃氢气/空气火焰形态,发现常压阶段火焰表面始终光滑,升压阶段形成了胞状火焰。Jiang等[8]实验分析了爆炸压力升高对胞状火焰结构的影响规律,指出爆炸压力升高初期,平均胞格面积快速减小,随着爆炸压力的继续升高,平均胞格面积趋于定值,且该定值与氢气掺混比无关。基于分形火焰理论,假设火焰完全湍流化,Nishimura等[3]对方程(1)进行了修正:

    dpdt=3(padia,maxp0)R2/3r1/3(pp0)1/γu[1(p0p)1/γupadia,maxppadia,maxp0]7/9SL (4)

    但该模型的预测精度强烈依赖于火焰失稳的临界半径,而且临界半径难以预先确定。

    虽然爆炸超压预测模型已有不少,但均未充分考虑火焰加速对爆炸超压的增强效应。基于此,本文揭示了定容燃烧阶段火焰失稳机理,通过引入褶皱因子,修正了光滑火焰模型,建立了耦合火焰不稳定的爆炸超压预测模型。

    图1给出了实验装置示意图。实验系统和实验方法说明详见文献[9-10],在本文仅做简单陈述。实验装置包括配气系统、高速纹影系统、压力采集系统、高压点火系统和时序控制系统。高速摄像机的拍摄频率为105 s−1。采用压电式压力传感器采集球形燃烧室内瞬态爆炸压力,压力数据使用横河数据采集仪记录,记录频率为105 s−1。高压点火器、高速摄像机、计算机均由时序控制器进行控制。

    图  1  实验装置图
    Figure  1.  Experimental apparatus

    忽略浮力不稳定,膨胀火焰失稳主要受热扩散不稳定和流体动力学不稳定控制。热扩散不稳定可用刘易斯数表征,流体动力学不稳定可用热膨胀比和火焰厚度表征。其中刘易斯数、热膨胀比和火焰厚度的计算公式如下[11-13]

    Le=DTDiM (5)
    σ=ρuρb (6)
    δ=2λρucpSL(TadT0)0.7 (7)

    式中:Le为刘易斯数,DT为热扩散系数,DiM为质量扩散系数,σ为热膨胀比,ρu为未燃气体密度,ρb为已燃产物密度,δ为火焰厚度,λ为热传导率,cp为定压比热容,SL为层流燃烧速度,Tad为绝热火焰温度,T0为燃料初始温度。

    考虑火焰前端未燃气体被绝热压缩,且假设燃烧反应瞬间达到平衡状态,进而计算整个爆炸过程中的刘易斯数、热膨胀比和火焰厚度。图2给出了甲烷/空气火焰形态特性和爆炸超压的耦合关系。本研究将爆炸压力低于初始压力10%的燃烧阶段视为定压阶段,爆炸压力高于初始压力10%的燃烧阶段视为定容阶段[14]。定压阶段,火焰表面整体光滑;定容阶段,火焰表面呈现典型的胞格结构。在整个爆炸过程中,刘易斯数变化不大,近似等于1,这意味着热扩散不稳定对火焰失稳或稳定均不起任何作用。热膨胀比和火焰厚度比值在常压阶段几乎不变,但随着爆炸压力的升高,热膨胀比和火焰厚度的比值快速升高。上述结果表明,对于当量比Φ=1.0的甲烷/空气爆炸,定容燃烧阶段火焰失稳可归因于增强的流体动力学不稳定。

    图  2  甲烷/空气火焰形态特性和爆炸超压的耦合关系(Φ=1.0)
    Figure  2.  Couplings of flame morphology and explosion pressure of stoichiometric methane/air mixture (Φ=1.0)

    图3给出了丙烷/空气火焰形态特性和爆炸超压的耦合关系。定压阶段和定容阶段的火焰形态迥异,胞状火焰仅在定容燃烧阶段形成。对于丙烷/空气火焰,尽管刘易斯数在整个爆炸过程中变化不大,但刘易斯数却始终大于1。这意味着热扩散不稳定能够对膨胀火焰起稳定效应。定压阶段,热膨胀比和火焰厚度比值变化不大,但在定容阶段,该比值迅速升高,这意味着流体动力学不稳定对火焰失稳效应逐渐增强。结合刘易斯的变化趋势可知,定容燃烧阶段,热扩散不稳定和流体动力学不稳定相互竞争,并最终导致火焰失稳。

    图  3  丙烷/空气火焰形态特性和爆炸超压的耦合关系(Φ=1.0)
    Figure  3.  Couplings of flame morphology and explosion pressure of stoichiometric propane/air mixture (Φ=1.0)

    基于2.1节可知,在增强的流体动力学不稳定作用下,胞状火焰会出现在定容燃烧阶段,这必然会增加火焰燃烧表面积,进而增强爆炸超压。因此,爆炸超压预测必须考虑定压阶段和定容阶段的火焰形态差异。通过引入褶皱因子ΞΔ,褶皱火焰模型建立如下:

    dpdt=3(padia,maxp0)R(pp0)1/γ[1(p0p)1/γpadia,maxppadia,maxp0]2/3(ΞΔSL) (8)

    式中:γ为未燃气体绝热系数,ΞΔ为褶皱因子。

    在增强的流体动力学不稳定作用下,胞状火焰呈分形特性,基于文献[15],采用等效面积和等效周长方法,可知爆炸过程中分形火焰的褶皱因子极限值是ΞΔmax=2.464 9。对于光滑火焰模型,褶皱因子始终等于1.0;对于湍流火焰模型,褶皱因子始终等于2.464 9;对于褶皱火焰模型,假设褶皱因子以指数形式从ΞΔ=1.0增长,且极限值是2.464 9,详细如下:

    ΞΔ=et,ΞΔ2.4649 (9)

    另外,本研究对比了绝热压缩和等温压缩对爆炸超压预测的影响规律。绝热压缩条件下,层流燃烧速度的计算公式如下[16-17]

    Φ=1.0的甲烷/空气:SL=0.36(TT0)1.612(pp0)0.374 (10)
    Φ=1.0的丙烷/空气:SL=0.345(TT0)1.93(pp0)0.19 (11)

    等温压缩条件下,整个爆炸过程中层流燃烧速度始终等于初始条件下层流燃烧速度,即Φ=1.0,对于甲烷/空气混合气体,SL=0.36 m/s;对于丙烷/空气混合气体,SL=0.345 m/s。

    图4给出了当量比Φ=1.0的丙烷/空气爆炸超压预测。可以发现,光滑火焰模型忽略了火焰不稳定对燃烧速率的增强作用,因此光滑火焰模型预测的爆炸压力远低于实验压力。对于光滑火焰模型,等温压缩预测的爆炸压力明显低于绝热压缩预测的爆炸压力,原因在于绝热压缩可增加层流燃烧速度。点火后,假设膨胀火焰完全湍流化,火焰褶皱程度被显著高估,进而造成理论预测压力偏高。充分考虑爆炸过程中火焰形态由光滑向胞状结构的转捩过程,假设褶皱因子以指数形式增加,褶皱火焰模型预测的爆炸压力和实验压力基本吻合。

    图  4  丙烷/空气爆炸超压预测(14 L燃烧腔室)
    Figure  4.  Explosion pressure prediction of stoichiometric propane/air mixture in a 14 L chamber

    图5给出了当量比Φ=1.0的甲烷/空气爆炸超压预测的尺度效应。燃烧腔室体积V=0.12 m3V=1.25 m3V=25.6 m3的压力数据摘自文献[18-20]。随着密闭腔室尺度的增加,湍流火焰模型均高估了爆炸超压,光滑火焰模型均低估了爆炸超压。另外,相比于丙烷/空气爆炸,甲烷/空气爆炸的绝热光滑火焰模型和等温光滑火焰模型的预测值相差不大。当燃烧腔室体积满足0.014 m3V≤1.25 m3时,实验压力处于褶皱火焰模型和绝热光滑模型计算值之间,这说明V≤1.25 m3时,褶皱火焰模型会高估火焰表面的褶皱程度。当燃烧腔室体积增大至25.6 m3时,流体动力学不稳定随着火焰尺度的增加而增强,火焰褶皱加剧,火焰传播速度加快,因此该腔室尺度下褶皱火焰模型能成功预测爆炸压力。

    图  5  甲烷/空气爆炸超压预测的尺度效应(Φ=1.0)
    Figure  5.  Scale effect of explosion pressure prediction of stoichiometric methane/air mixture (Φ=1.0)

    为预测密闭腔室内爆炸压力行为特性,本文首先揭示了火焰不稳定和爆炸超压的耦合机制,之后通过向光滑火焰模型引入褶皱因子,建立了考虑火焰不稳定的褶皱火焰模型和湍流火焰模型,最后揭示了尺度效应对爆炸超压预测的影响规律。结论如下:

    (1)随着爆炸压力的升高,流体动力学不稳定增强,火焰失稳加剧,且在定容燃烧阶段,胞状火焰形成;

    (2)相比于实验压力,光滑火焰模型忽略了火焰不稳定,进而导致理论预测值偏低,且绝热压缩下理论预测压力高于等温压缩下预测值;湍流火焰模型会高估膨胀火焰褶皱程度,理论计算值明显高于实验压力;

    (3)通过改变气体种类和密闭腔室体积,褶皱火焰模型在一定程度上可以成功再现甲烷/空气和丙烷/空气的爆炸压力。

  • 图  1  实验装置图

    Figure  1.  Experimental apparatus

    图  2  甲烷/空气火焰形态特性和爆炸超压的耦合关系(Φ=1.0)

    Figure  2.  Couplings of flame morphology and explosion pressure of stoichiometric methane/air mixture (Φ=1.0)

    图  3  丙烷/空气火焰形态特性和爆炸超压的耦合关系(Φ=1.0)

    Figure  3.  Couplings of flame morphology and explosion pressure of stoichiometric propane/air mixture (Φ=1.0)

    图  4  丙烷/空气爆炸超压预测(14 L燃烧腔室)

    Figure  4.  Explosion pressure prediction of stoichiometric propane/air mixture in a 14 L chamber

    图  5  甲烷/空气爆炸超压预测的尺度效应(Φ=1.0)

    Figure  5.  Scale effect of explosion pressure prediction of stoichiometric methane/air mixture (Φ=1.0)

  • [1] KWON O C, ROZENCHAN G, LAW C K. Cellular instabilities and self-acceleration of outwardly propagating spherical flames [J]. Proceedings of the Combustion Institute, 2002, 29(2): 1775–1783. DOI: 10.1016/S1540-7489(02)80215-2.
    [2] WU F J, JOMAAS G, LAW C K. An experimental investigation on self-acceleration of cellular spherical flames [J]. Proceedings of the Combustion Institute, 2013, 34(1): 937–945. DOI: 10.1016/j.proci.2012.05.068.
    [3] NISHIMURA I, MOGI T, DOBASHI R. Simple method for predicting pressure behavior during gas explosions in confined spaces considering flame instabilities [J]. Journal of Loss Prevention in the Process Industries, 2013, 26(2): 351–354. DOI: 10.1016/j.jlp.2011.08.009.
    [4] DAHOE A E, ZEVENBERGEN J F, LEMKOWITZ S M, et al. Dust explosions in spherical vessels: the role of flame thickness in the validity of the ‘cube-root law’ [J]. Journal of Loss Prevention in the Process Industries, 1996, 9(1): 33–44. DOI: 10.1016/0950-4230(95)00054-2.
    [5] BRADLEY D, MITCHESON A. Mathematical solutions for explosions in spherical vessels [J]. Combustion and Flame, 1976, 26: 201–217. DOI: 10.1016/0010-2180(76)90072-9.
    [6] LAUTKASKI R. Pressure rise in confined gas explosions [DB/OL]. (2005-12-30). https://www.vtt.fi/inf/julkaisut/muut/2005/PRO1_P1026_05.pdf.
    [7] KUZNETSOV M, KOBELT S, GRUNE J, et al. Flammability limits and laminar flame speed of hydrogen-air mixtures at sub-atmospheric pressures [J]. International Journal of Hydrogen Energy, 2012, 37(22): 17580–17588. DOI: 10.1016/j.ijhydene.2012.05.049.
    [8] JIANG Y H, LI G X, LI F S, et al. Experimental investigation of correlation between cellular structure of the flame front and pressure [J]. Fuel, 2017, 199: 65–75. DOI: 10.1016/j.fuel.2017.02.036.
    [9] LI Y C, BI M S, ZHANG S L, et al. Dynamic couplings of hydrogen/air flame morphology and explosion pressure evolution in the spherical chamber [J]. International Journal of Hydrogen Energy, 2018, 43(4): 2503–2513. DOI: 10.1016/j.ijhydene.2017.12.044.
    [10] LI Y C, BI M S, LI B, et al. Effects of hydrogen and initial pressure on flame characteristics and explosion pressure of methane/hydrogen fuels [J]. Fuel, 2018, 233: 269–282. DOI: 10.1016/j.fuel.2018.06.042.
    [11] LI H M, LI G X, SUN Z, et al. Experimental investigation on laminar burning velocities and flame intrinsic instabilities of lean and stoichiometric H2/CO/air mixtures at reduced, normal and elevated pressures [J]. Fuel, 2014, 135: 279–291. DOI: 10.1016/j.fuel.2014.06.074.
    [12] KIM W K, MOGI T, KUWANA K, et al. Self-similar propagation of expanding spherical flames in large scale gas explosions [J]. Proceedings of the Combustion Institute, 2015, 35(2): 2051–2058. DOI: 10.1016/j.proci.2014.08.023.
    [13] SUN Z Y, LIU F S, BAO X C, et al. Research on cellular instabilities in outwardly propagating spherical hydrogen-air flames [J]. International Journal of Hydrogen Energy, 2012, 37(9): 7889–7899. DOI: 10.1016/j.ijhydene.2012.02.011.
    [14] KATSUMI T, AIDA T, AIBA K, et al. Outward propagation velocity and acceleration characteristics in hydrogen-air deflagration [J]. International Journal of Hydrogen Energy, 2017, 42(11): 7360–7365. DOI: 10.1016/j.ijhydene.2016.06.165.
    [15] LI Y C, BI M S, GAO W. Theoretical pressure prediction of confined hydrogen explosion considering flame instabilities [J]. Journal of Loss Prevention in the Process Industries, 2019, 57: 320–326. DOI: 10.1016/j.jlp.2019.01.001.
    [16] GU X J, HAQ M Z, LAWES M, et al. Laminar burning velocity and Markstein lengths of methane-air mixtures [J]. Combustion and Flame, 2000, 121(1−2): 41–58. DOI: 10.1016/S0010-2180(99)00142-X.
    [17] METGHALCHI M, KECK J C. Laminar burning velocity of propane-air mixtures at high temperature and pressure [J]. Combustion and Flame, 1980, 38: 143–154. DOI: 10.1016/0010-2180(80)90046-2.
    [18] CASHDOLLAR K L, ZLOCHOWER I A, GREEN G M, et al. Flammability of methane, propane, and hydrogen gases [J]. Journal of Loss Prevention in the Process Industries, 2000, 13(3−5): 327–340. DOI: 10.1016/s0950-4230(99)00037-6.
    [19] GIERAS M, KLEMENS R. Experimental studies of explosions of methane-air mixtures in a constant volume chamber [J]. Combustion Science and Technology, 2009, 181(4): 641–653. DOI: 10.1080/00102200802665102.
    [20] MITTAL M. Explosion pressure measurement of methane-air mixtures in different sizes of confinement [J]. Journal of Loss Prevention in the Process Industries, 2017, 46: 200–208. DOI: 10.1016/j.jlp.2017.02.022.
  • 期刊类型引用(5)

    1. 贾进章,张先如,王枫潇. NH_4H_2PO_4-KHCO_3混合粉体对管网甲烷爆炸的抑制特性. 爆破器材. 2024(01): 43-50 . 百度学术
    2. 杨克,王跃胜,纪虹,郑凯,邢志祥,毕海普,蒋军成. 管道中纳米二氧化硅、二氧化碳和七氟丙烷气固混合物抑制瓦斯爆炸的试验研究. 安全与环境工程. 2023(05): 28-36 . 百度学术
    3. 钟飞翔,郑立刚,马鸿雁,杜德朋,王玺,潘荣锟. CH_4/O_2/CO_2预混体系爆炸动力学研究. 爆炸与冲击. 2022(01): 4-17 . 本站查看
    4. 杨克,王辰升,纪虹,郑凯,邢志祥,毕海普,蒋军成. 聚多巴胺包覆混合粉体抑制甲烷爆炸的实验研究. 化工学报. 2022(09): 4245-4254 . 百度学术
    5. 秦毅,陈小伟,黄维. 密闭空间可燃气体爆炸超压预测. 爆炸与冲击. 2020(03): 42-53 . 本站查看

    其他类型引用(3)

  • 加载中
图(5)
计量
  • 文章访问数:  5215
  • HTML全文浏览量:  1558
  • PDF下载量:  75
  • 被引次数: 8
出版历程
  • 收稿日期:  2019-01-04
  • 修回日期:  2019-02-21
  • 网络出版日期:  2020-01-25
  • 刊出日期:  2020-01-01

目录

/

返回文章
返回