Exponential convergence of a generalized FEM for heterogeneous reaction-diffusion equations

TL;DR

This paper introduces a generalized finite element method for heterogeneous reaction-diffusion equations that achieves exponential convergence and uniform accuracy with respect to the perturbation parameter, verified through numerical experiments.

Contribution

The paper develops a new GFE method with exponential decay error estimates for reaction-diffusion equations, including cases with large oversampling sizes, improving robustness and efficiency.

Findings

Exponential decay of local approximation errors with respect to oversampling size and degrees of freedom.

Uniform convergence of the method with respect to the perturbation parameter .

Numerical results confirming theoretical error estimates and convergence properties.

Abstract

A generalized finite element method is proposed for solving a heterogeneous reaction-diffusion equation with a singular perturbation parameter $\varepsilon$, based on locally approximating the solution on each subdomain by solution of a local reaction-diffusion equation and eigenfunctions of a local eigenproblem. These local problems are posed on some domains slightly larger than the subdomains with oversampling size $\delta^{\ast}$. The method is formulated at the continuous level as a direct discretization of the continuous problem and at the discrete level as a coarse-space approximation for its standard FE discretizations. Exponential decay rates for local approximation errors with respect to $\delta^{\ast}/\varepsilon$ and $\delta^{\ast}/h$ (at the discrete level with $h$ denoting the fine FE mesh size) and with the local degrees of freedom are established. In particular, it is shown that the method at the continuous level converges uniformly with respect to $\varepsilon$ in the standard $H^{1}$ norm, and that if the oversampling size is relatively large with respect to $\varepsilon$ and $h$ (at the discrete level), the solutions of the local reaction-diffusion equations provide good local approximations for the solution and thus the local eigenfunctions are not needed. Numerical results are provided to verify the theoretical results.