一种非定常NS方程并行求解设计.pdfVIP

中国]：程热物理学会 热机气动热力学 第Lf．=二届年会 编号：062006 一种非定常N-S方程并行求解设计 李雪松 徐建中 中国科学院工程热物理研究所 北京：100080 TEL．-010—mail：1 ixs78@sina．com 摘要： 为，解决汁算流体力学(CFD)中定常非定常计算与越来越人的计算量，并行计算已成为一 种现实有效的选择。论文首先研究了一种并行区域分解策略，该策略简单而商效，但需要算法配合． 为此，采用了一种与并行完全兼容的隐式方法DP．LUR方法。通过双时Il；lJ步长法。将DP．LUR方法延 伸应用到非定常计算中面不改变其原有的性质．最后分折了并行编程中的主要难点，提出解决方法， 即采用中间数据分离节点下标|j处理，并给出了并行程序的总体结构． 关键词：并行非定常隐式DP·LUR方法双时I．日J步长法 1 前言 叶轮机械计算流体动力学(CFD)技术正在由三维粘性定常计算向非定常计算方向 发展。一方面，从严格意义上说，真实湍流都是非定常的，这一点体现在湍流高级模拟 方法直接数值模拟(DNS)与大涡模拟(LES)中。另一方面，由于动静叶排间的相对运 动与相互干扰，以及每个叶片排流场的周向非均匀性，叶轮机械内部流场本质上是非定 常的，为了准确模拟叶轮机械多级流动，采用雷诺时均方法(RANS)也需要非定常计算。 然而，!{F定常计算的计算鼍远远大丁．定常汁算，这限制了叶轮机械非定常计算的应用。 近年来，随着计算机速度的进一步提高，特别是高速局域网技术的成熟，组建高经 济、高性能的机群系统已经成为可能。因此，并行计算也已经成为计算加速的一个重要 手段，越来越得到重视。而隐式算法则是定常计算有效的加速收敛的方法之一。冈此， 隐式与并行理应成为非定常计算的重要手段。 对于隐式算法，一类著名的分解方法是ADI方法【l21。ADI方法在二维计算时效果非 ，常好，f日众所周知，ADI在三维计算时存在明显的缺陷，如三维ADI方法需要更多的分 的方法：LU冈子分解方法pl。与ADI方法的比较显示LU方法更有吸引力【4】。随后，LU 、并行编程模式包括芡享内存模式与分布式内存模式。芡享内存模式可以I主tFortran或 C编译器自动完成，较为简单。但这一模式效率低、适用面窄。分布式内存模式编程较 难，但适用面广，可以应用独立于机器的通信库，如PVM、MPI等，因而可以使得同一 程序，运行于不同的讨+算平台上。分布式内存模式需要算法的配合。显式算法具有内在 的并行性，并行与串行运行结果町以完全一致。由于边界难以处理，大规模并行计算条 件F，隐式算法的实施复杂程度与收敛速度都受剑很大影响【9】，甚至影响剑计算稳定性。 因此，采用良好的并行逻辑以及与之匹配的非定常并行隐式算法足非常关键的。 652 2主导方程 下的张量表达式如卜．：． 堕+堡：翌 (1) a西考J d专J 可将方程(1)展开如下： (’7I 挈+箕+罂+罂：娑+掣+等 —(2) ．．=_-一‘一0-一=：一+一．．-一 at a毛arl a《 a芎 ℃q 蚕‘ Q=帕 舯矧耋：：嚣：端 H=JUIF+毛p+‘lH、 pw 其中9=【p∥户VpE]7为守恒变量，F、G、H为对流通量，F’、G’、H’ 为粘性通量，其具体形式可以参见一般的计算流体力学讲义。 3并行策略 奉文将采取的=}1：行策略为一