非结构网格Euler+并行算法研讨.pdfVIP

非结构网格Euler 并行算法研究 周磊 朱培烨 蔡江鹏 中国航空计算技术研究所 摘 要 对于非结构网格上的 Euler 流动的并行计算问题，本文针对同一网格，对比了坐标递 归二分法，惯性递归二分法和多级递归图二分法的区域分割结果。另一方面，对比了粗细两种 非结构网格方法的并行计算结果，得到它们的计算效率变化规律，分析这些变化规律的成因， 对实际工程中并行加速比的定义提出了自己的看法。 关键词 并行计算，加速比 1. 引言 为对飞行器外流绕流流场进行数值模拟，人们常常采用非结构网格来方便地对复杂的几何 外形进行空间离散。注意与结构化网格相比，非结构网格计算需要大得多的内存空间来存储单 元间不规律的拓扑连结关系，同时计算量也增加了很多。 所以，非结构网格上的流场数值模拟问题常常是计算规模庞大的工程问题，对此类问题进 行并行计算研究正是当前计算流体力学领域的热点问题,见[1][3] 。 本文通过对一个简单标模的非结构网格 Euler 流动的并行计算，讨论大规模工程计算所涉 及的区域分裂及加速比方面的问题。 2. 数值方法 我们采用四步龙格库塔有限体积方法对三维M6 机翼的定常绕流流场进行 Euler 数值模拟。 采用了局部时间步长等加速收敛的方法。对计算区域采用四面体单元的非结构网格进行空间离 散。 计算区域见下图： 图 1 计算区域示意图 图2 机翼附近局部放大图 3. 分区计算方法 分区算法的思想是将总体计算网格分成若干个块子计 算网格，每块计算网格单独求解，在每个计算时间步后各子 网格在其相邻子网格的边界处的解通过消息传递机制 1 (massage-passing )进行数据交换。区域分解的基本原则是，尽量使得每个子计算网格包含相同数 量的计算单元。 并行计算平台是 12 台 SUN 工作站以局域网络相联接，配以 MPI 并行通讯软件库构成，每 台工作站内存为 128 兆。 我们采用粗细两套网格（粗网格 5.6 万网格点，细网格 9.3 万网格点）来研究不同计算规模 的问题的并行计算效果。 我们分别采用了坐标递归二分法，惯性递归二分法（由坐标递归二分法改进而来）（见[4]） 和美国明尼苏达大学 George Karpis 等人提出的多级图分区理论及其相应的 Metis 软件进行负载平 衡分区，并比较了这三种分区方法对于同一问题（分割同一套密网格）的分区结果。 图3 总体网格分区图 图4 局部网格分区图 上面是用 Metis 程序包将细网格的 12 块子网格进行分区后的结果图，左图是全貌图，而右 图是左图中小黑框部分的局部放大图，注意：由于网格线消隐的作用，有些子区域看不到。经过 统计，上述分区是比较平均的（每块网格的单元数离理想平均单元数的最大偏差在 3%之内） 4. 结果与分析 首先，粗网格 12 块的并行计算的结果与单块计算结果一致，细网格 12 块的并行计算的结果 与 2 块并行计算结果一致。且并行计算的区域块数与收敛速度无关。 得到下面的各种非结构网格分区方法的计算结果对比表： rcb rib metis 平衡偏差 0.1845636 0.1205518 3.564 1856E-02 通讯量 31653 26748 1660 1 通讯次数 110 100 84 在上表中，RCB 表示递归坐标二分法，RIB 表示递归