摘要: 对常规异权值Monte Carlo(MC)方法进行改进,基于强核函数思想,通过颗粒群单重遍历即可求得强核函数最大值,采用接受-拒绝法随机搜寻凝并对,并利用搜寻过程中拒绝和接受的所有凝并对的信息来估计凝并事件的等待时间(时间步长),从而避免颗粒群的双重遍历,以提高MC的效率.对典型工况的模拟结果显示该快速方法计算代价仅为O(Ns),能够显著提高计算效率,同时保持足够的计算精度,较好地协调计算代价与计算精度之间的矛盾.
中图分类号:
徐祖伟, 赵海波, 刘昕, 郑楚光. 模拟颗粒凝并过程的快速Monte Carlo方法[J]. 计算物理, 2013, 30(1): 89-97.
XU Zuwei, ZHAO Haibo, LIU Xin, ZHENG Chuguang. Fast Monte Carlo Method for Particle Coagulation Dynamics[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2013, 30(1): 89-97.