０　引　言

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

目前，国内外对烟气脱硫吸收塔进行大量研究，主要采用实验方法，如研究塔的阻力特性、液滴运动速度沿塔高变化和ＴＣＡ塔内温度场分布等，这些研究对指导工业应用具有重要意义，但其结果往往只针对特定的设备或结构，具有较大的局限性．随着计算机技术的迅速发展，计算流体力学（ＣｏｍｐｕｔａｔｉｏｎａｌＦｌｕｉｄＤｙｎａｍｉｃ，ＣＦＤ）已成为研究三维流动的重要方法：周山明等［４］利用ＦＬＵＥＮＴ计算空塔和喷淋状态下的塔热态流场，结果表明脱硫吸收塔入口处流场变化最剧烈、压降损失最大，并根据计算结果改造来流烟道；孙克勤等采用混合网格和随机颗粒生成模型对烟气脱硫吸收塔的热态流场进行数值模拟；郭瑞堂等采用ＦＬＵＥＮＴ结合非稳态反应传质－反应理论对湿法脱硫液柱冲击塔内的流场和ＳＯ_２的吸收进行数值模拟．

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

１　基于ＲＡＮＳ求解器的ＣＦＤ数值模拟
方法
１．１　控制方程

时均的不可压缩连续性方程和ＮＳ方程
（ＲＡＮＳ方程）如下：

１．２　湍流模型和多相流模型

ＲＮＧｋ_-ε湍流模型提供针对低雷诺数有效黏性的微分解析式，具备数值稳定性好、求解压力梯度精确以及工程实用等优点，因此本文的数值计算采用ＲＮＧｋ_-ε湍流模型．多相流模型采用欧拉模型．

１．３　边界条件
（１）入口边界条件：采用速度入口边界条件Ｖ_ｉｎ＝Ｖ_０．
（２）出口边界条件：采用出流边界条件．
（３）物面条件：满足壁面黏附条件，壁面处流体速度与运动边界速度相同．

１．４　数值离散和求解
（１）时间项的离散：采用直接１阶隐式离散．
（２）空间项的离散：扩散项以中心差分格式进行差分，对流项采用２阶迎风格式．
采用ＳＩＭＰＬＥ法处理压力－速度耦合问题，离散方程以Ｇａｕｓｓ-Ｓｅｉｄｅｌ迭代法求解．

２　烟气脱硫数值模拟

数值模拟对象为某个用于烟气脱硫的脱硫吸收塔，配有喷枪喷射浆液用于烟气脱硫，其脱硫过程涉及浆液对烟气中有害气体的吸收、浆液中Ｃａ（ＯＨ）_２与烟气中硫化物的化学反应以及浆液的蒸发．考虑到具体计算的时间问题以及实际问题的复杂程度，本文作相应简化，不考虑浆液中Ｃａ（ＯＨ）_２与烟气中硫化物的化学反应以及实际脱硫过程中的传热蒸发．在研究浆液湿壁问题时，本文从烟气及浆液的流
动角度（速度分布）进行细致的分析研究，考虑到实际传热蒸发对烟气湿壁具有
很好的抑制作用，因此本文的分析结论偏于保守可靠．

２．１　数值模拟对象

吸收塔入口处烟气速度为６ｍ／ｓ，喷枪喷射浆液流速度为２５ｍ／ｓ．为便于分
析，建立固连于塔体的坐标系，并约定：吸收塔对称面所在的面为ｘＯｙ面，ｚ轴垂直于ｘＯｙ面并满足右手法则，吸收塔模型和坐标系见图１．采用非结构网格对吸收塔内计算域进行网格划分，网格数量约为２００万个．

图１　吸收塔模型和坐标系

２．２　流场分析

２．２．１　脱硫吸收塔内部流场分析

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

图２　对称面内烟气速度云图 图３　对称面内浆液速度云图

图４　吸收塔内部流线示意

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

（１）浆液喷射轮廓大于烟气，接近实际喷枪，喷射角度近６０°，参与烟气流动，在脱硫吸收塔内形成涡旋．

（２）３组喷枪所在区域的浆液与空气流场完全一致，选取对称面进行流线分析，绘出浆液空气流线图．因ｔｅｃｐｌｏｔ中流线是起点式绘制，而喷枪出口轮廓的起点皆为分离点，故起点处流线也各个分离，喷枪所在位置处速度云图和流线示意见图５．

图５　喷枪所在位置处速度云图和流线示意

（３）导流板处速度分布均匀且未形成涡旋，表明设计的导流板导流效果显著．
综上所述，从脱硫吸收塔整体及喷枪局部速度云图和流线可知，ＦＬＵＥＮＴ具备较好的模拟脱硫吸收塔在脱硫过程中烟气和浆液流动的能力．
　　
２．３　湿壁情况分析

进一步选取较小速度比例给出壁面浆液速度云图，以分析可能的浆液湿壁现象，图６～９为中间喷枪和左右２个喷枪在ｙＯｚ面内和相应ｘＯｚ面内的速度云图．由图６和７可知，左右２个喷枪喷出的浆液流贴近近壁面．进一步截取浆液速度剖面进行分析，在近最大轮廓面及以下１．５ｍ处截取剖面并精细显示速度云图，见图８和９．可知，左右２个喷枪喷射的浆液流速度在近壁面处达到３～６ｍ／ｓ，而从整体流线图看出塔内存在涡旋，考虑到实际喷枪出流速度大于２５ｍ／ｓ时将会使火焰喷射轮廓更大进而射到塔壁上，因此该脱硫吸收塔的设计方案有可能出现湿壁现象．建议将左右喷枪挪向中部，并收缩喷枪所在管道直径，增加喷枪数量．

３　结　论

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