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. © Springer-Verlag Berlin Heidelberg 2009.
CITATION STYLE
El Haddad, R., Lécot, C., & Venkiteswaran, G. (2009). Quasi-Monte Carlo simulation of diffusion in a spatially nonhomogeneous medium. In Monte Carlo and Quasi-Monte Carlo Methods 2008 (pp. 339–354). Springer Verlag. https://doi.org/10.1007/978-3-642-04107-5_21
Mendeley helps you to discover research relevant for your work.