气体动力学直接模拟MonteCarlo的高效GPU并行计算-计算物理.PDFVIP

第３２卷 第２期 计　 　 算　 　 物　 　 理 Ｖｏｌ．３２，Ｎｏ．２ ２０１５年３月 　 ＣＨＩＮＥＳＥ ＪＯＵＲＮＡＬ ＯＦ ＣＯＭＰＵＴＡＴＩＯＮＡＬ ＰＨＹＳＩＣＳ　 Ｍａｒ． ，２０１５ 文章编号：１００１⁃２４６Ｘ（２０１５）０２⁃０１６９⁃０８ 气体动力学直接模拟Ｍｏｎｔｅ Ｃａｒｌｏ 的 高效 ＧＰＵ并行计算 贺永翔，　 刘　 昕，　 赵海波∗ （华中科技大学煤燃烧国家重点实验室，武汉　 ４３００７４） 摘　 要：实现了基于计算统一设备架构（ＣＵＤＡ）的直接模拟Ｍｏｎｔｅ Ｃａｒｌｏ（ＤＳＭＣ）并行算法，改进了原有多图形处理 器（ＧＰＵ）数据之间传输并行算法，数值模拟计算二维Ｃｏｕｅｔｔｅ流和二维顶盖驱动方腔流，定量比较了ＣＰＵ、单 ＧＰＵ 和多ＧＰＵ并行计算的结果和计算时间．结果表明单ＧＰＵ并行计算相对ＣＰＵ计算的加速效果可以达到１０ ～３０倍， 双ＧＰＵ并行计算加速效果可以达到４０～６０倍，多ＧＰＵ并行计算的加速效率接近１００％，且计算精度能够得到良好 保证． 关键词：ＤＳＭＣ；ＣＵＤＡ；图形处理器；计算效率 中图分类号：Ｏ３５６ 文献标志码：Ａ ０　 引言 克努森数（Ｋｎ＝λ／ Ｌ）是描述流体流动从连续到稀薄介质过渡的一个无量纲参数，其大小为气体的平均 自由程与系统的特征长度的比值．当Ｋｎ数大于０１时，连续介质的质量、动量和能量方程中剪切应力和热量 ［１］ 不能由低阶的宏观量来表征，基于连续性假设的连续性方程和Ｎａｖｉｅｒ⁃Ｓｔｏｋｅｓ方程失效 ，只能通过求解 Ｂｏｌｔｚｍａｎｎ方程来得到气体速度分布函数．Ｂｏｌｔｚｍａｎｎ方程是一个多维多变量积分微分方程，给理论和数值求 解带来了很大的困难，目前只有在少数理想情况下才得到其解析解．针对这种情况，Ｂｉｒｄ 提出了直接模拟 ［２］ ［３－４］ Ｍｏｎｔｅ Ｃａｒｌｏ（ＤＳＭＣ）方法 ，并已被证明收敛于Ｂｏｌｔｚｍａｎｎ方程 ．ＤＳＭＣ方法具有更容易实现复杂而真实 的分子模型（比如ＶＨＳ、ＶＳＳ、ＧＨＳ等）、考虑化学反应和辐射的气体流动、模拟三维复杂流场等优势． 目前该 方法在超音速绕流飞行器，化学合成薄膜，微尺度／ 纳米尺度的流动中有着广泛的应用． ＤＳＭＣ方法源于分子动力学．分子动力学方法的基本思想是用大量的模拟分子运动代替真实流体运动， 且认为当两个运动分子间的距离小于一定值时将发生碰撞，碰撞后分子运动速度按经典力学的规律计算得 到．该方法通过模拟分子的运动轨迹来判断分子间是否发生碰撞，模拟所需的计算机时正比于模拟分子数Ｎ 的平方，因此通常分子动力学方法所追踪的分子数不能太大．为了克服这个缺陷，Ｂｉｒｄ提出了几率论方法判 断分子间的碰撞．认为在充分短的时间间隔内，模拟分子的运动与碰撞是解耦的．即在单个计算时间步长内， 模拟分子作匀速直线运动，之后从几率的角度出发处理模拟分子之间、以及模拟分子与壁面的碰撞问题． ＤＳＭＣ方法是通过大量模拟分子的运动与碰撞实现对真实气体流动的模拟，每个模拟分子代表相当数目的 真实分子．该方法的重要思想是碰撞与解耦，正确的碰撞取样方法是模拟能否还原真实流动过程的关键． ［５］ Ｂｉｒｄ提出了“时间计数器”（ｔｉｍｅ⁃ｃｏｕｎｔｅｒ，ＴＣ）方法 ，其本质是一种事件驱动的方法，主要思想是通过接受拒 绝法选择成功的碰撞对，根据成功选择碰撞对的核函数来估计碰撞等待时间，直到累积时间超过时间步长 Δｔ．但当网格中分子数不够大时，偶然中选的分子碰撞对核函数很小使碰撞频率畸变，导致误差．此外，沈青 ［６］ ［７］ 则建议了随机取样频率（ｒａｎｄｏｍｌｙ ｓａｍｐｌ