A parameter robust numerical method for a two dimensional reaction-diffusion problem.

In this paper a singularly perturbed reaction-diffusion partial differential equation in two space dimensions is examined. By means of an appropriate decomposition, we describe the asymptotic behaviour of the solution of problems of this kind. A central finite difference scheme is constructed for this problem which involves an appropriate Shishkin mesh. We prove that the numerical approximations are almost second order uniformly convergent (in the maximum norm) with respect to the singular perturbation parameter. Some numerical experiments are given that illustrate in practice the theoretical order of convergence established for the numerical method

The Construction of Finite Difference Approximations to Ordinary Differential Equations

Finite difference approximations of the form Σ^(si)_(i=-rj)d_(j,i)u_(j+i)=Σ^(mj)_(i=1) e_(j,if)(z_(j,i)) for the numerical solution of linear nth order ordinary differential equations are analyzed. The order of these approximations is shown to be at least r_j + s_j + m_j - n, and higher for certain special choices of the points Z_(j,i). Similar approximations to initial or boundary conditions are also considered and the stability of the resulting schemes is investigated

A parameter uniform fitted mesh method for a weakly coupled system of two singularly perturbed convection-diffusion equations

In this paper, a boundary value problem for a singularly perturbed linear system of two second order ordinary differential equations of convection- diffusion type is considered on the interval [0, 1]. The components of the solution of this system exhibit boundary layers at 0. A numerical method composed of an upwind finite difference scheme applied on a piecewise uniform Shishkin mesh is suggested to solve the problem. The method is proved to be first order convergent in the maximum norm uniformly in the perturbation parameters. Numerical examples are provided in support of the theory