《计算流体力学》2-迭代与时间相关法.docVIP

§3. 定常问题的迭代解法 1）在实际的计算中，网格点的数量非常多，需要求解大型的代数方程组。为此，可直接调用编程语言（C/C++、FORTRAN、Matlab等）的库来求解，也可使用下面讨论的迭代解法。 Laplace方程 的差分近似为 将它改写成 或 其中 当 时， 。 给定迭代的初始近似 ，最简单的迭代公式为 （ ） 称为简单迭代法。 2）简单迭代为给定迭代的初始近似 。 对 ： 在流场内的每一个网格点（内点）上用迭代公式 进行计算。OABC、OE、ED上，可直接将边界条件指定的流函数值代入。 在边界 CD上的每一个网格点（不包括端点C、D）处，用公式 计算。 对所有网格点，计算 。 若 ，迭代收敛。 加 ，返回步骤3。 3））时，上一个网格点（）处的流函数值已经从 更新为 ；而网格点（）位于上一行，其流函数值也已经从 更新为 。迭代算法可以充分利用这些最新得到的结果，这样就构成了高斯－赛德尔迭代法 如果迭代计算的过程中对网格点的扫描是按其他顺序进行的，则上式右边会有所不同。但高斯－赛德尔迭代法的基本思路总是不变的，即：充分利用已经得到的最新结果。 采用高斯－赛德尔迭代法，通常能够加快迭代的收敛，但并不总是如此。甚至可以找到反例，对某种特定的代数方程组，用简单迭代法收敛，用高斯－赛德尔迭代法反而不收敛。 4） ，即 将它与 的加权平均作为迭代