A Fast Stratified Sampling Simulation of Coagulation Processes

We develop a new version of the direct simulation Monte Carlo method for coagulation processes governed by homogeneous Smoluchowsky equations. The method is based on a subdivision of the set of particle pairs into classes, and on an efficient algorithm for sampling from a discrete distribution, the so-called Walker's alias method. The efficiency of the new method is drastically increased compared to the conventional methods, especially when the coagulation kernel is strongly varying. The method is applied to solving a problem of islands formation on a surface due to a diffusion controlled coagulation

[1]  Karl Sabelfeld,et al.  Stochastic Lagrangian models and algorithms for spatially inhomogeneous Smoluchowski equation , 2003, Math. Comput. Simul..

[2]  Wolfgang Wagner,et al.  An Efficient Stochastic Algorithm for Studying Coagulation Dynamics and Gelation Phenomena , 2000, SIAM J. Sci. Comput..

[3]  Daniel T. Gillespie,et al.  The Stochastic Coalescence Model for Cloud Droplet Growth. , 1972 .

[4]  Karl K. Sabelfeld,et al.  Monte Carlo simulation of the coagulation processes governed by Smoluchowski equation with random coefficients , 1997, Monte Carlo Methods Appl..

[5]  CoalescenceDavid J. Aldous.  Stochastic Coalescence , 1998 .

[6]  Daniel T. Gillespie,et al.  An Exact Method for Numerically Simulating the Stochastic Coalescence Process in a Cloud , 1975 .

[7]  J. Marqusee.  Dynamics of late stage phase separations in two dimensions , 1984 .

[8]  F. Family,et al.  Rate-equation approach to island capture zones and size distributions in epitaxial growth. , 2001, Physical review letters.

[9]  Karl K. Sabelfeld,et al.  Stochastic algorithms for solving Smolouchovsky coagulation equation and applications to aerosol growth simulation. , 1996, Monte Carlo Methods Appl..

[10]  Flavius Guias.  A Monte Carlo approach to the Smoluchowski equations , 1997, Monte Carlo Methods Appl..

[11]  Markus Kraft,et al.  An efficient stochastic algorithm for simulating Nano-particle dynamics , 2002 .

[12]  Kurt Liffman,et al.  A direct simulation Monte-Carlo method for cluster coagulation , 1992 .

[13]  Robert F. Tichy,et al.  Quasi-Monte Carlo Methods for Numerical Integration: Comparison of Different Low Discrepancy Sequences , 1996, Monte Carlo Methods Appl..

[14]  Alastair J. Walker,et al.  An Efficient Method for Generating Discrete Random Variables with General Distributions , 1977, TOMS.