Accelerated fitted operator finite difference method for singularly perturbed parabolic reaction-diffusion problems

Document Type : Research Paper

Authors

1 Department of Mathematics, College of Natural science, Jimma University, Ethiopia.

2 Institut De Mathematiques et de sciences physiques, Universit D’Abomey Calavi, Benin.

Abstract

This paper deals with the numerical treatment of singularly perturbed parabolic reaction-diffusion initial boundary value problems. Introducing a fitting parameter into the asymptotic solution and applying average finite difference approximation, a fitted operator finite difference method is developed for solving the problem. To accelerate the rate of convergence of the method, Richardson extrapolation technique is applied. The consistency and stability of the proposed method have been established very well to ensure the convergence of the method. Numerical experimentation is carried out on some model problems and both the results are presented in tables and graphs. The numerical results are compared with findings of some methods existing in the literature and found to be more accurate. Generally, the formulated method is consistent, stable, and more accurate than some methods existing in the literature for solving singularly perturbed parabolic reaction-diffusion initial boundary value problems.

Keywords