计算流体力学在烟气脱硫中的应用

摘要：为优化脱硫吸收塔内流场分布，提高脱硫效率，降低脱硫投资和运行成本，提出基于RANS方程和多相流模型的烟气脱硫数值模拟方法，基于FLUENT采用多相流模型针对吸收塔内烟气脱硫过程进行模拟.通过对吸收塔烟气和浆液区域进行建模及网格划分，并对内部和喷枪局部处流场进行数值模拟，分析得到浆液和烟气在吸收塔内的流动规律.模拟结果验证该方法的可行性.通过对可能的浆液湿壁现象进行分析，给出吸收塔和喷枪设计的建议.

关键词：脱硫吸收塔；烟气脱硫；多相流；浆液；湿壁

中图分类号：TK229.6；TB115.1 文献标志码：B

作者简介： 祁江涛(1982―)，男，湖北京山人，工程师，硕士，研究方向为流体力学非线性数值计算，（Email)

Application of computational fluid dynamics in

flue gas desulfurization

QI Jiangtao1, LI Yinghua1, WU Haolun2

(1. Research Department for Whole Performance, China Ship Scientific Rearch Center, Wuxi 214082, Jiangsu, China;

2. Research and Development Department, Wuxi Huaxing Electric Power Environment Protection Engineering Co., Ltd.,

Wuxi 214082, Jiangsu, China)

Abstract： To optimize the flow field distribution in desulfurization absorption tower, improve desulfurization efficiency, and decrease desulfurization investment and operation cost, a numerical simulation method of flue gas desulfurization based on RANS equations and multiphase flow model is proposed, and the process of flue gas desulfurization in absorption tower is simulated based on multiphase flow model by FLUENT. The areas of flue gas and serosity in absorption tower are modeled and meshed, the flow field of the area inside absorption tower and near jet guns is simulated, and the flow pattern of serosity and flue gas inside absorption tower is obtained. The simulation results show the feasibility of the method. The suggestions on the design of absorption tower and jet guns are proposed by the analysis on wetted wall phenomenon.

Key words： desulfurization absorption tower; flue gas desulfurization; multiphase flow; serosity; wetted wall

0引言

烟气脱硫是目前能大规模控制燃煤造成SO2污染最为有效的方法之一，而石灰石―石膏湿法脱硫技术以其脱硫效率高、吸收剂来源丰富、成本低廉、技术成熟和运行可靠等优点获得广泛应用.从气液两相流体力学和化学反应动力学的观点看，脱硫吸收塔内流体流动的目的是强化气液两相的混合和质量传递、延长气液两相在塔内的接触时间、增大气液两相的接触面积并尽量减小吸收塔的阻力.合理的塔内流场分布对提高脱硫效率、降低脱硫投资和运行成本都具有重要意义.

目前，国内外对烟气脱硫吸收塔进行大量研究,主要采用实验方法，如研究塔的阻力特性、液滴运动速度沿塔高变化［1］和TCA塔内温度场分布［2］等,这些研究对指导工业应用具有重要意义,但其结果往往只针对特定的设备或结构,具有较大的局限性.随着计算机技术的迅速发展,计算流体力学(Computational Fluid Dynamic，CFD) ［3］已成为研究三维流动的重要方法：周山明等［4］利用FLUENT计算空塔和喷淋状态下的塔热态流场，结果表明脱硫吸收塔入口处流场变化最剧烈、压降损失最大，并根据计算结果改造来流烟道；孙克勤等［5］采用混合网格和随机颗粒生成模型对烟气脱硫吸收塔的热态流场进行数值模拟；郭瑞堂等［6］采用FLUENT结合非稳态反应传质-反应理论对湿法脱硫液柱冲击塔内的流场和SO2的吸收进行数值模拟.

本文尝试应用FLUENT对某脱硫吸收塔内烟气脱硫过程进行初步数值模拟，通过对内部流场进行分析验证本文模拟的合理性，进而对脱硫过程中脱硫吸收塔内是否存在湿壁现象进行深入分析研究.

1基于RANS求解器的CFD数值模拟方法1.1控制方程

时均的不可压缩连续性方程和NS方程（RANS方程）如下：uixi=0(1)uit+ujuixj=－1ρpxi+氮xjνuixj－u′iu′j(2)

式中：ui,uj和p都为时均量；氮xj(－u′iu′j)为雷诺应力项.

1.2湍流模型和多相流模型

RNG k拨磐牧髂P吞峁攵缘屠着凳有效黏性的微分解析式，具备数值稳定性好、求解压力梯度精确以及工程实用等优点，因此本文的数值计算采用RNG k拨磐牧髂P.多相流模型采用欧拉模型.

1.3边界条件

(1)入口边界条件：采用速度入口边界条件Vin=V0.

(2)出口边界条件：采用出流边界条件.

(3)物面条件：满足壁面黏附条件，壁面处流体速度与运动边界速度相同.

1.4数值离散和求解

(1)时间项的离散：采用直接1阶隐式离散.

(2)空间项的离散：扩散项以中心差分格式进行差分，对流项采用2阶迎风格式.

采用SIMPLE法处理压力-速度耦合问题，离散方程以GaussSeidel迭代法求解.

2烟气脱硫数值模拟

数值模拟对象为某个用于烟气脱硫的脱硫吸收塔，配有喷枪喷射浆液用于烟气脱硫，其脱硫过程涉及浆液对烟气中有害气体的吸收、浆液中Ca(OH)2与烟气中硫化物的化学反应以及浆液的蒸发.考虑到具体计算的时间问题以及实际问题的复杂程度，本文作相应简化，不考虑浆液中Ca(OH)2与烟气中硫化物的化学反应以及实际脱硫过程中的传热蒸发.在研究浆液湿壁问题时，本文从烟气及浆液的流动角度（速度分布）进行细致的分析研究，考虑到实

图 1吸收塔模型和

坐标系

Fig.1Absorptiontower

model and

coordinate system际传热蒸发对烟气湿壁具有很好的抑制作用，因此本文的分析结论偏于保守可靠.

2.1数值模拟对象

吸收塔入口处烟气速度为6 m/s，喷枪喷射浆液流速度为25 m/s.为便于分析，建立固连于塔体的坐标系，并约定：吸收塔对称面所在的面为xOy面，z轴垂直于xOy面并满足右手法则，吸收塔模型和坐标系见图1.采用非结构网格对吸收塔内计算域进行网格划分，网格数量约为200万个.

2.2流场分析

2.2.1脱硫吸收塔内部流场分析

(1)烟气入口处速度分布均匀，稳定在6 m/s左右；导流板处烟气分布较为均匀；烟气在进入喷嘴时，由于与高速浆液进行动量交换，速度迅速增加到20 m/s以上，并显示出如喷枪一样的火焰状喷射轨迹，脱硫吸收塔对称面内烟气速度云图见图2；脱硫吸收塔内广大区域速度较小，根据连续性方程，出口处速度应较大，从图2中亦能看到剧增的出口速度.(2)喷枪处浆液速度云图见图3，显示出与实际情况相符的火焰状轮廓，在进入脱硫吸收塔内后，由于与烟气混合在一起，脱硫吸收塔内的浆液速度分布几乎与烟气一致.(3)进一步将脱硫吸收塔内流线示意绘出，见图4，可知脱硫吸收塔内烟气流动非常复杂，烟气在刚进入脱硫吸收塔内时流动均匀；而后进入主塔体时形成一对反向旋转的涡对，左边的很大，右边的相比之下较小但强度很强，并分别向上和向下卷曲延伸.导流板处及其上表面均无涡旋，只在趋于主塔体处形成较强的一次涡和二次分离涡.

图 2对称面内烟气速度云图

Fig.2Gas velocity contour in symmetric plane

图 3对称面内浆液速度

云图

Fig.3Serosity velocity

contour in symmetric

plane图 4吸收塔内部流线示意

Fig.4Streamline inside absorption tower

2.2.2三喷枪及导流板局部流场分析

(1)浆液喷射轮廓大于烟气，接近实际喷枪，喷射角度近60°，参与烟气流动，在脱硫吸收塔内形成涡旋.

(2)3组喷枪所在区域的浆液与空气流场完全一致，选取对称面进行流线分析，绘出浆液空气流线图.因tecplot中流线是起点式绘制，而喷枪出口轮廓的起点皆为分离点，故起点处流线也各个分离，喷枪所在位置处速度云图和流线示意见图5.

(a)z=0处烟气速度云图(b)z=0处浆液速度云图(c)z=2.1处烟气速度云图(d)z=2.1处浆液速度云图(e)z=-2.1处烟气速度云图(f)z=-2.1处浆液速度云图(g)z=0处烟气流线(h)z=0处浆液流线图 5喷枪所在位置处速度云图和流线示意

Fig.5Velocity contour and streamline at jet gun location

(3)导流板处速度分布均匀且未形成涡旋，表明设计的导流板导流效果显著.

综上所述，从脱硫吸收塔整体及喷枪局部速度云图和流线可知，FLUENT具备较好的模拟脱硫吸收塔在脱硫过程中烟气和浆液流动的能力.

2.3湿壁情况分析

进一步选取较小速度比例给出壁面浆液速度云图，以分析可能的浆液湿壁现象，图6~9为中间喷枪和左右2个喷枪在yOz面内和相应xOz面内的速度云图.由图6和7可知，左右2个喷枪喷出的浆液流贴近近壁面.进一步截取浆液速度剖面进行分析，在近最大轮廓面及以下1.5 m处截取剖面并精细显示速度云图，见图8和9.可知，左右2个喷枪喷射的浆液流速度在近壁面处达到3~6 m/s，而从整体流线图看出塔内存在涡旋，考虑到实际喷枪出流速度大于25 m/s时将会使火焰喷射轮廓更大进而射到塔壁上，因此该脱硫吸收塔的设计方案有可能出现湿壁现象.建议将左右喷枪挪向中部，并收缩喷枪所在管道直径，增加喷枪数量.

图 62个喷枪所在yOz面浆液速度云图

Fig.6Serosity velocity

contour in yOz

plane of two jet

guns图 7第3喷枪所在yOz面处浆液速度云图

Fig.7Serosity velocity

contour in yOz

plane of the third

jet guns图 8近最大轮廓面处

浆液速度云图

Fig.8Serosity velocity

contour near

maximum

profile plane

图 9近最大轮廓面1.5 m

处浆液速度云图

Fig.9Serosity velocity

contour at location

that is 1.5 m from

maximum profile

plane3结论

基于商业软件FLUENT采用多相流模型针对脱硫吸收塔内烟气脱硫过程进行初步模拟，建立基于RANS方程的烟气脱硫多相流数值模拟方法.通过对脱硫吸收塔内部和喷枪局部处流场进行分析，模拟得到浆液和烟气在脱硫吸收塔内的流动规律.从模拟结果看，本文建立的方法可行.通过对可能的浆液湿壁现象进行分析，认为此种脱硫吸收塔的设计方案有可能出现湿壁现象，建议将左右喷枪挪向中部，并且收缩喷枪所在管道直径，增加喷枪数量.下一步工作将在给定的喷枪速度范围内进行系列计算，综合比较给定湿壁影响区域，为施工设计提供参考；同时，也考虑在进一步的深入计算中加入传热蒸发模型.

参考文献：

［1］黄跃. 蒸汽透瓶中自发凝结及流动特性的理论和试验研究［D］. 西安: 西安交通大学, 1987.

［2］李亮, 丰镇平, 李国君. 一维喷管中存在自发凝结的跨音速湿蒸汽两相流动数值模拟［J］. 工程热物理学报, 2001, 22(6): 703705.

LI Liang, FENG Zhenping, LI Guojun. Numerical simulation of two phase wet steam transonic flow with spontaneous condensation in a 1D nozzle［J］. J Eng Thermophysics, 2001, 22(6): 703705.

［3］YOUNG J B. The spontaneous condensation of steam in supersonic nozzles［J］. Phys Chem Hydrodynamics, 1982, 228(2): 243274.

［4］周山明, 金保升, 仲兆平, 等. 大型烟气脱硫塔的流体动力学模拟及优化设计［J］. 东南大学学报:自然科学版, 2005, 35(1): 105110.

ZHOU Shanming, JIN Baosheng, ZHONG Zhaoping, et al. Flow simulation for large FGD scrubbers and optimization design［J］. J Southeast Univ: Nat Sci, 2005, 35(1): 105110.

［5］孙克勤, 周山明, 仲兆平, 等. 大型烟气脱硫塔的流体动力学模拟及优化设计［J］. 热能动力工程. 2005, 20(3): 270274.

SUN Keqin, ZHOU Shanming, ZHONG Zhaoping, et al. Hydrodynamic simulation and optimized design of a large sized flue gas desulfurizer［J］. J Eng Therm Energy & Power, 2005, 20(3): 270274.

［6］郭瑞堂, 高翔, 王君, 等. 液柱塔内流场和SO2吸收的CFD模拟和优化［J］. 浙江大学学报: 工学版, 2007, 41(3): 494498.

GUO Ruitang, GAO Xiang, WANG Jun, et al. CFD simulation and optimization of flow field and SO2 absorption in impinging stream scrubber［J］. J Zhejiang Univ: Eng Sci, 2007, 41(3): 494498.