This paper deals with numerical treatment of singularly perturbed parabolic differential equations having large time delay. The highest order derivative term in the equation is multiplied by a perturbation parameter ε , taking arbitrary value in the interval 0 , 1 . For small values of ε , solution of the problem exhibits an exponential boundary layer on the right side of the spatial domain. The properties and bounds of the solution and its derivatives are discussed. The considered singularly perturbed time delay problem is solved using the Crank-Nicolson method in temporal discretization and exponentially fitted operator finite difference method in spatial discretization. The stability of the scheme is investigated and analysed using comparison principle and solution bound. The uniform convergence of the scheme is discussed and proven. The formulated scheme converges uniformly with linear order of convergence. The theoretical analysis of the scheme is validated by considering numerical test examples for different values of ε .

This paper deals with solution methods for singularly perturbed delay differential equations having delay on the convection and reaction terms. The considered problem exhibits an exponential boundary layer on the left or right side of the domain. The terms with the delay are treated using Taylor?s series approximation and the resulting singularly perturbed boundary value problem is solved using a specially designed exponentially finite difference method. The stability of the scheme is analysed and investigated using a comparison principle and solution bound. The formulated scheme converges uniformly with linear order of convergence. The theoretical findings are validated using three numerical test examples.