Diffusion in a nonhomogeneous medium: quasi-random walk on a lattice

Abstract : We are interested in Monte Carlo (MC) methods for solving the diffusion equation: in the case of a constant diffusion coefficient, the solution is approximated by using particles and in every time step, a constant stepsize is added to or substracted from the coordinates of each particle with equal probability. For a spatially dependent diffusion coefficient, the naive extension of the previous method using a spatially variable stepsize introduces a systematic error: particles migrate in the directions of decreasing diffusivity. A correction of stepsizes and stepping probabilities has recently been proposed and the numerical tests have given satisfactory results. In this paper, we describe a quasi-Monte Carlo (QMC) method for solving the diffusion equation in a spatially nonhomogeneous medium: we replace the random samples in the corrected MC scheme by low-discrepancy point sets. In order to make a proper use of the better uniformity of these point sets, the particles are reordered according to their successive coordinates at each time step. We illustrate the method with numerical examples: in dimensions 1 and 2, we show that the QMC approach leads to improved accuracy when compared with the original MC method using the same number of particles.
Document type :
Journal articles
Liste complète des métadonnées

Cited literature [18 references]  Display  Hide  Download

Contributor : Christian Lecot <>
Submitted on : Wednesday, February 19, 2014 - 10:19:32 PM
Last modification on : Wednesday, April 17, 2019 - 2:30:21 PM
Document(s) archivé(s) le : Monday, May 19, 2014 - 2:20:23 PM


Files produced by the author(s)


  • HAL Id : hal-00949605, version 1



Rami El Haddad, Christian Lécot, Gopalakrishnan Venkiteswaran. Diffusion in a nonhomogeneous medium: quasi-random walk on a lattice. Monte Carlo Methods and Applications, De Gruyter, 2010, 16, pp.211-230. ⟨hal-00949605⟩



Record views


Files downloads