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

A splitting uniformly convergent method for one-dimensional parabolic singularly perturbed convection-diffusion systems

In this paper we deal with solving robustly and efficiently one-dimensional linear parabolic singularly perturbed systems of convection-diffusion type, where the diffusion parameters can be different at each equation and even they can have different orders of magnitude. The numerical algorithm combines the classical upwind finite difference scheme to discretize in space and the fractional implicit Euler method together with an appropriate splitting by components to discretize in time. We prove that if the spatial discretization is defined on an adequate piecewise uniform Shishkin mesh, the fully discrete scheme is uniformly convergent of first order in time and of almost first order in space. The technique used to discretize in time produces only tridiagonal linear systems to be solved at each time level; thus, from the computational cost point of view, the method we propose is more efficient than other numerical algorithms which have been used for these problems. Numerical results for several test problems are shown, which corroborate in practice both the uniform convergence and the efficiency of the algorithm

An almost third order finite difference scheme for singularly perturbed reaction–diffusion systems

AbstractThis paper addresses the numerical approximation of solutions to coupled systems of singularly perturbed reaction–diffusion problems. In particular a hybrid finite difference scheme of HODIE type is constructed on a piecewise uniform Shishkin mesh. It is proved that the numerical scheme satisfies a discrete maximum principle and also that it is third order (except for a logarithmic factor) uniformly convergent, even for the case in which the diffusion parameter associated with each equation of the system has a different order of magnitude. Numerical examples supporting the theory are given

An asymptotic-numerical hybrid method for singularly perturbed system of two-point reaction-diffusion boundary-value problems

This article focuses on the numerical approximate solution of singularly perturbed systems of secondorder reaction-diffusion two-point boundary-value problems for ordinary differential equations. To handle these types of problems, a numerical-asymptotic hybrid method has been used. In this hybrid approach, an efficient asymptotic method, the so-called successive complementary expansion method (SCEM) is employed first, and then a numerical method based on finite differences is applied to approximate the solution of corresponding singularly perturbed reactiondiffusion systems. Two illustrative examples are provided to demonstrate the efficiency, robustness, and easy applicability of the present method with convergence propertiesNo sponso