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

Abstract : 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.
Document type :
Conference papers
P. L'Ecuyer, A. Owen. Eighth international Conference on Monte Carlo and Quasi-Monte Carlo Methods in Scientific Computing, Jul 2008, Montréal, Canada. Springer, pp.339-354, 2009
Liste complète des métadonnées

https://hal.archives-ouvertes.fr/hal-00392298
Contributor : Christian Lecot <>
Submitted on : Saturday, June 6, 2009 - 5:59:22 PM
Last modification on : Thursday, January 11, 2018 - 6:12:26 AM

Identifiers

  • HAL Id : hal-00392298, version 1

Collections

Citation

Rami El Haddad, Christian Lécot, Gopalakrishnan Venkiteswaran. Quasi-Monte Carlo simulation of diffusion in a spatially nonhomogeneous medium. P. L'Ecuyer, A. Owen. Eighth international Conference on Monte Carlo and Quasi-Monte Carlo Methods in Scientific Computing, Jul 2008, Montréal, Canada. Springer, pp.339-354, 2009. 〈hal-00392298〉

Share

Metrics

Record views

160