A parameter-robust numerical method for a system of reaction-diffusion equations in two dimensions

摘要

A system of M(>= 2) coupled singularly perturbed linear reaction-diffusion equations is considered on the unit square. Under certain hypotheses on the coupling, a maximum principle is established for the differential operator. The relationship between compatibility conditions at the corners of the square and the smoothness of the solution on the closed domain is fully described. A decomposition of the solution of the system is constructed. A finite-difference method for the solution of the system on a Shishkin mesh is presented, and it is proved that the computed solution is second-order accurate (up to a logarithmic factor). Numerical results are given to support this result and to investigate the effect of weaker compatibility assumptions on the data. (c) 2007 Wiley Periodicals. Inc.