High accuracy multigrid solution of the 3D convection-diffusion equation

Murli M. Gupta, Jun Zhang

Research output: Contribution to journalArticlepeer-review

85 Scopus citations


We present an explicit fourth-order compact finite difference scheme for approximating the three-dimensional (3D) convection-diffusion equation with variable coefficients. This 19-point formula is defined on a uniform cubic grid. Fourier smoothing analysis is performed to show that the smoothing factor of certain relaxation techniques used with the scheme is smaller than 1. We design a parallelization-oriented multigrid method for fast solution of the resulting linear system using a four-color Gauss-Seidel relaxation technique for robustness and efficiency, and a scaled residual injection operator to reduce the cost of multigrid inter-grid transfer operator. Numerical experiments on a 16 processor vector computer are used to test the high accuracy of the discretization scheme as well as the fast convergence and the parallelization or vectorization efficiency of the solution method. Several test problems are solved and highly accurate solutions of the 3D convection-diffusion equations are obtained for small to medium values of the grid Reynolds number. Effects of using different residual projection operators are compared on both vector and serial computers.

Original languageEnglish
Pages (from-to)249-274
Number of pages26
JournalApplied Mathematics and Computation
Issue number2
StatePublished - Jul 15 2000

Bibliographical note

Funding Information:
This research was partially supported by a grant (DMS970001P) from Pittsburg Supercomputing Center.

ASJC Scopus subject areas

  • Computational Mathematics
  • Applied Mathematics


Dive into the research topics of 'High accuracy multigrid solution of the 3D convection-diffusion equation'. Together they form a unique fingerprint.

Cite this