Stochastic diffusion processes on Cartesian meshes

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS

卷:294

Article

Meinecke, Lina¹  Lotstedt, Per¹ 

[1] Uppsala Univ, Dept Informat Technol, Div Comp Sci, SE-75105 Uppsala, Sweden

关键词: Stochastic simulation;    Diffusion;    Cartesian mesh;   

DOI  :  10.1016/j.cam.2015.07.035

来源: Elsevier

【 摘 要 】

Diffusion of molecules is simulated stochastically by letting them jump between voxels in a Cartesian mesh. The jump coefficients are first derived using finite difference, finite element, and finite volume approximations of the Laplacian on the mesh. An alternative is to let the first exit time for a molecule in random walk in a voxel define the jump coefficient. Such coefficients have the advantage of always being non-negative. These four different ways of obtaining the diffusion propensities are compared theoretically and in numerical experiments. A finite difference and a finite volume approximation generate the most accurate coefficients. (C) 2015 Elsevier B.V. All rights reserved.

【 授权许可】

Free