Show simple item record

Files in this item


Item metadata

dc.contributor.authorFrittelli, Massimo
dc.contributor.authorMadzvamuse, Anotide
dc.contributor.authorSgura, Ivonne
dc.contributor.authorVenkataraman, Chandrasekhar
dc.identifier.citationFrittelli , M , Madzvamuse , A , Sgura , I & Venkataraman , C 2017 , ' Preserving invariance properties of reaction-diffusion systems on stationary surfaces ' , IMA Journal of Numerical Analysis , vol. Advance articles .
dc.identifier.otherPURE: 250911781
dc.identifier.otherPURE UUID: 42b806ce-990d-43a5-8504-82db8e470aba
dc.identifier.otherScopus: 85054132947
dc.identifier.otherWOS: 000491255100008
dc.descriptionThis work (AM, CV) is partly supported by the EPSRC grant number EP/J016780/1 and the Leverhulme Trust Research Project Grant (RPG-2014-149). The authors (MF, AM, IS CV) would like to 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] supported by EPSRC Grant Number EP/K032208/1. AM acknowledges funding from the European Union Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 642866 and was partially supported by a grant from the Simons Foundation. AM is a Royal Society Wolfson Research Merit Award Holder funded generously by the Wolfson Foundation.en
dc.description.abstractWe propose and analyse a lumped surface finite element method for the numerical approximation of reaction–diffusion systems on stationary compact surfaces in ℝ3. The proposed method preserves the invariant regions of the continuous problem under discretization and, in the special case of scalar equations, it preserves the maximum principle. On the application of a fully discrete scheme using the implicit–explicit Euler method in time, we prove that invariant regions of the continuous problem are preserved (i) at the spatially discrete level with no restriction on the meshsize and (ii) at the fully discrete level under a timestep restriction. We further prove optimal error bounds for the semidiscrete and fully discrete methods, that is, the convergence rates are quadratic in the meshsize and linear in the timestep. Numerical experiments are provided to support the theoretical findings. We provide examples in which, in the absence of lumping, the numerical solution violates the invariant region leading to blow-up.
dc.relation.ispartofIMA Journal of Numerical Analysisen
dc.rights© The authors 2017. Published by Oxford University Press on behalf of the Institute of Mathematics and its Applications. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.en
dc.subjectSurface finite elementsen
dc.subjectMass lumpingen
dc.subjectInvariant regionen
dc.subjectMaximum principleen
dc.subjectHeat equationen
dc.subjectPattern formationen
dc.subjectQA Mathematicsen
dc.subjectQC Physicsen
dc.titlePreserving invariance properties of reaction-diffusion systems on stationary surfacesen
dc.typeJournal articleen
dc.description.versionPublisher PDFen
dc.contributor.institutionUniversity of St Andrews. Applied Mathematicsen
dc.description.statusPeer revieweden

This item appears in the following Collection(s)

Show simple item record