Numerical simulation of a class of fractional subdiffusion equations via the alternating direction implicit method

In this article, a new numerical technique is proposed for solving the two-dimensional time fractional subdiffusion equation with nonhomogeneous terms. After a transformation of the original problem, standard central difference approximation is used for the spatial discretization. For the time step, a new fractional alternating direction implicit (FADI) scheme based on the L₁ approximation is considered. This FADI scheme is constructed by adding a small term, so it is different from standard FADI methods. The solvability, unconditional stability and H¹ norm convergence are proved. Finally, numerical examples show the effectiveness and accuracy of our proposed method.