Show simple item record

Files in this item


Item metadata

dc.contributor.authorFrittelli, Massimo
dc.contributor.authorMadzvamuse, Anotido
dc.contributor.authorSgura, Ivonne
dc.contributor.authorVenkataraman, Chandrasekhar
dc.identifier.citationFrittelli , M , Madzvamuse , A , Sgura , I & Venkataraman , C 2018 , ' Numerical preservation of velocity induced invariant regions for reaction-diffusion systems on evolving surfaces ' , Journal of Scientific Computing , vol. First Online .
dc.identifier.otherPURE: 253026619
dc.identifier.otherPURE UUID: b4bb042b-acb6-46c1-b31b-eb244be8287d
dc.identifier.otherScopus: 85047938572
dc.identifier.otherWOS: 000446594600013
dc.descriptionThe 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. MF’s and IS’s research work has been performed under the auspices of the Italian National Group for Scientific Calculus (GNCS-INdAM). This work (AM, CV) is partly supported by the EPSRC grant number EP/J016780/1 and the Leverhulme Trust Research Project Grant (RPG-2014-149). AM acknowledges funding from the European Union Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No 642866.en
dc.description.abstractWe propose and analyse a finite element method with mass lumping (LESFEM) for the numerical approximation of reaction–diffusion systems (RDSs) on surfaces in ℝ3 that evolve under a given velocity field. A fully-discrete method based on the implicit–explicit (IMEX) Euler time-discretisation is formulated and dilation rates which act as indicators of the surface evolution are introduced. Under the assumption that the mesh preserves the Delaunay regularity under evolution, we prove a sufficient condition, that depends on the dilation rates, for the existence of invariant regions (i) at the spatially discrete level with no restriction on the mesh size and (ii) at the fully-discrete level under a timestep restriction that depends on the kinetics, only. In the specific case of the linear heat equation, we prove a semi- and a fully-discrete maximum principle. For the well-known activator-depleted and Thomas reaction–diffusion models we prove the existence of a family of rectangles in the phase space that are invariant only under specific growth laws. Two numerical examples are provided to computationally demonstrate (i) the discrete maximum principle and optimal convergence for the heat equation on a linearly growing sphere and (ii) the existence of an invariant region for the LESFEM–IMEX Euler discretisation of a RDS on a logistically growing surface.
dc.relation.ispartofJournal of Scientific Computingen
dc.rights© The Author(s) 2018. This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.en
dc.subjectEvolving surfaceen
dc.subjectDilation rateen
dc.subjectHeat equationen
dc.subjectMaximum principleen
dc.subjectInvariant regionen
dc.subjectQA75 Electronic computers. Computer scienceen
dc.subjectQC Physicsen
dc.subjectQH301 Biologyen
dc.titleNumerical preservation of velocity induced invariant regions for reaction-diffusion systems on evolving 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