Frittelli, Massimo, Madzvamuse, Anotida, Sgura, Ivonne and Venktaraman, Chandrasekhar (2017) Preserving invariance properties of reaction–diffusion systems on stationary surfaces. IMA Journal of Numerical Analysis, 2017. pp. 1-36. ISSN 0272-4979
![]() |
PDF
- Published Version
Available under License Creative Commons Attribution. Download (649kB) |
Abstract
We propose and analyse a lumped surface finite element method for the numerical approximation of reaction–diffusion systems on stationary compact surfaces in R3. 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.
Item Type: | Article |
---|---|
Schools and Departments: | School of Mathematical and Physical Sciences > Mathematics |
Subjects: | Q Science > QA Mathematics |
Depositing User: | Billy Wichaidit |
Date Deposited: | 31 Oct 2017 09:47 |
Last Modified: | 22 Feb 2018 12:05 |
URI: | http://srodev.sussex.ac.uk/id/eprint/70777 |
View download statistics for this item
📧 Request an update