Diffusion in a nonhomogeneous medium: quasi-random walk on a lattice - Université Grenoble Alpes Accéder directement au contenu
Article Dans Une Revue Monte Carlo Methods and Applications Année : 2010

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


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.
Fichier principal
Vignette du fichier
ElHaddadMCM09.pdf (886.96 Ko) Télécharger le fichier
Origine Fichiers produits par l'(les) auteur(s)

Dates et versions

hal-00949605 , version 1 (19-02-2014)


  • 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, 2010, 16, pp.211-230. ⟨hal-00949605⟩
172 Consultations
390 Téléchargements


Gmail Mastodon Facebook X LinkedIn More