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

油气爆炸过程火焰燃烧模式的实验估计

张培理 杜扬

引用本文:
Citation:

油气爆炸过程火焰燃烧模式的实验估计

    作者简介: 张培理(1985-),男,博士研究生,zpl612323@163.com;
  • 基金项目: 国家自然科学基金项目 51276195
    后勤工程学院青年基金项目 YQ16-420802

  • 中图分类号: O389

Experimental estimation of the combustion regime in the oil-gas explosion process

  • CLC number: O389

  • 摘要: 首先分析讨论了油气爆炸过程中火焰燃烧模式的估计方法,然后在激波管内进行了低、中、高3次不同初始油气浓度条件下的油气爆炸实验,通过实验数据分别计算出了低、中、高初始油气浓度条件下油气爆炸在初期、中期和后期的丹姆克尔数和湍流雷诺数,最后依靠丹姆克尔数-湍流雷诺数图对低、中、高初始油气浓度条件下油气爆炸初期、中期和后期的火焰燃烧模式进行了定量估计。结果表明:低、中、高初始油气浓度条件下激波管油气爆炸过程初期、中期和后期的火焰燃烧模式均为漩涡内小火焰模式。
  • 图 1  湍流预混火焰的3种模式随DaRel0的分布

    Figure 1.  Distribution of the three turbulence premixed flame models based on values of Da and Rel0

    图 2  实验装置布置示意图

    Figure 2.  Arrangement of the experimental equipments

    图 3  流场速度和爆炸超压随时间的变化曲线

    Figure 3.  Variation curves of gas velocity and pressure vs. time

    图 4  3组电偶采集到的流场温度随时间的变化曲线

    Figure 4.  Temperature vs. time curves acquired by the three thermocouples

    图 5  油气爆炸过程火焰燃烧模式分布

    Figure 5.  Distribution of the flame regimes of oil-gas explosion process

    图 6  油气爆炸在148、152和156 ms时的火焰高速摄影照片

    Figure 6.  High speed photos of the flame when the time was 148, 152 and 156 ms

    表 1  各时刻激波管内气体流速、压力和已燃气体温度实验数据(ϕ=1)

    Table 1.  Flow velocity, pressure and burned gas temperature in the shock tube at different times (ϕ=1)

    t/ms v/(m·s-1) vrms* p/Pa Tmax/K
    50 55.82 34.62 98 530 1 127.39
    150 72.76 51.56 340 800 1 203.62
    250 29.14 7.94 486 190 1 324.05
    下载: 导出CSV

    表 2  参数BMB2ϕM取值

    Table 2.  Value of BM, B2 and ϕM

    燃料 ϕM BM/(cm·s-1) B2/(cm·s-1)
    甲醇 1.11 36.92 -140.51
    丙烷 1.08 34.22 -138.65
    异辛烷 1.13 26.32 -84.72
    RMFD-303 1.13 27.58 -78.34
    下载: 导出CSV

    表 3  各时刻的DaRel0计算值(ϕ=1)

    Table 3.  Calculated values of Da & Rel0 at different times (ϕ=1)

    t/ms Da Rel0
    50 0.19 4 603
    150 0.26 31 366
    250 21 4 012
    下载: 导出CSV

    表 4  各时刻激波管内气体流速、压力和已燃气体温度实验数据(ϕ=0.72)

    Table 4.  Flow velocity, pressure and burned gas temperature in the shock tube at different times (ϕ=0.72)

    t/ms v/(m·s-1) vrms* p/Pa Tmax/K
    50 65.13 45.45 115 120 1 158.42
    150 -5.74 25.42 322 410 1 201.77
    250 11.04 8.64 357 380 1 238.51
    下载: 导出CSV

    表 5  各时刻激波管内气体流速、压力和已燃气体温度实验数据(ϕ=1.28)

    Table 5.  Flow velocity, pressure and burned gas temperature in the shock tube at different times (ϕ=1.28)

    t/ms v/(m·s-1) vrms* p/Pa Tmax/K
    50 29.69 7.10 11 510 1 251.63
    150 86.91 64.32 227 780 1 296.15
    250 12.64 9.95 504 310 1 357.36
    下载: 导出CSV

    表 6  各时刻的DaRel0计算值

    Table 6.  Calculated values of Da & Rel0 at different times

    t/ms Da Rel0
    ϕ=0.72 ϕ=1.28 ϕ=0.72 ϕ=1.28
    50 0.03 1.31 6 767 93
    150 0.10 0.15 9 993 15 838
    250 0.31 1.76 3 589 4 970
    下载: 导出CSV
  • [1] Steinberg A M, Driscoll J F. Straining and wrinkling processes during turbulence-premixed flame interaction measured using temporally-resolved diagnostics[J]. Combustion and Flame, 2009, 156(12):2285-2306. doi: 10.1016/j.combustflame.2009.06.024
    [2] Shin D H, Lieuwen T. Flame wrinkle destruction processes in harmonically forced, laminar premixed flames[J]. Combustion and Flame, 2012, 159(11):3312-3322. doi: 10.1016/j.combustflame.2012.06.015
    [3] Yi Y, Geng L, Jing G. Experimental study on the fractal characteristic of methane explosion flame[J]. Safety Science, 2012, 50(4):679-683. doi: 10.1016/j.ssci.2011.08.052
    [4] Won S H, Windom B, Jiang B, et al. The role of low temperature fuel chemistry on turbulent flame propagation[J]. Combustion and Flame, 2014, 161(2):475-483. doi: 10.1016/j.combustflame.2013.08.027
    [5] Zhang M, Wang J, Wu J, et al. Flame front structure of turbulent premixed flames of syngas oxyfuel mixtures[J]. International Journal of Hydrogen Energy, 2014, 39(10):5176-5185. doi: 10.1016/j.ijhydene.2014.01.038
    [6] Mukaiyama K, Shibayama S, Kuwana K. Fractal structures of hydrodynamically unstable and diffusive-thermally unstable flames[J]. Combustion and Flame, 2013, 160(11):2471-2475. doi: 10.1016/j.combustflame.2013.05.017
    [7] Williams F A. Asymptotic methods in turbulent combustions[J]. AIAA Journal, 1986, 24:867-875. doi: 10.2514/3.9361
    [8] Abraham J, Williams F A, Bracco F V. A discussion of turbulent flame structure in premixed charges[R]. Warrendale, Pennsylvania, USA: Society of automotive engineers, 1985.
    [9] Zhang P, 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
    [10] Yang D, Li Z P, Hong O Y. Effects of humidity, temperature and slow oxidation reactions on the occurrence of gasoline-air explosions[J]. Journal of Fire Protection Engineering, 2013, 23(3):226-238. doi: 10.1177/1042391513486464
    [11] 杨辉, 崔鑫, 郑昕等.管道中湍流强度及湍流积分尺度随时间的变化研究[J].安全与环境工程, 2013, 20(4):102-104. doi: 10.3969/j.issn.1671-1556.2013.04.023
    Yang H, Cui X, Zheng X, et al.Investigation of variation of turbulence intensity and turbulent integral scale with time in tube[J]. Safety and Environmental Engineering, 2013, 20(4):102-104. doi: 10.3969/j.issn.1671-1556.2013.04.023
    [12] Sak C, Liu R, Ting D S-K, et al. The role of turbulence length scale and turbulence intensity on forced convection from a heated horizontal circular cylinder[J]. Experimental Thermal and Fluid Science, 2007, 31(4):279-289. doi: 10.1016/j.expthermflusci.2006.04.007
    [13] Metghalchi M, Keck J C. Burning velocities of mixtures of air with methanol, isooctane, and indolene at high pressure and temperature[J]. Combustion and Flame, 1982, 48:191-210. doi: 10.1016/0010-2180(82)90127-4
    [14] Spadling D B. Combustion and Mass Transfer[M]. New York: Pergamon, 1979:59.
    [15] Stephen R T. An introduction to combustion:concepts and application[M]. New York: The McGraw-Hill Companies, 2000:appendix C.
    [16] Ballal D R, Lefebvre A H. The structure and propagation of turbulent flames[J]. Proceedings of the Royoal Society of London, 1979, 344(1637):217-234.
    [17] Ballal D R, Lefebvre A H. The structure of a premixed turbulent flame[J]. Proceedings of the Royoal Society of London, 1979, 367(1730):253-280.
    [18] Poinsot T, Candel S, Trouvé A. Applications of direct numerical simulation to premixed turbulent combustion[J]. Progress in Energy and Combustion Science, 1995, 21(6):531-576. doi: 10.1016/0360-1285(95)00011-9
  • [1] 杜扬王世茂齐圣王波李阳超李国庆 . 油气在顶部含弱约束结构受限空间内的爆炸特性. 爆炸与冲击, 2017, 37(1): 53-60. doi: 10.11883/1001-1455(2017)01-0053-08
    [2] 张志彪黄风雷 . 内部爆炸加载下柱壳的环周分裂数. 爆炸与冲击, 2015, 35(5): 763-767. doi: 10.11883/1001-1455(2015)05-0763-05
    [3] 张社荣王高辉王超孙博 . 水下爆炸冲击荷载作用下混凝土重力坝的破坏模式. 爆炸与冲击, 2012, 32(5): 501-507. doi: 10.11883/1001-1455(2012)05-0501-07
    [4] 张启灵李端有李波 . 水下爆炸冲击作用下重力坝的损伤发展及破坏模式. 爆炸与冲击, 2012, 32(6): 609-615. doi: 10.11883/1001-1455(2012)06-0609-07
    [5] 朱锡牟金磊王恒张振华 . 水下爆炸载荷作用下加筋板的毁伤模式. 爆炸与冲击, 2010, 30(3): 225-231. doi: 10.11883/1001-1455(2010)03-0225-07
    [6] 牟金磊朱锡张振华王恒 . 水下爆炸载荷作用下加筋板的毁伤模式. 爆炸与冲击, 2009, 29(5): 457-462. doi: 10.11883/1001-1455(2009)05-0457-06
    [7] 侯海量朱锡梅志远 . 舱内爆炸载荷及舱室板架结构的失效模式分析. 爆炸与冲击, 2007, 27(2): 151-158. doi: 10.11883/1001-1455(2007)02-0151-08
    [8] 曹卫国徐森梁济元高伟潘峰饶国宁 . 煤粉尘爆炸过程中火焰的传播特性. 爆炸与冲击, 2014, 34(5): 586-593. doi: 10.11883/1001-1455(2014)05-0586-08
    [9] 喻健良闫兴清 . 硅酸铝棉对火焰速度和爆炸超压的抑制作用. 爆炸与冲击, 2013, 33(4): 363-368. doi: 10.11883/1001-1455(2013)04-0363-06
    [10] 高伟阿部俊太郎荣建忠土桥律 . 气流特征对水平长管内石松子粉尘爆炸火焰结构的影响. 爆炸与冲击, 2015, 35(3): 372-379. doi: 10.11883/1001-1455-(2015)03-0372-08
    [11] 薛冰马宏昊沈兆武余勇 . 爆炸容器内小药量实验动态标定压力传感器. 爆炸与冲击, 2015, 35(3): 437-441. doi: 10.11883/1001-1455(2015)03-0437-05
    [12] 朱跃进董刚刘怡昕范宝春蒋华 . 激波诱导火焰变形、混合和燃烧的数值研究. 爆炸与冲击, 2013, 33(4): 430-437. doi: 10.11883/1001-1455(2013)04-0430-08
    [13] 杜扬李国庆吴松林张培理周毅齐圣王世茂 . T型分支管道对油气爆炸强度的影响. 爆炸与冲击, 2015, 35(5): 729-734. doi: 10.11883/1001-1455(2015)05-0729-06
    [14] 李国庆杜扬齐圣王世茂李蒙李润 . 连续圆孔障碍物对油气泄压爆炸火焰特性影响大涡模拟. 爆炸与冲击, 2018, 38(6): 1286-1394. doi: 10.11883/bzycj-2017-0215
    [15] 高洪泉卢芳云王少龙颜澎罗永锋袁伟胡建 . 具有爆炸反应特性装置在密闭容器内爆炸威力的评价方法. 爆炸与冲击, 2011, 31(3): 306-310. doi: 10.11883/1001-1455(2011)03-0306-05
    [16] 赵志红郭建春 . 层内爆炸压裂岩石破碎颗粒尺寸的预测模型. 爆炸与冲击, 2011, 31(6): 669-672. doi: 10.11883/1001-1455(2011)06-0669-04
    [17] 范进徐大立任新见 . 接触爆炸荷载作用下带孔防护结构内冲击波的传播规律. 爆炸与冲击, 2014, 34(6): 658-666. doi: 10.11883/1001-1455(2014)06-0658-09
    [18] 王启睿张晓忠孔福利张福明 . 多级穿廊结构坑道口部内爆炸冲击波传播规律的实验研究. 爆炸与冲击, 2011, 31(5): 449-454. doi: 10.11883/1001-1455(2011)05-0449-06
    [19] 罗振敏王涛程方明宋钰吴慷 . 小尺寸管道内二氧化碳抑制甲烷爆炸效果的实验及数值模拟. 爆炸与冲击, 2015, 35(3): 393-400. doi: 10.11883/1001-1455-(2015)03-0393-08
    [20] 陈小伟张方举梁斌谢若泽徐艾明 . A3钢钝头弹撞击45钢板破坏模式的试验研究. 爆炸与冲击, 2006, 26(3): 199-207. doi: 10.11883/1001-1455(2006)03-0199-09
  • 加载中
图(6)表(6)
计量
  • 文章访问数:  426
  • HTML全文浏览量:  3
  • PDF下载量:  409
  • 被引次数: 0
出版历程
  • 收稿日期:  2014-11-10
  • 录用日期:  2015-03-15
  • 刊出日期:  2016-09-25

油气爆炸过程火焰燃烧模式的实验估计

    作者简介:张培理(1985-),男,博士研究生,zpl612323@163.com
  • 后勤工程学院军事供油工程系,重庆 401311
基金项目:  国家自然科学基金项目 51276195后勤工程学院青年基金项目 YQ16-420802

摘要: 首先分析讨论了油气爆炸过程中火焰燃烧模式的估计方法,然后在激波管内进行了低、中、高3次不同初始油气浓度条件下的油气爆炸实验,通过实验数据分别计算出了低、中、高初始油气浓度条件下油气爆炸在初期、中期和后期的丹姆克尔数和湍流雷诺数,最后依靠丹姆克尔数-湍流雷诺数图对低、中、高初始油气浓度条件下油气爆炸初期、中期和后期的火焰燃烧模式进行了定量估计。结果表明:低、中、高初始油气浓度条件下激波管油气爆炸过程初期、中期和后期的火焰燃烧模式均为漩涡内小火焰模式。

English Abstract

  • 油气(包括天然气、瓦斯和原油、成品油挥发物等)属于典型的易燃易爆气体,稍有不慎就可能引发火灾爆炸安全事故,造成重大人员伤亡和巨额经济损失。近年来,随着国民经济与社会生产的快速发展,受限空间油气火灾爆炸事故时有发生,如2007年我国西南某油库发生的油库爆炸事故;2013年11月,山东青岛原油输油管线发生爆炸事故;2014年7月,台湾高雄发生燃气爆炸事故,以上事故均带来严重的人员伤亡和财产损失。为了避免此类事故的发生,减少人员伤亡与财产损失,尚需加强对油气爆炸机理的研究,从而提出有效的油气爆炸防治措施。

    受限空间油气爆炸火焰属于典型的湍流预混火焰,湍流火焰的燃烧模式与油气爆炸过程的化学反应过程、湍流流动以及传热传质过程密切相关。一般来说,不同的油气爆炸火焰燃烧模式会导致不同的油气爆炸特性,如爆炸超压、超压上升速率、火焰速度等,究其内在的原因之一就是湍流会使火焰前锋面发生褶皱和扭曲[1],进而对油气爆炸过程的传热传质[2-3],甚至化学反应历程[4-5]产生影响。因此,研究受限空间油气爆炸火焰的燃烧模式对于进一步认识油气爆炸过程的机理和规律、构建油气爆炸过程数学模型具有十分重要的意义。

    • 对湍流火焰燃烧模式的研究一直是可燃气体燃烧和爆炸领域的热点和难点[6]。F.A.Williams[7]和J.Abraham等[8]根据湍流柯尔莫格洛夫微尺度lk、湍流积分尺度l0和层流火焰厚度δL等参数构成的判据,将湍流预混火焰分为3种模式:褶皱层流火焰模式、漩涡内小火焰模式和分布反应模式。具体方法是先计算湍流预混火焰的丹姆克尔数Da和湍流雷诺数Rel0,然后根据反映湍流预混火焰3种模式的湍流雷诺数-丹姆克尔数图(如图 1所示),来确定湍流火焰的燃烧模式。

      图  1  湍流预混火焰的3种模式随DaRel0的分布

      Figure 1.  Distribution of the three turbulence premixed flame models based on values of Da and Rel0

      图 1中,粗实线(lk/δL=1)上方是满足威廉斯-克里莫夫判据的褶皱层流火焰存在的条件,粗实线(l0/δL=1)下方是满足丹姆克尔判据的分布反应模式存在的条件,两条粗实线之间是满足漩涡内小火焰模式的分布条件。因此只要计算出油气爆炸过程的丹姆克尔数Da和湍流雷诺数Rel0,根据图 1就可以判断油气爆炸火焰燃烧的模式。

      丹姆克尔数Da在燃烧中是一个很重要的量纲一参数,定义为:

      $Da=\frac{{{\tau }_{\text{flow}}}}{{{\tau }_{\text{chem}}}}=\frac{{{l}_{0}}/v_{\text{rms}}^{*}}{{{\delta }_{\text{L}}}/{{S}_{\text{L}}}}=\left( \frac{{{l}_{0}}}{{{\delta }_{\text{L}}}} \right)\left( \frac{{{S}_{\text{L}}}}{v_{\text{rms}}^{*}} \right) $

      式中:τflow为流动特征时间,τchem为化学特征时间,vrms*为湍流脉动速度,SL为层流火焰速度。l0/SL为火焰几何尺度比,SL/vrms*为火焰相对湍流强度。

      湍流雷诺数Rel0的定义为:

      $ R{{e}_{{{l}_{0}}}}=\frac{\rho v_{\text{rms}}^{*}{{l}_{0}}}{\mu } $

      式中:ρ为燃烧火焰的密度,μ可燃气体的为动力粘度。

      从上面的分析可以看出,要想判断油气爆炸火焰的模式,必须得到油气爆炸过程的湍流脉动速度vrms*,湍流积分尺度l0,层流火焰厚度δL和层流火焰速度SL等参数,显然在实验室要想直接获得这些参数是非常困难的。

      本文拟通过激波管中的油气爆炸实验,测量油气爆炸过程中已燃气体和未燃气体的温度、平均流速、爆炸压力等数据来估计油气爆炸过得的丹姆克尔数Da和湍流雷诺数Rel0,进而估计油气爆炸过程的火焰模式。

    • 本文中的实验装置主要包括激波管、真空循环泵、油气雾化装置、数据采集系统、点火系统、高速摄影仪、电脑等。数据采集系统,主要包括压力采集系统、浓度采集仪、温度采集设备和流场速度设备。实验装置布置示意图如图 2所示。

      图  2  实验装置布置示意图

      Figure 2.  Arrangement of the experimental equipments

      激波管尺寸为200 mm×200 mm×6 100 mm,其中包括一段长300 mm的玻璃观察段。油气雾化装置的结构和工作原理见参考文献[9-10]。

      压力采集系统主要由压力传感器、采集卡和电脑组成。沿激波管共布置了3个压力传感器,实验中以这3个压力传感器采集到数据的平均值作为激波管内的压力值。浓度采集系统主要由GXH-1050型红外分析仪和NHA-502型汽车尾气分析仪组成。实验中油气浓度由GXH-1050型红外分析仪测试,而其他气体组分,如O2、CO、CO2等,均由NHA-502型汽车尾气分析仪测量。

      实验中未燃气体和已燃气体的温度由热电偶、采集卡和电脑组成的温度采集设备测量。根据实验需要,沿激波管共设置了3组热电偶,分别布置在激波管50、305和560 cm处。

      油气爆炸过程中湍流的脉动速度的测量较为困难,本文中拟先用6162型高温热线风速仪来记录油气爆炸过程中的瞬时气体流速,并根据整个油气爆炸持续的时间得到平均流速,然后根据任一时刻的瞬时速度和平均流速得到湍流脉动速度。该热线风速仪探头设置在激波管中间位置,可以每隔5~10 ms测一次气体的瞬时流速,并储存在计算机内。

      不同的初始条件会导致不同的实验结果。为了使实验结果具有代表性,选取低(1.20%)、中(1.65%)、高(2.10%)3种初始油气体积分数进行实验,经计算,其油气当量比ϕ(按异辛烷计算)分别为0.72、1.00和1.28。其他初始条件设置如下:初始温度为常温,压力为常压,点火能设为5 J。

      本文先以ϕ=1.00(即油气初始体积分数为1.65%)时的实验结果为例,给出火焰燃烧模式估计的详细步骤与结果;然后再列出初始油气体积分数分别为1.20%和2.10%时的实验数据,并分别估计其燃烧模式;最后根据3次燃烧模式计算的结果,分析和讨论激波管内油气爆炸火焰的燃烧模式。

    • 图 3为初始油气体积分数为1.65%时(ϕ=1),流场速度v和爆炸超压p随时间变化的曲线。由图 3可以看出整个油气爆炸持续了约288 ms,在此过程中激波管内部的气体速度在-60~100 m/s之间不断振荡波动,其平均流速v≈21.20 m/s。激波管内气体流速任一时刻的脉动量vrms*可以由瞬时流速v(t)和平均流速v计算:

      $ v_{\text{rms}}^{*}=\left| v\left( t \right)-\bar{v} \right| $

      图  3  流场速度和爆炸超压随时间的变化曲线

      Figure 3.  Variation curves of gas velocity and pressure vs. time

      图 4为初始油气体积分数为1.65%时(ϕ=1),3组热电偶采集到的流场温度T随时间的变化曲线。从图 4可以看出,未燃气体的温度Tu约为300 K,1#、2#和3#热电偶采集到的已燃气体的最大温度Tmax分别约为1 127.39、1 203.62和1 324.05 K。已燃气体的最大温度正是火焰经过热电偶时的燃烧气体的温度,火焰离开热电偶后,热电偶处的气体温度由于剧烈的热损失开始下降,之后温度趋于平稳,此时的温度即为燃烧后的气体温度为Tb。可以认为热电偶采集到的最大温度Tmax即为油气爆炸已燃气体的瞬时温度。由于1#、2#和3#热电偶分别设置在激波管的前端、中部和尾部,因此其各自采集到的最大温度代表了油气爆炸初期、充分发展和后期的燃烧气体的温度。

      图  4  3组电偶采集到的流场温度随时间的变化曲线

      Figure 4.  Temperature vs. time curves acquired by the three thermocouples

      本文中拟分别对油气爆炸初期(50 ms)、充分发展(150 ms)和后期(250 ms)的火焰燃烧模式进行估计。从图 3图 4可以得到ϕ=1,t=50、150和250 ms时激波管内气体流速、压力和已燃气体温度等数据,如表 1所示。

      t/ms v/(m·s-1) vrms* p/Pa Tmax/K
      50 55.82 34.62 98 530 1 127.39
      150 72.76 51.56 340 800 1 203.62
      250 29.14 7.94 486 190 1 324.05

      表 1  各时刻激波管内气体流速、压力和已燃气体温度实验数据(ϕ=1)

      Table 1.  Flow velocity, pressure and burned gas temperature in the shock tube at different times (ϕ=1)

    • 湍流积分尺度代表了湍流涡的平均尺度,根据杨辉等[11]和C.Sak等[11]的研究,管道或坑道内湍流的平均积分尺度约为管道或坑道横截面积的1/20~1/5。为了便于估算,设积分尺度为激波管横截面尺寸的1/10,即l0=0.02 m。根据式(1),先计算50 ms时的丹姆克尔数。

      其中流动特征时间为τflow=l0/vrms*=0.000 6 s。

      化学特征时间τchem=δL/SL,根据文献[13]:

      $ {{S}_{\text{L}}}={{S}_{\text{L, ref}}}{{\left( \frac{{{T}_{\text{u}}}}{{{T}_{\text{u, ref}}}} \right)}^{\gamma }}{{\left( \frac{p}{{{p}_{\text{ref}}}} \right)}^{\beta }} $

      $ {{S}_{\text{L, ref}}}={{B}_{\text{M}}}+{{B}_{2}}{{\left( \phi -{{\phi }_{\text{M}}} \right)}^{2}} $

      式中:Tu, ref=298 K;pref为标准大气压;BMB2ϕM是由燃料类型确定的常数,由表 2给出;温度指数γ和压力指数β是当量比ϕ的函数,其经验公式如下[13]

      $ \left\{ \begin{align} & \gamma =2.18-0.8\left( \phi -1 \right)\text{ } \\ & \beta =-0.16+0.22\left( \phi -1 \right) \\ \end{align} \right. $

      燃料 ϕM BM/(cm·s-1) B2/(cm·s-1)
      甲醇 1.11 36.92 -140.51
      丙烷 1.08 34.22 -138.65
      异辛烷 1.13 26.32 -84.72
      RMFD-303 1.13 27.58 -78.34

      表 2  参数BMB2ϕM取值

      Table 2.  Value of BM, B2 and ϕM

      根据表 2,选异辛烷的BM,B2ϕM参考值来估算SL。根据式5(b),SL, ref=24.89 cm/s;油气体积分数为1.65%时当量比ϕ=1,则根据式(6),γβ为2.18和-0.16。因此根据表 1SL=0.253 7 m/s。

      层流火焰厚度的计算采用D.B.Spalding的计算式[14]

      $ {{\rm{ \mathsf{ δ} }}_{\rm{L}}} = \frac{{2\alpha }}{{{S_{\rm{L}}}}} $

      式中:α为热扩散率,用油气爆炸50 ms时平均温度0.5(Tb+Tu)=713.7 K下空气的物性[15]来近似,并进行压力修正,得α=α0(pref/p)=1.04×10-4 m2/s,α0为空气的热扩散率。因此,δL=8.20×10-4 m。由此可得,τchem=3.23×10-3 s。根据式(1),Da=0.19。

      油气混合物的动力粘度μ仍以空气在Tb下的动力粘度μb来近似,根据文献[15],μb=4.56×10-5 N·s/m2,则湍流雷诺数

      $R{e_{{l_0}}} = \frac{{\rho v_{{\rm{rms}}}^*{l_0}}}{\mu } = \frac{{\left( {p/R{T_{\rm{b}}}} \right)v_{{\rm{rms}}}^*{l_0}}}{{{\mu _b}}} = 4\;603 $

      式中:R为理想气体常数。

      同理,计算出油气爆炸150和250 ms时的丹姆克尔数Da和湍流雷诺数Rel0,如表 3所示。

      t/ms Da Rel0
      50 0.19 4 603
      150 0.26 31 366
      250 21 4 012

      表 3  各时刻的DaRel0计算值(ϕ=1)

      Table 3.  Calculated values of Da & Rel0 at different times (ϕ=1)

      表 4~5分别是油气当量比为0.72和1.28时激波管内50、150和250 ms时刻气体流速、压力和已燃气体温度的实验数据。根据上文的方法和步骤,当油气当量比分别为0.72和1.28时50、150和250 ms时刻的丹姆克尔数Da和湍流雷诺数Rel0计算结果如表 6所示。

      t/ms v/(m·s-1) vrms* p/Pa Tmax/K
      50 65.13 45.45 115 120 1 158.42
      150 -5.74 25.42 322 410 1 201.77
      250 11.04 8.64 357 380 1 238.51

      表 4  各时刻激波管内气体流速、压力和已燃气体温度实验数据(ϕ=0.72)

      Table 4.  Flow velocity, pressure and burned gas temperature in the shock tube at different times (ϕ=0.72)

      t/ms v/(m·s-1) vrms* p/Pa Tmax/K
      50 29.69 7.10 11 510 1 251.63
      150 86.91 64.32 227 780 1 296.15
      250 12.64 9.95 504 310 1 357.36

      表 5  各时刻激波管内气体流速、压力和已燃气体温度实验数据(ϕ=1.28)

      Table 5.  Flow velocity, pressure and burned gas temperature in the shock tube at different times (ϕ=1.28)

      t/ms Da Rel0
      ϕ=0.72 ϕ=1.28 ϕ=0.72 ϕ=1.28
      50 0.03 1.31 6 767 93
      150 0.10 0.15 9 993 15 838
      250 0.31 1.76 3 589 4 970

      表 6  各时刻的DaRel0计算值

      Table 6.  Calculated values of Da & Rel0 at different times

      表 3表 6中的数据,绘入丹姆克尔数-湍流雷诺数图中,如图 5所示。由图 5可以看出,油气当量比分别为0.72, 1.00和1.28时,油气爆炸50、150和250 ms时的火焰燃烧模式均落入了漩涡内小火焰模式区域。也就是说,低中高初始油气浓度下,油气爆炸初期、充分发展和后期的火焰燃烧模式均为漩涡小火焰模式。这一结论与D.R.Ballal等[16-17]对预混丙烷火焰燃烧模式的实验研究结论一致,他们通过实验发现大部分受限流动中的预混丙烷火焰燃烧模式与漩涡内小火焰模式相关。

      图  5  油气爆炸过程火焰燃烧模式分布

      Figure 5.  Distribution of the flame regimes of oil-gas explosion process

      图 6是当初始油气体积分数为1.65%时油气爆炸在148、152和156 ms时刻火焰的高速摄影照片。由图 6可以清楚地看出油气爆炸火焰前锋产生了明显的褶皱,并且随着油气爆炸的发展,火焰前锋后生成了许多小火焰点。这些小火焰点很可能与火焰前锋褶皱的不断加强破碎过程有关。

      图  6  油气爆炸在148、152和156 ms时的火焰高速摄影照片

      Figure 6.  High speed photos of the flame when the time was 148, 152 and 156 ms

      在漩涡内小火焰模式下,油气爆炸火焰的燃烧支持漩涡破碎模型的理论[18],即认为燃烧速度取决于未燃气体破碎成更小微团的速度,由于不断地破碎,使得未燃混合物与已燃热烟气之间有足够的界面进行反应[4],这表明,并不是化学反应的速率决定着燃烧速度,而是湍流混合速度控制着燃烧过程。这一点在对受限空间油气爆炸火焰燃烧的数值分析模型进行选择时至关重要,该结果为受限空间油气爆炸火焰燃烧的数值分析模型的选择提供了依据。

    • 本文首先根据实验先得到了低、中、高3种初始油气浓度条件下油气爆炸过程的气体流速、压力以及火焰温度。通过对实验数据的分析和处理,分别确定了低、中、高3种初始油气浓度条件下油气爆炸过程初期(50 ms)、中期(150 ms)和后期(250 ms)时的气体脉动速度、压力、已燃气体温度和未燃气体温度,最后根据丹姆克尔数和湍流雷诺数对低、中、高3种初始油气浓度条件下油气爆炸过程中的初期、中期和后期火焰燃烧模型分别进行了定量估计。结果表明,低、中、高3种初始油气浓度条件下,激波管油气爆炸过程初期、中期和后期的火焰燃烧模式均为漩涡内小火焰模式;在漩涡内小火焰模式下,油气爆炸过程中火焰的燃烧支持漩涡破碎模型理论,即燃烧速度取决于未燃气体破碎成更小微团的速度,这也表明并不是化学反应的速率决定着燃烧速度,而是湍流混合速度控制着燃烧过程。本文的研究方法和结论可为进一步研究油气爆炸火焰燃烧模式及其数值分析模型的选择提供的参考。

参考文献 (18)

目录

    /

    返回文章
    返回