A computational approach for 2D elliptic singularly perturbed weakly coupled systems of convection–diffusion type with multiple scales and parameters in the diffusion and the convection terms

Abstract

In this work, we consider the efficient resolution of a 2D elliptic singularly perturbed weakly coupled system of convection–diffusion type, which has small parameters at both the diffusion and the convection terms. We assume that the diffusion parameters can be distinct at each equation of the system, and also, they can have different orders of magnitude, but the convection parameter is the same at both equations of the system. Then, in general, overlapping regular or parabolic boundary layers can appear in the exact solution. The continuous problem is approximated by using the standard upwind finite difference scheme, which is constructed on a special piecewise uniform Shishkin mesh. Then, the numerical scheme is an almost first‐order uniformly convergent method with respect to all the perturbation parameters. Some numerical results, obtained with the numerical algorithm for one test problem, are presented, which show the good performance of the proposed numerical method and also corroborate in practice the theoretical results.