Abstract In this paper, we consider a time-dependent diffusion problem with two-sided Riemann-Liouville fractional derivatives.By introducing a fractional-order flux as auxiliary variable, we establish the saddle-point camo party invitations variational formulation, based on which we employ a locally conservative mixed finite element method to approximate the unknown function, its derivative and the fractional flux in space and use the backward Euler scheme to discrete essie well nested energy the time derivative, and thus propose a fully discrete expanded mixed finite element procedure.We prove the well-posedness and the optimal order error estimates of the proposed procedure for a sufficiently smooth solution.
Numerical experiments are presented to confirm our theoretical findings.