Double fast algorithm for solving time-space fractional diffusion problems with spectral fractional Laplacian

Abstract

This paper presents an efficient and concise double fast algorithm to solve high dimensional time-space fractional diffusion problems with spectral fractional Laplacian. We first establish semi-discrete scheme of time-space fractional diffusion equation, which uses linear finite element or fourth-order compact difference method combining with matrix transfer technique to approximate spectral fractional Laplacian. Then we introduce a fast time-stepping L1 scheme for time discretization. The proposed scheme can exactly evaluate fractional power of matrix and perform matrix-vector multiplication at per time level by using discrete sine transform, which doesn't need to resort to any iteration method and can significantly reduce computation cost and memory requirement. Further, we address stability and convergence analyses of full discrete scheme based on fast time-stepping L1 scheme on graded time mesh. Our analysis shows that the choice of graded mesh factor $\omega =(2-\alpha )/\alpha$ shall give an optimal temporal convergence $\mathcal{O}(N^{-(2-\alpha )})$ with $N$ denoting the number of time mesh. Finally, ample numerical examples are delivered to illustrate our theoretical analysis and the efficiency of the suggested scheme.

References in corpus (6)

• Nonlocal diffusion and applications
• Numerical Methods for Fractional Diffusion
• Boundary problems for the fractional and tempered fractional operators
• A comparative study on nonlocal diffusion operators related to the fractional Laplacian
• Analysis of numerical methods for spectral fractional elliptic equations based on the best uniform rational approximation
• Numerical approximation for fractional diffusion equation forced by a tempered fractional Gaussian noise