第13章_计算流体力学CFD(5)

,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,第,13,章 计算流体力学,CFD(5),6,计算流体力学的基本方法,6.1 Lax-,Wendroff,方法,Lax-,Wendroff,方法,Lax-,Wendroff,方法是一种显式有限差分方法，适合于推进求解。,二维时间推进网格,Lax-,Wendroff,方法,Lax-,Wendroff,方法在时间和空间上都具有二阶精度。,二维时间推进网格,Lax-,Wendroff,方法,非定常二维无粘流（欧拉方程非守恒形式）：,Lax-,Wendroff,方法,Lax-,Wendroff,显式推进求解,(,沿时间方向进行泰勒级数展开,),：,Lax-,Wendroff,方法,空间导数采用中心差分：,Lax-,Wendroff,方法,(,),求对时间,t,的二阶导数：,Lax-,Wendroff,方法,Lax-,Wendroff,方法,Lax-,Wendroff,方法,Lax-,Wendroff,显式推进求解 ：,6.2,MacCormack,方法,MacCormack,方法,MacCormack,方法在时间和空间上都具有二阶精度。,MacCormack,方法是一种显式有限差分方法，适合于推进求解。,MacCormack,方法比,Lax-,Wendroff,方法应用起来更简单。,MacCormack,方法,校正步,预估步,MacCormack,方法,预估步：空间导数用向前差分计算。,MacCormack,方法,预估步：空间导数用向前差分计算。,预估值：,MacCormack,方法,校正步：空间导数用向后差分计算。,MacCormack,方法,MacCormack,方法,在,MacCormack,方法中，预估步用向前差分，校正步用向后差分；也可以预估步用向后差分，校正步用向前差分。或者在时间推进解法的相继两个时间步中轮流使用这两种办法。,6.3,粘性流动、守恒形式和空间推进,6.3.1,粘性流动,粘性流动,粘性流动的控制方程是,N-S,方程。,对定常流动，,N-S,方程的数学性质更多地表现为椭圆型的，不能采用,Lax-,Wendroff,方法和,MacCormack,方法求解。,对非定常流动，可以采用,Lax-,Wendroff,方法或,MacCormack,方法求解,N-S,方程。,6.3.2,守恒形式,守恒形式,可以采用,Lax-,Wendroff,方法或,MacCormack,方法求解,U,的分量在各时间步的值。,非定常守恒形式欧拉方程（二维）：,6.3.3,空间推进,空间推进,定常守恒型二维欧拉方程：,对于亚声速流动，上述方程是椭圆型的，所有空间推进方法都不适用，,MacCormack,方法也不适用。,空间推进,对于超声速流动，上述方程是双曲型的，空间推进方法适用，,MacCormack,方法也适用。,定常守恒型二维欧拉方程：,空间推进,MacCormack,方法：,定常守恒型二维欧拉方程：,空间推进,预测步,：（向前差分）,预估值：,空间推进,预估值：,空间推进,校正步,：（向后差分）,6.4,松弛法及其在低速无粘流动中的应用,松弛法及其在低速无粘流动中的应用,松弛法特别适合于求解椭圆型偏微分方程，常被用来求解无粘亚声速的低速流动。,松弛法及其在低速无粘流动中的应用,考虑无粘不可压流体的二维无旋流动，控制方程为,Laplace,方程：,松弛法及其在低速无粘流动中的应用,松弛法是一种迭代法,上标,n,和,n+1,表示迭代次数,松弛法及其在低速无粘流动中的应用,松弛法是一种迭代法,松弛法及其在低速无粘流动中的应用,松弛法是一种迭代法,松弛法及其在低速无粘流动中的应用,松弛法是一种迭代法,从左至右扫描,松弛法及其在低速无粘流动中的应用,松弛法是一种迭代法,当所有网格点处的 都小于一个预定的值时，迭代收敛。,松弛法及其在低速无粘流动中的应用,运用逐次松弛法可加快收敛的过程。,从左至右扫描,从下至上扫描,松弛法及其在低速无粘流动中的应用,运用逐次松弛法可加快收敛的过程。,是松弛因子，如果,1,，叫做逐次超松弛法；如果,1,，叫做逐次低松弛法。,松弛法及其在低速无粘流动中的应用,运用逐次松弛法可加快收敛的过程。,选取合适的值，可以减少迭代次数，从而减少计算时间。在某些问题中，迭代次数可减少到原来的,1/30,6.5,数值耗散、色散及人工粘性,数值耗散、色散及人工粘性,一维波动方程：,差分方程：,截断误差：,数值耗散、色散及人工粘性,差分方程：,泰勒级数展开：,数值耗散、色散及人工粘性,差分方程：,将泰勒级数展开代入差分方程得：,数值耗散、色散及人工粘性,差分方程：,将泰勒级数展开代入差分方程得：,数值耗散、色散及人工粘性,差分方程：,等号右边将对,t,的偏导数转化为对,x,的偏导数得：,数值耗散、色散及人工粘性,差分方程：,偏微分方程（修正方程）：,一维波动方程（偏微分方程）：,数值耗散、色散及人工粘性,差分方程：,一维波动方程（偏微分方程）：,差分方程的精确解是上述一维波动方程的数值解（含误差）,数值耗散、色散及人工粘性,差分方程：,差分方程的精确解是上述修正方程的精确解（不含误差）,偏微分方程（修正方程）：,数值耗散、色散及人工粘性,差分方程：,偏微分方程（修正方程）：,一维波动方程（偏微分方程）：,数值耗散、色散及人工粘性,修正方程等号右端的项是截断误差，如果截断误差的主项是偶数阶导数，数值解将主要表现出耗散行为；如果主项是奇数阶导数，数值解将主要表现出色散行为。,偏微分方程（修正方程）：,数值耗散、色散及人工粘性,等号右端的偶数阶导数项起数值耗散的作用，奇数阶导数项起数值色散的作用。,偏微分方程（修正方程）：,数值耗散、色散及人工粘性,数值耗散的作用很象物理粘性，二阶导数项前的系数被称为人工粘性。,偏微分方程（修正方程）：,数值耗散、色散及人工粘性,数值耗散的影响会将波抹平,数值耗散、色散及人工粘性,色散导致波的不同相位在传播中产生畸变，表现为波前和波后出现振荡。,数值耗散、色散及人工粘性,尽管人工粘性降低了解的精度，但通常有助于提高解的稳定性。,偏微分方程（修正方程）：,6.6,交替方向隐式,(ADI),方法,交替方向隐式,(ADI),方法,考虑二维热传导方程：,等号右端有五个未知量，不能得到三对角方程组，不能采用托马斯算法（追赶法）求解。,采用,Crank-Nicolson,方法（隐式）：,交替方向隐式,(ADI),方法,考虑二维热传导方程：,第一步：,时间步长为 ，空间导数采用中心差分，只对,x,的导数采用隐式处理。,交替方向隐式,(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,压力修正公式,压力修正公式,压力修正公式为：,压力修正公式,压力修正公式为：,上述压力修正公式具有椭圆型的性质，可以用松弛法数值求解。,在不可压流场中，压力的扰动将会传遍整个流场，这与上述方程的椭圆型性质相吻合。,压力修正公式,压力修正公式为：,压力修正公式是压力修正 的泊松方程的中心差分表达式。,上述泊松方程中的二阶偏导数用中心差分替代。,式中：,压力修正公式,压力修正 的泊松方程,(,为椭圆型,),：,d,相当于一个质量源项。,6.7.5,数值方法：,SIMPLE,方法,数值方法：,SIMPLE,方法,SIMPLE,是,Semi-implicit method for pressure-linked equation (,压力耦合方程的半隐式算法,),的缩写。,数值方法：,SIMPLE,方法,SIMPLE,算法的步骤如下：,1,）在右图所示的交错网格上分别给出,数值方法：,SIMPLE,方法,SIMPLE,算法的步骤如下：,2,）求出,采用动量方程求解。,数值方法：,SIMPLE,方法,2,） 的求法：,X,方向的动量方程：,数值方法：,SIMPLE,方法,2,） 的求法：,在,a,点：,在,b,点：,数值方法：,SIMPLE,方法,X,方向的动量方程：,差分方程：,数值方法：,SIMPLE,方法,差分方程：,X,方向的动量方程：,数值方法：,SIMPLE,方法,2,） 的求法：,数值方法：,SIMPLE,方法,2,） 的求法：,数值方法：,SIMPLE,方法,2,） 的求法：,Y,方向的动量方程：,数值方法：,SIMPLE,方法,2,） 的求法：,在,c,点：,在,d,点：,数值方法：,SIMPLE,方法,Y,方向的动量方程：,差分方程：,数值方法：,SIMPLE,方法,2,） 的求法：,数值方法：,SIMPLE,方法,3,）将 和 代入压力修正公式，在所有内部网格点上求解,数值方法：,SIMPLE,方法,SIMPLE,算法的步骤如下：,4,）在所有内部网格点上计算,数值方法：,SIMPLE,方法,SIMPLE,算法的步骤如下：,5,）将 作为新的,，重复步骤,(2),至步骤,(5),，直到收敛。收敛的合理标准是质量源项,d,趋于零。,数值方法：,SIMPLE,方法,对于某些应用，压力修正公式会发散，而不是收敛，此时，可采用低松弛：,为低松弛因子，建议取为,0.8,6.7.6,压力修正法的边界条件,压力修正法的边界条件,对不可压粘性流动，如果给定下列边界条件，则物理问题是唯一确定的：,压力修正法的边界条件,1,）在入流边界上，,p,和,v,给定，,u,是变化的。,为零,给定，并保持不变,压力修正法的边界条件,2,）在出流边界上，,p,给定，,u,和,v,是变化的。,为零,压力修正法的边界条件,3,）在壁面上，给定粘性无滑移条件，于是壁面速度为零。,压力修正法的边界条件,3,）在壁面上：,在壁面附近， 很小，假设 ，则有,压力修正法的边界条件,3,）在壁面上：,压力修正法的边界条件,压力修正方程具有椭圆型性质，在计算区域的整个边界上都给定了关于压力的边界条件。,6.8,用于,CFD,的计算机绘图技术,6.8.1,xy,图,6.8.2,等值线图,6.8.3,向量图和流线图,向量图和流线图,高超声速飞行器表面上的三维向量图和流线图,6.8.4,网格图,机翼绕流计算的三维网格图,6.8.5,组合图,