IDEAS home Printed from https://ideas.repec.org/a/eee/matcom/v47y1998i2p153-163.html
   My bibliography  Save this article

Simulation of diffusion using quasi-random walk methods

Author

Listed:
  • Coulibaly, Ibrahim
  • Lécot, Christian

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.

Suggested Citation

  • Coulibaly, Ibrahim & Lécot, Christian, 1998. "Simulation of diffusion using quasi-random walk methods," Mathematics and Computers in Simulation (MATCOM), Elsevier, vol. 47(2), pages 153-163.
    • Handle: RePEc:eee:matcom:v:47:y:1998:i:2:p:153-163
    DOI: 10.1016/S0378-4754(98)00100-1
    as

    Download full text from publisher

    File URL: http://www.sciencedirect.com/science/article/pii/S0378475498001001
    Download Restriction: Full text for ScienceDirect subscribers only
    File URL: https://libkey.io/10.1016/S0378-4754(98)00100-1?utm_source=ideas
    LibKey link: if access is restricted and if your library uses this service, LibKey will redirect you to where you can use your library subscription to access this item
    ---><---

    As the access to this document is restricted, you may want to search for a different version of it.

    Citations

    Citations are extracted by the CitEc Project, subscribe to its RSS feed for this item.
    as


    Cited by:

    1. Lécot, Christian & Schmid, Wolfgang Ch., 2001. "Particle approximation of convection–diffusion equations," Mathematics and Computers in Simulation (MATCOM), Elsevier, vol. 55(1), pages 123-130.
    2. El Haddad Rami & Lécot Christian & Venkiteswaran Gopalakrishnan, 2010. "Diffusion in a nonhomogeneous medium: quasi-random walk on a lattice," Monte Carlo Methods and Applications, De Gruyter, vol. 16(3-4), pages 211-230, January.

    Corrections

    All material on this site has been provided by the respective publishers and authors. You can help correct errors and omissions. When requesting a correction, please mention this item's handle: RePEc:eee:matcom:v:47:y:1998:i:2:p:153-163. See general information about how to correct material in RePEc.

    If you have authored this item and are not yet registered with RePEc, we encourage you to do it here. This allows to link your profile to this item. It also allows you to accept potential citations to this item that we are uncertain about.

    We have no bibliographic references for this item. You can help adding them by using this form .

    If you know of missing items citing this one, you can help us creating those links by adding the relevant references in the same way as above, for each refering item. If you are a registered author of this item, you may also want to check the "citations" tab in your RePEc Author Service profile, as there may be some citations waiting for confirmation.

    For technical questions regarding this item, or to correct its authors, title, abstract, bibliographic or download information, contact: Catherine Liu (email available below). General contact details of provider: http://www.journals.elsevier.com/mathematics-and-computers-in-simulation/ .

    Please note that corrections may take a couple of weeks to filter through the various RePEc services.

    IDEAS is a RePEc service. RePEc uses bibliographic data supplied by the respective publishers.