A novel and simple spectral method for nonlocal PDEs with the fractional Laplacian

Abstract

We propose a novel and simple spectral method based on the semi-discrete Fourier transforms to discretize the fractional Laplacian $(-\Delta {)}^{\frac{\alpha }{2}}$. Numerical analysis and experiments are provided to study its performance. Our method has the same symbol $|\mathbf{\xi }{|}^{\alpha }$ as the fractional Laplacian $(-\Delta {)}^{\frac{\alpha }{2}}$ at the discrete level, and thus it can be viewed as the exact discrete analogue of the fractional Laplacian. This {\it unique feature} distinguishes our method from other existing methods for the fractional Laplacian. Note that our method is different from the Fourier pseudospectral methods in the literature which are usually limited to periodic boundary conditions (see Remark \ref{remark0}). Numerical analysis shows that our method can achieve a spectral accuracy. The stability and convergence of our method in solving the fractional Poisson equations were analyzed. Our scheme yields a multilevel Toeplitz stiffness matrix, and thus fast algorithms can be developed for efficient matrix-vector multiplications. The computational complexity is $\mathcal{O}(2N\log (2N))$, and the memory storage is $\mathcal{O}(N)$ with $N$ the total number of points. Extensive numerical experiments verify our analytical results and demonstrate the effectiveness of our method in solving various problems.

References in corpus (10)

• Ten equivalent definitions of the fractional Laplace operator
• Fractional Laplacian in Bounded Domains
• A fractional Laplace equation: regularity of solutions and Finite Element approximations
• A short FE implementation for a 2d homogeneous Dirichlet problem of a 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
• Accurate numerical methods for two and three dimensional integral fractional Laplacian with applications
• A comparative study on nonlocal diffusion operators related to the fractional Laplacian
• Computing the ground and first excited states of the fractional Schrodinger equation in an infinite potential well
• A unified meshfree pseudospectral method for solving both classical and fractional PDEs
• Numerical approximations for the tempered fractional Laplacian: Error analysis and applications

Cited by in corpus (1)

• Approximating and preconditioning the stiffness matrix in the GoFD approximation of the fractional Laplacian