计算物理 ›› 2013, Vol. 30 ›› Issue (1): 89-97.

• 论文 • 上一篇    下一篇

模拟颗粒凝并过程的快速Monte Carlo方法

徐祖伟, 赵海波, 刘昕, 郑楚光   

  1. 华中科技大学 煤燃烧国家重点实验室, 武汉 430074
  • 收稿日期:2012-04-13 修回日期:2012-08-06 出版日期:2013-01-25 发布日期:2013-01-25
  • 通讯作者: 赵海波(1977-),男,湖南宁乡,教授,博士,从事多相流系统的实验、模型及模拟研究,E-mail:klinsmannzhb@163.com
  • 作者简介:徐祖伟(1986-),男,河南,主要从事多相流热物理研究,E-mail:xu_zuwei@163.com
  • 基金资助:
    国家自然科学基金(50876037);教育部新世纪优秀人才支持计划(NCET-09-0395);多相复杂系统国家重点实验室开放课题(MPCS-2011-D-02)资助项目

Fast Monte Carlo Method for Particle Coagulation Dynamics

XU Zuwei, ZHAO Haibo, LIU Xin, ZHENG Chuguang   

  1. State Key Laboratory of Coal Combustion, Huazhong University of Science and Technology, Wuhan 430074, China
  • Received:2012-04-13 Revised:2012-08-06 Online:2013-01-25 Published:2013-01-25

摘要: 对常规异权值Monte Carlo(MC)方法进行改进,基于强核函数思想,通过颗粒群单重遍历即可求得强核函数最大值,采用接受-拒绝法随机搜寻凝并对,并利用搜寻过程中拒绝和接受的所有凝并对的信息来估计凝并事件的等待时间(时间步长),从而避免颗粒群的双重遍历,以提高MC的效率.对典型工况的模拟结果显示该快速方法计算代价仅为O(Ns),能够显著提高计算效率,同时保持足够的计算精度,较好地协调计算代价与计算精度之间的矛盾.

关键词: 颗粒群平衡模拟, 凝并动力学, Monte Carlo方法, 异数目权值, 粒颗粒尺度分布

Abstract: We propose a fast random simulation strategy based on differentially weighted MC. The strategy improves computation efficiency significantly, and guarantees enough calculation accuracy, thus coordinates contradiction between computation cost and computation accuracy. The main idea is based on majorant kernel. It is possible to transfer a traditional coagulation kernel to a majorant kernel through splitting and amplifying slightly. The maximum of majornant kernel is obtained by single looping over all simulation particles. The maximum majornant kernel is used to approximate the maximum coagulation kernel in particle population, and is further used to search coagulation particle pairs randomly with acceptance-rejection method. The waiting time (time-step) for a coagulation event is calculated by summing coagulation kerenls of particle pairs involved in acceptance/rejection processes.. Double looping in normal Monte Carlo simulation is avoided and computation efficiency is improved greatly.

Key words: population balance modeling, coagulation, Monte Carlo method, differentially weighted simulation particles, particle size distribution

中图分类号: