{"title":"Stochastic models of advection-diffusion in layered media","authors":"Elliot J. Carr","doi":"arxiv-2409.08447","DOIUrl":null,"url":null,"abstract":"Mathematically modelling diffusive and advective transport of particles in\nheterogeneous layered media is important to many applications in computational,\nbiological and medical physics. While deterministic continuum models of such\ntransport processes are well established, they fail to account for randomness\ninherent in many problems and are valid only for a large number of particles.\nTo address this, this paper derives a suite of equivalent stochastic\n(discrete-time discrete-space random walk) models for several standard\ncontinuum (partial differential equation) models of diffusion and\nadvection-diffusion across a fully- or semi-permeable interface. Our approach\ninvolves discretising the continuum model in space and time to yield a Markov\nchain, which governs the transition probabilities between spatial lattice sites\nduring each time step. Discretisation in space is carried out using a standard\nfinite volume method while two options are considered for discretisation in\ntime. A simple forward Euler discretisation yields a stochastic model taking\nthe form of a local (nearest-neighbour) random walk with simple analytical\nexpressions for the transition probabilities while an exact exponential\ndiscretisation yields a non-local random walk with transition probabilities\ndefined numerically via a matrix exponential. Constraints on the size of the\nspatial and/or temporal steps are provided for each option to ensure the\ntransition probabilities are non-negative. MATLAB code comparing the stochastic\nand continuum models is available on GitHub\n(https://github.com/elliotcarr/Carr2024c) with simulation results demonstrating\ngood agreement for several example problems.","PeriodicalId":501369,"journal":{"name":"arXiv - PHYS - Computational Physics","volume":"21 1","pages":""},"PeriodicalIF":0.0000,"publicationDate":"2024-09-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"arXiv - PHYS - Computational Physics","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/arxiv-2409.08447","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0
Abstract
Mathematically modelling diffusive and advective transport of particles in
heterogeneous layered media is important to many applications in computational,
biological and medical physics. While deterministic continuum models of such
transport processes are well established, they fail to account for randomness
inherent in many problems and are valid only for a large number of particles.
To address this, this paper derives a suite of equivalent stochastic
(discrete-time discrete-space random walk) models for several standard
continuum (partial differential equation) models of diffusion and
advection-diffusion across a fully- or semi-permeable interface. Our approach
involves discretising the continuum model in space and time to yield a Markov
chain, which governs the transition probabilities between spatial lattice sites
during each time step. Discretisation in space is carried out using a standard
finite volume method while two options are considered for discretisation in
time. A simple forward Euler discretisation yields a stochastic model taking
the form of a local (nearest-neighbour) random walk with simple analytical
expressions for the transition probabilities while an exact exponential
discretisation yields a non-local random walk with transition probabilities
defined numerically via a matrix exponential. Constraints on the size of the
spatial and/or temporal steps are provided for each option to ensure the
transition probabilities are non-negative. MATLAB code comparing the stochastic
and continuum models is available on GitHub
(https://github.com/elliotcarr/Carr2024c) with simulation results demonstrating
good agreement for several example problems.