Constant-number Monte Carlo simulation of population balances

A Monte Carlo method for the simulation of growth processes is presented, in which the number of particles is kept constant, regardless of whether the actual process results in a net loss (as in coagulation) or net increase (as in fragmentation) of particles. General expressions are derived for the inter-event time, number concentration, and expected average size as a function of time. The method is applied to two coagulation models, constant kernel and Brownian coagulation, and it is shown to provide an accurate description of the dynamical evolution of systems undergoing coagulation. This algorithm allows for indefinitely long simulations. The error scales as the inverse square root of the number of particles in the simulation and grows logarithmically in time.

[1]  Y. Tang,et al.  Monte Carlo Simulation of Agglomeration and Grinding , 1997 .

[2]  D. Ramkrishna.  The Status of Population Balances , 1985 .

[3]  Themis Matsoukas,et al.  Dynamics of aerosol agglomerate formation , 1991 .

[4]  Margaritis Kostoglou,et al.  Evaluation of Zero Order Methods for Simulating Particle Coagulation , 1994 .

[5]  P. A. Daniels,et al.  The mass distribution of an accreting system of dust particles – a Monte Carlo simulation , 1981 .

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

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

[8]  S. Friedlander,et al.  The self-preserving particle size distribution for coagulation by brownian motion☆ , 1966 .

[9]  F. Leyvraz.  Rate Equation Approach to Aggregation Phenomena , 1986 .

[10]  Doraiswami Ramkrishna,et al.  Analysis of population balance—IV: The precise connection between Monte Carlo simulation and population balances , 1981 .

[11]  Octave Levenspiel,et al.  A Monte Carlo treatment for reacting and coalescing dispersed phase systems , 1965 .

[12]  Robert M. Ziff,et al.  Coagulation equations with gelation , 1983 .

[13]  Doraiswami Ramkrishna,et al.  Simulation of particulate systems using the concept of the interval of quiescence , 1977 .

[14]  Monte Carlo results for random coagulation , 1985 .

[15]  P. Meakin.  The Growth of Fractal Aggregates , 1987 .

[16]  Sotiris E. Pratsinis,et al.  Time-Lag for Attainment of the Self-Preserving Particle Size Distribution by Coagulation , 1994 .

[17]  J. P. Gooch,et al.  Monte Carlo simulation of size-enlargement mechanisms in crystallization , 1996 .

[18]  W. T. Scott.  Poisson Statistics in Distributions of Coalescing Droplets , 1967 .

[19]  D'arcy W. Thompson.  On growth and form i , 1943 .

[20]  W. H. Weinberg,et al.  Theoretical foundations of dynamical Monte Carlo simulations , 1991 .

[21]  P. Das.  Monte Carlo simulation of drop breakage on the basis of drop volume , 1996 .