基于区域分解并行SIMPLER算法研究.pdf

第14卷3期 应用基础与工程科学学报 Vol，14，No．3 2006 2006年9月 OFBASICSCIENCEANDENGINEERINGSeptember JOURNAL 文章编号：1005-0930(2006)03-0316-08中图分类号：V211 文献标识码：A 肖曼玉， 欧阳洁， 李永刚 (西北工业大学应用数学系，陕西西安710072) 摘要：在同位网格上建立了基于重叠型区域分解的并行SIMPLER算法，给出了 区域划分和数据交换的实施细则．并在集群系统下通过对顶盖驱动方腔流进行 模拟验证了算法的可靠性．数值结果表明，加速比成线性增加，随着计算规模的 增大，并行效率可达到95％以上． 关键词：区域分解；不可压流；SIMPLER算法：并行计算 流体流动的数值模拟广泛应用于气象、航天、机械、采矿等自然研究和工程计算的各 for Linked MethodPressure Revised)算法是目前 个领域．SIMPLER(Semi—Implicit Equations 在流动、传热和燃烧等领域内关于流场求解较为成熟的方法，已经在计算流体力学和数值 传热学界得到了广泛的应用． 近年来，由于计算流体力学研究问题规模的不断增大，现有计算机资源渐已无法满足 其日益增长的需求．随着高性能并行计算与通信的迅速发展，针对流体流动数值模拟的并 国内外备受关注的热点和前沿课题．鞠立力等…对非定常不可压N-s方程压力修正投影 法的并行算法进行了研究，重点构造了压力Poisson方程快速求解的并行算法．李树民 等[21找到一种适合压力修正算法的内边界耦合条件，初步实现了不可压N—S方程的隐式 并行．尤其值得注意的是，近年来发展迅速的区域分解算法因其适用于构造并行算法的特 点而逐渐被推广到计算流体力学领域．J 孔渗流的模拟． 本文基于区域分解原理进行了同位网格下并行SIMPLER算法的研究，给出了该算法 实施的具体步骤以及实现数据交换与负载平衡的措施，并通过数值试验验证了算法的可 靠性，考察了算法的并行效能。 1 N—S方程离散及基于同位网格的SIMPLER算法 考虑如下二维不可压缩流的N-S方程 V·H=0 (1) 塑+M·V (2) tOp7一比h=H Rg 收稿日期：2005-06—16；修订日期：2006-07-03 基金项目：国家自然科学基金重大项目；陕西省自然科学基金(2005A16) 作者简介：肖曼玉(1980一)，女，博士研究生． 万方数据 No．3 肖曼玉等：基于区域分解的并行SIMPLER算法研究 317 式中球=(“，秽)表示所求速度向量，凰表示雷诺数，P代表压力． 1．1 交错网格与同位网格 若基于图1所示的交错网格系统H o对区域进行离散，速度u、口与压力P将分别存储 于三套网格．并行计算时，将导致程序设计的繁琐以及数据结构和信息传递设计等方面的 诸多不便．在此，我们采用如图2所示的同位网格∞o离散区域，速度“、秽和压力P采用相 同的控制容积．同时，运用动量插值法∞J，通过对控制容积界面流速的离散，引人相邻两 点的压差来保证压力与速度的耦合． _ - ● 一 } ● ● ●