Higher order global solution and normalized flux for singularly perturbed reaction-diffusion problems (Q979301)

The authors consider a one dimensional singularly perturbed reaction-diffusion equation. A modified Shishkin mesh is introduced and a higher order compact finite difference solution on this mesh is presented. Piecewise cubic interpolants for both the exact solution and discrete solution are formulated. Thanks to the modfied Shishkin mesh, the authors prove that the convergence is uniform in the sense that the convergence accuracy is the same for any value of the diffusion parameter \(\varepsilon\). More precise, the convergence order analysis contains two principal results. The first result states that the method is almost of fourth convergence order. The second result states that the normalized flux of the piecewise cubic interpolant of the discrete solution approximates the normalized flux of the piecewise cubic interpolant of the exact solution by order three, almost everywhere, and by order four at the mid-points of the mesh. The theoretical results are confirmed by numerical examples.