数值模拟报告.doc

密闭空间内甲烷气体的爆炸数值模拟 陈婷婷2220130040 季晓林 2120130301 摘要：密闭受限空间中可燃气体的爆炸研究对于石油及天然气工业的安全生产具有重要意义。以RNGK-ε湍流模型基础,建立了可燃气体单步化学反应湍流爆炸模型, 以有限体积法求解爆炸流动及反应控制方程,从而对二维受限空间中可燃气体爆炸的过程及规律进行了数值模拟,模拟结果与实验数据有着较好的吻合性。所做的工作为受限空间中可燃气体爆炸特性及规律的进一步研究及工业防爆抑爆技术的工艺实施、系统设计和关键参数计算提供了理论依据。 关键词：湍流模型 二维受限空间 有限体积法 数值模拟 1引言 工业上由气体爆炸引起的事故屡见不鲜，往往造成重大的经济损失和严重的人员伤亡。气体爆炸常发生在密闭空间或者起始阶段发生在密闭空间，因为泄露的可燃气体在封闭场所容易形成可燃气云。密闭空间的长度与直径比( L /D) 较大时，火焰传播与压力变化过程更加复杂。因此研究和分析L /D 较大的容器内气体爆炸具有非常重要的实际意义。 本文对内径100 mm，长3000 mm(L/D=)的容器内甲烷－空气爆炸过程进行了数值模拟。数值模型采用RNG k- ε方法计算湍流，采用基于梯度方法对燃烧过程建立模型。模型中考虑了甲烷浓度、湍流、温度和压力等因素对燃烧速率的影响。模拟计算得到的爆炸最大压力与爆炸时间与实验结果吻合，验证了数值结果的有效性。数值结果揭示了爆炸过程中火焰的传播规律、流场特性。 2 数学模型及数值方法 基于实验研究和分析，管道中的可燃气体爆炸是典型的湍流爆炸。其本质是一带压力波的高湍流度、高反应速率的燃烧过程。爆炸过程不仅存在一般湍流燃烧的影响因素，还有爆炸过程所特有的高反应速率特性以及压力波的传播、压力波与火焰的正反馈机制，基于此，本文将建立管道中可燃气体爆炸过程的湍流爆炸理论分析模型。鉴于湍流爆炸模型的复杂性, 将该模型分为描述流场流动的湍流模型和基于该湍流流场的爆炸燃烧模型进行研究。采用基于总能方程的RNG k- ε湍流模型模拟爆炸过程中的湍流流动。 2.1 可燃气体爆炸流场湍流模型的建立 以甲烷气体为爆炸介质, 其爆炸单步化学反应为: CH4+2O2= CO2+ 2H2O+ Q， Q= 880 kJ/ mol 为反应释放的能量。 对于单步不可逆化学反应，采用阿仑尼乌斯（Ahrrenius） 兼顾湍流模拟精度和计算量, 本文采用两方程湍流模型为流场湍流模型。由于RNG k- ε模型具有严格统计分析得到的系数和有效粘度的微分格式的特点, 使其在广泛的流动范围内比标准k- ε模型精确, 并且RNG k- ε 模型在快速应变流的计算中也具有一致认可的良好精度 。因此, 以RNG k- ε模型作为爆炸流场的湍流模型。采用总能形式的能量方程能较好地描述高速流动, 而且在求解时可以利用矢通分裂方法提高压力计算的精度, 因此, 采用总能E= e+ ( 1/ 2) ui uj 作为能量的度量。 气体爆燃数学上为理想气体由于热导致膨胀流动 基本方程： 质量守恒方程： 能量守恒方程： 动量守恒方程： 湍流方程： 湍流动能 耗散率 燃料质量分数： 2.2 物理模型 如图1 所示, 受限空间为封闭空间, 其直径为100 mm，长度为3 m，整个空间内充满了甲烷- 空气的混合物。在管道封闭端有一高温启爆源, 以直接诱导爆炸。 图1 模型几何结构 采用与x , y 坐标平行的正交网格对整个计算域进行网格划分。 边界条件: 壁面采用无滑移、无传热的边界条件。采用壁面函数法补充壁面区域的流动条件。由于非平衡壁面函数法考虑了压力梯度的影响和边界区域k 和ε的产生与消耗, 在有压力梯度存在的流场模拟中仍能得出正确的模拟结果, 因此采用非平衡壁面函数法处理壁面边界条件。 初始条件: 启爆源T0= 1.5kK, 高温火源诱导爆炸其他区域: p 0= 0, T0= 300 K, 其中压力为101. 325 kPa惰化区的惰气按不同工况设为不同的值。( FVM)求解所建立的理论模型, 各控制方程中的粘性通量的处理直接利用中心差分格式: φf = ( 1/ 2)( φ0+ φ1) + ( 1/ 2) (Λφ0* rρ0+Λφ0* rρ1) , 中心差分格式离散粘通量可以使粘通量的计算保持二阶精度。 无粘通量处理较为复杂, 利用矢通分裂法进行处理， Ff=1/2( FR+ FL) —1/2Γ| A | ( QR- QL) 其中: | A | = M| | M- 1 3 数值结果分析 如表一所示 时间 （ms）））））））））））））—2768K之间。 图2为不同时刻管道内压力分布图,可以清晰地看出管道内