方腔自然对流换热问题的LBM并行计算.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：083440 方腔自然对流换热问题的 LBM 并行计算 赵婷，杨茉*，李佳明，卢玫 （上海理工大学动力工程学院，上海，200093） (Tel: 021 Email：zhaoting_842@163.com) 摘要：本文利用消息传递接口库实施数据传递，将 SPMD 编程模式的并行算法与 LBM 相结合，提出 了LBM 的两种并行计算方法，在小型局域网下实现了对二维方腔自然对流换热问题的LBM 并行模拟。 数值结果表明，所建立的并行计算环境和区域分解并行算法能够得到正确和收敛的数值结果。对并行 算法作出的性能分析表明，所给出的区域分解并行算法得到了明显的加速效率。 关键词：并行计算，LBM，区域分解，方腔自然对流 0 前言 格子玻尔兹曼方法(Lattice-Boltzmann method, LBM) 由格子气自动机（Lattice Gas Automata ，LGA ）发展而来，从微观层次描述流体内部流动，是近年来涌现出来的一种 [1] 新型的能够与现代高性能计算机紧密耦合的本质并行算法 ，在流体动力学及其相关领 [2],[3] [4] 域得到了成功的应用 。通过多尺度展开能够从LBM 的演化方程得到宏观方程 。同 时，LBM具有算法简单，空间、时间完全离散及粒子相互作用的局域性等优点使其适合 分布式并行。消息传递接口（MPI ）是国内外在高性能计算中最广泛使用的并行编程环 境，已成为一种并行程序的标准[5] 。而消息传递并行性的开发中基于区域分解形式的单 控制流多数据流（Simple Program Multiple Date，SPMD）编程是最普遍且并行效率最高 的形式。 本文将 SPMD 的编程模式和LBM 相结合，通过 MPI 在 Windows XP PC 机组成的 小型局域网环境下搭建了并行计算平台，研究基于 LBM 的高效计算机并行算法，给出 了 LBM 区域分解并行方法实施的具体步骤以及实现数据交换的措施，并通过对二维方 腔自然对流换热问题的并行模拟实验比较，考察和验证所给出的并行算法的并行性能和 可靠性。 1计算模型及LBM模型 1.1 物理模型 图1示出了本文要计算的二维方腔内自然对流换热问题物理模型。方腔上下底面绝 热，左边为高温壁面T ，右边为低温壁面T ，方腔高度为H ，方腔内充满空气，Pr =0.7， h c β ⋅ ΔT ⋅ g ⋅ H 3 Ra 。 ν D 资助项目：国家自然科学基金项目(No.；上海市重点学科建设项目资助(No. J50501) *通讯作者：上海军工路 516 号，200093，Tel: 021Email: yangm66@ g 图 1 计算模型 1.2 数学模型 LBM计算选用模型为GUO等[6]提出的TD2G9模型，即二维空间中分布函数和速度等 分解为图2所示的9个方向，c 表达式为： i