(7.3.1)--6.5-松弛迭代法讲义.pdfVIP

《计算流体力学基础》讲义 第六章 简单的计算流体力学技术 6.5 松弛迭代法 松弛技术是一种适用于椭圆型差分方程的求解技术。在空气动力学中，低速、亚声速无 黏流动的控制方程为椭圆型偏微分方程，因此可以运用松弛迭代方法来进行流场求解。它可 以是显式的，也可以是隐式的。本课程仅讨论显式松弛技术，也称为点迭代技术。 假设求解的问题是二维无黏无旋不可压流动，流场速度势函数为 ，其控制方程是  Laplace 方程(6.31) ： 2 2     2 + 2 0 (6.31) x y 假设在均匀网格上求解该方程，那么其二阶中心差分格式如式(6.32)所示。  −2 +  −2 + i+1, j i , j i−1, j i , j+1 i , j i , j−1 2 + 2 0 (6.32) (x ) (y ) 假如求解域如图6.4 所示，边界由20 个离散点组成，内点共有15 个。速度势1 ~ 20 在边界上的值1 ~ 20 均为已知的，而在其他15 个内点上的值都是未知的。式(6.32)可以 在任何一个内点使用，这样便可以得到一个具有 15 个线性代数方程系统，该系统共包含15 个未知数。如何求解此线性方程组系统。通常有三类求解方案。第一种是采用标准 Cramer 法则，然而，其需要的运算次数是相当大的。就图示的例子来说，需要求解一个 15×15 的 行列式的值。对于真实问题，网格点数成百上千万，使用Cramer 法则求解该类问题是行不 通的。第二种较为合理一些的方法是高斯消去法。不过，最简单的方法是松弛技术。 图 6.4 松弛迭代法示意图 n 假设方程(6.32) 中4 个变量是已知的，并且认为它们的值是第 次迭代步上的值，同时 - 1 - 《计算流体力学基础》讲义 假定只有一个变量是未知的，而且它是当前第n +1步的值。选择 作为未知量，由式(6.32) i, j 得： 2 2 n +n n +n  n+1 (x ) (y ) i+1, j i−1, j i , j+1 i , j−1 i, j 2 2  2 =+ 2  (6.33) 2(x) +2(y )  (x) (y ) 