Spectral methods for two-dimensional space and time fractional Bloch-Torrey equations

Abstract

In this paper, we consider the numerical approximation of the space and time fractional Bloch-Torrey equations. A fully discrete spectral scheme based on a finite difference method in the time direction and a Galerkin-Legendre spectral method in the space direction is developed. In order to reduce the amount of computation, an alternating direction implicit (ADI) spectral scheme is proposed. Then the stability and convergence analysis are rigorously established. Finally, numerical results are presented to support our theoretical analysis.