A fitted numerical method for singularly perturbed parabolic reaction-diffusion problems

This paper treats a time-dependent singularly perturbed reaction-diffusion problem. We semidiscretize the problem in time by means of the classical backward Euler method. We develop a fitted operator finite difference method (FOFDM) to solve the resulting set of linear problems (one at each time level). We prove that the underlying fitted operator satisfies the maximum principle. This result is then used in the error analysis of the FOFDM. The method is shown to be first order convergent in time and second order convergent in space, uniformly with respect to the perturbation parameter. We test the method on several numerical examples to confirm our theoretical findings.