Lumped finite elements for reaction-cross-diffusion systems on stationary surfaces

Massimo Frittelli, Anotida Madzvamuse, Ivonne Sgura, Chandrasekhar Venkataraman

Research output: Contribution to journalArticlepeer-review

16 Citations (Scopus)
1 Downloads (Pure)

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.
Original languageEnglish
Pages (from-to)3008-3023
JournalComputers and Mathematics with Applications
Volume74
Issue number12
Early online date30 Aug 2017
DOIs
Publication statusPublished - 15 Dec 2017

Keywords

  • Surface finite elements
  • Mass lumping
  • Invariant region
  • Reaction-cross-diffusion
  • Convergence analysis
  • Pattern formation
  • Rosenzweig-MacArthur

Fingerprint

Dive into the research topics of 'Lumped finite elements for reaction-cross-diffusion systems on stationary surfaces'. Together they form a unique fingerprint.

Cite this