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 (Q1675441)

The authors propose a uniformly convergent numerical scheme for singularly perturbed systems of weakly coupled convection-diffusion two-point boundary value problems (BVPs) with discontinuous convection coefficients and source terms. The solution of this system of BVPs exhibit interior layers. To discretize the domain piecewise-uniform Shishkin meshes are used. The differential equations are approximated by the cubic spline scheme in the interior layer regions, where the meshes are fine, and by the classical finite difference scheme in the outer regions, where the meshes are coarse. At the discontinuity point a one-sided second-order difference scheme is used. The proposed hybrid scheme is almost second-order convergent with a logarithmic term. Theoretical error estimates are obtained, and numerical experiments are carried out.