Error expansion for an upwind scheme applied to a two-dimensional convection-diffusion problem

Research output: Contribution to journalArticlepeer-review

Abstract

We consider a singularly perturbed convection-diffusion problem in a rectangular domain. It is solved numerically using a first-order upwind finite-difference scheme on a tensor-product piecewise-uniform Shishkin mesh with O(N) mesh points in each coordinate direction. It is known [G. I. Shishkin, Grid Approximations of Singularly Perturbed Elliptic and Parabolic Equations, Russian Academy of Sciences, Ural Branch, Ekaterinburg, Russia, 1992 (in Russian)] that the error is almost-first-order accurate in the maximum norm. We decompose the error into a sum of continuous almost-first-order terms and the almost-second-order residual under the assumption ε < CN-1, where ε is the singular perturbation parameter and C is a constant. This error expansion is applied to obtain maximum-norm error estimates for the Richardson extrapolation technique and derive bounds on the errors in approximating the derivatives of the true solution by divided differences of the computed solution. The analysis uses a decomposition of the true solution requiring fewer compatibility conditions than in earlier publications. Numerical results are presented that support our theoretical results.

Original languageEnglish
Pages (from-to)1851-1869
Number of pages19
JournalSIAM Journal on Numerical Analysis
Volume41
Issue number5
DOIs
Publication statusPublished - 2003
Externally publishedYes

Keywords

  • Approximation of derivatives
  • Convection-diffusion
  • Error expansion
  • Richardson extrapolation
  • Shishkin mesh
  • Singular perturbation
  • Upwind scheme

Fingerprint

Dive into the research topics of 'Error expansion for an upwind scheme applied to a two-dimensional convection-diffusion problem'. Together they form a unique fingerprint.

Cite this