A characterisation of oscillations in the discrete two-dimensional convection-diffusion equation

Howard C. Elman, Alison Ramage

Research output: Contribution to journalArticle

12 Citations (Scopus)

Abstract

It is well known that discrete solutions to the convection-diffusion equation contain nonphysical oscillations when boundary layers are present but not resolved by the discretisation. However, except for one-dimensional problems, there is little analysis of this phenomenon. In this paper, we present an analysis of the two-dimensional problem with constant flow aligned with the grid, based on a Fourier decomposition of the discrete solution. For Galerkin bilinear finite element discretisations, we derive closed form expressions for the Fourier coefficients, showing them to be weighted sums of certain functions which are oscillatory when the mesh Péclet number is large. The oscillatory functions are determined as solutions to a set of three-term recurrences, and the weights are determined by the boundary conditions. These expressions are then used to characterise the oscillations of the discrete solution in terms of the mesh Péclet number and boundary conditions of the problem.

Original languageEnglish
Pages (from-to)263-288
Number of pages26
JournalMathematics of Computation
Volume72
Issue number241
Early online date4 Dec 2001
DOIs
Publication statusPublished - 1 Jan 2003

Fingerprint

Convection-diffusion Equation
Oscillation
Mesh
Boundary conditions
Finite Element Discretization
Fourier coefficients
Weighted Sums
Galerkin
Recurrence
Boundary Layer
Closed-form
Discretization
Grid
Decompose
Boundary layers
Term
Decomposition
Convection

Keywords

  • computation mathematics
  • convection-diffusion equations
  • oscillations

Cite this

@article{9d316f3e0f1c401e9aceaa1c2bfa33e6,
title = "A characterisation of oscillations in the discrete two-dimensional convection-diffusion equation",
abstract = "It is well known that discrete solutions to the convection-diffusion equation contain nonphysical oscillations when boundary layers are present but not resolved by the discretisation. However, except for one-dimensional problems, there is little analysis of this phenomenon. In this paper, we present an analysis of the two-dimensional problem with constant flow aligned with the grid, based on a Fourier decomposition of the discrete solution. For Galerkin bilinear finite element discretisations, we derive closed form expressions for the Fourier coefficients, showing them to be weighted sums of certain functions which are oscillatory when the mesh P{\'e}clet number is large. The oscillatory functions are determined as solutions to a set of three-term recurrences, and the weights are determined by the boundary conditions. These expressions are then used to characterise the oscillations of the discrete solution in terms of the mesh P{\'e}clet number and boundary conditions of the problem.",
keywords = "computation mathematics, convection-diffusion equations, oscillations",
author = "Elman, {Howard C.} and Alison Ramage",
year = "2003",
month = "1",
day = "1",
doi = "10.1090/S0025-5718-01-01392-8",
language = "English",
volume = "72",
pages = "263--288",
journal = "Mathematics of Computation",
issn = "0025-5718",
number = "241",

}

A characterisation of oscillations in the discrete two-dimensional convection-diffusion equation. / Elman, Howard C.; Ramage, Alison.

In: Mathematics of Computation, Vol. 72, No. 241, 01.01.2003, p. 263-288.

Research output: Contribution to journalArticle

TY - JOUR

T1 - A characterisation of oscillations in the discrete two-dimensional convection-diffusion equation

AU - Elman, Howard C.

AU - Ramage, Alison

PY - 2003/1/1

Y1 - 2003/1/1

N2 - It is well known that discrete solutions to the convection-diffusion equation contain nonphysical oscillations when boundary layers are present but not resolved by the discretisation. However, except for one-dimensional problems, there is little analysis of this phenomenon. In this paper, we present an analysis of the two-dimensional problem with constant flow aligned with the grid, based on a Fourier decomposition of the discrete solution. For Galerkin bilinear finite element discretisations, we derive closed form expressions for the Fourier coefficients, showing them to be weighted sums of certain functions which are oscillatory when the mesh Péclet number is large. The oscillatory functions are determined as solutions to a set of three-term recurrences, and the weights are determined by the boundary conditions. These expressions are then used to characterise the oscillations of the discrete solution in terms of the mesh Péclet number and boundary conditions of the problem.

AB - It is well known that discrete solutions to the convection-diffusion equation contain nonphysical oscillations when boundary layers are present but not resolved by the discretisation. However, except for one-dimensional problems, there is little analysis of this phenomenon. In this paper, we present an analysis of the two-dimensional problem with constant flow aligned with the grid, based on a Fourier decomposition of the discrete solution. For Galerkin bilinear finite element discretisations, we derive closed form expressions for the Fourier coefficients, showing them to be weighted sums of certain functions which are oscillatory when the mesh Péclet number is large. The oscillatory functions are determined as solutions to a set of three-term recurrences, and the weights are determined by the boundary conditions. These expressions are then used to characterise the oscillations of the discrete solution in terms of the mesh Péclet number and boundary conditions of the problem.

KW - computation mathematics

KW - convection-diffusion equations

KW - oscillations

UR - http://www.scopus.com/inward/record.url?scp=85009726896&partnerID=8YFLogxK

U2 - 10.1090/S0025-5718-01-01392-8

DO - 10.1090/S0025-5718-01-01392-8

M3 - Article

VL - 72

SP - 263

EP - 288

JO - Mathematics of Computation

JF - Mathematics of Computation

SN - 0025-5718

IS - 241

ER -