A second order numerical method for a class of parameterized singular perturbation problems on adaptive grid

Abstract

A nonlinear singularly perturbed boundary value problem depending on a parameter is considered. First, we solve the problem using the backward Euler finite difference scheme on an adaptive grid. The adaptive grid is a special nonuniform mesh generated through equidistribution principle by a positive monitor function depending on the solution. The behavior of the solution, the stability and the error estimates are discussed. Then, the Richardson extrapolation technique is applied to improve the accuracy of the computed solution associated to the backward Euler scheme. The proofs of the uniform convergence for the backward Euler scheme and the Richardson extrapolation are carried out. Numerical experiments validate the theoretical estimates and indicates that the estimates are sharp.

Keywords

• singular perturbation
• layer adapted mesh
• parameterized problems
• boundary layer
• richardson extrapolation