A higher order numerical scheme for singularly perturbed parabolic turning point problems exhibiting twin boundary layers

In this article, a parameter-uniform numerical method is presented to solve one-dimensional singularly perturbed parabolic convection-diffusion multiple turning point problems exhibiting two exponential boundary layers. We study the asymptotic behaviour of the solution and its partial derivatives. The problem is discretized using the implicit Euler method for time discretization on a uniform mesh and a hybrid scheme for spatial discretization on a generalized Shishkin mesh. The scheme is shown to be ε-uniformly convergent of order one in time direction and order two in spatial direction upto a logarithmic factor. Numerical experiments are conducted to validate the theoretical results. Comparison is done with the upwind scheme on a uniform mesh as well as on the standard Shishkin mesh to demonstrate the higher order accuracy of the proposed scheme on a generalized Shishkin mesh.