Fourier analysis of multigrid for a model two-dimensional convection-diffusion equation

Howard C. Elman, Alison Ramage

Research output: Contribution to journalArticle

3 Citations (Scopus)

Abstract

We present a Fourier analysis of multigrid for the two-dimensional discrete convection-diffusion equation. For constant coefficient problems with grid-aligned flow and semi-periodic boundary conditions, we show that the two-grid iteration matrix can be reduced via a set of orthogonal transformations to a matrix containing individual 4×4 blocks. This enables a trivial computation of the norm of the iteration matrix demonstrating rapid convergence in the case of both small and large mesh Peclet numbers, where the streamline-diffusion discretisation is used in the latter case. We also demonstrate that these results are strongly correlated with the properties of the iteration matrix arising from Dirichlet boundary conditions.

Original languageEnglish
Pages (from-to)283-306
Number of pages24
JournalBIT Numerical Mathematics
Volume46
Issue number2
DOIs
Publication statusPublished - 1 Jun 2006

Fingerprint

Fourier analysis
Convection-diffusion Equation
Fourier Analysis
Iteration
Boundary conditions
Streamline Diffusion
Grid
Orthogonal Transformation
Peclet number
Discrete Equations
Periodic Boundary Conditions
Model
Dirichlet Boundary Conditions
Trivial
Discretization
Mesh
Norm
Convection
Coefficient
Demonstrate

Keywords

  • convection-diffusion
  • convergence
  • multigrid

Cite this

@article{4a9e5299aa1941d8a3f132e731236540,
title = "Fourier analysis of multigrid for a model two-dimensional convection-diffusion equation",
abstract = "We present a Fourier analysis of multigrid for the two-dimensional discrete convection-diffusion equation. For constant coefficient problems with grid-aligned flow and semi-periodic boundary conditions, we show that the two-grid iteration matrix can be reduced via a set of orthogonal transformations to a matrix containing individual 4×4 blocks. This enables a trivial computation of the norm of the iteration matrix demonstrating rapid convergence in the case of both small and large mesh Peclet numbers, where the streamline-diffusion discretisation is used in the latter case. We also demonstrate that these results are strongly correlated with the properties of the iteration matrix arising from Dirichlet boundary conditions.",
keywords = "convection-diffusion, convergence, multigrid",
author = "Elman, {Howard C.} and Alison Ramage",
year = "2006",
month = "6",
day = "1",
doi = "10.1007/s10543-006-0052-5",
language = "English",
volume = "46",
pages = "283--306",
journal = "BIT Numerical Mathematics",
issn = "0006-3835",
number = "2",

}

Fourier analysis of multigrid for a model two-dimensional convection-diffusion equation. / Elman, Howard C.; Ramage, Alison.

In: BIT Numerical Mathematics, Vol. 46, No. 2, 01.06.2006, p. 283-306.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Fourier analysis of multigrid for a model two-dimensional convection-diffusion equation

AU - Elman, Howard C.

AU - Ramage, Alison

PY - 2006/6/1

Y1 - 2006/6/1

N2 - We present a Fourier analysis of multigrid for the two-dimensional discrete convection-diffusion equation. For constant coefficient problems with grid-aligned flow and semi-periodic boundary conditions, we show that the two-grid iteration matrix can be reduced via a set of orthogonal transformations to a matrix containing individual 4×4 blocks. This enables a trivial computation of the norm of the iteration matrix demonstrating rapid convergence in the case of both small and large mesh Peclet numbers, where the streamline-diffusion discretisation is used in the latter case. We also demonstrate that these results are strongly correlated with the properties of the iteration matrix arising from Dirichlet boundary conditions.

AB - We present a Fourier analysis of multigrid for the two-dimensional discrete convection-diffusion equation. For constant coefficient problems with grid-aligned flow and semi-periodic boundary conditions, we show that the two-grid iteration matrix can be reduced via a set of orthogonal transformations to a matrix containing individual 4×4 blocks. This enables a trivial computation of the norm of the iteration matrix demonstrating rapid convergence in the case of both small and large mesh Peclet numbers, where the streamline-diffusion discretisation is used in the latter case. We also demonstrate that these results are strongly correlated with the properties of the iteration matrix arising from Dirichlet boundary conditions.

KW - convection-diffusion

KW - convergence

KW - multigrid

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

U2 - 10.1007/s10543-006-0052-5

DO - 10.1007/s10543-006-0052-5

M3 - Article

AN - SCOPUS:33745297120

VL - 46

SP - 283

EP - 306

JO - BIT Numerical Mathematics

JF - BIT Numerical Mathematics

SN - 0006-3835

IS - 2

ER -