Lumped finite elements for reaction-cross-diffusion systems on stationary surfaces
Abstract
We consider a lumped surface finite element method (LSFEM) for the spatial approximation of reaction-diffusion equations on closed compact surfaces in R3 in the presence of cross-diffusion. We provide a fully-discrete scheme by applying the implicit-explicit (IMEX) Euler method. We provide sufficient conditions for the existence of polytopal invariant regions for the numerical solution after spatial and full discretisations. Furthermore, we prove optimal error bounds for the semi- and fully-discrete methods, that is the convergence rates are quadratic in the meshsize and linear in the timestep. To support our theoretical findings, we provide two numerical tests. The first test confirms that in the absence of lumping numerical solutions violate the invariant region leading to blow-up due to the nature of the kinetics. The second experiment is an example of Turing pattern formation in the presence of cross-diffusion on the sphere.
Citation
Frittelli , M , Madzvamuse , A , Sgura , I & Venkataraman , C 2017 , ' Lumped finite elements for reaction-cross-diffusion systems on stationary surfaces ' , Computers and Mathematics with Applications , vol. 74 , no. 12 , pp. 3008-3023 . https://doi.org/10.1016/j.camwa.2017.07.044
Publication
Computers and Mathematics with Applications
Status
Peer reviewed
ISSN
0898-1221Type
Journal article
Rights
© 2017 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).
Description
All the authors (AM, IS, CV, MF) thank the Isaac Newton Institute for Mathematical Sciences for its hospitality during the programme (Coupling Geometric PDEs with Physics for Cell Morphology, Motility and Pattern Formation; EPSRC EP/K032208/1). This work (AM) has received funding from the European Union's Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 642866. AM and CV acknowledge support from the Engineering and Physical Sciences Research Council (EP/J016780/1) on Modelling, analysis and simulation of spatial patterning on evolving biological surfaces and the Leverhulme Trust Research Project Grant (RPG-2014-149) on Unravelling new mathematics for 3D cell migration. AM was partially supported by a fellowship from the Simons Foundation. AM is a Royal Society Wolfson Research Merit Award Holder, generously funded by the Wolfson Foundation.Collections
Items in the St Andrews Research Repository are protected by copyright, with all rights reserved, unless otherwise indicated.