CFD计算中某些并行方法的讨论研究.pdfVIP

CFD计算中某些并行方法的讨论 OFD计算中某些并行方法的讨论 刘鑫1郑忠华2赵慧勇2杨顺华2陆林生1 14083 1江南计算技术研究所无锡2 2中国空气动力研究与发展中心绵阳621000 摘要：本文主要讨论了复杂组合体的CFD数值模拟中的某些并行方法。首先介绍单区域分解 策略需要考虑并行效率、通信优化等问题，并提出适合多区域分解的静态负载平衡算法。在重 叠网格和拼接网格技术方面，主要介绍了找重方法和插值方法。另外，本文还讨论了隐式方法 求解线化方程时不完全的Lu分解可能带来的解的扰动问题以及当前隐式方法带来的串并行结 果不一致的问题，并给出了解决方法。实验结果表明上述讨论的并行方法是有效的。 1引言 本文主要讨论复杂组合体的CFD并行计算，控制方程是考虑化学非平衡反应的三维NS方 程： 二+——+——+——=—』士—二+—二+．Y (11) 警+娑+警+警：婴+孥+孚+s (1．1) dt dx dy dZ d)c d))dZ 式中：Q=p，pu，pv，pw，肛，胪，y是守恒变量；F、G、E是无粘通量向量； E、G，、E。是粘通量向量；S是化学反应源项。湍流模型采用k一￡二方程模型，数值计算 方法采用隐式有限体积法u1，采用消息传递模型编写了MPI程序。本文着重讨论这一课题的并 行计算方法，并重点在区域分解方法、重叠网格和拼接网格处理方法以及隐式求解方法等方面 进行一些探讨。 本文首先讨论了单区域分解策略需要考虑并行效率、通信优化、并行可扩展性以及工程实 现可行性等问题，并提出适合多区域分解的静态负载平衡算法：一维优先的规则分块算法和基 于自动重分块的不规则分块算法，分别处理中小规模和大规模并行的问题。在重叠网格找重方 法上，主要介绍了线性搜索方法、八叉树搜索方法以及基于k-d树的搜索方法；在拼接网格插 值方法上，介绍了基于型函数的双线性插值方法。另外，本文讨论了隐式方法求解线化方程时 不完全的LU分解可能带来的解的扰动问题，并提出牛顿法和预条件子共轭梯度法来修正线化 方程的解，提高解的稳定性。在线化方程并行求解时，如果虚边界用上一时刻的值，不是严格 的隐式方法，因此引入流水线并行和迭代法尝试解决弱隐式方法带来的串并行结果不一致的问 题。实验结果表明上述讨论的并行方法是有效的。 2区域分解方法 2．1概述 一般地，我们将相同的指令集作用于数据集所有元素的并行计算方法称为数据并行方法 旧1。区域分解是数据并行方法的核心，是在科学计算中常用的一种并行方法。由于在航空航天 CFD数值计算中大多采用Euler坐标系，在计算坐标系下利用结构化网格离散化后，求解空间 一般由一个(称为单区域)或多个(称为多区域)网格单元规则分布的区域组成，物理量在空 间坐标方向的依赖关系往往呈现出边缘依赖性，使区域分解方法(在CFD中常称分块方法)在 航空航天CFD得到了广泛应用蠊…。 按域划分方式区域分解方法允重迭区域和非重迭区域两种分解方法。设求解区域Q分解为： Q=IIQ； (2．1) 若Q，nQ，=0(1≠，)，唰称非重迭区域分解方法，否则是重迭区域分解方法。按虚边 界处理方法划分，区域分解还可分为通信型和类物理边界两种分解方法。对于通信型分解方法， 当前区域虚边界的物理量由相邻块的物理量通信得到。对于类物理边界分解方法，当前区域虚 387 第十一届全国激波与激波管学术会议 边界的物理量类似于真实边界一样给定边界条件，各个区域可以划分不同类型网格，独立求解， 可获得粒度粗并行，求解方法也灵活，但是要研究新的计算方法，例如Schware迭代技术，而 且如果虚边界条件给定不合适，可能导致收敛于非物理解。由于本文的研究课题没有采用这种 方法，对此方法不进一步讨论，只对单区域和多区域分解方法进行进一步的讨论。 2．2单区域分解的某些基本原则 单区域分解是