A quasi–Monte Carlo scheme for Smoluchowski’s coagulation equation

Lécot, Christian; Wagner, Wolfgang

doi:10.1090/S0025-5718-04-01627-8

A quasi–Monte Carlo scheme for Smoluchowski’s coagulation equation
HTML articles powered by AMS MathViewer

by Christian Lécot and Wolfgang Wagner PDF

Math. Comp. 73 (2004), 1953-1966 Request permission

Abstract:

This paper analyzes a Monte Carlo algorithm for solving Smoluchowski’s coagulation equation. A finite number of particles approximates the initial mass distribution. Time is discretized and random numbers are used to move the particles according to the coagulation dynamics. Convergence is proved when quasi-random numbers are utilized and if the particles are relabeled according to mass in every time step. The results of some numerical experiments show that the error of the new algorithm is smaller than the error of a standard Monte Carlo algorithm using pseudo-random numbers without reordering the particles.

References

David J. Aldous, Deterministic and stochastic models for coalescence (aggregation and coagulation): a review of the mean-field theory for probabilists, Bernoulli 5 (1999), no. 1, 3–48. MR 1673235, DOI 10.2307/3318611
Hans Babovsky, On a Monte Carlo scheme for Smoluchowski’s coagulation equation, Monte Carlo Methods Appl. 5 (1999), no. 1, 1–18. MR 1684990, DOI 10.1515/mcma.1999.5.1.1
Ye. R. Domilovskiy, A. A. Lushnikov, and V. N. Piskurov, Monte Carlo simulation of coagulation processes; Russian transl. in Izv. Acad. Sci. USSR Atmospher. Ocean. Phys. 15 (1979), no. 2. MR 580838
Andreas Eibeck and Wolfgang Wagner, An efficient stochastic algorithm for studying coagulation dynamics and gelation phenomena, SIAM J. Sci. Comput. 22 (2000), no. 3, 802–821. MR 1784826, DOI 10.1137/S1064827599353488
K.-T. Fang, F.J. Hickernell, and H. Niederreiter (Eds.), Monte Carlo and Quasi–Monte Carlo Methods 2000, Springer, Berlin, 2002.
Henri Faure, Discrépance de suites associées à un système de numération (en dimension $s$), Acta Arith. 41 (1982), no. 4, 337–351 (French). MR 677547, DOI 10.4064/aa-41-4-337-351
A.L. Garcia, C. Van den Broeck, M. Aertsens, and R. Serneels, A Monte Carlo simulation of coagulation, Phys. A 143 (1987), 535–546.
D.N. Gillespie, An exact method for numerically simulating the stochastic coalescence process in a cloud, J. Atmospheric Sci. 32 (1975), 1977–1989.
Ole J. Heilmann, Analytical solutions of Smoluchowski’s coagulation equation, J. Phys. A 25 (1992), no. 13, 3763–3771. MR 1172075, DOI 10.1088/0305-4470/25/13/024
Christian Lécot, Error bounds for quasi-Monte Carlo integration with nets, Math. Comp. 65 (1996), no. 213, 179–187. MR 1325870, DOI 10.1090/S0025-5718-96-00690-4
K. Liffman, A direct simulation Monte Carlo method for cluster coagulation, J. Comput. Phys. 100 (1992), 116–127.
Harald Niederreiter, Point sets and sequences with small discrepancy, Monatsh. Math. 104 (1987), no. 4, 273–337. MR 918037, DOI 10.1007/BF01294651
Harald Niederreiter, Random number generation and quasi-Monte Carlo methods, CBMS-NSF Regional Conference Series in Applied Mathematics, vol. 63, Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 1992. MR 1172997, DOI 10.1137/1.9781611970081
K. K. Sabelfeld, S. V. Rogasinsky, A. A. Kolodko, and A. I. Levykin, Stochastic algorithms for solving Smolouchovsky coagulation equation and applications to aerosol growth simulation, Monte Carlo Methods Appl. 2 (1996), no. 1, 41–87. MR 1395316, DOI 10.1515/mcma.1996.2.1.41
M. v. Smoluchowski, Versuch einer mathematischen Theorie der Koagulationskinetik kolloider Lösungen, Z. Phys. Chem. 92 (1916), 129–168.
M. v. Smoluchowski, Drei Vorträge über Diffusion, Brownsche Molekularbewegung und Koagulation von Kolloidteilchen, Phys. Z. 17 (1916), 557–571, 585–599.
J.L. Spouge, Monte Carlo results for random coagulation, J. Colloid Interface Sci. 107 (1985), 38–43.
S. C. Zaremba, Some applications of multidimensional integration by parts, Ann. Polon. Math. 21 (1968), 85–96. MR 235731, DOI 10.4064/ap-21-1-85-96