适用于混合网格的改进雅可比迭代法及其应用.pdfVIP

2016 年 3 月 March 2016 北京航空航天大学学报 Journal of Beijing University of AeronautÎcs and AstronautÍcs Vo l. 42 NO.3 第42 卷第3 期 http:// bhxb. buaa. edu. cn jbuaa@ buaa. edu. cn DOI: 10. 13700/j. bh. 1001-5965.2015.0197 适用于混合网格的改进雅可比迭代法及其应用 黄宇1 ，阎超1 ，\袁武2 (1.北京航空航天大学航空科学与工程学院，北京 100083; 2. 中国科学院计算机网络信息中心超级计算中心，北京 100190) 摘 要: LU-SGS 因有较高的鲁棒性和小的内存需求而得到广泛应用，然而用于混合 网格计算前需要进行网格排序和分组来实现算法并行;此外，LU-SGS 格式收敛效率不高。针 对这些缺点，本文提出了一种改进的适用于复杂混合网格的雅可比迭代方法，无需网格排序和 分组就可实现算法的并行化，且有较快的收敛速度。该方法编程实现简单，易于采用 OpenMP 实现并行。算例研究表明，相比于 LU-SGS 格式，在各来流条件下，本文提出的方法收敛速度 更快，鲁棒性好，并行和串行结果一致，且内存需求增加很少。 关 键 词:混合网格;隐式算法;并行计算; OpenMP; 雅可比迭代;网格排序 中图分类号; V221 + .3; TB553 文献标识码:A 文章编号: 1001-5965 (2016 )03 -0551-11 计算流体力学(Computatíonal Fluíd Dynam- 上的基于共享内存 OpenMP 并行模式的 ícs ,CFD) 已得到广泛应用。随着需要模拟的外 GMRES + LU-SGS 隐式格式。为了并行化，该格 式事先也需以某种方法对网格进行分组，并将网 形越来越复杂，混合网格由于其能充分利用非结 格在处理器多个进程间划分。但该方法改变了 构网格的生成优势和对复杂外形的强大适应能 串行程序执行顺序，造成串行、并行结果不一 力，已得到越来越多的使用。然而，由于混合网格 致，且实现复杂。此外当采用动网格计算时，网 单元间的无序性和复杂的数据结构，相比于结构 网格其计算效率较低且存储需求较大[IJ 。另一 格排序和分组需反复进行，为求解带来了额外 的时间硝耗。 方面，显式格式如 Runge-Kutta 法等虽实现简单， 为了克服上述问题，本文提出了一种无需网 计算量小，然而在高雷诺数流动计算中，由于稳定 格排序和分组的改进雅可比迭代方法。将每个单 性限制，时间步长不能过大，计算效率不高。由于 元同其相邻单元作为独立的块处理，并行地直接 上述原因，为了在混合网格上提高计算效率，隐式 方法如 Gauss-Seídel 、 LU-SGS , BLU -SGS 和 GMRES 求解各块相应的方程组;为了使算法并行化，抛弃 等格式在近年来得到广泛的研究和应用[1 川。 了 LU-SGS 中的前扫描和后扫描模式，依次扫描