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

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.