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
Additional Information
- Christian Lécot
- Affiliation: Laboratoire de Mathématiques, Université de Savoie, Campus scientifique, 73376 Le Bourget-du-Lac cedex, France
- Email: Christian.Lecot@univ-savoie.fr
- Wolfgang Wagner
- Affiliation: Weierstraß-Institut für Angewandte Analysis und Stochastik, Mohrenstraße 39, 10117 Berlin, Germany
- Email: wagner@wias-berlin.de
- Received by editor(s): November 11, 2002
- Received by editor(s) in revised form: March 14, 2003
- Published electronically: January 5, 2004
- Additional Notes: Computation was supported by the Centre Grenoblois de Calcul Vectoriel du Commissariat à l’Énergie Atomique, France
- © Copyright 2004 American Mathematical Society
- Journal: Math. Comp. 73 (2004), 1953-1966
- MSC (2000): Primary 65C05; Secondary 70-08, 82C80
- DOI: https://doi.org/10.1090/S0025-5718-04-01627-8
- MathSciNet review: 2059745