北航计算流体力学第11课.docVIP

用于求解Laplace方程的网格系统 , , 式中 求解三对角方程的追赶法简介：（以方程组为例） （1） （2） （3） （4） （5） 由式（1）得： （即） 代入式（2）得： （即） 继续上面的工作，可得： 最后求得： 然后回代，依次求出从而完成了求解工作。 以下是一段求解任意三对角方程组的Fortran程序 Subroutine TRIDAG SUBROUTINE TRIDAG (AA,BB,CC,FF,N1,N) C C******************************************************************************************* C SOLUTION OF A TRIDIAGONAL SYSTEM OF N-N1+1 EQUATIONS OF THE FORM C C AA(K)*X(K-1) + BB(K) * X(K) + CC(K) * X(K+1) = FF(K) K=N1, ?,N C C K RANGING FROM N1 TO N C THE SOLUTION X(K) IS STORED IN FF(K) C AA(N1) AND CC(N) ARE NOT USED C AA,BB,CC,FF ARE VECTORS WITH DIMENSION N, TO BE SPECIFIED IN THE C CALLING PROGRAM C******************************************************************************************* C DIMENSION AA(1),BB(1),CC(1),FF(1) BB(N1)=1./BB(N1) AA(N1)=FF(N1)*BB(N1) N2=N1+1 N1N=N1+N DO 10 K=N2,N K1=K-1 CC(K1)=CC(K1)*BB(K1) BB(K) =BB(K)-AA(K)*CC(K1) BB(K) =1./BB(K) AA(K) =(FF(K)-AA(K)*AA(K1))*BB(K) 10 CONTINUE C C BACK SUBSTITUTION C FF(N)=AA(N) DO 20 K1=N2,N K=N1N-K1 FF(K)=AA(K)-CC(K)*FF(K+1) 20 CONTINUE RETURN END Laplace方程 （6-1） 五点格式： （6-2） 其精度为。 设，五点格式写成如下形式 （6-3） 式中，；。 迭代法框图 一．点迭代法 1．Jacob迭代法（简单迭代法） 将式（6-3）写成 （6-4） 2．Gauss-Siedel点迭代 充分利用新值，可将式（6-4）写成 （6-5） Gauss-Siedel迭代比Jacob迭代快一倍。 二．线迭代法 1．Jacob线迭代法（简单线迭代法） （6-6） 构成三对角方程，可用追赶法求解。 2．Gauss-Siedel线迭代 （6-7） 一般都采用Gauss-Siedel迭代。 以上是行迭代，同样可采用列迭代 （6-8） 交替方向隐式方法(Alternating Direction Implicit Method (ADI) ) 也可以先用列迭代后用行迭代。 超松驰法（Successive Over-Relaxation Method） 1．点超松驰法（PSOR） 式中，称为松驰因子 合并以上两式，得到 （6-10） 2．线超松驰法（LSOR） 将行迭代（6-7）写成 合并以上两式，整理得： 行迭代： （6-11） 列迭代： （6-12） 3．加速ADI法（Accelerated ADI——AADI法） （6-13） 几种方法的大致迭代次数 加速方式 迭代方式 Jacobi Gauss-Siedel SOR 点迭代法 2000 1000 100 线迭代法 1000 500