Fast Fourier-like Mapped Chebyshev Spectral-Galerkin Methods for PDEs with Integral Fractional Laplacian in Unbounded Domains

by   Changtao Sheng, et al.

In this paper, we propose a fast spectral-Galerkin method for solving PDEs involving integral fractional Laplacian in R^d, which is built upon two essential components: (i) the Dunford-Taylor formulation of the fractional Laplacian; and (ii) Fourier-like bi-orthogonal mapped Chebyshev functions (MCFs) as basis functions. As a result, the fractional Laplacian can be fully diagonalised, and the complexity of solving an elliptic fractional PDE is quasi-optimal, i.e., O((N_2N)^d) with N being the number of modes in each spatial direction. Ample numerical tests for various decaying exact solutions show that the convergence of the fast solver perfectly matches the order of theoretical error estimates. With a suitable time-discretization, the fast solver can be directly applied to a large class of nonlinear fractional PDEs. As an example, we solve the fractional nonlinear Schrödinger equation by using the fourth-order time-splitting method together with the proposed MCF-spectral-Galerkin method.



There are no comments yet.


page 20


Generalised Hermite functions and their applications in Spectral Approximations

In 1939, G. Szegö first introduced a family of generalised Hermite polyn...

Highly accurate operator factorization methods for the integral fractional Laplacian and its generalization

In this paper, we propose a new class of operator factorization methods ...

Computing solution landscape of nonlinear space-fractional problems via fast approximation algorithm

The nonlinear space-fractional problems often allow multiple stationary ...

A universal solution scheme for fractional and classical PDEs

We propose a unified meshless method to solve classical and fractional P...

A unified meshfree pseudospectral method for solving both classical and fractional PDEs

In this paper, we propose a meshfree method based on the Gaussian radial...

Numerical Approximation of the Fractional Laplacian on R Using Orthogonal Families

In this paper, using well-known complex variable techniques, we compute ...

A contour method for time-fractional PDEs and an application to fractional viscoelastic beam equations

We develop a rapid and accurate contour method for the solution of time-...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.