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

作者：徐祖伟; 赵海波; 刘昕; 郑楚光

来源：计算物理, 2013, 30(01): 89-97.

DOI：10.19596/j.cnki.1001-246x.2013.01.012

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

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

• 单位
  华中科技大学; 煤燃烧国家重点实验室