三维湍流Rayleigh-Bénard热对流高效并行直接求解方法.pdfVIP

物理 学报 ActaPhys．Sin．Vo1．64，No．15(2015) 154702 三维湍流Rayleigh—B6nard热对流的高效并行 直接求解方法术 张义招 包芸十 (中山大学力学系，广州 510275) (2014年 l1月3日收到；2015年3月4日收到修改稿) 高Ra数Rayleigh—B6nard热对流的湍流特性研究是当前国际上的一个热门研究课题，DNS模拟计算是 研究该课题的重要手段之一．当计算规模增大而网格数巨大时计算工作难 以实现，高Ra湍流热对流 的数值 模拟研究面临重大挑战．本文创建了大规模高效并行计算的三维湍流热对流直接求解方法．采用FFT变换解 耦压力泊松方程，将其变换成沿Z方向上的块三对角方程组，并利用块三对角方程的MPI与0penMP联立的 大规模高效并行近似解求解方案，创建了可以高效并行计算的热对流直接求解方法．通过对该方法并行效率 的验证计算，证明新的直接求解并行计算方法具有很好的并行效率和计算时效．三维窄方腔热对流的计算结 果表明，本文方法计算的三维热对流特性是合理的．本文创建的可大规模高效并行计算的三维湍流热对流直 接求解方法，也很可能是关于计算流体力学不可压NS方程大规模高效并行计算在特殊情况中计算技术上的 一 个突破． 关键词：泊松方程直接求解，并行计算，Rayleigh—Bdnard热对流，DNS PACS：47．27．te DOI：10．7498／aps．64．154702 的流动信息，但是DNS数值模拟需要花费大量的 1 引 言 计算资源．随着计算机计算能力的不断提高，国际 上这方面的研究在近几年得到发展．Ahlers等在对 热对流广泛存在于天体、太阳、地球地幔、大 湍流热对流研究的综述 中指出，相对于实验研究， 气、海洋环流等 自然界运动，以及核反应堆堆芯冷 RB热对流 的DNS模拟研究存在 的主要困难是需 却系统、化工产品生产、电子元件设计等工程应 要巨大计算工作量 [4]．近几年，热对流的数值模拟 用领域 中．热对流研究对天体 内部运动规律、全 研究正在 向高Ra数的RB对流计算进发．宽高比 球气候变化、海洋环流、地磁变化、反应堆设计 为1／2的高Ra数计算结果较早获得 5[】．Shishkina 等 问题具有重要意义．作为一种典型的热对流运 和Wagner模拟计算热对流，使用了限体积法计算 动，Rayleigh-B6nard热对流 (RB对流)是一种 由底 部水平平板加热流体后所形成的位于两个平行平 了Ra≤2X100的情况，而对Ra=2×10加的情 况则 引入了大涡模型 6[，．Stevens等采用DNS计 板之间流体 自然对流过程．由于实验和分析的简便 性，使得RB对流成为热对流研究领域 中最为典型 算了三维圆柱体 内的RB热对流 其 中F=0．5和 的流动现象 [1]-近两年最新的实验一直努力增加 Pr=0．7，Ra数从2×106到2×1011，所用最大网 Ra数的范围，已做到了10 ≤Ra≤1015，并得到 格数为 (No×A ×人 )1081×351×1301[引． 了Nu(Ra1的标度关