化工流体力学 第六章.pptVIP

鼓泡塔中汽液两相流的数值计算： 鼓泡塔是汽液两相进行质量、动量和能量传递的重要设备，广泛应用于化工、石油及生物化工等 领域中。 双流体模型:鼓泡塔一般为圆柱形结构，气体从底部进入鼓泡床后，产生垂直向上运动的气泡，一 般情况下气泡在塔内分布不均匀，可用双流体模型（欧拉-欧拉模型）来描述汽液两相流动。作出 假设：气泡与液体可视为互相穿透的连续介质，两种流体以一定几率出现，用气含率与液含率表 征，汽液两相流的基本方程由连续性方程与动量方程组成。对于含有Np相的流体其连续方程如下： 式中右端项描述为由α相到β的质量传递。r α代表α相的体积分数。忽略相间质量传递与源相： 动量方程为： 方程右端描述了所有作用在控制体上α相流体上的力，总压力项，黏性应力项以及重力与相间相 互作用项Fα,i。只考虑阻力的影响，可根据Clift等所得的公式： 式中，阻力系数CD 0.44,；气泡直径db 4mm为常数； 为气泡相对于流体的滑移速度； 为流体流动速度。 对于连续相可用单相标准湍流模型k-e模型来模拟汽液两相流中连续相中的湍流现象。其湍流动能 k与湍流耗散率e方程如下： k-e模型中的常数可查表得。汽液流中气泡诱导的湍流对流动有一定影响，可通过改进湍流动能与 湍流耗散率产生项Gα体现，该项描述相间的相互作用。如下式： 式中， 1.44与 1.92。 初始边界条件： 垂直上升汽液两相流在鼓泡塔内的流动物理模型和求解区域如下图所示，设定入口及出口边界条 件，划分网格。在求解时应考虑与单相流的不同点，从汽液两个连续方程出发，应用SIMPLE方 法，进行压力校正和含气率求解。 结果与讨论： 鼓泡塔内液体速度矢量与气含率等值线分布图如下所示，图a为t 700s时塔内速度矢量与气含率的 瞬时分布，可知塔内存在涡旋回流区。图b为历时1400s时，塔内速度矢量与气含率平均分布，其 速度矢量平均值呈抛物线分布，在塔轴线附近上升的速度值最大，在壁面附近沿壁面向下的速度 值最大。含气率的平均值在轴线附近最大。 下图为气体进口速度分别在0.15cm/s，0.5cm/s，1.0cm/s和2.0cm/s时，塔内h 1.6m处液体相应速度 沿径向分布的实验与数值计算结果比较，可知当速度在0.15cm/s，1.0cm/s与0.5cm/s时，径向速度 实验与数值计算结果较为吻合，而当气体进口速度较大（2cm/s）时，两者相差较大。 下图为气体进口速度分别在0.15cm/s，0.5cm/s，1.0cm/s和2.0cm/s时，塔内h 1.6m处含气率沿径向 的分布的实验与数值计算结果比较。由图可知，其实验值的分布较计算值陡峭，轴线附近含气率 的实验值比计算值大而壁面附近的实验值比计算值小。这表明双流体模型能有效预测鼓泡塔中气 液两相流流动的主要参数。 * * 第六章 计算流体力学 本章要点： 1. 流体力学主要方程与模型 2.数值方法（有限元法） 3.节典型流场（绕流，管流）的数值计算 4.节鼓泡反应器中流场的数值模拟 概 述 计算流体力学采用数值模拟来研究流体运动规律，具有耗费少，时间短，重复性好等优点。 解决问题思路：遵循以下步骤（1）进行简化以确定合理的物理模型。 （2）根据所确定物理模型建立相应的数学模型。 （3）对数学模型中各偏微分方程组选用合适的方法进行离散成代数 方程组，用恰当的格式进行求解。 （4）对所得的数值结果进行验证。 6.1 流体力学主要方程与模型 不可压N-S方程 本节主要介绍不可压N-S方程的统一形式，并推导出涡流函数方程。并简要介绍湍流模型。 在现代计算流体力学里，N-S方程写成如下统一形式： 式中各项代表意义与第三章相同。 涡流函数方程 对于二维不可压黏性流动，忽略体积力，其直角坐标下无量纲方程形式： 动量方程中有压力项，为了消去压力项，引入流函数ψ： 这时，连续方程自动满足。将上式带入下二维涡量的定义式中： 可得泊松方程： 涡量传递方程为： 上两式即为涡流函数方程的N-S 湍流模型： 第三章所述的雷诺方程不封闭，因此必须补充方程，所补充的方程有零方程、一方程和二方程模型， 还有代数应力模型（ASM）和雷诺应力方程模型（DSM）等。 单方程模型：湍流动能方程形式如下： 式中 为湍流产生项，σk为湍流脉动动能的普朗特数，常取为1，cD为常数， 湍流的长度标尺l可由经验公式得到。单方程模型考虑到了湍流动能的对流与扩散，湍流粘度联系到整 个流场。 双方程模型： 湍流耗散率ε常用形式如下： 湍流动能κ方程可改为如下形式： 模型中各经验常数取值如下：cμ 0.9，c1 1.44，c2 1.92，σk 1.0，σε 1.3，σT 1.0。 Κ和ε是表征湍流发生、发展与耗散的重要特征量