Numerical approximations for the tempered fractional Laplacian: Error analysis and applications

Abstract

In this paper, we propose an accurate finite difference method to discretize the $d$-dimensional (for $d\ge 1$) tempered integral fractional Laplacian and apply it to study the tempered effects on the solution of problems arising in various applications. Compared to other existing methods, our method has higher accuracy and simpler implementation. Our numerical method has an accuracy of $O(h^{\epsilon })$, for $u\in C^{0,\alpha +\epsilon }(\bar{\Omega })$ if $\alpha <1$ (or $u\in C^{1,\alpha -1+\epsilon }(\bar{\Omega })$ if $\alpha \ge 1$) with $\epsilon >0$, suggesting the minimum consistency conditions. The accuracy can be improved to $O(h^2)$, for $u\in C^{2,\alpha +\epsilon }(\bar{\Omega })$ if $\alpha <1$ (or $u\in C^{3,\alpha -1+\epsilon }(\bar{\Omega })$ if $\alpha \ge 1$). Numerical experiments confirm our analytical results and provide insights in solving the tempered fractional Poisson problem. It suggests that to achieve the second order of accuracy, our method only requires the solution $u\in C^{1,1}(\bar{\Omega })$ for any $0<\alpha <2$. Moreover, if the solution of tempered fractional Poisson problems satisfies $u\in C^{p,s}(\bar{\Omega })$ for $p=0,1$ and $0<s\le 1$, our method has the accuracy of $O(h^{p+s})$. Since our method yields a (multilevel) Toeplitz stiffness matrix, one can design fast algorithms via the fast Fourier transform for efficient simulations. Finally, we apply it together with fast algorithms to study the tempered effects on the solutions of various tempered fractional PDEs, including the Allen-Cahn equation and Gray-Scott equations.

arXiv comment21 pages, 11 figures, 3 tables

References in corpus (10)

• Fractional Quantum Mechanics and Levy Path Integrals
• Complex Patterns in a Simple System
• A fractional Laplace equation: regularity of solutions and Finite Element approximations
• Fractional Diffusion Equation for a Power-Law-Truncated Levy Process
• Accurate numerical methods for two and three dimensional integral fractional Laplacian with applications
• Computing the ground and first excited states of the fractional Schrodinger equation in an infinite potential well
• A Riesz basis Galerkin method for the tempered fractional Laplacian
• Finite difference schemes for the tempered fractional Laplacian
• Rational spectral methods for PDEs involving fractional Laplacian in unbounded domains
• Algorithm implementation and numerical analysis for the two-dimensional tempered fractional Laplacian

Cited by in corpus (4)

• A unified meshfree pseudospectral method for solving both classical and fractional PDEs
• Tackling the Curse of Dimensionality in Fractional and Tempered Fractional PDEs with Physics-Informed Neural Networks
• A novel and simple spectral method for nonlocal PDEs with the fractional Laplacian
• A universal solution scheme for fractional and classical PDEs