CFD计算过程和算法.pdf

CFD计计算算过过程程和和算算法法

各位同学，⼤家好，我是七师兄，今天我们来学习Airpak⾼级班的第四节的内容。《计算流体

⼒学的计算过程和算法》

CFD数值模拟⼀般遵循以下⼏个步骤:

(1)建⽴所研究问题的物理模型,再将其抽象成为数学、⼒学模型。然后确定要分析的⼏何体

的空间影响区域。

(2)建⽴整个⼏何形体与其空间影响区域,即计算区域的CAD模型,将⼏何体的外表⾯和整个计

算区域进⾏空间⽹格划分。⽹格的稀疏以及⽹格单元的形状都会对的计算产⽣很⼤的影响。不

同的算法格式为保证计算的稳定性和计算效率,⼀般对⽹格的要求也不同。

(3)加⼊求解所需要的初始条件,⼊⼝与出⼝处的边界条件⼀般为速度、压⼒条件。

(4)选择适当的算法,设定具体的控制求解过程和精度的⼀些条件,对所需分析的问题进⾏求

解,并保存数据⽂件结果。

(5)选择合适的后处理器(PostProcessor)读取计算结果⽂件,分析并显⽰出来。

以上这些步骤构成了CFD数值模拟的全过程。其中数学模型的建⽴我们会在后⾯的程中

讲到。

2.在运⽤CFD⽅法对⼀些实际问题进⾏模拟时,常常需要设置⼯作环境、边界条件和动量守

恒和选择算法等,特别是算法的选择对模拟的效率及其正确性有很⼤影响,需要特别重视。

区域离散化是利⽤⼀组有限个离散的点来代替原来连续的空间。实施过程是把所计算的区

域划分成许多互不重叠的⼦区域,确定每个⼦区域的节点位置和该节点所代表的控制体积。节点

是需要求解未知物理量的⼏何位置、控制体积、应⽤控制⽅程或守恒定律区域,的最⼩⼏何单

位。

我们来看下，这个图⽚，左边是⼀个房间，⼀个三维模型，他是⼀个连续空间，我们知道

连续空间，它⾥⾯⽆数个节点，为了⽅便研究这个三维空间，我们将它⽹格划分成有限个细⼩

的单元，这样的话，就⽅便计算机进⾏迭代计算了。

常⽤的离散化⽅法：有限差分法、有限单元法和有限体积法

1）有限差分法：是数值解法中最经典的⽅法,它是将求解区域划分为差分⽹格,⽤有限个⽹

格节点代替连续的求解域,然后将偏微分⽅程(控制⽅程)的导数⽤差商代替,推导出含有离散点上

有限个未知数的差分⽅程组。

这种⽅法产⽣和发展得⽐较早,也⽐较成熟,较多⽤于求解双曲线和地物线型问题。⽤它求解

边界条件较为复杂,尤其是求解椭圆型问题时,不如有限元法或有限体积法⽅便。

构造差分的⽅法有多种形式,⽬前主要采⽤的是奉勒级数展开⽅法,其基本的差分表达式主要

有四种形式,⼀阶向前差分,⼀阶向后差分,⼀阶中⼼差分和⼆阶中⼼差分，其中前两种形式为⼀阶

计算精度,后两种形式为⼆阶计算精度,通过对时间和空间⼏种不同差分形式的组合,可以组合成不

同的差分计算格式。

2）有限单元法

有限单元法是将⼀个连续的求解域任意分成适当形状的许多微⼩单元,并于各⼩单元分⽚构

造插值函数,然后根据极值原理(变分或加权余量法)将问题的控制⽅程转化为有单元上的有限元⽅

程,把总体的极值作为各单元极值之和,即将局部单元合成，形成嵌⼊了指定边界条件的代数⽅程

组,求解该⽅程组就得到各节点上待求的函数恒.

有限单元求解的速度⽐有限差分法和有限体积法慢,在商⽤CFD软件中应⽤并不泛,⽬前常

⽤的商⽤CFD软件中,只有FIDAP采⽤的是有单限元法。

3）有限体积法

有限体积法⼜称为控制体积法,是将计算区域划分为⽹格,并使每个⽹络周围有⼀个互不重复

的控制体积,将待解的微分⽅程对每个控制体积积分,从⽽得到⼀组离散⽅程。

其中的未知数是⽹格节点上的因变量,⼦域法加离散是有限体积法的基本思想。有限体积法

的基本思路易于理解,并能得出直接的物理解释。

离散⽅程的物理意义是因变量在有限⼤⼩的控制体积中的守恒原理,如同微分⽅表⽰因变量

在⽆限⼩的控制体积中的守恒原理⼀样。

有限体积法得出的离散⽅程要求因变量的积分守恒对任意⼀组控制集体都得到满⾜,对整个

计算区域⾃然也得到满⾜,这是有限体积法的优点。

有⼀些离散⽅法,如有限差分法,仅当⽹格极其细密时,离散⽅程才满⾜积分守恒；⽽有限体积

法即使在粗⽹格情况下,也显⽰出准确的积分守恒。

就离散⽅法⽽⾔,有限体积法可视作有限单元法和有限差分法的中间产物。

有限差分法较直观,理论成熟,精度可选,但对于不规则区域的处理较为烦琐,虽然⽹格⽣成可以

使有限差分法应⽤于不规则区域,但对于区域的连续性等