A fourth-order extrapolated compact difference method for time-fractional convection-reaction-diffusion equations with spatially variable coefficients

This paper is concerned with numerical methods for a class of time-fractional convection-reaction-diffusion equations. The convection and reaction coefficients of the equation may be spatially variable. Based on the weighted and shifted Grünwald-Letnikov formula for the time-fractional derivative and a compact finite difference approximation for the spatial derivative, we establish an unconditionally stable compact difference method. The local truncation error and the solvability of the resulting scheme are discussed in detail. The stability of the method and its convergence of third-order in time and fourth-order in space are rigorously proved by the discrete energy method. Combining this method with a Richardson extrapolation, we present an extrapolated compact difference method which is fourth-order accurate in both time and space. A rigorous proof for the convergence of the extrapolation method is given. Numerical results confirm our theoretical analysis, and demonstrate the accuracy of the compact difference method and the effectiveness of the extrapolated compact difference method.