A Riesz basis Galerkin method for the tempered fractional Laplacian

TL;DR

This paper introduces a Riesz basis Galerkin method for solving equations involving the tempered fractional Laplacian, providing theoretical analysis and efficient computational techniques for boundary value problems.

Contribution

The work develops a new Galerkin finite element approach using Riesz bases for the tempered fractional Laplacian, including convergence analysis and efficient matrix generation methods.

Findings

Proved well-posedness of the Galerkin formulation.

Established convergence of the finite element methods.

Validated methods through numerical experiments.

Abstract

The fractional Laplacian $\Delta^{\beta/2}$ is the generator of $\beta$-stable L\'evy process, which is the scaling limit of the L\'evy fight. Due to the divergence of the second moment of the jump length of the L\'evy fight it is not appropriate as a physical model in many practical applications. However, using a parameter $\lambda$ to exponentially temper the isotropic power law measure of the jump length leads to the tempered L\'evy fight, which has finite second moment. For short time the tempered L\'evy fight exhibits the dynamics of L\'evy fight while after sufficiently long time it turns to normal diffusion. The generator of tempered $\beta$-stable L\'evy process is the tempered fractional Laplacian $(\Delta+\lambda)^{\beta/2}$ [W.H. Deng, B.Y. Li, W.Y. Tian, and P.W. Zhang, Multiscale Model. Simul., in press, 2017]. In the current work, we present new computational methods for the tempered fractional Laplacian equation, including the cases with the homogeneous and nonhomogeneous generalized Dirichlet type boundary conditions. We prove the well-posedness of the Galerkin weak formulation and provide convergence analysis of the single scaling B-spline and multiscale Riesz bases finite element methods. We propose a technique for efficiently generating the entries of the dense stiffness matrix and for solving the resulting algebraic equation by preconditioning. We also present several numerical experiments to verify the theoretical results.