Uniformly convergent numerical method for solving modified Burgers’ equations on a non-uniform mesh

Abstract

In this paper, we consider the one-dimensional modified Burgers’ equation in the finite domain. This type of problem arises in the field of sonic boom and explosions theory. At the high Reynolds’ number there is a boundary layer in the right side of the domain. From the numerical point of view, one of the difficulties in dealing with this problem is that even smooth initial data can give rise to solution varying regions, i.e., boundary layer regions. To tackle this situation, we propose a numerical method on non-uniform mesh of Shishkin type, which works well at high as well as low Reynolds number. The proposed numerical method comprises of Euler implicit and upwind f inite difference scheme. First we discretize in the temporal direction by means of Euler implicit method which yields the set of ordinary differential equations at each time level. The resulting set of differential equations are approximated by upwind scheme on Shishkin mesh. The proposed method has been shown to be parameter uniform and of almost first order accurate in the space and time. An extensive amount of analysis has been carried out in order to prove parameter uniform convergence of the method. some test examples have been solved to verify the theoretical results.