Fast High-Order Difference Scheme for the Modified Anomalous Subdiffusion Equation Based on Fast Discrete Sine Transform

The modified anomalous subdiffusion equation plays an important role in the modeling of the processes that become less anomalous as time evolves. In this paper, we consider the efficient difference scheme for solving such time-fractional equation in two space dimensions. By using the modified L1 method and the compact difference operator with fast discrete sine transform technique, we develop a fast Crank-Nicolson compact difference scheme which is proved to be stable with the accuracy of O τ min 1 + α , 1 + β + h 4 . Here, α and β are the fractional orders which both range from 0 to 1, and τ and h are, respectively, the temporal and spatial stepsizes. We also consider the method of adding correction terms to efficiently deal with the nonsmooth problems. Numerical examples are provided to verify the effectiveness of the proposed scheme.