Abstract
Monte-Carlo (MC) method is widely adopted to take into account general dynamic equation (GDE) for particle coagulation, however popular MC method has high computation cost and statistical fatigue. A new Multi-Monte-Carlo (MMC) method, which has characteristics of time-driven MC method, constant number method and constant volume method, was promoted to solve GDE for coagulation. Firstly MMC method was described in details, including the introduction of weighted fictitious particle, the scheme of MMC method, the setting of time step, the judgment of the occurrence of coagulation event, the choice of coagulation partner and the consequential treatment of coagulation event. Secondly MMC method was validated by five special coagulation cases in which analytical solutions exist. The good agreement between the simulation results of MMC method and analytical solutions shows MMC method conserves high computation precision and has low computation cost. Lastly the different influence of different kinds of coagulation kernel on the process of coagulation was analyzed: constant coagulation kernel and Brownian coagulation kernel in continuum regime affect small particles much more than linear and quadratic coagulation kernel, whereas affect big particles much less than linear and quadratic coagulation kernel.
Similar content being viewed by others
References
Tucker W G. An overview of PM 2.5 sources and control strategies[J].Fuel Processing Technology, 2000,65(1): 379–392.
Meng Z, Dabdub D, Seinfeld J H. Size-resolved and chemically resolyed model of atmospheric aerosol dynamics[J].Journal of Geophysical Research, 1998,103 (3): 3419–3435.
Debry E, Sportisse B, Jourdain B. A stochastic approach for the numerical simulation of the general dynamics equation for aerosols[J].Journal of Computational Physics, 2003,184 (2): 649–669.
Liffman K. A direct simulation Monte Carlo method for cluster coagulation[J].Journal of Computational Physics, 1992,100 (1): 116–127.
Kruis F E, Maisels A, Fissan H. Direct simulation Monte Carlo method for particle coagulation and aggregation[J].AICHE Journal, 2000,46 (9): 1735–1742.
Smith M, Matsoukas T. Constant-number Monte Carlo simulation of population balances[J].Chemical Engineering Science, 1998,53 (9): 1777–1786.
Lee K, Matsoukas T. Simultaneous coagulation and break-up using constant-N Monte Carlo [J].Powder Technology, 2000,110 (1/2): 82–89.
Lin Y, Lee K, Matsoukas T. Solution of the population balance equation using constantnumber Monte Carlo[J].Chemical Engineering Science, 2002,57 (12): 2241–2252.
Nanbu K. Direct simulation scheme derived from the Boltzmann equation I—Monocomponent gases[J].Journal of the Physical Society of Japan, 1980,49 (11): 2042–2049.
Nanbu K, Yonemura S. Weighted particles in coulomb collision simulations based on the theory of a cumulative scattering angle[J].Journal of Computational Physics, 1998,145 (2): 639–654.
Zhao Haibo, Zheng Chuguang, Xu Minghou. Multi-Monte Carlo method for particle coagulation: description and validation[J].Applied Mathematics and Computation, 2005,168 (2).
Friedlander S K, Wang C S. The self-preserving particle size distribution for coagulation by Brownian motion[J].Journal of Colloid and Interface Science, 1966,22 (1): 126–132.
Lee K W, Lee Y J, Han D S. The log-normal size distribution theory for Brownian coagulation in the low Knudsen number regime[J].Journal of Colloid and Interface Science, 1997,188 (2): 486–492.
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by LI Jia-chun
Project supported by the National Key Basic Research and Development Program of China (No. 2002CB211602); the National Natural Science Foundation of China (No. 90410017)
Rights and permissions
About this article
Cite this article
Hai-bo, Z., Chu-guang, Z. & Ming-hou, X. Multi-monte-carlo method for general dynamic equation considering particle coagulation. Appl Math Mech 26, 953–962 (2005). https://doi.org/10.1007/BF02464246
Received:
Revised:
Issue Date:
DOI: https://doi.org/10.1007/BF02464246
Key words
- numerical solution
- population balance equation
- fictitious particle
- particle size distribution
- constant kernel
- computation cost
- computation precision