Simulation of diffusion using quasi-random walk methods

Ibrahim Coulibaly and Christian Lécot

Mathematics and Computers in Simulation (MATCOM), 1998, vol. 47, issue 2, 153-163

Abstract: We present a method for computer simulation of diffusion. The method uses quasi-random walk of particles. We consider a pure initial value problem for a simple diffusion equation in s space dimensions. We introduce s spatial steps Δxi. A semi-discrete approximation to the equation is obtained by replacing the spatial derivatives with finite differences. N particles are sampled from the initial distribution. The time interval is partitioned into subintervals of length Δt. The discretization in time is obtained by resorting to the forward Euler method. In every time step the particle movement is regarded as an approximate integration is s+1 dimensions. A quasi-Monte Carlo estimate for the integral is obtained by using a (0, s+1)-sequence. A key element in successfully applying the low discrepancy sequence is a technique involving renumbering the particles at each time step. We prove that the computed solution converges to the solution of the semi-discrete equation as N→∞ and Δt→0. We present numerical tests which show that random walk results are improved with quasi random sequences and renumbering.

Keywords: Diffusion; Monte Carlo methods; Random walk; Discrepancy; Error bounds (search for similar items in EconPapers)
Date: 1998
References: View complete reference list from CitEc
Citations: View citations in EconPapers (2)

Downloads: (external link)
http://www.sciencedirect.com/science/article/pii/S0378475498001001
Full text for ScienceDirect subscribers only

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:eee:matcom:v:47:y:1998:i:2:p:153-163

DOI: 10.1016/S0378-4754(98)00100-1

Access Statistics for this article

Mathematics and Computers in Simulation (MATCOM) is currently edited by Robert Beauwens

More articles in Mathematics and Computers in Simulation (MATCOM) from Elsevier
Bibliographic data for series maintained by Catherine Liu ().