多松弛时间格子BoItzmann方法在GPU上的实现.pdfVIP

第 28卷 第 3期 计算机与应用化学 、，01．28．NO．3 2011年 3月 28f1 ComputersandAppliedChemistry March 2011 多松弛时间格子 BoItzmann方法在GPU上的实现 张云 一一，Y-4,伟 ，葛蔚 ，杨朝合 (1．中国石油大学重质油加]：国家重点实验室，山东，青岛，266555： 2．中国科学院过程工程研究所多相复杂系统 国家重点实验室，北京，100190) 摘要：近年来，随着统一计算设备构架 (CUDA)的出现，高端 图形处理器 (GPU)在 图像处理、计算流体力学等科学计算领 域的应用得到了快速发展 。属于介观数值方法的格子Boltzmann方法 (LBM)是 1种新的计算流体力学 (CFD)方法，具有算 法简单 、能处理复杂边界条件 、压力能够直接求解等优势 ，在 多相流 、湍流、渗流等领域得到了广泛应用。LBM 由于具有 内在 的并行性，特别适合在 GPU上计算 。采用多松弛时问模型 (MRT)的LBM，受松弛因子的影响较小并且数值稳定性较好 。本 文实现了MRT-LBM在基于CUDA 的GPU上的计算，并通过计算流体力学经典算例——二维方腔流来验证计算的正确性。在 雷诺数 Re=[10，10]之间，计算了多达 26种雷诺数的算例，并将 Re=l0，4×10，10，2x10，5x10，7．5x10算例对应的主涡中 心坐标与文献 中结果进行 了对 比。计算结果与文献数值实验符合较好，从而验证了算法实现 的正确忡 ，并显示出MRT-LBM 具 有更优的数值稳定性。奉文还分析 了在 GPU上 MRT-LBM 的计算性 能并与CPU的计算进行了比较，结果表明，GPU可以极大 地加快 MRT-LBM 的计算，NVIDIATeslaC2050相对于单核 IntelXeon5430CPU的加速 比约为60侪 。 关键词：格子 Boltzmann方法；多松弛时间模型；方腔流；GPU 中图分类号 ：TQ015．9；TP391．9；06-39 文献标识码：A 文章编号：1001—4160(2011)03—265—269 1 引言 Boltzmann方法(MI己T-LBM)在GPU上计算的实现方法。 日前 ，高端的图形处理器 (graphicsprocessingunit， 2 算法实现 GPU)计算性能可以达到每秒万亿次级别，这相当于 1 2．1 MRT-LBM 的理论基础 个高性能计算集群系统，远远大于主流CPU的计算能力， 在这里我们采用二维D2Q9模型f6进行说 明，格子 已被用于图形处理以外的通用计算领域 ，如数据处理， 科学计算等。自从NVIDIA公司于 2007年推出针对通用 Boltzmann方程可以写为， 计算 GPU (generalpurposeGPU， GPGPU)的编程平台 Z(x+eidt，t+6t)一 (x，t)=Q ， (1) 其 中， 是t时刻在位置 x处速度为e的分布 函数 ， CUDA (ComputeUnifiedDeviceArchitecture) 以来 ， i=0，1…，8，具体分布如图 1所示，Q，为碰撞函数。在 GPGPU得到了快速发展_lJ。 MRT-LBM 中碰撞函数为加l， 格子 Boltzmann方法 (1atticeBoltzmannmethod， LBM)直接从离散模型出发，应用质量、动量和能量守 Q=一M