Numerical Approximation of the Fractional Rayleigh–Stokes Problem Arising in a Generalised Maxwell Fluid

This paper presents a numerical technique to approximate the Rayleigh–Stokes model for a generalised Maxwell fluid formulated in the Riemann–Liouville sense. The proposed method consists of two stages. First, the time discretization of the problem is accomplished by using the finite difference. Second, the space discretization is obtained by means of the predictor–corrector method. The unconditional stability result and convergence analysis are analysed theoretically. Numerical examples are provided to verify the feasibility and accuracy of the proposed method.