Report 1997-03

hp FEM for Reaction-Diffusion Equations. I: Robust Exponential Convergence

J.M. Melenk and C. Schwab

Abstract: A singularly perturbed reaction-diffusion equation in two dimensions is considered. We assume analyticity of the input data, i.e., the boundary of the domain is an analytic curve and the right hand side is analytic. We show that the hp version of the finite element method leads to robust exponential convergence provided that one layer of needle elements of width $O(p \varepsilon)$ is inserted near the domain boundary, that is, the rate of convergence is $O({ exp} (-b p))$ and independent of the perturbation parameter $\varepsilon$. We also show that the Spectral Element Method based on the use of a Gauss-Lobatto quadrature rule of order O(p) for the evaluation of the stiffness matrix and the load vector retains the exponential rate of convergence.

Keywords: boundary layer, singularly perturbed problem, asymptotic expansions, error bounds

Paper: Available as PDF (529 KB) or as hardcopy to order reports@sam.math.ethz.ch.

Publishing information: SIAM J. Num. Anal., Vol. 35, (1998), p. 1520-1557