In this paper, a uniformly convergent numerical scheme is designed for solving singularly perturbed reaction-diffusion problems. The problem is converted to an equivalent weak form and then a Galerkin finite element method is used on a piecewise uniform Shishkin mesh with linear basis functions. The convergence of the developed scheme is proved and it is shown to be almost fourth order uniformly convergent in the maximum norm. To exhibit the applicability of the scheme, model examples are considered and solved for different values of a singular perturbation parameter ε and mesh elements. The proposed scheme approximates the exact solution very well.

Keywords: Finite Element, Fitted Mesh, Parameter Uniform, Singularly Perturbed.