瓦斯爆炸运动火焰生产压力波的数值模拟.pdf

1994-2010 China Academic Journal Electronic Publishing House. All rights reserved. 收稿日期 20041116 基金项目国家自然科学基金重点项目50134040 作者简介吴 兵19672 , 男,山西省阳泉市人,副教授,工学博士,从事安全工程方面的研究 1 E-mail wbelcyvip. sina. com 3Donald N H C, Henry E P. M athematical study of propagating flame and its aerodynam ics in a coalm ine passway. Report of Investiga2 tion7908, Pisttsburgh M ining and Safety Research Center, Pisttsburgh, Califonia, U SA, 1974. 第34卷 第4期 中国矿业大学学报 Vol . 34 No. 4 2005年7月 Journalof China U niversity ofM ining 2.华北科技学院 煤矿安全技术培训中心,河北 燕郊 101601 摘要从三维N2S方程出发,用TVD格式,对瓦斯爆炸过程中火焰产生压力波的过程进行了数 值模拟.在此基础上,模拟了氢氧燃烧驱动的破膜过程以及破膜前后压缩波、 稀疏波对火焰阵面 的影响 1 同时,也研究了瓦斯爆炸过程中,压力波、 火焰与障碍物的相互作用 1 数值模拟结果与 理论分析吻合较好,从而进一步验证了该程序能处理含有化学反应和复杂管道的预混可燃气体 爆炸问题. 关键词瓦斯爆炸;数值模拟;火焰;压力波;障碍物 中图分类号 TD 712 文献标识码 A Numerical Si mulation of Shock W ave Induced byMoving Flame for Gas Explosion WU Bing 1, ZHAN G L i2cong1, XU Jing2de2 1. Schoolof Resources and Safety Engineering, China U niversity ofM ining 2. The Training Centre of CoalM ine Safety Technology, North China Institute of Science and Technology, Yanjiao, Hebei 101601, China Abstract Three di mensionalN2S equations and TVD scheme were used to establish the numerical model that is applied to si mulate the pressure wave induced by flame during the gas explosion. On the basis of thismodel, the process of the diaphragm bursting drived by hydrogen2air m ixture gas deflagration was si mulated.The influence of the compressed wave and rare factionalwave on the flame wave was analyzed. A t the same ti me, the interaction among shock wave, flame wave and barrier was studied.The results show that the numerical si mulation is in accord w ith the theory analysis .Therefore, this numerical si mulation program can be used to analyze the pre2 m ixed com2 bustible gas explosion that includes chem ical reaction in the tube w ith complex conditions . Key words gas explosion; numerical si mulation; flame; shock wave; barrier 矿井瓦斯爆炸后,火焰锋面沿巷道传播超前于 火焰锋面的压力波,使煤尘从沉积状态变为飞扬状 态,形成煤尘云3.在一定条件下,煤尘云会引起二 次爆炸.从本质上看爆炸是火焰生成压力波的过 程,搞清火焰传播生成压力波的物理机制具有十分 重要的意义[122].在实际研究中,主要考虑爆炸过程 中压缩波、 稀疏波、 火焰锋面和障碍物之间的相互 作用[324].本文采用氢氧燃烧驱动引燃激波管中按 1994-2010 China Academic Journal Electronic Publishing House. All rights reserved. 当量比预混的甲烷和空气混合气体 1 其工作原理 是点火后,火焰在管道中自左向右传播,在火焰的 前方诱导出压力波,压力波在膜片上反射,当膜片 处的压力大于破膜压力时,膜片破裂.实践证明,预 混气体燃烧导致的爆炸均与空间受限有关[5],如井 下巷道壁面、 矿车、 风门和密闭墙等.本文在对障碍 物的影响研究的同时,对风门和密闭墙也进行了相 关分析.可以把风门和密闭墙看作膜,则理论上,破 膜前,火焰面诱导的压力波反射后与火焰锋面相互 作用,使得火焰传播速度减小甚至滞止;但由于压 缩波后温度、 压力的增大,对燃烧又起到促进作用. 破膜后,稀疏波会大大加快火焰化学反应的传播 速度,但因稀疏波降低了甲烷混合物的温度、 压力 等参数,会导致燃烧速率下降,甚至造成火焰熄灭. 因此,如何使得上述过程匹配,保证稀疏波既加速 火焰锋面,又不至于造成火焰熄灭是我们研究问题 的关键所在. 1 数学模型 对于三维非定常、 黏性、 可压缩预混气体爆炸 问题可以建立如下数学模型[627] 5Qg 5t 5Eg 5x 5Fg 5y 5Gg 5z 5Egy 5x 5Fgy 5y 5Ggy 5z S, 其中S为气相反应源项 S - w1- ∑w 1kpk -w2- ∑w 2kph -wns- ∑w nskpk 0 0 0 wsQs - Ρij 5ui 5xj - Λi Θ 2 5Θ 5xj 5p 5xj - Θ Ε - C1 Ε k Ρij ui 5xi pi- C2ΘΕ 2 k . 数值计算时,采用隐式方法处理化学反应源 项,采用LU分解和迎风TVD格式分别处理对流 项的隐式和显式部分,扩散项采用中心差分. 2 计算结果及分析 初始条件初始压力p0为1个大气压,初始温 度T0为300 K,初始的马赫数为0;初始气体混合 物的温度、 压力和浓度场沿径向均匀分布;初始浓 度场分为两部分,即当1≤IX≤4181时,为按当 量比42混合的H2 ,O 2和N2,即H2O2激波点 火段;当4181≤IX≤M X时,为按当量比9. 5 混合的CH4 ,O 2和N2,其中IX, 41, 81和M X为无 量纲距离. 边界条件反应装置是左端封闭、 右端开口、 边 界光滑的方形钢管如图1所示,管子的抗压强度足 够大.沿壁面,气体取无滑移速度条件,温度取绝热 条件,气相压力法向梯度取为零.出口采用外推条 件,膜设为固壁条件.方管的长径比为15∶1,断面 为40 mm40 mm. H2 A IRCH4 A IR 图1 试验装置 Fig. 1 Experi mental device 由于壁面附近,流场结构变化剧烈,采用指数 加密方式处理,另外,为了减少计算量,取管道的一 半作为计算域图2. 图2 计算域网格 Fig. 2 The mesh of calculation 为了研究火焰和压力波的情况,我们选择方管 的一个面IY 21,IZ 11进行分析. 图3所示为nn 500时压力、 温度和流线图, 其中nn为无量刚时间步数.从图3中可见,在左端 面点火后,火焰阵面自左向右传播并诱导出压力 波,其超压为2个大气压左右. 图3 点火示意图 Fig. 3 Schematic diagram of ignition 波的实质是扰动的传播图4.点火后,波阵 424 中国矿业大学学报 第34卷 1994-2010 China Academic Journal Electronic Publishing House. All rights reserved. 面成为已受扰动影响的质点和未受扰动影响的质 点之间的 “分界面”,随着燃烧的不断进行,球形波 变成平面波. 图4 平面波示意图 Fig. 4 Schematic diagram of plane wave 从图5可见,压力波到达膜片后,在膜片上反 射为左行传播的压力波,随着波的不断传播,与滞 后的火焰面相互作用,产生压缩波和透射波,左行 的压力波使火焰速度减小,甚至滞止.当膜片处的 压力大于破膜压力时,膜片瞬时破裂,同时产生右 行激波和左行的稀疏波.稀疏波和后面的火焰面又 相互作用,使火焰加速. 图5 破膜过程参数变化示意图 Fig. 5 Schematic diagram of parameter change for the diaphragm bursting 火焰与障碍物相互作用的情况如图6所示,由 图可见,随着火焰与障碍物的不断接近,火焰锋面 开始变形,在障碍物的左侧由于壁面的反射使火焰 滞后,在障碍物的上面,特别是左上角处,火焰的传 播速度较快[8].随着时间的推移,火焰阵面从障碍 物上方越过,同时火焰阵面被拉长,在障碍物的两 侧附近能看到湍流漩涡的存在. 由于化学反应的延迟性,质点在激波压缩后, 经过一段诱导时间才能点火,紧随激波后面的是化 学反应诱导区[9].从图7中可见破膜后甲烷并没 有立即在膜处被点燃,而是经过一段距离后才被氢 氧燃烧的火焰点燃,且在障碍物的左上角最为剧 烈. 图6 火焰与障碍物作用示意图 Fig. 6 Schematic diagram interaction between flame and barrier 图7 CO组分变化图 Fig. 7 D iagram of species change for CO 3 结 论 1 瓦斯爆炸后,火焰阵面附近区域,障碍物左 上角和管子封闭端附近区域温度变化较为陡峭;火 焰逐渐接近障碍物时,温度上升很快,但由于管壁 的吸收,温度开始下降.总的来说,在火焰传播的管 道中设置的障碍物对气相火焰具有加速的作用,加 速的机理障碍物诱导的湍流区对燃烧过程的正反 524第4期 吴 兵等瓦斯爆炸运动火焰生成压力波的数值模拟 1994-2010 China Academic Journal Electronic Publishing House. All rights reserved. 馈造成的. 2 当火焰遇到障碍物时,火焰面发生变形,燃 烧面扩大,燃烧速度加快,火焰对前驱冲击波提供 的能量突然增加,使得前驱冲击波的压力和速度等 参数突然上升,出现火焰正的激励效应.另外由于 管道的截面突然减小,气流速度增大,湍流度增加. 火焰传播过程中的湍流效应是产生冲击波的主要 因素,爆炸火焰的传播速度的大小直接影响着爆炸 冲击波的生成和加强程度. 3 破膜后压力变化氢氧燃烧破膜后如果后 面没有甲烷化学反应的话,相当于氢氧产生的冲击 波在空气中传播,膜片后面产生左行稀疏波,传到 左壁面,反射的稀疏波会透过膜片传播过去,压力 逐渐的衰减的.如果后面有甲烷化学反应,会使冲 击波加强.另外破膜需要一个过程,同时损耗能量, 破膜后由于与前面的未扰动气体有压差,波速又会 增大,破膜压力的衰减是很小的.破膜压力越高,火 焰在管内传播的时间越长. 参考文献 [1] 余立新,孙文超,吴承康.半开口管道中的氢 空气火 焰加速和压力发展过程[J ].工程热物理学报, 2001, 225 6372640. Yu L X, SunW C, W u C K. Flame acceleration and pressure development of H22A ir in a sem i2open tube [J ]. Journalof Engineering Thermophysics, 2001, 22 5 6372640. [2] 高泰荫,黄军涛. CH42O2混合气中爆燃爆震转捩的数 值模拟[J ].爆炸与冲击, 1998, l184 3232330. Gao T Y, Huang J T. N umerical si mulation of defla2 gration to detonation transitionDDTof CH42O2gas m ixture [J ]1Explosion and Shock W aves, 1998, l18 4 3232330. [3] 林伯泉.瓦斯爆炸动力学特征参数的测定及其分析 [J ].煤炭学报, 2002, 272 1642167. L in B Q1The measurement and analysis of dynam ics feature parameter in gas explosion[J ]. Journalof Chi2 na Coal Society, 2002, 272 1642167. [4] 吴 兵.矿井半封闭空间瓦斯爆燃过程热动力学研究 [D ].北京中国矿业大学资源与安全工程学院, 2003. [5] 菅从光,林柏泉,翟 成.瓦斯爆炸过程中爆炸波的 传播规律[J ].中国矿业大学学报, 2003, 324 3632 366. Jian C G, L in B Q , ZhaiC. Structure variation law of explosion wave in gas explosion[J ]. Journal of China U niversity of M ining Technology, 2003, 32 4 3632366. [6] Sm irnovN N , Panfilov I I . Deflagration to detonation transitionin combustible gasm ixtures[J ]. Combustion and Flame, 1995101 912100. [7] U lrich B, M artin S. N umerical si mulation of prem ix2 ed combustion processes in closed tubes[J ]. Combus2 tion and Flame, 1998114 3972419. [8] Rankin D D, M ccann M A. Overpressures from non2 detonating,baffleac2celerated turbulentflamesin tubes[J ].Combustion and Flame, 2000 120 5042 514. [9] Clifford L J, M ilneA M. N umericalmodeling of che2 m istry and gas dynam ics during shock2induced ethy2 lene combustion [J ].Combustion and Flame,1996 104 3112327. 责任编辑 陈其泰 624 中国矿业大学学报 第34卷