A uniformly convergent scheme to solve two-dimensional parabolic singularly perturbed systems of reaction-diffusion type with multiple diffusion parameters

Abstract

In this work, we deal with solving two-dimensional parabolic singularly perturbed systems of reaction-diffusion type where the diffusion parameters at each equation of the system can be small and of different scale. In such case, in general, overlapping boundary layers appear at the boundary of the spatial domain and, because of this, special meshes are required to resolve them. The numerical scheme combines the central difference scheme to discretize in space and the fractional implicit Euler method together with a splitting by components to discretize in time. If the fully discrete scheme is defined on an adequate piecewise uniform Shishkin mesh in space then it is uniformly convergent of first order in time and of almost second order in space. Some numerical results illustrate the theoretical results. © 2020 John Wiley & Sons, Ltd.