Abstract:
This work aims to construct an efficient and highly accurate numerical method to address the time singularity at involved in a class of time-fractional parabolic integro-partial differential equations in one and two dimensions. The L2- scheme is used to discretize the time-fractional operator, whereas a modified version of the composite trapezoidal approximation is employed to discretize the Volterra operator in time. Subsequently, it helps to convert the proposed model into a second-order boundary value problem in a semi-discrete form. The multi-dimensional Haar wavelets are then used for grid adaptation and efficient computations for the two-dimensional problem, whereas the standard second-order approximations are employed to approximate the spatial derivatives for the one-dimensional case. The stability analysis is carried out on an adaptive mesh in time. The convergence analysis leads to accurate solution in the space-time domain for the one-dimensional problem having time singularity based on the norm for a suitable choice of the grading parameter. Furthermore, it provides accurate solution for the two-dimensional problem having unbounded time derivative at . The analysis also highlights a higher order accuracy for a sufficiently smooth solution resides in even if the mesh is discretized uniformly. The truncation error estimates for the time-fractional operator, integral operator, and spatial derivatives are presented. In addition, we have examined the impact of various parameters on the robustness and accuracy of the proposed method. Numerous tests are performed on several examples in support of the theoretical analysis. The advancement of the proposed methodology is demonstrated through the application of the time-fractional Fokker-Planck equation and the fractional-order viscoelastic dynamics having weakly singular kernels. It also confirms the superiority of the proposed method compared with existing approaches available in the literature.