Fast fourier-like mapped chebyshev spectral-galerkin methods for pdes with integral fractional laplacian in unbounded domains

Changtao Sheng, Jie Shen, Tao Tang, Li-Lian Wang, Huifang Yuan

Research output: Contribution to journalJournal articlepeer-review

33 Citations (Scopus)

Abstract

In this paper, we propose a fast spectral-Galerkin method for solving PDEs involving an integral fractional Laplacian in Rd, which is built upon two essential components: (i) the Dunford- Taylor formulation of the fractional Laplacian; and (ii) Fourier-like biorthogonal mapped Chebyshev functions (MCFs) as basis functions. As a result, the fractional Laplacian can be fully diagonalized, and the complexity of solving an elliptic fractional PDE is quasi-optimal, i.e., O((N log2 N)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.

Original languageEnglish
Pages (from-to)2435-2464
Number of pages30
JournalSIAM Journal on Numerical Analysis
Volume58
Issue number5
DOIs
Publication statusPublished - 2020

Scopus Subject Areas

  • Numerical Analysis
  • Computational Mathematics
  • Applied Mathematics

User-Defined Keywords

  • Biorthogonal basis functions
  • Dunford-Taylor formula
  • Integral fractional Laplacian
  • Mapped Chebyshev functions
  • Nonlocal/singular operators

Fingerprint

Dive into the research topics of 'Fast fourier-like mapped chebyshev spectral-galerkin methods for pdes with integral fractional laplacian in unbounded domains'. Together they form a unique fingerprint.

Cite this