Abstract
A fitted tension spline numerical scheme for a singularly perturbed parabolic problem (SPPP) with time delay is proposed. The presence of a small parameter ε as a multiple of the diffusion term leads to the suddenly changing behaviors of the solution in the boundary layer region. This results in a challenging duty to solve the problem analytically. Classical numerical methods cause spurious nonphysical oscillations unless an unacceptable number of mesh points is considered, which requires a large computational cost. To overcome this drawback, a numerical method comprising the backward Euler scheme in the time direction and the fitted spline scheme in the space direction on uniform meshes is proposed. To establish the stability and uniform convergence of the proposed method, an extensive amount of analysis is carried out. Three numerical examples are considered to validate the efficiency and applicability of the proposed scheme. It is proved that the proposed scheme is uniformly convergent of order one in both space and time. Further, the boundary layer behaviors of the solutions are given graphically.