A parameter-uniform second order numerical method for a weakly coupled system of singularly perturbed convection–diffusion equations with discontinuous convection coefficients and source terms

Abstract

In this article, a parameter-uniform hybrid numerical method is presented to solve a weakly coupled system of two singularly perturbed convection–diffusion equations with discontinuous convection coefficients and source terms. Due to these discontinuities, interior layers appear in the solution of the problem considered. The hybrid numerical method uses the standard finite difference scheme in the coarse mesh region and the cubic spline difference scheme in the fine mesh region which is constructed on piecewise-uniform Shishkin mesh. Second order one sided difference approximations are used at the point of discontinuity. Error analysis is carried out and the method ensures that the parameter-uniform convergence of almost the second order. Numerical results are provided to validate the theoretical results.