A numerical method for a system of singularly perturbed reaction-diffusion equations (Q1612394)

It is well-known that standard numerical methods are not appropriate for singularly perturbed problems. The authors concentrate on the case of a Dirichlet problem for a system of two singularly perturbed reaction-diffusion problems, with diffusion coefficients \(\varepsilon_1\) and \(\varepsilon_2\) such that \(\varepsilon_1=\varepsilon\) and \(\varepsilon_2=1\). A numerical parameter-uniform method whose solutions converge pointwise at all points of the domain independently of the singular perturbation parameters is constructed and analysed. The authors establish essentially second-order parameter-uniform convergence. Numerical results are presented, which illustrate the theoretical results. This paper extends the results for piecewise-uniform meshes applied to a single one-dimensional singularly perturbed reaction-diffusion problem [\textit{J. J. H. Miller, E. O'Riordan} and \textit{G. I. Shishkin}, Fitted mesh methods for the singular perturbation reaction diffusion problem, in: Proceedings of the Fifth International Colloquium on Numerical Analysis, 1996, Plovdiv, Bulgaria (Academic Publications, 1997), pp. 99-105] to a system of two coupled reaction-diffusion problems.