颗粒凝并动力学Monte Carlo方法的高效GPU并行计算.pdfVIP

2014 年 第 59 卷 第 14 期： 1358 ~ 1368 《中国科学》杂志社 论 文 SCIENCE CHINA PRESS 颗粒凝并动力学Monte Carlo 方法的高效GPU 并行计算 * 赵海波 , 徐祖伟, 刘昕, 史家伟, 郑楚光 华中科技大学煤燃烧国家重点实验室, 武汉 430074 * 联系人, E-mail: klinsmannzhb@163.com 2013-10-15 收稿, 2014-01-21 接受, 2014-04-01 网络版发表 国家自然科学基金、教育部新世纪优秀人才支持项目(NCET-09-0395)和多相复杂系统国家重点实验室开放课题 (MPCS-2011-D-02)资助 摘要 Monte Carlo(MC) 方法作为一种求解颗粒群平衡方程(PBE) 的有效方法(PBMC), 由于 关键词 它对多维问题的适应性、符合实际颗粒动力学特征的离散和随机本质、程序结构相对简单、 颗粒群平衡模拟 易于编程实现等优点受到人们持久、普遍的关注. 但在涉及到颗粒凝并问题时, 常规的 凝并 PBMC 方法计算代价较高, 与模拟颗粒数目的平方成正比, 限制了其工程应用. 并行计算技 随机模拟 并行计算 术的快速发展, 特别是近年来NVIDIA 公司提出的计算统一设备架构(CUDA)为PBMC 的快 CUDA 速高效模拟提供了一个良好的平台. 本文在CUDA 平台上实现了颗粒凝并动力学PBMC 的图 计算效率 形处理器(GPU)并行计算(分别实现了累计概率法和接受-拒绝法选择凝并对)及中央处理器 (CPU) 的协同处理, 与目前广泛运行于 CPU 的串行计算相比, 取得了精确的计算结果和非常 明显的加速, 计算代价仅与颗粒数目成正比, 在当前主流 GPU/CPU 设备上能够达到上百倍 的加速比. 颗粒群平衡模拟(PBM)研究的是离散单元(如 环芳香烃的生成长大; 燃烧源可吸入颗粒物的生 多相流中离散的固体颗粒、液滴或气泡等, 下文均 成、生长和迁徙) 、大气环境(如室内、室外气溶胶 用颗粒表示)在各种动力学事件(包括碰撞、凝并(或 的生成、团聚等微物理化学过程; 自然环境中 团聚、凝聚) 、冷凝/蒸发、成核、沉积等)主导下颗 雨、雾、冰雹等的形成和降落) 、化学工程(如流化 粒群的内部变量(通常为单组分颗粒的尺度, 也可 床、精馏塔等过程设备中团聚体和气泡的破碎和聚 包括颗粒中各组分质量、颗粒表面积、孔隙率、分 并; 乳化液滴颗粒的聚合和表面化学反应; 结晶过 形维数等) 的分布函数的时间和空间演变过程. 程中晶体的生成和生长; 纳米粉体制备过程中的 PBM 既可以作为一种单独的数值模拟手段来研究 快速化学反应、凝并、摩擦破碎和干燥; 胶体的团 离散系统的动力学演变过程, 也可以与多相流模 聚、微混合及结构重建现象; 异类液相萃取过程中 型耦合起来以考虑离散相-离散相、离散相-连续相 液滴颗粒的聚并、破碎、相间质量交换和相分离); 间的相互作用( 即四向耦合的多相湍流模型), 包括 生物医学工程(如细胞的生成、生长和分裂以及内 离散单元的碰撞、聚并、破碎以及相变等关键热点, 部蛋白质、DNA 和RNA 的再分布)等领域. 从而推动了多相流研究向更为科学的定量描述和 在 PBM 中, 离散系统被视为 1 个服从统计学 精确控制的转变.