封闭管道内气体燃烧火焰传播特性的研究.docVIP

封闭管道内气体燃烧火焰传播特性的研究 　　【摘 要】 为了研究可燃气体在管道内的燃烧规律，选取甲烷与空气预混气体为可燃气体，建立内置障碍物情况不同的封闭管道模型，运用CFD计算流体动力学方法，对管道内气体燃烧过程进行了数值模拟。监测了火焰面的形状、位置、速度的变化;分析了挡板位置、数量以及阻塞率对火焰面传播过程的影响。研究表明：火焰面传播速度随挡板数量、阻塞率的增加而增大，挡板和点火区的距离与管径的比值在0.5～3范围内时，火焰面传播速度达到最大值。研究为工程结构的防爆设计提供依据。 　　【关键词】 封闭管道 气体燃烧 障碍物 计算流体动力学 火焰传播 　　1 引言 　　2014年6月3日，重庆市万圣经济开发区砚石台煤矿发生重大瓦斯爆炸事故，造成22人死亡;8月1日，台湾高雄前镇区发生丙烯气体泄漏并爆炸，造成28人死亡，302人受伤。无论是工业生产，还是日常生活中，可燃气体安全问题已引起人们的高度关注。 　　气体爆炸问题的研究方法，主要有解析法，实验方法和数值方法。解析法中，常用的解析模型有等温模型、绝热模型和一般模型[1]。模型对实际情况进行了较大的简化，分析存在一定的误差。相比实验方法，数值方法成本低，安全性高，已经被许多研究人员采用。采用数值方法主要有集中参数方法（LMP）和计算流体动力学方法（CFD）。CFD方法能够反映变量在空间的分布，可以适应形状复杂的情况，得到了广泛的应用。运用CFD方法，王志荣等[2]模拟了连通容器内气体的爆炸，结果表明连通管后的容器内产生预压缩，湍流导致爆炸压力增加;李帆等[3]模拟了光滑管道内气体的爆炸，得到了不同时刻管内压力、湍流动能的分布;毕明树等[4]对管道内可燃气爆炸进行了一维数值模拟。 　　实际上，输气管道在使用过程中，可燃气中含有的杂质将在管道内沉积、聚集，导致管道流通面积减小。在长径比L/D较大（L/D5）且内部有较密集的隔板和设备等障碍物的密封体内，气体燃烧过程更加复杂[1、5-6]。为更好的对管道及其围护结构进行抗爆设计，有必要对不同阻塞率下管道内火焰传播特性进行研究。 　　本文利用计算流体动力学软件Fluent对内径为20mm，长300mm（L/D=15）的管道内甲烷与空气预混气体燃烧过程进行数值模拟。 　　2 数值模型 　　建立内径为20mm，长为300mm的管道二维模型，网格大小为1mm×1mm，共约6000个网格。四周边界条件为固定壁面，无滑移。由于气体爆炸发生时间较短，故在模拟过程中忽略热损失，假设壁面为绝热。管道内充满甲烷与空气的混合气体，甲烷含量为10%，混合气体温度为300K。化学反应方程式为： 　　以管道左端面中心为圆点，设置半径为3mm的半圆形点火区，点火区温度为1800K。如图1中，红色区域所示。 　　本文共建立7个管道模型。0号为内壁光滑管道模型，管内无挡板。其他模型，分别在管道内距离点火区l处，设置1道或2道挡板，对管道形成一定的阻塞。管道截面如图1所示。具体情况见表1。 　　3 计算结果与分析 　　3.1 内壁光滑管道内气云燃烧模拟 　　内壁光滑管道内气体模拟燃烧过程中，火焰面位置及形状变化如图2所示。点火初期，火焰面近似成半圆形向外扩散，由于侧壁温度与环境温度相同，为300K，火焰面沿轴向延伸，距离壁面较近处还未燃烧。0.5ms时刻，已燃区域逐渐饱满，火焰面成平面，化学反应主要集中于已燃区域右侧。 　　由于湍流作用，火焰面有一定的厚度，继续向前传播。2ms时刻，火焰面形成指向已燃区域的尖端，形成“郁金香”火焰面，火焰发展为明显的湍流，火焰面厚度进一步增加。此后，火焰面大致维持该形状，逐渐向右传播，至所有气体燃烧完毕。 　　每隔0.1ms在FLUENT中监测管道中心线上火焰面中心位置，得到火焰面位置与时间的关系曲线。对该曲线进行拟合并求导，可得火焰面中心的速度与时间关系。火焰面的位置、速度与时间关系如图3所示。 　　点火初期，点火温度较高，由阿累尼乌斯方程[1]可知，化学反应速率较大，火焰面以较快速度向前传播，随着燃烧的进行，火焰面速度逐渐减小。当火焰面以稳定的“郁金香”形向前传播时，一方面，因为监测点的局限性以及湍流的存在，造成监测数据不能准确反映整个火焰面的位置，引起一定的误差;另一方面，由于受到侧壁及管道右端面反射的超压影响，传播速度呈现一定的波动性。 　　3.2 阻塞率对火焰传播速度的影响 　　分别监测0、1、2号模型火焰面中心位置。记录火焰面中心速度与时间关系，如图4所示。在有挡板的管道中，火焰面传播速度在壁面反射超压的作用下依然有波动，最大速度与平均速度均远大于光滑管道。挡板对火焰面传播有很大影响。 　　由于挡板的存在，挡板左侧形成燃烧腔，腔内超压值较大，挡板处形成超压出口，造成气体喷