Numerical approximations for fractional elliptic equations via the method of semigroups

Abstract

We provide a novel approach to the numerical solution of the family of nonlocal elliptic equations $(-\Delta {)}^{s}u=f$ in $\Omega$, subject to some homogeneous boundary conditions $\mathcal{B}(u)=0$ on $\partial \Omega$, where $s\in (0,1)$, $\Omega \subset {\mathbb{R}}^n$ is a bounded domain, and $(-\Delta {)}^s$ is the spectral fractional Laplacian associated to $\mathcal{B}$ on $\partial \Omega$. We use the solution representation $(-\Delta {)}^{-s}f$ together with its singular integral expression given by the method of semigroups. By combining finite element discretizations for the heat semigroup with monotone quadratures for the singular integral we obtain accurate numerical solutions. Roughly speaking, given a datum $f$ in a suitable fractional Sobolev space of order $r\ge 0$ and the discretization parameter $h>0$, our numerical scheme converges as $O(h^{r+2s})$, providing super quadratic convergence rates up to $O(h^4)$ for sufficiently regular data, or simply $O(h^{2s})$ for merely $f\in L^2(\Omega )$. We also extend the proposed framework to the case of nonhomogeneous boundary conditions and support our results with some illustrative numerical tests.

arXiv comment26 pages, 5 figures and 2 tables. To appear in ESAIM Mathematical Modelling and Numerical Analysis

References in corpus (7)

• Numerical solving the boundary value problem for fractional powers of elliptic operators
• What Is the Fractional Laplacian?
• Regularity theory and high order numerical methods for the (1d)-Fractional Laplacian
• A PDE approach to fractional diffusion: a posteriori error analysis
• Fractional Laplacians in bounded domains: Killed, reflected, censored and taboo Lévy flights
• Derivation of the nonlocal pressure form of the fractional porous medium equation in the hydrological setting
• Finite difference methods for fractional Laplacians

Cited by in corpus (1)

• Approximate representation of the solutions of fractional elliptical BVP through the solution of parabolic IVP