×

\(hp\) FEM for reaction-diffusion equations. I: Robust exponential convergence. (English) Zbl 0972.65093

Summary: 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(\text{exp} (-b p))\) and independent of the perturbation parameter \(\varepsilon\).
Additionally, we show that the use of numerical quadrature for the evaluation of the stiffness matrix and the load vector retains the exponential rate of convergence. In particular, the spectral element method based on the use of a Gauss–Lobatto quadrature rule with \((p+1)\times (p+1)\) points yields robust exponential convergence.

MSC:

65N30 Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs
35B25 Singular perturbations in context of PDEs
35J65 Nonlinear boundary value problems for linear elliptic equations
65N35 Spectral, collocation and related methods for boundary value problems involving PDEs
65N12 Stability and convergence of numerical methods for boundary value problems involving PDEs
PDFBibTeX XMLCite
Full Text: DOI