Adaptive two-regime method: application to front propagation

Abstract

The Adaptive Two-Regime Method (ATRM) is developed for hybrid (multiscale) stochastic simulation of reaction-diffusion problems. It efficiently couples detailed Brownian dynamics simulations with coarser lattice-based models. The ATRM is a generalization of the previously developed Two-Regime Method [Flegg et al, Journal of the Royal Society Interface, 2012] to multiscale problems which require a dynamic selection of regions where detailed Brownian dynamics simulation is used. Typical applications include a front propagation or spatio-temporal oscillations. In this paper, the ATRM is used for an in-depth study of front propagation in a stochastic reaction-diffusion system which has its mean-field model given in terms of the Fisher equation [Fisher, Annals of Eugenics, 1937]. It exhibits a travelling reaction front which is sensitive to stochastic fluctuations at the leading edge of the wavefront. Previous studies into stochastic effects on the Fisher wave propagation speed have focused on lattice-based models, but there has been limited progress using off-lattice (Brownian dynamics) models, which suffer due to their high computational cost, particularly at the high molecular numbers that are necessary to approach the Fisher mean-field model. By modelling only the wavefront itself with the off-lattice model, it is shown that the ATRM leads to the same Fisher wave results as purely off-lattice models, but at a fraction of the computational cost. The error analysis of the ATRM is also presented for a morphogen gradient model.

12 pages, 13 figures, submitted to Journal of Chemical Physics

References in corpus (2)

• A practical guide to stochastic simulations of reaction-diffusion processes
• Fundamental Limits to Position Determination by Concentration Gradients

Cited by in corpus (11)

• Hybrid approaches for multiple-species stochastic reaction-diffusion models
• Efficient Reactive Brownian Dynamics
• A general approximation for the dynamics of quantitative traits
• A Monte Carlo simulation for kinetic chemotaxis models: an application to the traveling population wave
• Particle-based and Meshless Methods with Aboria
• Efficient parameter sensitivity computation for spatially-extended reaction networks
• Optimisation of Simulations of Stochastic Processes by Removal of Opposing Reactions
• Spatial-Stochastic Simulation of Reaction-Diffusion Systems
• Efficient simulation techniques for biochemical reaction networks
• A variational-level-set based partitioning method for block-structured meshes
• The pseudo-compartment method for coupling PDE and compartment-based models of diffusion