A computational method for the coupled solution of reaction-diffusion equations on evolving domains and manifolds : application to a model of cell migration and chemotaxis

MacDonald, G. and MacKenzie, J.A. and Nolan, M. and Insall, R.H. (2016) A computational method for the coupled solution of reaction-diffusion equations on evolving domains and manifolds : application to a model of cell migration and chemotaxis. Journal of Computational Physics, 309. 207–226. ISSN 0021-9991 (https://doi.org/10.1016/j.jcp.2015.12.038)

[thumbnail of MacDonald-etal-JOCP-2015-A-computational-method-for-the-coupled-solution-of-reaction-diffusion-equations-on-evolving]
Preview
Text. Filename: MacDonald_etal_JOCP_2015_A_computational_method_for_the_coupled_solution_of_reaction_diffusion_equations_on_evolving.pdf
Final Published Version
License: Creative Commons Attribution 4.0 logo

Download (3MB)| Preview

Abstract

In this paper, we devise a moving mesh finite element method for the approximate solution of coupled bulk-surface reaction-diffusion equations on an evolving two dimensional domain. Fundamental to the success of the method is the robust generation of bulk and surface meshes. For this purpose, we use a novel moving mesh partial differential equation (MMPDE) approach. The developed method is applied to model problems with known analytical solutions; these experiments indicate second-order spatial and temporal accuracy. Coupled bulk-surface problems occur frequently in many areas; in particular, in the modelling of eukaryotic cell migration and chemotaxis. We apply the method to a model of the two-way interaction of a migrating cell in a chemotactic field, where the bulk region corresponds to the extracellular region and the surface to the cell membrane.