Srinivasan Natesan1, S. Gowrisankar2
CMES-Computer Modeling in Engineering & Sciences, Vol.88, No.4, pp. 245-268, 2012, DOI:10.3970/cmes.2012.088.245
Abstract In this article, we propose a parameter-uniform computational technique to solve singularly perturbed parabolic initial-boundary-value problems exhibiting parabolic layers. The domain is discretized with a uniform mesh on the time direction and a nonuniform mesh obtained via equidistribution of a monitor function for the spatial variable. The numerical scheme consists of the implicit-Euler scheme for the time derivative and the classical central difference scheme for the spatial derivative. Truncation error, and stability analysis are carried out. Error estimates are derived, and numerical examples are presented. More >