分离器数值模拟 fluent

7.数值模拟与结果分析

7.1数值计算方法简介

计算流体力学作为流体力学研究中的一门新兴分支，正在工业和科研领域内发挥越来越重要的作用。将CFD工具运用到分离机械的研究中，也成为工程技术人员改进设计、提高效率的有效手段，是CFD应用的前沿。一些成熟的算法，模型也以商业软件的形式出现在工程及科研领域。相比研究单位自行开发的计算程序，商业计算软件一般具有以下特点：

①通用性广。由于商业软件面向的用户对象广泛，处理的实际问题多种多样，因此其覆盖的应用范围要尽可能广。

②计算稳定性好。多数软件经过不同研究领域内的算例测试，对不同类型的问题具有较好的适应能力。

③使用方便，商业软件经过不同友好的用户界面，方便用户的使用。

④一般商业软件也存在一些明显的不足，例如：算法相对陈旧，不能紧跟CFD研究领域内的最新成果；与不同行业内的实际要求存在一定的距离，难以将各研究单位已有的研究成果结合到商业软件中。这在一定程度上限制了商业软件在工程实际中的应用。

FLUENT是由美国FLUENT公司于1983推出的CFD软件。它是继PHOENICS软件之后的第二个投放市场的基于有限体积法的软件。FLUENT是目前功能最全面、适用性最广、国内使用最广泛的CFD软件之一。本文运用fluent软件对离心式分离器的内流场进行分析计算，fluent公司是享誉世界的最大计算流体力学软件供应商，fluent软件能够精确地模拟无粘流、层流、湍流、化学反应、多相流等复杂的流动现象。应用领域包括：航空航天、汽车设计、生物医药、化学处理、石油天然气、发电系统电子半导体、涡轮设计、HVAC、玻璃加工等。FLUENT 具有精度高，收敛快，稳定性好等特点。

Gambit是前置处理器，能针对及其复杂的几何外形生成三维四面体，六面体的非结构化网格及混合网格。该模块还具有方便的网络检查功能，对网络单元体积、扭曲率、长细比等影响收敛和稳定的参数进行统计并生成报告。

7.2 计算流体力学基础

在流体力学的研究中，常用的方法有理论研究方法、数值计算方法和实验研究方法。理论研究方法的特点是：能够清晰、普遍地揭示出流动的内在规律，但该方

法目前只局限于少数比较简单的理论模型，而且需要研究者具有较高的理论素养和数

学功底。实验研究方法的特点就是结果可靠，但其局限性在于相似准则不能全部满足、

尺寸限制、边界影响等，同时，实验研究需要场地、仪器设备和大量的经费，研究周

期也比较长。数值计算方法所需要的时间和费用都较少，并且具有较高的精度，目前

在流体力学的研究中扮演着越来越重要的角色。数值计算方法要求对问题的物理特性

有足够的了解，并能建立较精确的描述方程组。

计算流体力学的求解过程大致可以分成下列若干步骤：

（1）建立基本守恒方程组

数值模拟的第一步是由流体力学、热力学、传热传质学、燃烧学及热等离子的基本原理出发，建立质量、动量、能量、组分、湍流特性的守恒方程组，如连续方程、动量方程、组分方程、湍能方程等。对湍流、多相流等，由不同的模拟理论出发，往往基本守恒方程组也不相同，因此，如何构造基本方程组，也是模拟理论的重要部分。FLUENT中所采用的都是比较成熟的模拟理论和数值模型。

（2）建立或选择模型或封闭方法

建立的基本方程往往是不封闭的，特别是湍流、甚至多相流、化学反应流更是如此。例如，动量方程中的脉动速度关联项（雷诺应力项），能量方程中的湍流导热项及辐射项，扩散方程中的扩散项及湍流反应项等都是未知的。解决这一问题，使方程组封闭，就是模拟理论的关键问题。FLUENT中已经预设了许多的物理模型，如湍流模型、两相流模型、湍流反应模型、辐射换热模型、污染物生成模型等。用户可以根据具体问题选择不同的模型，也可以自己通过实验事实或物理概念的基本假设来构造各个过程的模型。

（3）确定初始与边界条件

数值模拟的第三步是必须按给定的几何形状和尺寸，由问题的物理特征出发，确定计算域并给定计算域的进出口，轴线（或对称面）及各壁面或自由面处条件。对湍流和多相流动还需要分别给出各相的各变量的时均值及脉动值的各初始条件与边界条件。正确给定边界条件是十分重要的。边界条件是否合理往往也是数值模拟成败的关键问题之一，初始条件是所研究对象在过程开始时候各个求解变量的空间分布情况。对于瞬态的非定常问题必须给定初始条件。对于定常问题，不需要初始条件。

（4）划分计算网格

采用数值方法求解控制方程时，都是想办法将控制方程在空间区域上进行离散，然后求解得到离散方程组，其本质就是把连续的空间变量用离散的网格点上的变量来近似，连

续的控制方程在离散后就成为所有网格点上变量非线性方程组。要想在空间区域上离散控制方程组，必须使用网格。现在已经发展出多种对各种区域进行离散以生成网格的方法，网格生成技术也成为CFD领域的一个独特分支。可以说，网格生成占据了整个CFD任务的70%以上的工作量。

（5）建立离散化方程

用数值方法求解偏微分方程组，必须将该方程组离散化，即把计算域内有限数量位置（网格节点或网格控制体中心点）上的因变量作为基本未知量来处理，从而建立一系列关于这组未知量的代数方程组，然后通过求解代数方程组来得到这些节点上的值。

对于所引入的因变量在节点之间的分布，假设及推导离散化方程的方法不同，形成了有限差分法，有限容积法，有限元法，或有限分析法等不同类型的离散化方法。FLUENT 使用的是有限容积法。在同一种离散化方法中，例如有限体积法，对方程中对流项所采用的离散格式不同，也将导致不同形式的离散方程，这种离散格式通常称为空间差分格式，FLUENT提供了多种离散格式供选择，如中心差分、一阶迎风格式、二阶迎风格式、QUICK 格式三阶MUCSL格式等，还有为可压缩流动中激波等间断捕捉设计的Roe格式，AUSM 类格式等。对于非定常问题还有涉及时间上的差分，FLUENT提供了隐式、显式的一阶和二阶的时间差分格式。

（6）制定求解方法

对离散完成的差分方程组已经有各种不同的求解方法。例如涡量—流函数算法、基于压力的压力—速度修正算法（SIMPLE系列算法），基于密度的耦合隐式或显式时间推进求解算法，矢通量分裂方法和通量差分分裂方法等。针对代数方程组的求解有三角矩阵法（追赶法）、逐线迭代、松弛高斯赛德尔迭代方法等。针对两相流和有反应的流动又有一些更专门的解法，如颗粒与流线的耦合PISC法，加速化学反应计算而设计的ISAT算法等。

FLUENT已经在求解器内设计了目前多数已经成熟的求解方法供用户选择，对算法及相应参数意义的详细了解有助于我们正确的使用FLUENT来设置所需的各种松弛因子、算法参数，提高计算效率。

（7）除了上述基本解法以外，还要针对具体问题的特点，研究一些计算方法的细节或称计算技巧。例如对于合理而经济的网格划分与安排，有时要选择随机过程的空间或时间而变化的网格系，以便不抹掉物理特征而又经济。又如对不规则形状边界的处理，松弛