Loading…

Stochastic diffusion processes on Cartesian meshes

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...

Full description

Saved in:
Bibliographic Details
Published in:Journal of computational and applied mathematics 2016-03, Vol.294, p.1-11
Main Authors: Meinecke, Lina, Lötstedt, Per
Format: Article
Language:English
Subjects:
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary: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.
ISSN:0377-0427
1879-1778
1879-1778
DOI:10.1016/j.cam.2015.07.035