第3章椭圆方程迭代法介绍.docVIP

椭圆型问题的差分法 §3－1 流体力学中的椭圆型问题 ·无旋流场中 速度势 （Laplace Eq.） ·二维不可压定常流动，利用涡－流函数表示： ·不可压分离流问题中，扰动压力场： ·定常的N－S方程求解问题 ·在网格自动生成中，求解椭圆型方程的网格生成方法 由于椭圆型方程的数学性质：求解域内部任何一点的解函数依赖于所有边界上的边界条件，因此从数值计算方法来看，就不能从一部分边界起步进行推进计算到另外的边界，这与发展方程的求解方法有很大的差别，椭圆型方程的数值求解方法，只能是在整个流场中进行迭代计算来求解。 §3－2 椭圆型问题的迭代法求解 （一）迭代法的基本概念 例：方程 （ Poisson 方程） 二维 差分离散 ………………..(*) 写成矩阵形式代数方程组为： …………………………………….. (1) 其中 一般地，对于线性方程组有，欲求未知函数的解矢量 若为非奇异矩阵，即：，则 由于是个阶数甚大的矩阵（非三对角），直接求解，或利用Gauss消去法求逆矩阵，计算量及所需计算机的内存都将十分巨大，所以在实际计算中不希望采用直接法求解。 迭代法的基本思想是：定义一个序列，当时，，从而得到方程(1)的解。 迭代法设法给出的迭代关系。（通常为计算方便，迭代法采取，使之简单） 若（即迭代关系式）与迭代步k无关，则称为平稳迭