《FLUENT 14.0超级学习手册》——1.2 计算流体力学（CFD）基础

本节书摘来自异步社区《FLUENT 14.0超级学习手册》一书中的第1章，第1.2节,作者： 唐家鹏 更多章节内容可以访问云栖社区“异步社区”公众号查看。

1.2 计算流体力学（CFD）基础

FLUENT 14.0超级学习手册
计算流体动力学（Computational Fluid Dynamics，CFD）是近代流体力学、数值数学和计算机科学结合的产物，是一门具有强大生命力的边缘科学。

1.2.1 CFD概述
CFD以电子计算机为工具，应用各种离散化的数学方法，对流体力学的各类问题进行数值实验、计算机模拟和分析研究，以解决各种实际问题。

计算流体力学和相关的计算传热学、计算燃烧学的原理是用数值方法求解非线性联立的质量、能量、组分、动量和自定义的标量的微分方程组，求解结果能预报流动、传热、传质、燃烧等过程的细节，并成为过程装置优化和放大定量设计的有力工具。计算流体力学的基本特征是数值模拟和计算机实验，它从基本物理定理出发，在很大程度上替代了耗资巨大的流体动力学实验设备，在科学研究和工程技术中产生巨大的影响。

计算流体力学是目前国际上一个强有力的研究领域，是进行传热、传质、动量传递及燃烧、多相流和化学反应研究的核心和重要技术，广泛应用于航天设计、汽车设计、生物医学工业、化工处理工业、涡轮机设计、半导体设计、HAVC&R 等诸多工程领域，板翅式换热器 设计是CFD技术应用的重要领域之一。

CFD在最近20 年中得到飞速的发展，除了计算机硬件工业的发展给它提供了坚实的物质基础外，还主要因为无论分析的方法或实验的方法都有较大的限制。例如，由于问题的复杂性，既无法作分析解，也因费用昂贵而无力进行实验确定，而CFD 的方法正具有成本低和能模拟较复杂或较理想的过程等优点。

经过一定考核的CFD软件可以拓宽实验研究的范围，减少成本昂贵的实验工作量。在给定的参数下用计算机对现象进行一次数值模拟相当于进行一次数值实验，历史上也曾有过首先由CFD数值模拟发现新现象而后由实验予以证实的例子。

CFD软件一般都能推出多种优化的物理模型，如定常和非定常流动、层流、紊流、不可压缩和可压缩流动、传热、化学反应等。对每一种物理问题的流动特点，都有适合它的数值解法，用户可选择显式或隐式差分格式，以期在计算速度、稳定性和精度等方面达到最佳。

CFD 软件之间可以方便地进行数值交换，并采用统一的前、后处理工具，这就省去了科研工作者在计算机方法、编程、前后处理等方面投入的重复、低效的劳动，而可以将主要精力和智慧用于物理问题本身的探索上。

1.2.2 CFD求解力学问题的过程
所有CFD问题的求解过程都可用图1-4表示。如果所求解的问题是瞬态问题，则可将图1-4的过程理解为一个时间步的计算过程，循环这一过程求解下个时间步的解。下面对各求解步骤进行简单介绍。

1．建立控制方程
建立控制方程是求解任何问题前都必须首先进行的。一般来讲，这一步是比较简单的。因为对于一般的流体流动而言，可直接写出其控制方程。假定没有热交换发生，则可直接将连续方程与动量方程作为控制方程使用。一般情况下，需要增加湍流方程。

2．确定边界条件和初始条件
初始条件与边界条件是控制方程有确定解的前提，控制方程与相应的初始条件、边界条件的组合构成对一个物理过程完整的数学描述。

初始条件是所研究对象在过程开始时刻各个求解变量的空间分布情况。对于瞬态问题，必须给定初始条件。对于稳态问题，不需要初始条件。

边界条件是在求解区域的边界上所求解的变量或其导数随地点和时间的变化规律。对于任何问题，都需要给定边界条件。

3．划分计算网格
采用数值方法求解控制方程时，都是想办法将控制方程在空间区域上进行离散，然后求解得到离散方程组。要想在空间域上离散控制方程，必须使用网格。现已发展出多种对各种区域进行离散以生成网格的方法，这些方法统称为网格生成技术。

不同的问题采用不同数值解法时，所需要的网格形式是有一定区别的，但生成网格的方法基本是一致的。目前网格分结构网格和非结构网格两大类。

简单地讲，结构网格在空间上比较规范，如对一个四边形区域，网格往往是成行成列分布的，行线和列线比较明显。而非结构网格在空间分布上没有明显的行线和列线。

对于二维问题，常用的网格单元有三角形和四边形等形式；对于三维问题，常用的网格单元有四面体、六面体、三菱体等形式。在整个计算域上，网格通过节点联系在一起。

目前各种CFD软件都配有专用的网格生成工具，如FLUENT 使用Gambit作为前处理软件。多数CFD软件可接收采用其他CAD或CFD/FEM软件产生的网格模型。例如，FLUENT可以接收ANSYS所生成的网格。

4．建立离散方程
对于在求解域内所建立的偏微分方程，理论上是有真解（或称精确解或解析解）的。但由于所处理问题自身的复杂性，一般很难获得方程的真解。因此就需要通过数值方法把计算域内有限数量位置（网格节点或网格中心点）上的因变量值当做基本未知量来处理，从而建立一组关于这些未知量的代数方程组，然后通过求解代数方程组来得到这些节点值，而计算域内其他位置上的值则根据节点位置上的值来确定。

由于所引入的应变量在节点之间的分数假设及推导离散化方程的方法不同，所以形成了有限差分法、有限元法、有限元体积法等不同类型的离散化方法。

对于瞬态问题，除了在空间域上的离散外，还要涉及在时间域上的离散。离散后，将要涉及使用何种时间积分方案的问题。

5．离散初始条件和边界条件
前面所给定的初始条件和边界条件是连续性的，如在静止壁面上速度为0，现在需要针对所生成的网格，将连续型的初始条件和边界条件转化为特定节点上的值，如静止壁面上共有90个节点，则这些节点上的速度值应均设为0。

商用CFD软件往往在前处理阶段完成网格划分后，直接在边界上指定初始条件和边界条件，然后由前处理软件自动将这些初始条件和边界条件按离散的方式分配到相应的节点上。

6．给定求解控制参数
在离散空间上建立了离散化的代数方程组，并施加离散化的初始条件和边界条件后，还需要给定流体的物理参数和湍流模型的经验系数等。此外，还要给定迭代计算的控制精度、瞬态问题的时间步长和输出频率等。

7．求解离散方程
进行上述设置后，生成了具有定解条件的代数方程组。对于这些方程组，数学上已有相应的解法，如线性方程组可采用Gauss消去法或Gauss-Seidel迭代法求解，而对于非线性方程组，可采用Newton-Raphson方法。

商用CFD软件往往提供多种不同的解法，以适应不同类型的问题。这部分内容属于求解器设置的范畴。

8．显示计算结果
通过上述求解过程得出了各计算节点上的解后，需要通过适当的手段将整个计算域上的结果表示出来ÿ