手机版
           

Fluent气体扩散仿真还原密闭车间H₂S泄漏扩散路径的完整建模

发布时间:2026-06-19   来源:科研学术网    
字号:

Fluent气体扩散仿真在这次脱硫车间泄漏事故复盘中,把”毒气怎么扩散的”这个问题从猜测变成了可追溯的浓度场演化记录。一台脱硫塔底部法兰泄漏硫化氢,泄漏速率约4.2 g/s,持续15分钟。事故未造成人员伤亡,但车间内三个区域的可燃气体探测器报警时序与直觉判断严重不符。Fluent气体扩散仿真的任务就是解释这个”不符”。

事故背景与仿真目标

事发车间尺寸24 m × 16 m × 8 m,内部布置脱硫塔、循环泵、控制柜及两套机械通风系统。顶部4台排风风机总风量12000 m³/h,底部2台补风风机总风量8000 m³/h。H₂S密度1.19 kg/m³(标况),IDLH(立即危及生命和健康浓度)为100 ppm,8小时容许浓度10 ppm。

事故时通风系统运行正常,但控制柜区域的探测器在第7分钟率先报警,而距离泄漏源更近的循环泵区域探测器在第9分钟才触发。运维团队的第一反应是”探测器坏了”。仿真需要回答:这个报警时序是否合理?

几何处理与网格划分

几何模型基于车间竣工图重建,保留脱硫塔、管道、泵体、控制柜、电缆桥架等主要障碍物,因为这些障碍物对流场结构的塑造作用不可忽略。通风口、门缝、泄漏法兰位置按实际尺寸建模。

网格策略认定混合网格更适合这类复杂几何。车间主体空间采用hexcore六面体核心网格,基础尺寸200 mm;障碍物表面及壁面附近生成5层边界层棱柱网格,首层厚度8 mm,增长比1.2,保证近壁y+在30-60范围。泄漏法兰口局部加密至10 mm。最终网格约156万单元,最大skewness 0.87,满足Fluent质量要求。

网格无关性验证用了三套网格:92万、156万、240万。监测点A(控制柜区域,距泄漏源水平距离11 m)的H₂S浓度在156万与240万网格下偏差小于3.5%,认定156万网格足够。

物理模型与边界条件设置

湍流模型选用RNG k-ε。理由是这个场景的流动以中等雷诺数湍流为主,存在回流和弯曲流线,RNG k-ε对这类流动的预测优于标准k-ε,计算成本又低于RSM。近壁处理采用标准壁面函数 [1]。

species transport模型启用,H₂S定义为第二相,空气为carrier。不启用化学反应——H₂S在空气中常温下不发生反应。Schmidt数取0.7,湍流Schmidt数0.7。

边界条件如下:

  • 泄漏口:质量流量入口,H₂S纯气体,4.2 g/s,温度25°C,水力直径15 mm
  • 排风口:压力出口,表压-120 Pa
  • 补风口:压力入口,表压0 Pa,温度25°C
  • 壁面:无滑移,绝热
  • 门缝:压力出口,表压0 Pa,模拟车间非完全密闭

泄漏口湍流参数按湍流强度5%、水力直径计算。

收敛判断与瞬态计算

这个仿真认定稳态先跑通流场、再用瞬态追踪扩散路径是更稳妥的策略。稳态计算迭代至连续性方程残差降至1×10⁻⁴以下,k和ε残差1×10⁻⁵,各速度分量残差1×10⁻⁴。同时监测排风口H₂S质量流量,待其趋于稳定后认定稳态收敛。

瞬态计算时间步长0.5 s,总时长900 s(15分钟)。每个时间步最大20次迭代。收敛判据为每个时间步残差降至1×10⁻⁴以下且监测点浓度连续两个时间步变化率小于1%。

时间步长选择经过敏感性分析:0.5 s与0.25 s在监测点A的浓度到达时间偏差小于4%,0.5 s足够。

扩散路径与报警时序验证

仿真结果直接回答了最初的疑问。H₂S从泄漏口喷出后,初始动量使其向上偏移约0.6 m,随后被顶部排风气流捕获,沿车间顶部向远端扩散。在到达控制柜区域前,气流遇到2.4 m高度的电缆桥架发生偏转,部分H₂S被”引导”向下沉降,在控制柜区域形成局部高浓度区。控制柜区域监测点A在第412秒(约6.9分钟)浓度达到100 ppm,与实际报警时间第7分钟高度吻合。

循环泵区域虽然距离泄漏源更近(4.5 m),但处于脱硫塔尾流屏蔽区,气流速度低,H₂S需要绕过塔体才能到达,浓度上升缓慢,第547秒(约9.1分钟)才达到100 ppm,同样与实际报警时序吻合。

差距不会说谎。仿真复现的报警时序证明探测器工作正常,”近源先报”的直觉判断在复杂流场中并不成立。Fluent气体扩散仿真揭示了一个被忽视的扩散路径:电缆桥架的导流效应使远端区域反而先达到危险浓度。

更值得警醒的是,车间内还存在一个通风死角区域(西北角,靠近控制柜后方),仿真显示该区域在第9分钟H₂S浓度达到730 ppm,是IDLH的7.3倍,而该区域未安装探测器。如果泄漏持续时间更长或人员进入该区域,后果将远比实际事故严重。

仿真结论提交后,企业在西北角通风死角增设2台探测器,并在电缆桥架下方加装导流板破坏偏转效应。改造后重新仿真验证,控制柜区域浓度达100 ppm的时间推迟至第11分钟,为应急响应争取了额外4分钟。

这次Fluent气体扩散仿真最有价值的输出不是浓度云图本身,而是它让运维团队放弃了对”探测器可能坏了”的执念,转而认真对待流场结构塑造扩散路径的事实。直觉在复杂湍流场中是靠不住的,仿真是把直觉换成物理的唯一途径。

图说天下

×
abaqus仿真
ansys仿真
comsol仿真
fluent仿真
力学仿真
多相流仿真
流体/流动仿真