Aspects of an adaptive finite element method for the fractional Laplacian: a priori and a posteriori error estimates, efficient implementation and multigrid solver

Abstract

We develop all of the components needed to construct an adaptive finite element code that can be used to approximate fractional partial differential equations, on non-trivial domains in $d\ge 1$ dimensions. Our main approach consists of taking tools that have been shown to be effective for adaptive boundary element methods and, where necessary, modifying them so that they can be applied to the fractional PDE case. Improved a priori error estimates are derived for the case of quasi-uniform meshes which are seen to deliver sub-optimal rates of convergence owing to the presence of singularities. Attention is then turned to the development of an a posteriori error estimate and error indicators which are suitable for driving an adaptive refinement procedure. We assume that the resulting refined meshes are locally quasi-uniform and develop efficient methods for the assembly of the resulting linear algebraic systems and their solution using iterative methods, including the multigrid method. The storage of the dense matrices along with efficient techniques for computing the dense matrix vector products needed for the iterative solution is also considered. The performance and efficiency of the resulting algorithm is illustrated for a variety of examples.

References in corpus (3)

• Fractional Laplacians on domains, a development of Hörmander's theory of mu-transmission pseudodifferential operators
• A short FE implementation for a 2d homogeneous Dirichlet problem of a Fractional Laplacian
• Towards an Efficient Finite Element Method for the Integral Fractional Laplacian on Polygonal Domains

Cited by in corpus (27)

• Numerical methods for nonlocal and fractional models
• Towards an Efficient Finite Element Method for the Integral Fractional Laplacian on Polygonal Domains
• What Is the Fractional Laplacian?
• A priori error estimates for the optimal control of the integral fractional Laplacian
• Hybrid Finite Element - Spectral Method for the Fractional Laplacian: Approximation Theory and Efficient Solver
• $\mathcal{H}$-matrix approximability of inverses of discretizations of the fractional Laplacian
• Space-time adaptive finite elements for nonlocal parabolic variational inequalities
• Local convergence of the FEM for the integral fractional Laplacian
• Finite element approximations of the nonhomogeneous fractional Dirichlet problem
• Petrov-Galerkin Method for Fully Distributed-Order Fractional Partial Differential Equations
• A cookbook for finite element methods for nonlocal problems, including quadrature rules and approximate Euclidean balls
• Optimal operator preconditioning for pseudodifferential boundary problems
• Exponential Convergence of hp FEM for the Integral Fractional Laplacian in Polygons
• Fractional Dissipative PDEs
• A grid-overlay finite difference method for the fractional Laplacian on arbitrary bounded domains
• A novel and simple spectral method for nonlocal PDEs with the fractional Laplacian
• A simple solver for the fractional Laplacian in multiple dimensions
• Efficient quadrature rules for finite element discretizations of nonlocal equations
• Plain convergence of adaptive algorithms without exploiting reliability and efficiency
• Fast Fourier-like Mapped Chebyshev Spectral-Galerkin Methods for PDEs with Integral Fractional Laplacian in Unbounded Domains
• An a posteriori error estimator for the spectral fractional power of the Laplacian
• On the stability of Scott-Zhang type operators and application to multilevel preconditioning in fractional diffusion
• Meshfree finite difference solution of homogeneous Dirichlet problems of the fractional Laplacian
• A fast solver for spectral element approximation applied to fractional differential equations using hierarchical matrix approximation
• Isogeometric Collocation Method for the Fractional Laplacian in the 2D Bounded Domain
• A variational scheme for hyperbolic obstacle problems
• Approximating and preconditioning the stiffness matrix in the GoFD approximation of the fractional Laplacian