Error Analysis of Nonlinear Time Fractional Mobile/Immobile Advection-Diffusion Equation with Weakly Singular Solutions

In this paper, a weighted and shifted Grünwald-Letnikov difference (WSGD) Legendre spectral method is proposed to solve the two-dimensional nonlinear time fractional mobile/immobile advection-dispersion equation. We introduce the correction method to deal with the singularity in time, and the stability and convergence analysis are proven. In the numerical implementation, a fast method is applied based on a globally uniform approximation of the trapezoidal rule for the integral on the real line to decrease the memory requirement and computational cost. The memory requirement and computational cost are O(Q) and O(QK), respectively, where K is the number of the final time step and Q is the number of quadrature points used in the trapezoidal rule. Some numerical experiments are given to confirm our theoretical analysis and the effectiveness of the presented methods.