Importance of the Voronoi domain partition for position-jump reaction-diffusion processes on non-uniform rectilinear lattices

C.A. Yates, R.E. Baker

Research output: Contribution to journalArticlepeer-review

6 Citations (Scopus)

Abstract

Position-jump processes are used for the mathematical modeling of spatially extended chemical and biological systems with increasing frequency. A large subset of the literature concerning such processes is concerned with modeling the effect of stochasticity on reaction-diffusion systems. Traditionally, computational domains have been divided into regular voxels. Molecules are assumed well mixed within each of these voxels and are allowed to react with other molecules within the same voxel or to jump to neighboring voxels with predefined transition rates. For a variety of reasons implementing position-jump processes on irregular grids is becoming increasingly important. However, it is not immediately clear what form an appropriate irregular partition of the domain should take if it is to allow the derivation of mean molecular concentrations that agree with a given partial differential equation for molecular concentrations. It has been demonstrated, in one dimension, that the Voronoi domain partition is the appropriate method with which to divide the computational domain. In this Brief Report, we investigate theoretically the propriety of the Voronoi domain partition as an appropriate method to partition domains for position-jump models in higher dimensions. We also provide simulations of diffusion processes in two dimensions in order to corroborate our results.
Original languageEnglish
Article number054701
Number of pages4
JournalPhysical Review E
Volume88
Issue number5
DOIs
Publication statusPublished - 2013

Keywords

  • cell migration, reaction diffusion, non-uniform domain, Voronoi domain partition, rectilinear meshes, density, concentration

Fingerprint Dive into the research topics of 'Importance of the Voronoi domain partition for position-jump reaction-diffusion processes on non-uniform rectilinear lattices'. Together they form a unique fingerprint.

Cite this