Monte Carlo simulation with fixed steplength for diffusion processes in nonhomogeneous media

JOURNAL OF COMPUTATIONAL PHYSICS

卷:239

Article

Hoyuelos, M.¹ 

[1] Univ Nacl Mar del Plata, Fac Ciencias Exactas & Nat, Inst Invest Fis Mar del Plata, CONICET UNMdP, RA-7600 Mar Del Plata, Buenos Aires, Argentina

关键词: Diffusion;    Monte Carlo;    Simulation;   

DOI  :  10.1016/j.jcp.2012.12.029

来源: Elsevier

【 摘 要 】

Monte Carlo simulation is one of the most important tools in the study of diffusion processes. For constant diffusion coefficients, an appropriate Gaussian distribution of particle's steplengths can generate exact results, when compared with integration of the diffusion equation. It is important to notice that the same method is completely erroneous when applied to non-homogeneous diffusion coefficients. A simple alternative, jumping at fixed steplengths with appropriate transition probabilities, produces correct results. Here, a model for diffusion of calcium ions in the neuromuscular junction of the crayfish is used as a test to compare Monte Carlo simulation with fixed and Gaussian steplength. (C) 2013 Elsevier Inc. All rights reserved.

【 授权许可】

Free