Higher order global solution and normalized flux for singularly perturbed reaction-diffusion problems |
| |
Authors: | R.K. Bawa C. Clavero |
| |
Affiliation: | Department of Computer Science, Punjabi University, Patiala 147 002, India CPS, Department of Applied Mathematics, University of Zaragoza, 50018 Zaragoza, Spain |
| |
Abstract: | In this paper a computational technique is proposed for obtaining a higher order global solution and global normalized flux of singularly perturbed reaction-diffusion two-point boundary-value problems. The HOC (higher order compact) finite difference scheme developed in Gracia et al. (2001) [4] and which is constructed on an appropriate piecewise uniform Shishkin mesh, has been considered to find an almost fourth order convergent solution at mesh points. Using these values, piecewise cubic interpolants based approximations for solution and normalized flux in whole domain have been defined. It has been shown that the global solution and the global normalized flux are also uniformly convergent. Moreover, for the global solution, the order of uniform convergence in the whole domain is optimal, i.e., it is the same as this one obtained at mesh points, whereas, for the global normalized flux, the uniform convergence is almost third order, except at midpoints of the mesh, where it is also almost fourth order. Theoretical error bounds have been provided along with some numerical examples, which corroborate the efficiency of the proposed technique to find good approximations to the global solution and the global normalized flux. |
| |
Keywords: | Singular perturbation problems Reaction-diffusion HOC scheme Shishkin mesh Global solution Global normalized flux |
本文献已被 ScienceDirect 等数据库收录! |
|