- Journal Home
- Volume 42 - 2024
- Volume 41 - 2023
- Volume 40 - 2022
- Volume 39 - 2021
- Volume 38 - 2020
- Volume 37 - 2019
- Volume 36 - 2018
- Volume 35 - 2017
- Volume 34 - 2016
- Volume 33 - 2015
- Volume 32 - 2014
- Volume 31 - 2013
- Volume 30 - 2012
- Volume 29 - 2011
- Volume 28 - 2010
- Volume 27 - 2009
- Volume 26 - 2008
- Volume 25 - 2007
- Volume 24 - 2006
- Volume 23 - 2005
- Volume 22 - 2004
- Volume 21 - 2003
- Volume 20 - 2002
- Volume 19 - 2001
- Volume 18 - 2000
- Volume 17 - 1999
- Volume 16 - 1998
- Volume 15 - 1997
- Volume 14 - 1996
- Volume 13 - 1995
- Volume 12 - 1994
- Volume 11 - 1993
- Volume 10 - 1992
- Volume 9 - 1991
- Volume 8 - 1990
- Volume 7 - 1989
- Volume 6 - 1988
- Volume 5 - 1987
- Volume 4 - 1986
- Volume 3 - 1985
- Volume 2 - 1984
- Volume 1 - 1983
Cited by
- BibTex
- RIS
- TXT
We consider the approximation of systems of reaction-diffusion equations, with the finite element method. The highest derivative in each equation is multiplied by a parameter $\varepsilon \in (0,1]$, and as $\varepsilon \rightarrow 0$ the solution of the system will contain boundary layers. We extend the analysis of the corresponding scalar problem from [Melenk, IMA J. Numer. Anal. 17(1997), pp. 577-601], to construct a finite element scheme which includes elements of size $\mathcal{O}(\varepsilon p)$ near the boundary, where $p$ is the degree of the approximating polynomials. We show that, under the assumption of analytic input data, the method yields exponential rates of convergence, independently of $\varepsilon $, when the error is measured in the energy norm associated with the problem. Numerical computations supporting the theory are also presented, which also show that the method yields robust exponential convergence rates when the error in the maximum norm is used.
We consider the approximation of systems of reaction-diffusion equations, with the finite element method. The highest derivative in each equation is multiplied by a parameter $\varepsilon \in (0,1]$, and as $\varepsilon \rightarrow 0$ the solution of the system will contain boundary layers. We extend the analysis of the corresponding scalar problem from [Melenk, IMA J. Numer. Anal. 17(1997), pp. 577-601], to construct a finite element scheme which includes elements of size $\mathcal{O}(\varepsilon p)$ near the boundary, where $p$ is the degree of the approximating polynomials. We show that, under the assumption of analytic input data, the method yields exponential rates of convergence, independently of $\varepsilon $, when the error is measured in the energy norm associated with the problem. Numerical computations supporting the theory are also presented, which also show that the method yields robust exponential convergence rates when the error in the maximum norm is used.