A uniformly convergent scheme on a nonuniform mesh for convection-diffusion parabolic problems

In this paper we construct a numerical method to solve one-dimensional time-dependent convection-diffusion problem with dominating convection term. We use the classical Euler implicit method for the time discretization and the simple upwind scheme on a special nonuniform mesh for the spatial discretization. We show that the resulting method is uniformly convergent with respect to the diffusion parameter. The main lines for the analysis of the uniform convergence carried out here can be used for the study of more general singular perturbation problems and also for more complicated numerical schemes. The numerical results show that, in practice, some of the theoretical compatibility conditions seem not necessary.

[1]  L. P. Shishkina,et al.  FITTED MESH METHODS FOR PROBLEMS WITH PARABOLIC BOUNDARY LAYERS , 1998 .

[2]  Martin Stynes,et al.  Numerical methods for time-dependent convection-diffusion equations , 1988 .

[3]  R. Kellogg,et al.  Analysis of some difference approximations for a singular perturbation problem without turning points , 1978 .

[4]  G. I. SHISHKIN,et al.  Grid approximation of singularly perturbed boundary value problems with convective terms , 1990 .

[5]  J. C. Jorge,et al.  A fractional step method on a special mesh for the resolution multidimensional evolutionary convection-diffusion problems , 1998 .

[6]  Lutz Tobiska,et al.  Numerical Methods for Singularly Perturbed Differential Equations , 1996 .

[7]  Martin Stynes,et al.  The midpoint upwind scheme , 1997 .

[8]  J. J. H Miller,et al.  Applications of advanced computational methods for boundary and interior layers , 1993 .

[9]  M. Stynes,et al.  Finite-element methods for singularly perturbed high-order elliptic two-point boundary value problems. II: convection-diffusion-type problems , 1995 .

[10]  Piet Hemker,et al.  ε-uniform schemes with high-order time-accuracy for parabolic singular perturbation problems , 2000 .

[11]  John J. H. Miller Fitted Numerical Methods for Singular Perturbation Problems: Error Estimates in the Maximum Norm for Linear Problems in One and Two Dimensions , 1996 .

[12]  Eugene O'Riordan,et al.  Special meshes for finite difference approximations to an advection-diffusion equation with parabolic layers , 1995 .

[13]  JohnM . Miller,et al.  Robust Computational Techniques for Boundary Layers , 2000 .

[14]  J. C. Jorge,et al.  An alternating direction scheme on a nonuniform mesh for reaction-diffusion parabolic problems , 2000 .