The hp finite element method for singularly perturbed systems of reaction-diffusion equations

  • Authors:

  • Christos Xenophontos

  • Affiliations:

  • Department of Mathematics and Statistics, University of Cyprus, Nicosia, Cyprus

  • Venue:
  • Neural, Parallel & Scientific Computations
  • Year:
  • 2008

Quantified Score

Hi-index 0.00

Visualization

Abstract

We consider the approximation of a coupled system of two singularly perturbed reaction-diffusion equations, with the finite element method. The solution to such problems contains boundary layers which overlap and interact, and the numerical approximation must take this into account in order for the resulting scheme to converge uniformly with respect to the singular perturbation parameters. We present results on a high order hpfinite element scheme which includes elements of size O(Ɛp) and O(µp) near the boundary, where Ɛ, µ are the singular perturbation parameters and pis the degree of the approximating polynomials. Under the assumption of analytic input data, the method yields exponential rates of convergence as p→ ∞, independently of Ɛ and µ. Numerical computations supporting the theory are also presented.