的求解过程

为了进行CFD计算,用户可借助商用软件来完成所需要的任务,也可自己直接编写计算程序。两种方法的基本工作过程是相同的,本节将介绍基本计算思路。无论是流动问题还是传热问题,无论是稳态问题还是瞬态问题,其求解过程都可用图9-1表示。

如果所求解的问题是瞬态问题,则可将上图的过程理解为一个时间步的计算过程,循环这一过程求解下个时间步的解。下面对各求解步骤作一简单介绍。

(1)建立控制方程

建立控制方程,是求解任何问题前都必须首先进行的。一般来讲,这一步是比较简单的,因为对于一般的流体流动而言,可直接写出其控制方程。例如,对于水流在水轮机内的流动分析问题,若假定没有热交换发生,则可直接将连续方程与动量方程作为控制方程使用。当然,由于水轮机内的流动大多是处于湍流范围,因此,一般情况下,需要增加湍流方程。

(2)确立初始条件和边界条件

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

图9-1　CFD求解过程

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

边界条件是在求解区域的边界上所求解的变量或其导数随地点和时间的变化规律。对于任何问题,都需要给定边界条件。例如,在锥管内的流动,在锥口进口断面上,我们可给定速度、压力沿半径方向的分布,而在管壁上,对速度取无滑移边界条件。

(3)划分计算网格

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

不同的问题采用不同数值解法时,所需要的网格形式是有一定区别的,但生成网格的基本方法是一致的。目前,网格分结构网格和非结构网格两大类。简单地讲,结构网格在空间上比较规范,如对一个四边形区域,网格往往是成行成列分布的,行线和列线比较明显。而对非结构网格在空间分布上没有明显的行线和列线。

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

目前各种CFD软件都配有专用的网格生成工具,多数CFD软件可接收其他CAD或CFD/FEM软件产生的网格模型。

(4)建立离散方程

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

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

(5)离散初始条件和边界条件

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

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

(6)给定求解控制参数

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

(7)求解离散方程

在进行了上述设置后,生成了具有定解条件的代数方程组。对于这些方程组,数学上已有相应的解法,如线性方程组可采用Gauss消去法或Gauss-Seidel迭代法求解,而对非线性方程组,可采用Newton-Raphson方法。在商用CFD软件中,往往提供多种不同的解法,以适应不同类型的问题。

(8)判断解的收敛性

对于稳态问题的解,或是瞬态问题在某个特定时间步上的解,往往需要通过多次迭代才能得到。有时,因网格形式或网格太小、对流项的离散插值格式等原因,可能导致解的发散。对于瞬态问题,若采用显示格式进行时间域上的积分,当时间步长过大时,也可能造成解的振荡或发散。因此,在迭代过程中,要对解的收敛性随时进行监视,并在系统达到指定精度后,结束迭代过程。

(9)显示和输出计算结果并进行后处理

通过上述求解过程得出了各计算节点上的解后,需要通过适当的手段将整个计算域上的结果表示出来。这时,我们可以采用线值图、矢量图、等值线图、流线图、云图等方式对计算结果进行表示。所谓线值图,是指在二维或三维空间上,将横坐标取为空间长度或时间历程,将纵坐标取为某一物理量,然后用光滑曲线或曲面在坐标系内绘制出某一物理量沿空间或时间的变化情况。矢量图是直接给出二维或三维空间里矢量(如速度)的方向及大小,一般用不同颜色和长度的箭头表示速度矢量。矢量图可以比较清楚地让用户发现其中存在的涡旋区。等值线图是用不同颜色的线条表示相等物理量(如温度)的一条线。流线图是用不同颜色线条表示质点运动轨迹。云图是使用渲染的方式,将流场某个截面上的物理量(如压力或温度)用连续变化的颜色块表示其分布。