Monte Carlo simulations of two-component drop growth by stochastic coalescence 1Universidad Autónoma de la Ciudad de México, México City, México
2Centro de Ciencias de la Atmósfera, Universidad Nacional Autónoma de México, México City, México
Received: 18 Jan 2008 – Published in Atmos. Chem. Phys. Discuss.: 16 Apr 2008 – Published: 18 Feb 2009Abstract. The evolution of two-dimensional drop distributions is simulated in this
study using a Monte Carlo method. The stochastic algorithm of Gillespie (1976) for chemical reactions
in the formulation proposed by Laurenzi et al. (2002) was used to simulate the kinetic behavior of the drop population.
Within this framework, species are defined as droplets of specific size and
aerosol composition. The performance of the algorithm was checked by a
comparison with the analytical solutions found by Lushnikov (1975) and
Golovin (1963) and with finite difference solutions of the two-component
kinetic collection equation obtained for the Golovin (sum) and hydrodynamic
kernels. Very good agreement was observed between the Monte Carlo
simulations and the analytical and numerical solutions. A simulation for
realistic initial conditions is presented for the hydrodynamic kernel. As
expected, the aerosol mass is shifted from small to large particles due to
collection process. This algorithm could be extended to incorporate various
properties of clouds such several crystals habits, different types of
soluble CCN, particle charging and drop breakup.
Citation: Alfonso, L., Raga, G. B., and Baumgardner, D.: Monte Carlo simulations of two-component drop growth by stochastic coalescence, Atmos. Chem. Phys., 9, 1241-1251, doi:10.5194/acp-9-1241-2009, 2009.