The blending region hybrid framework for the simulation of stochastic reaction-diffusion processes

Christian A Yates, Adam George, Armand Jordana, Cameron A. Smith, Andrew B. Duncan, Konstantinos C Zygalakis

Research output: Contribution to journalArticlepeer-review

Abstract

The simulation of stochastic reaction-diffusion systems using fine-grained representations can become computationally prohibitive when particle numbers become large. If particle numbers are sufficiently high then it may be possible to ignore stochastic fluctuations and use a more efficient coarse-grained simulation approach. Nevertheless, for multiscale systems which exhibit significant spatial variation in concentration, a coarse-grained approach may not be appropriate throughout the simulation domain. Such scenarios suggest a hybrid paradigm in which a computationally cheap, coarse-grained model is coupled to a more expensive, but more detailed fine-grained model enabling the accurate simulation of the fine-scale dynamics at a reasonable computational cost. In this paper, in order to couple two representations of reaction-diffusion at distinct spatial scales, we allow them to overlap in a "blending region". Both modelling paradigms provide a valid representation of the particle density in this region. From one end of the blending region to the other, control of the implementation of diffusion is passed from one modelling paradigm to another through the use of complementary "blending functions" which scale up or down the contribution of each model to the overall diffusion. We establish the reliability of our novel hybrid paradigm by demonstrating its simulation on four exemplar reaction-diffusion scenarios.
Original languageEnglish
Article number20200563
Number of pages36
JournalJournal of the Royal Society, Interface
Volume17
Issue number171
Early online date21 Oct 2020
DOIs
Publication statusPublished - 28 Oct 2020

Fingerprint Dive into the research topics of 'The blending region hybrid framework for the simulation of stochastic reaction-diffusion processes'. Together they form a unique fingerprint.

Cite this