Hybrid Finite Element - Spectral Method for the Fractional Laplacian: Approximation Theory and Efficient Solver

Abstract

A numerical scheme is presented for approximating fractional order Poisson problems in two and three dimensions. The scheme is based on reformulating the original problem posed over $\Omega$ on the extruded domain $\mathcal{C}=\Omega \times [0,\infty )$ following Caffarelli and Silvestre (2007). The resulting degenerate elliptic integer order PDE is then approximated using a hybrid FEM-spectral scheme. Finite elements are used in the direction parallel to the problem domain $\Omega$, and an appropriate spectral method is used in the extruded direction. The spectral part of the scheme requires that we approximate the true eigenvalues of the integer order Laplacian over $\Omega$. We derive an a priori error estimate which takes account of the error arising from using an approximation in place of the true eigenvalues. We further present a strategy for choosing approximations of the eigenvalues based on Weyl's law and finite element discretizations of the eigenvalue problem. The system of linear algebraic equations arising from the hybrid FEM-spectral scheme is decomposed into blocks which can be solved effectively using standard iterative solvers such as multigrid and conjugate gradient. Numerical examples in two and three dimensions show that the approach is quasi-optimal in terms of complexity.

References in corpus (4)

• An extension problem related to the fractional Laplacian
• Aspects of an adaptive finite element method for the fractional Laplacian: a priori and a posteriori error estimates, efficient implementation and multigrid solver
• Towards an Efficient Finite Element Method for the Integral Fractional Laplacian on Polygonal Domains
• Regularity of spectral fractional Dirichlet and Neumann problems

Cited by in corpus (11)

• Numerical methods for nonlocal and fractional models
• What Is the Fractional Laplacian?
• A unified meshfree pseudospectral method for solving both classical and fractional PDEs
• Fast implicit difference schemes for time-space fractional diffusion equations with the integral fractional Laplacian
• Local convergence of the FEM for the integral fractional Laplacian
• Exponential Convergence of $hp$ FEM for Spectral Fractional Diffusion in Polygons
• A simple solver for the fractional Laplacian in multiple dimensions
• An a posteriori error estimator for the spectral fractional power of the Laplacian
• Certified reduced basis methods for fractional Laplace equations via extension
• A multigrid-reduction-in-time solver with a new two-level convergence for unsteady fractional Laplacian problems
• A Direct Sampling Method for the Inversion of the Radon Transform