CHINESE JOURNAL OF COMPUTATIONAL PHYSICS ›› 2013, Vol. 30 ›› Issue (1): 89-97.

Previous Articles     Next Articles

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

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

CLC Number: