Quasi-Monte Carlo simulation of diffusion in a spatially nonhomogeneous medium - Archive ouverte HAL Accéder directement au contenu
Communication Dans Un Congrès Année : 2009

Quasi-Monte Carlo simulation of diffusion in a spatially nonhomogeneous medium

Résumé

We propose and test a quasi-Monte Carlo (QMC) method for solving the diffusion equation in the spatially nonhomogeneous case. For a constant diffusion coefficient, the Monte Carlo (MC) method is a valuable tool for simulating the equation: the solution is approximated by using particles and in every time step the displacement of each particle is drawn from a Gaussian distribution with constant variance. But for a spatially dependent diffusion coefficient, the straightforward extension using a spatially variable variance leads to biased results. A correction to the Gaussian steplength was recently proposed and provides satisfactory results. In the present work, we devise a QMC variant of this corrected MC scheme. We present the results of some numerical experiments showing that our QMC algorithm converges better than the corresponding MC method for the same number of particles.
Fichier non déposé

Dates et versions

hal-00392298 , version 1 (06-06-2009)

Identifiants

  • HAL Id : hal-00392298 , version 1

Citer

Rami El Haddad, Christian Lécot, Gopalakrishnan Venkiteswaran. Quasi-Monte Carlo simulation of diffusion in a spatially nonhomogeneous medium. Eighth international Conference on Monte Carlo and Quasi-Monte Carlo Methods in Scientific Computing, Jul 2008, Montréal, Canada. pp.339-354. ⟨hal-00392298⟩
110 Consultations
0 Téléchargements

Partager

Gmail Facebook X LinkedIn More