摘要：超细微颗粒物的尺度分布由通用动力学方程(GDE)描述,该文发展了求解同时考虑凝并和破碎的GDE的一种多重Monte Carlo(MMC)算法,该算法可以抛弃当前MonteCarlo算法中广泛采用的'子系统'的概念,引入了加权的'虚拟颗粒'的概念,且基于时间驱动,模拟过程中保持虚拟颗粒数目和计算区域体积不变.首先详细描述了MMC算法,包括时间步长的确定方法、同时处理凝并和破碎的方案、凝并和破碎事件发生的判断准则、凝并伙伴的确定、处理凝并和破碎事件的后果时保持虚拟颗粒数目恒定的方法等;然后利用MMC算法对几种特殊工况进行了数值模拟,模拟结果与理论值符合很好,验证了MMC算法的计算精度,且该算法具有较小的计算代价.

关键词：多重蒙特卡罗算法 凝并 破碎 通用动力学方程 超细微颗粒 

单位：华中科技大学; 煤燃烧国家重点实验室; 湖北省; 武汉市; 430074