求解考虑颗粒凝并的通用动力学方程的多重Monte Carlo算法

摘要

MonteCarlo(MC)方法被广泛用于通用动力学方程的求解,然而普通MC方法的计算代价较高而计算精度不稳定·　提出一种新的多重MonteCarlo(MMC)算法来求解GDE,该算法同时具有基于时间驱动MC方法、常数目法和常体积法的特点·　首先详细介绍了该算法,包括加权虚拟颗粒的引入,MMC算法的计算流程,时间步长的设置,颗粒是否发生凝并事件的判断,凝并伙伴的寻找,凝并事件的后果处理·　然后利用MMC算法对存在理论分析解的5种特殊工况进行数值求解,模拟结果与理论解符合很好,证明MMC算法具有良好的计算精度和较低的计算代价·　最后分析了不同类型的凝并核对于凝并过程的影响,常凝并核和连续区布朗凝并核对小颗粒影响大一些。