非结构同位网格SIMPLE类算法收敛性能比较$V-空气动力学学报.PDF

第 卷 第 期 空 气 动 力 学 学 报 ， !! #$% !! ’$ 年 月 !#! !$%’()!*+! ,+)+! ， !(() (* +,- !(() !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 文章编号： （ ） (!./01/!. !(() (0(!/*0(2 非结构同位网格,+*-.$ 类算法收敛性能比较 1 ! 赖锡军 ，汪德 ，傅源方 （ 中国科学院南京地理与湖海研究所，江苏南京 ； 1 !1(((/ 河海大学环境科学与工程学院，江苏南京 ； 中国科学院力学所，北京 ） ! !1((*/ 1(((/( 摘 要：以非结构网格的 算法为基础，将算法扩展为 算法。利用 角的斜方腔流动计算成果，分 +34567 +345678 (9 析了非结构同位网格的+34567: +345678 算法的收敛性能；比较了因网格的非正交而引入的非正交项的取舍对该 算法收敛性能的影响；并采用显式校正步法对+345678 算法进行了显式校正。比较表明，在非结构同位网格+340 5678 算法中可忽略非正交项，但有必要对压力作亚松弛。显式校正步法可显著地加速在非结构网格上求解’ ; + 方程的收敛性能，而且在不同的松弛因子组合下，均有较好的收敛速率。 关键词：非结构同位网格； 算法；收敛性能；显式校正步法 +345678 中图分类号： 文献标识码： #!11 用显式校正步法（ ）对 / 引 言 ,N-%=I=B I$DD,IB=$? B,- +345678 算法进行了显式校正。 在采用原始变量法求解粘性不可压缩流体时，需 0 数值方法 要求解压力场，常用联立求解法；压力校正法；压力 5$=$? 方程法；人为压缩法。其中压力校正法自 0 1 0 控制方程 1*@! 年由5ABA?CAD 与+-A%E=?F 提出+34567 算法以来， 不可压缩流动的连续方程和动量