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

DOI10.1137/17M1144696zbMATH Open1402.65148arXiv1709.01639OpenAlexW2751463447WikidataQ129451279 ScholiaQ129451279MaRDI QIDQ4580288FDOQ4580288

Mark Ainsworth, Christian Glusa

Publication date: 14 August 2018

Published in: SIAM Journal on Scientific Computing (Search for Journal in Brave)

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 $mathcalC=Omegaimes[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.

Full work available at URL: https://arxiv.org/abs/1709.01639

zbMATH Keywords

eigenvaluefractional Laplaciannumerical approximationfinite element-spectral methodinteger Laplacian

Mathematics Subject Classification ID

Complexity and performance of numerical algorithms (65Y20) Iterative numerical methods for linear systems (65F10) A priori estimates in context of PDEs (35B45) Estimates of eigenvalues in context of PDEs (35P15) Fractional partial differential equations (35R11) Numerical methods for eigenvalue problems for boundary value problems involving PDEs (65N25) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30) Spectral, collocation and related methods for boundary value problems involving PDEs (65N35)

Cites Work

• Title not available (Why is that?)
• Title not available (Why is that?)
• Toward the optimal preconditioned eigensolver: Locally optimal block preconditioned conjugate gradient method
• A PDE approach to fractional diffusion in general domains: a priori error analysis
• Numerical approximation of fractional powers of elliptic operators
• Theory and practice of finite elements.
• On the spectrum of two different fractional operators
• Finite element approximation of eigenvalue problems
• Extension Problem and Harnack's Inequality for Some Fractional Operators
• Finite Element Interpolation of Nonsmooth Functions Satisfying Boundary Conditions
• An Extension Problem Related to the Fractional Laplacian
• On the Schrödinger equation and the eigenvalue problem
• Second term of the spectral asymptotic expansion of the Laplace-Beltrami operator on manifolds with boundary
• Regularity of spectral fractional Dirichlet and Neumann problems
• Multilevel methods for nonuniformly elliptic operators and fractional diffusion
• Aspects of an adaptive finite element method for the fractional Laplacian: a priori and a posteriori error estimates, efficient implementation and multigrid solver
• Tomorrow’s Dynamics
• Hybrid Finite Element--Spectral Method for the Fractional Laplacian: Approximation Theory and Efficient Solver

Cited In (32)

• Efficient spectral methods for PDEs with spectral fractional Laplacian
• An efficient and accurate numerical method for the spectral fractional Laplacian equation
• Local Convergence of the FEM for the Integral Fractional Laplacian
• Model reduction for fractional elliptic problems using Kato's formula
• A Fast Solver for the Fractional Helmholtz Equation
• The finite element method for fractional diffusion with spectral fractional Laplacian
• A multigrid-reduction-in-time solver with a new two-level convergence for unsteady fractional Laplacian problems
• What is the fractional Laplacian? A comparative review with new results
• Reduced Basis Methods for Fractional Laplace Equations via Extension
• A Quadrature Scheme for Steady-State Diffusion Equations Involving Fractional Power of Regularly Accretive Operator
• A Unified Meshfree Pseudospectral Method for Solving Both Classical and Fractional PDEs
• A unified view of some numerical methods for fractional diffusion
• Tensor product method for fast solution of optimal control problems with fractional multidimensional Laplacian in constraints
• Solving fractional Laplacian visco-acoustic wave equations on complex-geometry domains using Grünwald-formula based radial basis collocation method
• IMEX Hermite--Galerkin Spectral Schemes with Adaptive Time Stepping for the Coupled Nonlocal Gordon-Type Systems in Multiple Dimensions
• High Order Conservative Finite Difference/Fourier Spectral Methods for Inviscid Surface Quasi-Geostrophic Flows
• A Hermite spectral method for fractional convection diffusion equations on unbounded domains
• A Probabilistic Scheme for Semilinear Nonlocal Diffusion Equations with Volume Constraints
• Exponential convergence of \textit{hp} FEM for spectral fractional diffusion in polygons
• A reduced basis method for fractional diffusion operators. II
• Hybrid Finite Element--Spectral Method for the Fractional Laplacian: Approximation Theory and Efficient Solver
• Iterative low-rank approximation solvers for the extension method for fractional diffusion
• An a posteriori error estimator for the spectral fractional power of the Laplacian
• Quasi-optimal convergence rate for an adaptive method for the integral fractional Laplacian
• Fast Fourier-like Mapped Chebyshev Spectral-Galerkin Methods for PDEs with Integral Fractional Laplacian in Unbounded Domains
• A reduced basis method for fractional diffusion operators. I
• Spectral, tensor and domain decomposition methods for fractional PDEs
• Hybrid fully spectral linearized scheme for time‐fractional evolutionary equations
• A Simple Solver for the Fractional Laplacian in Multiple Dimensions
• Dispersion analysis of the gradient weighted finite element method for acoustic problems in one, two, and three dimensions
• The numerical solution of fractional integral equations via orthogonal polynomials in fractional powers
• Numerical methods for nonlocal and fractional models

Uses Software

• lobpcg.m

Recommendations

• Exponentially accurate spectral and spectral element methods for fractional ODEs 👍 👎
• Efficient spectral methods for PDEs with spectral fractional Laplacian 👍 👎
• A fast solver for spectral elements applied to fractional differential equations using hierarchical matrix approximation 👍 👎
• An efficient and accurate numerical method for the spectral fractional Laplacian equation 👍 👎
• The spectral collocation method for efficiently solving PDEs with fractional Laplacian 👍 👎
• A Hybrid Spectral Element Method for Fractional Two-Point Boundary Value Problems 👍 👎
• Finite element approximation for the fractional eigenvalue problem 👍 👎
• Hybrid fully spectral linearized scheme for time‐fractional evolutionary equations 👍 👎
• A Petrov-Galerkin spectral element method for fractional elliptic problems 👍 👎
• The finite element method for fractional diffusion with spectral fractional Laplacian 👍 👎

This page was built for publication: Hybrid Finite Element--Spectral Method for the Fractional Laplacian: Approximation Theory and Efficient Solver

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q4580288)