A numerical approach for a two-parameter singularly perturbed weakly-coupled system of 2-D elliptic convection–reaction–diffusion PDEs

In this work, we consider the numerical approximation of a two dimensional elliptic singularly perturbed weakly-coupled system of convection–reaction–diffusion type, which has two different parameters affecting the diffusion and the convection terms, respectively. The solution of such problems has, in general, exponential boundary layers as well as corner layers. To solve the continuous problem, we construct a numerical method which uses a finite difference scheme defined on an appropriate layer-adapted Bakhvalov–Shishkin mesh. Then, the numerical scheme is a first order uniformly convergent method with respect both convection and diffusion parameters. Numerical results obtained with the algorithm for some test problems are presented, which show the best performance of the proposed method, and they also corroborate in practice the theoretical analysis.