第13章计算流体力学cfd5.ppt

第13章 计算流体力学CFD(5) 数值方法：SIMPLE方法 SIMPLE是Semi-implicit method for pressure-linked equation (压力耦合方程的半隐式算法)的缩写。 数值方法：SIMPLE方法 SIMPLE算法的步骤如下： 1）在右图所示的交错网格上分别给出 数值方法：SIMPLE方法 SIMPLE算法的步骤如下： 2）求出 采用动量方程求解。 数值方法：SIMPLE方法 2） 的求法： X方向的动量方程： 数值方法：SIMPLE方法 2） 的求法： 在a点： 在b点： 交替方向隐式(ADI)方法 第一步： 简化为三对角形式 交替方向隐式(ADI)方法 第一步： 对每一个固定的j，对所有的i联立形成方程组。 对不同的j，重复上述过程。 交替方向隐式(ADI)方法 考虑二维热传导方程： 第二步： 时间步长为 ，空间导数采用中心差分，只对y的导数采用隐式处理。 交替方向隐式(ADI)方法 第二步： 简化为三对角形式 交替方向隐式(ADI)方法 第二步： 对每一个固定的i，对所有的j联立形成方程组。 对不同的i，重复上述过程。 交替方向隐式(ADI)方法 两步结束之后，T在时间方向上推进了一个时间步长?t. 考虑二维热传导方程： 推进过程只涉及三对角方程组。 交替方向隐式(ADI)方法 第一步，差分方程的x方向是隐式的。 考虑二维热传导方程： 所以这种方法叫交替方向隐式方法(Alternating Direction Implicit, ADI) 第二步，差分方程的y方向是隐式的。 交替方向隐式(ADI)方法 考虑二维热传导方程： ADI格式对t,x,y都是二阶精度的 截断误差为： 6.7 压力修正法及其在不可压粘性流动中的应用 压力修正法及其在不可压粘性流动中的应用 不可压无粘流动受椭圆型偏微分方程控制（不可压欧拉方程），松弛法是求解椭圆型问题经典的数值方法，本质上是一个迭代过程。 压力修正法及其在不可压粘性流动中的应用 不可压粘性流动的控制方程是不可压的N－S方程，这个方程具有椭圆型和抛物型的混合特性，松弛法不是特别适用。 压力修正法及其在不可压粘性流动中的应用 压力修正法也是一种迭代过程，在不可压N－S方程的数值求解中得到了广泛的应用。 6.7.1 不可压N－S方程 不可压N－S方程 假设?=常数，?=常数，可压缩N－S方程转化为不可压N－S方程： 上述四个方程封闭，含 四个未知数。 6.7.2 交错网格的应用 交错网格的应用 二维不可压流体的连续性方程为： 中心差分格式为： 右上角是u的值，左下角是v的值 速度会出现右图的棋盘式分布 交错网格的应用 右上角是u的值，左下角是v的值 可压流动中不会发生右图的问题，因为连续性方程中包含了密度对时间和空间的变化。 在可压缩流动中，右图速度的棋盘分布经过一个时间步就会被抹平。 交错网格的应用 二维不可压流体压力梯度采用中心差分： 棋盘式的离散压力分布 压力会出现右图的棋盘式分布 交错网格的应用 在交错网格上使用中心差分就不会出现速度和压力的棋盘式分布问题。 交错网格 交错网格的应用 在(i-1,j), (i,j), (i+1,j), (i,j+1),(i,j-1)等图中的实心原点上计算压力 交错网格 交错网格的应用 在(i-1/2,j), (i+1/2,j)等图中的空心原点上计算u 交错网格 在(i,j-1/2), (i,j+1/2)等图中的空心原点上计算v 交错网格的应用 连续性方程在网格点(i,j)的中心差分表达式为： 交错网格 6.7.3 压力修正法的基本原理 压力修正法的基本原理 压力修正法本质上是一种迭代法，思路如下： 1）迭代开始时，先给定压力的初始近似p* 2）用p*的值从动量方程中求解u,v,w,得到与p*有关的u*,v*,w* 压力修正法的基本原理 压力修正法本质上是一种迭代法，思路如下： 修正后的压力为 3）将u*,v*,w*代入连续性方程，它们不一定满足连续性方程。用连续性方程构造压力的修正量 ,加到p*上，使速度场满足连续性方程。 修正后的速度为 速度修正量 可以从 得到。 压力修正法的基本原理 压力修正法本质上是一种迭代法，思路如下： 4） 用步骤3)中修正后的压力做为新的p*，回到步骤2)。重复这个过程，直到速度场满足连续性方程为止。 这样就得到修正好了的流场。 6.7.4 压力修正公式 压力修正公式 压力修正公式为： 压力修正公式 压力修正公式为： 上述压力修正公式具有椭圆型的性质，可以用松弛法数值求解。