高等计算流体力学讲义(8).docVIP

PAGE1

PAGE20

MACROBUTTONMTEditEquationSection高等计算流体力学讲义（8）

第三章不可压缩流动的数值方法

§1基本方程及其性质

一、基本方程

考虑不可压缩NS方程：

（1）

（2）

其中粘性应力为，

（3）

如果粘性系数为常数，

（4）

经无量纲化，常粘性系数不可压缩NS方程可以写为：

，

其中为运动粘性系数。NS方程也可以写为无量纲化形式

其中已经吸收到p中（代表）。不可压缩方程的边界条件为：

固体壁面：，

进口条件：，

出口条件：。

不可压缩方程中的压力场可以相差任一常数而对速度场无影响，所以压力场只是在相差任意常数的条件下是确定的。为了确定全场压力值，还应指定流场中某一点的压力。

二、不可压N－S方程的特点：

方程为二阶偏微分方程，二阶项中包含参数（粘性系数）。

边界层、分离、湍流…

方程是非线性的，表现为对流项。

对一维问题，非线性项为。假定u的波数为k的Fourier分量为

（5）

则：。即振幅由；波数由。也就是说，振幅呈现非线性变化，且可以产生高频成分。粘性的作用，使得解的结构进一步复杂化，考虑模型方程

把（5）式带入模型方程，得

可见，雷诺数越大，或频率越低（流动结构的尺度越大），振幅衰减越慢。

综上所述：由于非线性的作用，会产生高频的流动结构；在大雷诺数的条件下，这些高频结构有较长的生命周期，并且与衰减缓慢的低频结构相互作用，使得流动表现出复杂的的非线性、多尺度特征。

没有压力对时间的偏导数。

由于方程表现出很强的椭圆形，不能利用比较成熟的发展型偏微分方程的数值求解的理论和方法，造成数值求解的困难。尤其是没有明显的计算压力的方程，严格地说，必须耦合求解动量方程和连续方程才能求出压力场和速度场。这种做法，计算量是非常巨大的。目前，不可压缩流动的数值方法还远不如可压缩流动的数值方法更成熟。

三、耦合求解方案

NS方程的一种可能的求解方案是动量方程和连续方程完全耦合求解，例如可以采用下面的方法。

差分算子的定义为：

表示离散变量；i，j分别为x，y方向的单位向量；，分别为x，y方向的空间步长。其中对流项显式处理，用二阶精度的Adams-Bashforth格式离散：

。

然而，这种方法需要处理非常庞大、形状不规则、刚性很强的稀疏矩阵，计算量非常巨大且不易收敛。所以很少使用。实用的求解不可压缩流动的数值方法都需要把连续方程和动量方程在一定程度上进行解耦。这种解耦处理方法可以使求解效率显著提高，但也由此产生了一系列的新问题。

§2MAC方法

（F.H.HarlowandJ.E.Welch,Numericalcalculationoftime-dependentviscousincompressibleflowoffluidwithfreesurface,Phys.Fluids8,2182(1965).）

一、基本方程

考虑不可压缩NS方程，

（1）

（2）

（3）

在上面的方程中，压力场相当于一种约束条件，使得由（2）、（3）求得的速度场满足连续方程（1）。联立求解（1）～（3）可以求得速度和压力，但计算量很大。为了高效率的计算压力，MAC方法中首先要导出压力满足的方程，称为压力波松方程。对上述方程进行变换，整理后得：

（4）

其中

考虑连续方程，有。但是，我们暂时保留与D有关的量。（4）式中左侧的第2、3项可以简化为：

，

这样，（4）式可以化为：

（5）

（5）式称为压力Poisson方程。经上述变换后，我们把N－S方程（1）、（2）、（3）式用（5）、（2）、（3）式代替。注意到：一般来说，（5）、（2）、（3）式与原始的N－S方程并不严格等价。也就是说，满足（5）、（2）、（3）式的解并不能严格保证。这就是为什么要在（4）、（5）式中保留与散度相关的量的原因。

二、求解步骤

下面介绍MAC方法的求解步骤：

i）求解压力Poisson方程（已知，求）

（5）式的时间导数项用前差离散：

我们希望，所以取。数值实验表明，这种处理方法可以使散度的绝对值保持在非常接近于零的水平，从而不影响数值解的有效性。此时（5）式写为：

（6）

对于二维问题，其中的Laplace算子可以用标准的5点中心格式离散，即

。

（6）式中的其他项用中心差分离散。（6）式称为离散压力泊松方程，可以用求解泊松方程的常规数值方法求