非结构网格数值计算并行处理的重排序算法研究.pdfVIP

． 第十一届全国激波与激波管学术会议 ————————————————————————————————————————————————————二．一- 一 非结构网格数值计算并行处理的重排序算法 吴颖川、王兰、乐嘉陵 中国空气动力研究与发展中心吸气式高超声速技术实验室，四川，绵阳，621000 摘要：本文介绍了非结构网格数值计算的并行处理技术，讨论了非结构网格并行处理时的重排序和分组算法。 0引言 非结构网格并行计算的分区处理是一个比较复杂的问题：给定P个处理器，整个计算区域 分为P个子区域(子网格)，每个处理器计算一个子网格。由于子网格之间的相互依赖性，处理 器之间需要进行通讯。对于frontalgrid generator[1]，节点按空间排列；Delauney generator[2]或者自适应网格[3]节点随机排列。很多网格生成器并没有对网格重新排序。 对非结构网格进行重新排序，可以减少不同处理器网格之间的相关性，提高并行效率。 ll 重排序方法包括ReverseCuthi Poole McKee(RCM)[4]、Gibbs 法和Farhat[7]的贪婪法。RCM和GFS用于减小并行网各块之间的通讯带宽；贪婪法按照相邻的 网格点进行块分配，每一块中的节点在内存中是连续的，这种方法基于结构网格的经验。 对于隐式方法，比如SGS和LUSGS等，非结构网格没有结构网格的扫描线，通过重新排序 可以确定非结构网格的扫描顺序。对于带预处理的Krylov字空间方法[8]，重排序对收敛性 影响很大。 本文给出了我们正在开发的高超声速非结构网格计算软件中采用的非结构网格的重排序 算法。算法分为排序(reorder)和分组(color)两步。排序基于带宽减小方法；分组的方法以 减少组内节点的相关性和提高高速缓存(cache)命中率为目的。 1带宽减小方法 这一类重排序方法使稀疏矩阵的非零项靠近对角项，用于矩阵预处理和网格分块。 Cuthill Poole 有几种减少通讯带宽的算法比如Reverse McKee(RCM)[4]和Gibbs 网格在内存中是连续的。在我们的程序中采用的是GPS方法，图1是40个非结构网格点排序前 和排序后的网格连接图。 2分组 分组的主要目的避免相连的节点有相同的标记。分组的原则是：一个组的节点与其他组 的节点相连，同一组内的节点之间没有连接关系。 算法可简单描述为： 1) 设定一个起始点，标记为1，设I=l； 2) 将所有没有标记的与标记为I的点相连的点标记为I+i，I=I+I； 3) 如果还有没有标记的点，重复第二步； 4) 检查每一个标记的节点，如果同一标记的节点之间有连接关系，则继续将他们分为 两个标记。 图2是对NACA0012算例的非结构网格进行分组标记后的示意图，相同颜色代表同一组的网 格点。图3是从不同起始点开始分组的示意图。 3算例 两个NACA0012翼型跨声速流动的典型算例 A．M．=0．8，攻角Q=1．25。。计算所得流场的压力云图和马赫数云图见图4(a)(b)。 366 非结构网格数值计算并行处理的重排序算法 自适应后的网格见图4(C)，翼型表面的压力系数分布与[10]的结果比较见图4(d)。 B M。=1．2，攻角a=7．0。。计算所得流场的压力云图和马赫数云图见图5(a)(b)。自 适应后的网格见图5(C)，翼型表面的压力系数分布与[10]的结果比较见图5(d)。 4结论 我们在正在开发的并行高超声速非结构计算软件上成功实现了重新排序算法， 重新排序的效果有待于用更多的算例进行验证。 参考文献 Ell and frontal fornode J．D．Muller．P．L．RoeH．Deconinck．May．1992．Aapproach