Abstract
In this paper, we construct and analyze a linearized finite difference/Galerkin–Legendre spectral scheme for the nonlinear multiterm Caputo time fractional-order reaction-diffusion equation with time delay and Riesz space fractional derivatives. The temporal fractional orders in the considered model are taken as 0 < β 0 < β 1 < β 2 < ⋯ < β m < 1 . The problem is first approximated by the L 1 difference method on the temporal direction, and then, the Galerkin–Legendre spectral method is applied on the spatial discretization. Armed by an appropriate form of discrete fractional Grönwall inequalities, the stability and convergence of the fully discrete scheme are investigated by discrete energy estimates. We show that the proposed method is stable and has a convergent order of 2 − β m in time and an exponential rate of convergence in space. We finally provide some numerical experiments to show the efficacy of the theoretical results.