管道内甲烷_空气预混爆炸燃烧的数值模拟fluent

速率的提高使反应释热速率增大 ,导致了燃烧产物的进一步
膨胀加速 ,压力增大 ;压力的增大又推动了火焰面前方可燃
·29 ·
气体的流动 ,导致了火焰面附近更高的湍流动能和燃烧速 率 ,在这种正反馈机制作用下 ,最终使火焰加速传播 。
图 1 不同时刻管道内温度分布
2. 2 不同时刻管道内压力分布
图 2 为不同时刻管道内压力分布图 ,可以清晰地看出管 道内气体爆燃过程中的压力变化情况 。
1. 3. 2 边界条件
在模拟气体爆炸过程时 ,二维管道结构需设置 4 种类型
的边界条件 ,即对称 、进口 、出口 、绝热壁面 。当管道为密闭
Fra Baidu bibliotek
空间时 ,可以不设进口和出口条件 。
对称边界条件 :
v =0,
9vx 9x
=
0
,
9Yi 9x
=
0
,
9T 9x
=
0
。
壁面边界条件 :
由于气体爆炸在瞬间即告完成 ,故在模拟时可以忽略热
·30 ·
工业安全与环保 Industrial Safety and Environmental
Protection
2A0p0r9il
年第 2009
35
卷第
4
期
氯乙烯聚合反应器爆炸危险性分析
郭尹亮
(东北大学资源与土木工程学院 沈阳 110004)
摘 要 通过故障树方法辨识出可能导致氯乙烯聚合反应器发生爆炸的基本事件 ,及其与聚合釜爆炸之间的逻辑关 系 。利用危险与可操作性研究方法分析了工艺中可能导致聚合釜爆炸的工艺偏离及其产生的原因和后果 ,明晰了工艺危险 性的主要控制因素 。针对聚氯乙烯生产工艺提出安全控制措施 。 关键词 聚氯乙烯 聚合釜 爆炸 工艺危险性
图 2 不同时刻管道内压力分布
对比图 1 和图 2 可以看出 ,不同时刻管道内最大压力集 中分布于火焰面前方一个很薄的区域内 ,形成了压力波 。
2. 3 不同时刻管道内速度分布
管道内预混可燃气体被点火源点火后 ,在最初的时间内 (5 ms) ,可燃气体由于点火源的加热而开始发生膨胀运动 , 此时最高速度为 18 m/ s 左右 。经过一段时间的运动后 (21. 4 ms) ,火焰速度增加 ,该时刻管道中部气体流动速度较快 ,气 体流动加快了火焰向前传播 ,顶部火焰传播速度增加较快 , 而壁面处气体流动较慢 ,火焰传播速度增加较慢 ,因此 ,火焰 在轴向拉伸较长 ,此时轴线最大速度可达 36. 3 m/ s。随着火 焰的继续发展 (74. 4 ms) ,受到密闭管道管壁的约束 ,火焰阵 面两侧气体出现反向流动 ,火焰面后端气体出现与火焰传播 方向的反向流动 ,而前端气体流动与火焰传播方向相同 ,从 而使火焰传播速度有所降低 ,此时管道内最大速度降为 13 m/ s。随着火焰传播的继续进行 (254. 4 ms) ,管道尾端管壁 对气体流动的障碍作用越来越大 ,火焰传播速度持续下降 , 管道轴线上的速度只有 5 m/ s。
P ρ
;
h
=
Yi hi ; hi = CP, jd T ; Tref = 298. 15
i
Tref
→
K; Yi 为组分 i 的质量分数 ; keff 为有效导热系数 ; J i 为组分 i
的扩散准热通量 ; Sh 为化学反应生成热 。
湍流方程 :
在标准 k - ε模型中 , k 和ε是 2 个基本未知量 , 与之相
0 前言
现代工业爆炸事故中 ,相当一部分是由管内悬浮流燃烧 引起的 。目前常用的解析模型有等温爆炸模型 、绝热爆炸模 型和一般模型[1] ,它们均假设整个容器内压力在任何时刻都 均匀分布 ,没有考虑流体动力学方程 ,也没有考虑化学反应 动力学问题 。随着计算机技术的发展 ,数值模拟方法能考虑 更多的问题 ,已成为求解爆炸这类复杂问题的有效手段 。在 此 ,利用描述爆燃过程的流体力学基本方程组和燃烧反应方 程 ,运用差分方法 ,对管道内可燃气体的爆燃进行数值模拟 。
图 3 不同时刻管道中火焰传播速度 - 距离分布
从图 3 中可以看到 ,火焰开始阶段速度增加较快 ,通过 管道 1/ 3 后 ,火焰开始减速 ,火焰速度在管道的 1/ 3 附近达 到最大值 ,继而火焰传播速度开始下降 ,这是因为管道为密 闭空间 ,火焰向右传播过程中 , 火焰面前方的压力波挤压火 焰面右边的未燃混合气体 ,在管道尾端受到管壁的阻碍 ,并 逐渐产生了与火焰传播方向相反的气流 ,从而降低了火焰的 传播速度 。
1. 2 控制方程
以 CFD (Computational Fluid Dynamics) 为基础的数值模拟 方法基本的求解思路 ,是建立控制方程组 ,这是求解任何问 题前都必须首先进行的 。由流体力学 、热力学 、传热传质学 、 燃烧学等基本原理出发 ,建立质量 、动量 、能量 、组分 、湍流特
碍和管端的约束作用 ,火焰改变方向与已燃区内的已燃气体
相作用 ,使已燃区内发生湍流 (10 ms) 。壁面的约束加剧了
湍流的产生 ,在湍流旋涡的作用下 ,火焰开始发生变形 、拉
伸 ,火焰进入加速传播阶段后 (20 - 40 ms) ,火焰面的变大使
更多的可燃气体参与反应 ,增大了燃烧反应速率 ,燃烧反应
定义源项 。
1. 3 边界条件
对于两端密封长度为 L = 790 mm、管道内体积分数 10 %
的甲烷2空气混合物爆炸 ,采用二维湍流燃烧模型 。
1. 3. 1 流场初始化条件
管道密闭 ,0 ≤x ≤L ,0 ≤r ≤R ,流场初始化条件如下 :
反应方程式 :CH4 + 2O2 CO2 + 2H2O
应的输运方程为 :
99ρt +
9 (ρkvx) 9x
=
9 9x
[
(μ
+
μi σk
)
9k 9r
]
+
Gk
+
Gb
-
ρε-
YM
+ Sk
(4)
3 项目资助 :山西省留学人员科研资助项目 。
9 (ρε) 9t
+
9 (ρεvr) 9r
=
9 9r
[
(μ+
μt σε
)
99εx ] +
ε C1ε k ( Gk +
C3εGb)
99ρt +
9 9xj
(ρuj
)
=0
(1)
动量守恒方程
:
9 9t
(ρui )
+
9 9xj
(ρuj ui )
=
-
9p 9xi
+
9τij 9xj
(2)
能量守恒方程 : 99ρt (ρE) + ·[ v (ρE +ρ) ] = ·( keff
∑→
T-
hiJ i) + Sh
(3)
i
T
∑ ∫ 式中 , E
=
h-
Numerical Simulation of Detonation and Combustion of Methane - air Premix in Pipeline LI Fan TAN Ying - xin
( School of Chemical Engineering and Environmental , North University of China Taiyuan 030051) Abstract In the actual production process , it is easy to cause the gaseous fuel and the air mixing , and to achieve the inflammable density in many situations. If there is a little ignition source existed , it possibly leads to combustion , explosion hazards , which could cause serious property damage and human casualties. In this paper , the methane is taken as the object of the study ; the air is selected as the oxidized gas. The numerical simulation to the mixture has been carried out . The conflagration pressure field , the temperature field and the density field changing data in the different time have been obtained. The theoretical basis and the data for guarding against explosions , the explosion sup2 pression for project have been provided. Keywords methane gas deflagration numerical simulation computational fluid dynamics
=0,
9p (0 , t) 9x
=
0,
9ρ(0 , t) 9x
=0,
9T(0 , t) 9x
=0,
9Yi (0 , t) 9x
= 0。
2 结果与分析
2. 1 不同时刻管道内火焰温度分布
图 1 为点火后不同时刻管道内火焰的温度场分布图 ,可
以较为直观地看到气体爆炸的发展变化过程 。
从图 1 可见 ,当火焰渐渐接近壁面后 ,由于受到壁面阻
摘 要 在实际生产过程中 ,有很多情况能使气体燃料与空气混合 ,达到可燃浓度 ,此时若有点火源存在 ,就可能酿成 燃烧 、爆炸灾害 ,造成严重的财产损失和人员伤亡 。选取甲烷作为研究对象 ,选用空气作为氧化性气体 ,对其进行数值模拟 。 获得了不同时刻爆燃的压力场 、温度场 、密度场的变化数据 ,为工程上防爆 、抑爆 、泄爆提供了理论基础和数据 。 关键词 甲烷气体 爆燃 数值模拟 计算流体动力学
相对于实验研究 ,数值模拟具有研究成本低 ,周期短 ,无 实验仪器干扰 ,能获得完整的数据 ,能将计算情况在计算机 屏幕形象地再现等独特的优点[2] 。 1 物理数学模型
1. 1 物理模型
建立 1 个合理的物理模型 ,是解决流体力学问题的重要 环节 。整个爆炸过程可以看作是由点燃开始 ,一层层新混合 气依次点火 ,并向未燃混合气传播的过程 。这样整个爆炸场 就形成了 3 个区域 :燃烧区 、火焰面 、未燃区 。
(6)
初始时刻 t ( t0) = 0 ;初始速度 v ( t0) = 0 。
已燃区 : T = 2 000 K, P0 = 101 325 Pa , YCH4 = 0 , YO2 = 0 , YCO2 = 0. 151 4 , YH2O = 0. 123 9 。
未燃区 : T = 300 K, P0 = 101 325 Pa , YCH4 = 0. 055 , YO2 = 0. 22 , YCO2 = 0 , YH2O = 0 。
损失 ,因 此 假 设 壁 面 为 绝 热 。v ( R , t )
= 0,
9ρ( R , t) 9x
= 0,
9T( R , t) 9x
= 0,
9Yi ( R , t) 9x
= 0;绝热壁面热流密度为零 ,即
qw
= 0 ,壁面无滑移 。
管道两端边界条件 :
x =0和
x = 800
mm ,0 ≤t ≤∞处 , v (0 , t)
性等守恒方程组 。这些方程构建成非线性偏微分方程组 ,不
用经典的分析法 ,而只能用数值方程求解 。单相层流动的基
本方程组已经很少存在争议 ,对湍流 ,由不同的模拟理论出
发 ,往往基本守恒方程组也不相同 ,因此 ,构造基本方程组 ,
就成为模拟理论的重要部分 。常用的爆燃气体动力学方程
组为[3 ] :
质量守恒方程 :
-
C2ερεk2 + Sε
(5)
式中 , Gk 为平均速度梯度下的湍流动能产生项 ; Gb 为浮力
湍流动能产生项 ,对不可压流体 , Gb = 0 ; YM 代表可压湍流中
的脉动膨胀贡献 , YM 、C1ε、C2ε和 C3ε为常数 ;σk 和σε 分别为
湍流动能 k 和耗散率ε对应的 Prandtl 数 ; Sk 和 Sε 为用户自
·28 ·
工业安全与环保 Industrial Safety and Environmental
Protection
2A0p0r9il
年第 2009
35
卷第
4
期
安全工程
管道内甲烷 - 空气预混爆炸燃烧的数值模拟 3
李帆 谭迎新
(中北大学化工与环境学院 太原 030051)