This paper provides an efficient numerical scheme for approximating the solution of the time-fractional parabolic equation on 2D curved domain. Here, the solution of the problem exhibits a weak singularity at the initial time \(t=0\) . The method is based on applying the L1 formula to approximate the Caputo time-fractional derivative and using the isoparametric finite element method to approximate the spatial direction. A fully discrete scheme is then constructed using numerical quadrature. Since \(\Omega _{h}\) differs from \(\Omega \) , the error estimates of the boundary terms on the region between \(\Omega _{h}\) and \(\Omega \) and the effect of numerical quadrature are both considered. Afterward, the stability analysis and optimal error estimates for the fully discrete scheme are proved in detail. Finally, some numerical experiments are presented to verify the theoretical results.
