湍球塔气体流动的数值模拟.doc

湍球塔气体流动的数值模拟

0前言
用湍球塔进行烟气脱硫，其脱硫效率在很大程度上取决于塔内发生的流体力学行为。通过前一阶段的实验发现，空塔气速分布，支撑板、挡板、漩流板的压
降及塔的总压降对湍球塔的高效连续运行有直接影响。这些参数是湍球塔最基本的特性参数，反映出塔板结构的合理性及操作过程中所需消耗的能量[1，2]。为
减少实验次数并更详尽地了解湍球塔内气体流动状况，很有必要对湍球塔内的气体流动进行数值模拟。

1方法简介
计算流体力学(简称CFD)是20世纪60年代伴随计算机技术迅速崛起的学科。CFD的应用使实验次数减少，节省了大量资金和时间，并能解决某些由于实验技术所限难以进行测量的问题，它是研究各种流体现象，设计、操作和研究各种流动系统和流动过程的有利工具口“。所以，尝试采用这种方法进行湍球塔内的气体流动模拟，分析塔内的气速分布及压力损失，为湍球塔的优化设计提供依据。
目前国外有很多发展成熟的商业CFD软件，这些软件一般包括3个主要部分：前处理器、解算器、后处理器。现采用Fluent6．0进行计算，它的解算器
采用完全的非结构化网格和控制体积法，适用于低速不可压流动、跨音速流动乃至可压缩性强的超音速和高超音速流动等各种复杂的流场，也完全适合于湍球塔内的气体流动模拟。

2几何建模与网格划分
在Fluent中，求解区域是用网格分割成有限个控制体(ControlVolumes，CVs)。同有限差分不同的是，网格为控制体积的边界，而不是计算节点。为保证守恒，CVs必须是不重叠的。因此网格生成质量对计算精度与稳定性影响极大，在几何形状复杂的区域上要生成好网格也是相当困难的。
现采用Gambit2．0进行几何建模与网格划分。该软件包含全面的几何建模能力和功能强大的网格划分工具，可划分出包含边界层等CFD特殊要求的高质量的网格。三维几何模型的建立及网格划分在整个模拟过程中是非常重要也是非常困难的一步。
图1为用Gambit2．0建立的湍球塔三维几何模型，基本上与实际实验装置在尺寸及结构上完全一致，仅湍球层塔板间距比实验装置中的间距要短些，实验装置中湍球层塔板间距为1．5 m，Gambit2．0建立的湍球塔几何模型中湍球层塔板间距为0．5 m。这主要是从减少生成网格数量和运算时间的角度考虑，且这样简化对不模拟湍球层小球的情况并没有影响。
由于实验装置为一塔设备，支撑板、挡板及漩流板的结构较复杂，且都存在极小的截面尺寸，这对网格的整体划分带来很多困难，图1所示的湍球塔模型总高2．85 m，塔径0．4m，最大截面尺寸为1 m，最小截面尺寸仅为0．004 m。为保证几何形状复杂且截面尺寸较小的几块塔板的网格质量，在尝试了多种方法后，利用Gambit2．0的网格自适应功能，对截面尺寸大的地方先用大尺度进行面网格划分，然后对整体用适合于最小截面尺寸处的小尺度进行划分，总共生成体的数量为206 550个。

图1 Gambit建立的湍球塔几何模型

3在FI uent中边界条件的设置及求解
为使问题简化，便于计算及边界条件的设置，对湍球塔内的气体流动作了以下几点假设：假设一：流体的各运动参数与时间无关，将流体流动看作定常流动；假设二：因为是冷态实验，所以认为整个模拟过程为等温过程；假设三：湍球塔内的气体流速较低，其Ma<0．3，所以在数值模拟时，将湍球塔内的气
体流动看作不可压缩流动。
在以上假设的基础上，F1uent中边界条件的设置主要涉及以下方面：
(1)定义所求解的模型方程。Fluent6．0中提供了以下湍流模型：Spalart—Allmaras模型、标准k—e模型、标准k一∞模型，但没有一个湍流模型对所有的问题是通用的。选择模型时主要依靠以下几点：流体是否可压、建立特殊的可行的问题、精度要求、计算机的能力及时间的限制。
Spalart-Allmams模型是设计用于航空领域的，主要是墙壁束缚流动，且已显示出很好的效果。在透平机械中的应用也愈加广泛。
标准k—e模型是最简单且完整的湍流模型，它是2个方程的模型，要解2个变量，即速度和长度。在Fluent中，标准k—e模型自从被Launder和Spalding
提出后，就变成工程流场计算中主要的工具了。适用范围广、经济、合理的精度，这就是它在工业流场和热交换模拟中有广泛应用的原因。这是个半经验的
公式，是从实验现象中总结出来的。
标准k一(I)模型是基于Wilcoxk一(1)模型，是为考虑低雷诺数、可压缩性和剪切流传播而修改的。wilcoxk一(1)模型预测了自由剪切流传播速率，像尾
流、混合流动、平板绕流、圆柱绕流和放射状喷射，因而可应用于墙壁束缚流动和自由剪切流动。
考虑湍球塔内的气体流动状况及以上的3点假设，本文选取了标准k—e模型。标准k—e模型的方程包括湍流动能方程k和扩散方程e：

方程(1)、(2)中G。表示由层流速度梯度而产生的湍流动能；C。是由浮力产生的湍流动能；y。由于在可压缩湍流中过渡的扩散产生的波动；C。c2、C，是常量；盯k和∥。是k方程和e方程的湍流Prandtl数；S。和．s。是用户定义的。
方程中的M，为湍流速度，由式(3)确定：

模型常量C1。=1．44，c2。=1．92，G=0．09，crk=1．0，矿。=1_3。Fluent中这些常量是从试验中得来的，包括空气、水的基本湍流。这些常量对于大多数情况是适用的，用户可在粘性模型面板中修改，本文在计算时均选取默认值。
(2)定义流动介质属性。实际操作时是用空气进行实验，所以在Fluent中把流体性质也定义为空气，空气的各物性参数取默认值。
(3)定义边界条件值。在Gambit中已定义好压力入口及压力出口，在Fluent中找到相应的面，根据不同的实验工况，定义出总压、静压及入口的坐标方向，Fluent本身会在初始化时自动计算出入口速度值。
(4)定义松弛因子和离散方程的格式。F1uent6．0采用压力校正法(SIMPLE法)作为其低速计算模块，同时也提供其他算法，如SIMPLER、SIMPLEC和PISO等。选取SIMPLE法，压力、动量、能量、密度方程的松弛因子均取默认值。
(5)定义收敛精度。Fluent6．0中默认的收敛精度为10～，经多次计算结果与实验结果对比，并考虑节省计算时间，取连续性方程的收敛精度为104。
(6)计算。把迭代次数设为1 000，计算到396步即收敛。

4计算结果
图2为空塔气速分布曲线，该工况进风管处静压为1 000Pa、全压为1 118 Pa，塔内气速数据是用烟道分析仪在湍球层顶部沿塔径方向的9个测点依次测得。图3为用Fluent在和实验工况完全相同的边界条件下计算出的塔内气速分布。

通过对比可看出，在塔径[一0．05，0．05]范围内，Fluent模拟数据与实验数据吻合较好。在[一O．20，一0．05]及[0．05，0．20]范围内误差较大，这主要是由实验装置中粗糙的璧面及开孔对靠近璧面处的气体流动有较大影响，而模型中的璧面是完全理想化的也没有开孔的影响。塔径中心处的速度波谷是由其上部的莲蓬头(见图1)造成的。
气体通过塔板的压降是塔板的主要流体力学特性，它不仅影响塔板的操作，还决定塔内轴向压力分布和全塔压降。表l给出不同工况下总压降、支撑板、挡板、漩流板实验值与模拟值的对比，通过对比可发现，最大误差为155．8％，最小误差仅为2％，4个工况的总压降平均误差为26．6％。最小的误差在挡板处，较大的误差都是发生在漩流板处，主要是由于档板结构较简单但旋流板结构非常复杂，旋流板处的气流变化较大，造成在数据测量时产生一定的误差。

参考文献：
[1]李爱民，楚华，魏励宏，等．湍球塔脱硫的冷态试验研究[J]．中国电力，2002，36(2)：57—59．
[2]楚 华，李爱民，魏励宏，等．脱硫湍球塔阻力特性试验及因次分析[J]．环境工程，2003，(增刊)：140—143．
[3]吴子牛．计算流体力学基本原理[M]．北京：科学出版社，2001．
[4][美]WF休斯，JA布赖顿．流体力学[M]．麦格劳一希尔教育出版集团，北京：科学出版社，2002．
[5]P JWitt，JHPerry，MPSchwarz．A numericalmodelforpredictingbubble formation in a 3D Fui dized bed[J]．Applied MathematicalModelling，1998，(22)：1071—1080．
[6]E 0lmos，c Gentric，ch Vial，e￡d．Nume“cal simulation of multi—phase flow in bubble olumnreactorS Innuence of bubble coalescenceand break—up[J]．Chemical Engineering Science，2001，56：6359—6365．