English

Numerical approximations for fractional elliptic equations via the method of semigroups

Numerical Analysis 2019-10-31 v3 Numerical Analysis Analysis of PDEs

Abstract

We provide a novel approach to the numerical solution of the family of nonlocal elliptic equations (Δ)su=f(-\Delta)^su=f in Ω\Omega, subject to some homogeneous boundary conditions B(u)=0\mathcal{B}(u)=0 on Ω\partial \Omega, where s(0,1)s\in(0,1), ΩRn\Omega\subset \mathbb{R}^n is a bounded domain, and (Δ)s(-\Delta)^s is the spectral fractional Laplacian associated to B\mathcal{B} on Ω\partial \Omega. We use the solution representation (Δ)sf(-\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 ff in a suitable fractional Sobolev space of order r0r\geq 0 and the discretization parameter h>0h>0, our numerical scheme converges as O(hr+2s)O(h^{r+2s}), providing super quadratic convergence rates up to O(h4)O(h^4) for sufficiently regular data, or simply O(h2s)O(h^{2s}) for merely fL2(Ω)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.

Keywords

Cite

@article{arxiv.1812.01518,
  title  = {Numerical approximations for fractional elliptic equations via the method of semigroups},
  author = {Nicole Cusimano and Félix del Teso and Luca Gerardo-Giorda},
  journal= {arXiv preprint arXiv:1812.01518},
  year   = {2019}
}

Comments

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