1、第51卷第4期2023 年8 月D0I:10.7631/issn.1000-2243.22465N,对甲烷-空气预混气体爆炸火焰传播稳定性的影响福州大学学报（自然科学版）Journal of Fuzhou University(Natural Science Edition)Vol.51 No.4Aug.2023文章编号：10 0 0-2 2 43(2 0 2 3)0 4-0 58 2-0 7王金贵，何伟杰，张成名，张苏1（1.福州大学环境与安全工程学院，福建福州350 10 8；2.32 50 4部队，江西上饶33410 0）摘要：基于自主搭建的1m球形爆炸实验系统，通过分析火焰结构与传播过
2、程来探究不同体积分数下的甲烷和氮气对火焰不稳定性的影响及其机制.研究结果表明，甲烷-空气预混气体火焰发展过程与甲烷体积分数直接相关，当甲烷体积分数接近爆炸极限，火焰明显上浮，浮力不稳定性是火焰失稳的主要因素当甲烷体积分数在7%13%范围内时，火焰总体呈规则球形，在热-质扩散因素和流体力学因素相互促进或制约作用下，火焰总体不稳定性呈现先迅速增强后缓慢减弱的趋势.通过改变氮气添加比，发现流体力学不稳定性对氮气添加比的变化更为敏感，当氮气添加比（体积分数）大于8%时，在流体力学不稳定性因素的主导影响下火焰稳定性迅速增强.关键词：甲烷爆炸；惰化；马克斯坦长度；火焰不稳定性中图分类号：X932Influ
3、ence of N,on the flame stability of methane-air premixed gasWANG Jingui,HE Wejie,ZHANG Chengming,ZHANG Su(1.College of Environment and Safety Engineering,Fuzhou University,Fuzhou,Fujian 350108,China;2.Unit 32504,Shangrao,Jiangxi 334100,China)Abstract:Based on the self-developed 1 m spherical explosi
4、on experimental system,investigating theeffect of methane and nitrogen on flame instability at different volume fractions and their mechanisms byanalyzing flame structure and propagation processes.The results show that the flame development processof methane-air premixed gas is directly related to t
5、he volume fraction of methane,and when the volumefraction of methane is close to the explosion limit,the flame obviously floats upward and buoyancy insta-bility is the main factor of flame destabilization.When the volume fraction of methane is in the rangefrom 7%to 13%,the flame has a regular spheri
6、cal shape in general,and the overall instability of theflame shows a trend of rapid enhancement and then slow weakening under the effect of heat-massdiffusion factors and hydrodynamic factors promoting or constraining each other.By changing the volumefraction of nitrogen,it was found that the hydrod
7、ynamic instability was more sensitive to the change ofvolume fraction of nitrogen,and when the volume fraction of nitrogen was greater than 8%,the flamestability was enhanced rapidly under the dominant influence of hydrodynamic instability factors.Keywords:methane explosion;inerting;Markstein length
8、;flame instability0引言甲烷作为当前较为低碳并广泛使用的清洁能源 1，其储存、运输和使用等环节依然面临爆炸的风险.甲烷的燃爆、抑爆研究大致分为两类，一类是测定火焰燃烧速度、爆炸超压、爆炸温度及其他基础参数探究气体配比、点火能量或障碍物等对爆炸特性的影响 2-5，另类主要是使用高速摄像、纹影摄像等技术研究各工况对火焰结构、火焰褶皱、爆炸波形态等火焰演化过程的影响 6-8 .Hu等 9 在定容燃烧弹内对甲烷-氢气-空气预混火焰的胞状不稳性进行实验研究，发现火焰不稳定性与初始压力呈正相关，而初始温度对火焰不稳定性影响较小.苗海燕等 10 在定容燃烧弹内研究惰性气体抑制下掺氢天然气-
9、空气预混气收稿日期：2 0 2 2-11-0 6通信作者：张苏（198 7-），博士，副教授，主要从事安全与应急管理领域研究，基金项目：福建省自然科学基金面上资助项目（2 0 2 0 J01505）文献标识码：A第4期体的火焰特性，结果表明：氢气的马克斯坦长度与惰性气体的添加比呈反比，而天然气的实验结果与氢气的相反，惰性气体的添加有利于天然气燃烧火焰的稳定.Prathap等 1研究CO，和N，对H,/CO预混气体燃烧速度和火焰不稳定性的影响，发现CO，对预混气体反应的抑制作用要强于N2.Kwon等12 1 则在高压条件下研究压力对于H,/O2/N，混合物火焰不稳定性的影响，发现初始压力的增加会
10、导致火焰厚度降低，即火焰的不稳定性增强.预混火焰的火焰不稳定性有3种：流体力学不稳定性、热-质扩散不稳定性和浮力不稳定性 13-15。以往研究大多将3种火焰不稳定性单独进行分析，并且浮力不稳定性的研究缺乏实验现象的支撑，本研究将利用中尺寸（1m）球形容器，通过高速摄像机记录火焰结构演化过程，更好地呈现出浮力对火焰结构的影响，研究不同甲烷体积分数及不同氮气配比下球形燃烧扩展火焰的失稳机理，分析3种火焰不稳定性的竞争及影响机制，研究成果对于完善甲烷抑爆机制具有重要意义.1实验系统实验系统由1m球形爆炸容器、高速摄影系统、数据采集系统、配气系统及点火系统5部分组成.实验系统如图1所示.CHN,王金贵
11、，等：N，对甲烷-空气预混气体爆炸火焰传播稳定性的影响真空压力表1m球形爆炸容器点火电极250mm泄压口coAirHD15-10型数显高压脉冲放电电源583压力传感器信号采集系统高速摄像机回电脑DDS信号发生器X真空泵图11m球形爆炸实验系统Fig.1 1 m spherical explosion test system为研究爆炸火焰结构，在1m球形爆炸容器前后两侧装有2 面直径2 50 mm、厚度35mm的光学石英玻璃视窗，为高速摄影机广角拍摄火焰提供光学通道.DDS信号发生器可同时触发高速摄像机及压力信号采集系统.配气系统包括真空泵、配气柱、气瓶及数字压力表等.点火系统由脉冲式点火器和铂
12、金点火电极组成，点火电极位于容器中心，两根电极间隙为3mm，点火电压为15kV.所有实验均在一个大气压和285K的初始温度下进行.2火焰稳定性理论及参数计算方法2.1马克斯坦长度在密闭容器中心点火，混合气火焰呈球形扩散.随着火焰的发展，火焰表面积逐渐增大，火焰锋面持续拉伸，导致表面开始出现胞状火焰引起燃烧速率的增加从而造成火焰的不稳定.在球型火焰尚未受到其他因素影响时，拉伸层流火焰传播速度S，和拉伸率K存在线性关系16 .即S。-S,=L,K式中：S。为未拉伸火焰速度；S，为拉伸层流火焰速度；L，为未燃烧预混气体的马克斯坦长度；K为拉伸率.从式（1)可知马克斯坦长度是拉伸层流火焰速度与拉伸率直
13、线的斜率相反数.马克斯坦长度表征了层流预混火焰对拉伸响应的敏感程度，反应了气体爆炸初期层流火焰的稳定性，马克斯坦长度越小，火焰越容易受火焰不稳定性影响.(1)http:/5842.2刘易斯数刘易斯数Le是混合物的热扩散率与贫反应物质量扩散率的比值.在富燃料区域，质量扩散率取氧气的质量扩散率；而在贫燃料区域，则取燃料的质量扩散率刘易斯数用于表征火焰的热-质扩散不稳定性，刘易斯数越大，火焰表面的扰动越强，火焰稳定性就越差.当Le1时，即热扩散系数 质量扩散系数.火焰向未燃区凸起部分，热量耗散大，火焰传播速度随即减小，凸起部分逐渐趋于平缓，使火焰褶皱面趋于稳定，当Le1时，即热扩散系数 1 3%7%
14、1 2%1 1%9%10%.从图4 可知，各甲烷体积分数火焰上浮率从高到低依次为：6%7%1 3%1 2%1 1%9%1 0%.综合图3 4 可知，上浮率随着甲烷体积分数的增大先减后升，而火焰速度（S，）随着甲烷体积分数的增大而先增后减.火焰速度越低，浮力不稳定性的影响越明显，如图5 所示.由图5 可知各甲烷体积分数火焰上浮率均大于0，说明各体积分数下均存在浮力不稳定性，但火焰速度较低时，浮力不稳定性的影响更明显，对比图3中各体积分数下的甲烷火焰上浮情况也可看出，甲烷接近爆炸极限（6%与1 3%）时，火焰速度小，火焰受浮力不稳定性作用增强而导致明显上浮.分析认为火焰受拉伸力、浮力和重力的耦合作
15、用，但在慢速燃烧扩散过程中，因热释放速率低，热损失大，浮力是导致火焰出现不稳定性的主要因素，因此火焰几何中心随着燃烧速率的降低而逐渐上移，呈现不规则形状.3.1.2甲烷体积分数对热-质扩散、D-L不稳定性的影响采用前述公式计算不同甲烷体积分数火焰的不稳定性参数，结果如表2 所示.图6 7 分别为马克斯坦长度(L)刘易斯数(Le)和火焰厚度(L)随着甲烷体积分数变化曲线.Tab.2Parameter calculation results(no inerting)参数6Ly/mm20.91Le0.887 600./mm从图3 可以看出，甲烷体积分数为7%1 3%范围内时，火焰总体为球形，对于球形
16、火焰可用马克斯坦长度衡量火焰总体不稳定性。观察图6 可以发现，甲烷体积分数小于1 0%时，甲烷燃烧球形火焰的马克斯坦长度随着甲烷体积分数的增大而快速减小，当甲烷体积分数大于1 0%时，马克斯坦长度也随之缓慢6%CH1.87%CH.9%CH.1.410%CH.-11%CH12%CH,1.013%CH0.60.2-0.2025 5075100125150175200R,/mm图4各甲烷体积分数上浮率变化曲线图Fig.4Variation curve of upwelling rate at eachmethane volume fractions2.52.01.S.w/1.51.00.50567 891011121314图5火焰传播速度和上浮率变化曲线图Fig.5Variation curves of flame propagationspeed and uplift rate表2 参数计算结果(无惰化)PcH,/%799.502.460.876 280.853 790.104 50.051 2http:/1.50一一火焰传播速度-上浮率1.251.000.750.500.250-0.25