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

《湍球塔气体流动的数值模拟》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、湍球塔气体流动的数值模拟0前言用湍球塔进行烟气脱硫，其脱硫效率在很大程度上取决于塔内发生的流体力学行为。通过前一阶段的实验发现，空塔气速分布，支撑板、挡板、漩流板的压降及塔的总压降对湍球塔的高效连续运行有直接影响。这些参数是湍球塔最基本的特性参数，反映出塔板结构的合理性及操作过程中所需消耗的能量[1，2]。为减少实验次数并更详尽地了解湍球塔内气体流动状况，很有必要对湍球塔内的气体流动进行数值模拟。1方法简介计算流体力学(简称CFD)是20世纪60年代伴随计算机技术迅速崛起的学科。CFD的应用使实验次数减少，节省了大量资金和时间，并能解决某些由于实验技术所限难以进行测量的问题，它是研

2、究各种流体现象，设计、操作和研究各种流动系统和流动过程的有利工具口“。所以，尝试采用这种方法进行湍球塔内的气体流动模拟，分析塔内的气速分布及压力损失，为湍球塔的优化设计提供依据。目前国外有很多发展成熟的商业CFD软件，这些软件一般包括3个主要部分：前处理器、解算器、后处理器。现采用Fluent6．0进行计算，它的解算器采用完全的非结构化网格和控制体积法，适用于低速不可压流动、跨音速流动乃至可压缩性强的超音速和高超音速流动等各种复杂的流场，也完全适合于湍球塔内的气体流动模拟。2几何建模与网格划分在Fluent中，求解区域是用网格分割成有限个控制体(ControlVolumes，CVs

3、)。同有限差分不同的是，网格为控制体积的边界，而不是计算节点。为保证守恒，CVs必须是不重叠的。因此网格生成质量对计算精度与稳定性影响极大，在几何形状复杂的区域上要生成好网格也是相当困难的。现采用Gambit2．0进行几何建模与网格划分。该软件包含全面的几何建模能力和功能强大的网格划分工具，可划分出包含边界层等CFD特殊要求的高质量的网格。三维几何模型的建立及网格划分在整个模拟过程中是非常重要也是非常困难的一步。图1为用Gambit2．0建立的湍球塔三维几何模型，基本上与实际实验装置在尺寸及结构上完全一致，仅湍球层塔板间距比实验装置中的间距要短些，实验装置中湍球层塔板间距为1．5m

4、，Gambit2．0建立的湍球塔几何模型中湍球层塔板间距为0．5m。这主要是从减少生成网格数量和运算时间的角度考虑，且这样简化对不模拟湍球层小球的情况并没有影响。由于实验装置为一塔设备，支撑板、挡板及漩流板的结构较复杂，且都存在极小的截面尺寸，这对网格的整体划分带来很多困难，图1所示的湍球塔模型总高2．85m，塔径0．4m，最大截面尺寸为1m，最小截面尺寸仅为0．004m。为保证几何形状复杂且截面尺寸较小的几块塔板的网格质量，在尝试了多种方法后，利用Gambit2．0的网格自适应功能，对截面尺寸大的地方先用大尺度进行面网格划分，然后对整体用适合于最小截面尺寸处的小尺度进行划分，总共

5、生成体的数量为206550个。图1Gambit建立的湍球塔几何模型3在FIuent中边界条件的设置及求解为使问题简化，便于计算及边界条件的设置，对湍球塔内的气体流动作了以下几点假设：假设一：流体的各运动参数与时间无关，将流体流动看作定常流动；假设二：因为是冷态实验，所以认为整个模拟过程为等温过程；假设三：湍球塔内的气体流速较低，其Ma<0．3，所以在数值模拟时，将湍球塔内的气体流动看作不可压缩流动。在以上假设的基础上，F1uent中边界条件的设置主要涉及以下方面：(1)定义所求解的模型方程。Fluent6．0中提供了以下湍流模型：Spalart—Allmaras模型、标准k—e模型

6、、标准k一∞模型，但没有一个湍流模型对所有的问题是通用的。选择模型时主要依靠以下几点：流体是否可压、建立特殊的可行的问题、精度要求、计算机的能力及时间的限制。Spalart-Allmams模型是设计用于航空领域的，主要是墙壁束缚流动，且已显示出很好的效果。在透平机械中的应用也愈加广泛。标准k—e模型是最简单且完整的湍流模型，它是2个方程的模型，要解2个变量，即速度和长度。在Fluent中，标准k—e模型自从被Launder和Spalding提出后，就变成工程流场计算中主要的工具了。适用范围广、经济、合理的精度，这就是它在工业流场和热交换模拟中有广泛应用的原因。这是个半经验的公式，是

7、从实验现象中总结出来的。标准k一(I)模型是基于Wilcoxk一(1)模型，是为考虑低雷诺数、可压缩性和剪切流传播而修改的。wilcoxk一(1)模型预测了自由剪切流传播速率，像尾流、混合流动、平板绕流、圆柱绕流和放射状喷射，因而可应用于墙壁束缚流动和自由剪切流动。考虑湍球塔内的气体流动状况及以上的3点假设，本文选取了标准k—e模型。标准k—e模型的方程包括湍流动能方程k和扩散方程e：方程(1)、(2)中G。表示由层流速度梯度而产生的湍流动能；C。是由浮力产生的湍流动