A novel higher-order numerical method for parabolic integro-fractional differential equations based on wavelets and $L2$-$1_σ$ scheme

This study aims to construct an efficient and highly accurate numerical method to solve a class of parabolic integro-fractional differential equations, which is based on wavelets and $L2$-$1_\sigma$ scheme. Specifically, the Haar wavelet decomposition is used for grid adaptation and efficient computations, while the high order $L2$-$1_\sigma$ scheme is considered to discretize the time-fractional operator. Second-order discretizations are used to approximate the spatial derivatives to solve the one-dimensional problem, while a repeated quadrature rule based on trapezoidal approximation is employed to discretize the integral operator. In contrast, we use the semi-discretization of the proposed two-dimensional model based on the $L2$-$1_\sigma$ scheme for the fractional operator and composite trapezoidal approximation for the integral part. The spatial derivatives are then approximated using two-dimensional Haar wavelets. In this study, we investigated theoretically and verified numerically the behavior of the proposed higher-order numerical methods. In particular, stability and convergence analyses are conducted. The obtained results are compared with those of some existing techniques through several graphs and tables, and it is shown that the proposed higher-order methods have better accuracy and produce less error compared to the $L1$ scheme in favor of fractional-order integro-partial differential equations.