Parameter uniform essentially first order convergence of a fitted mesh method for a class of parabolic singularly perturbed Robin problem for a system of reaction-diffusion equations

In this paper, a class of linear parabolic systems of singularly perturbed second order differential equations of reaction-diffusion type with initial and Robin boundary conditions is considered. The components of the solution $\vec u$ of this system exhibit parabolic boundary layers with sublayers. A numerical method composed of a classical finite difference scheme on a piecewise uniform Shishkin mesh is suggested. This method is proved to be first order convergent in time and essentially first order convergent in the space variable in the maximum norm uniformly in the perturbation parameters