A high order convergent numerical method for singularly perturbed time dependent problems using mesh equidistribution

Abstract

The purpose of this paper is to introduce a high order convergent numerical method for singularly perturbed time dependent problems using mesh equidistribution. The discretization is based on the backward Euler scheme in time and a high order non-monotone scheme in space. In time direction we consider a uniform mesh, while in spatial direction we construct an adaptive mesh through equidistribution of a monitor function involving appropriate power of the solution’s second derivative. The method is analysed in two steps, splitting the time and space discretization errors. We establish that the method is uniformly convergent with optimal order having order one in time and order four in space. Further, we use the Richardson extrapolation technique for improving the order of convergence from one to two in time. Numerical experiments are presented to confirm the theoretically proven convergence result.