代数多重网格求解器及其在不可压粘笥流动数值模拟中的应用.pdfVIP

代数多重网格求解器及其在不可压粘性 流动数值模拟中的应用 陈兵 李玉成 赖国璋 (大连理工大学。大连116024) 摘要 本文在代数多重网格法中采用较新的方法构造插值算子和祖网格方程组系数矩阵。构造了求 w线性代数方超组的快这有效，适用范围广的代数多重网格法求解器.并将该求解器用于不可压粘性流 体运动的\-S方程分步解法中压力Poisson方程的求解.用本方法进行了二维方腔拖带流和近底水平皿 柱绕流问题的数值模拟，与前人和物模实验结果相符合，检验了该求解5适应不规则网格的能力，表明 本文的代数多重网格法求解器是高效快速的。尤其适合网格结点多的大规模问题的求解。 关”代”多重1148t·不可压”性流体.Poisson71气 、 前言 在让鱼园蛀助尝 (CFO)中，控制流体运动的偏微分荃本方程通常用有限元及有限 差分等方法进行离散，用离散的偏微分方程的解来近似原米的偏微分方程的解。线性偏 微分方程离散后形成线性代数方程组:而非线性偏微分方程一般经过线住化处理，形成 一系列线性代数方程组。通过迭代过程求解。因此，求解线性代数方程组是计算流体动 力学研究中重要的环节。有限元和有限差分等方法形成的线性代数方程组的系数矩阵， 多为大型的正定、对角占优的稀疏矩阵。对这种矩阵的求解可用 Gauss消元法等直接解 法，也可用各种迭代方法。如Jacobi迭代、Gauss-Seidel迭代等。直接解法需要的存储 量大，要存储带宽内的所有元素，而迭代法只存储系数矩阵的非零元素，可大大节省存 储量。前人研究表明。线性代数方程组的阶数N和求解该方程组的计算量 T与求解空间 的解析度Q、空间维数d的关系可表示如下 (Webster,1996) N一Q0 (1) T一Qar (2) 对于Gauss消元法，夕等于3，由此可见，对于空间结点数较多或三维情况下 (ct-3)，用 直接解法是不可取的.最理想的情况是61等于1，此时求解的计算量与空间解析度为线性 关系，即计算量随着方程组中方程个数的增加而线性增大。一般的点松弛迭代方法，例 如Jacobi或SOR点迭代格式中，某一点上解的校正是建立在相邻的点上的信息，误差校 正信息的传播是逐点推进的，因此不可能通过点松弛法的一次扫描对整个求解空间进行 误差的校正，从而计算量与空间离散点数不成线性关系。Webster(1996)得出点松弛方 .280 . ! 1 月 J 马 法中Q可表示为 夕二1+21d (3) d为空间维数。正是由于点松弛法误差校正推进的特点，使它对网格尺度波长的误差消减 很快，而对长波误差消减较慢，即它对消减误差的高频分量效率较高。而对误差的低频 分量效率较低。为了克服普通点松弛迭代法的这个缺点，多重网格法是较好的选择。简 单地说，多重网格法在最初的细网格的墓础上，通过选抒一部分结点构成新的粗网格， 并基于此粗网格构造新的线性代数方程组，此方程组比原有的方程组规模缩减，但性质 接近。显然，用此粗网格方程组结合适当的点松弛迭代法，可以有效消减大于原来的细 网格尺度的