Computers and Mathematics with Applications, cilt.128, ss.108-120, 2022 (SCI-Expanded, Scopus)
We consider the multi-term time fractional diffusion equation on a bounded convex domain. We analyze the L1 method on a graded mesh in time to compensate for the weak singularity of the solution and a weak Galerkin finite element method in space discretization. The stability analyses are presented for both semi-discrete and fully-discrete schemes and we prove that two schemes are unconditionally stable. Error estimates in L2-norm and a discrete H1 equivalent norm for both schemes are rigorously derived. Further we discuss the optimal spatial order error estimates in L2-norm. Finally, we give some numerical experiments to show the efficiency of the proposed method.