Analysis of stationary iterative methods for the discrete convection-diffusion equation with a 9-point compact scheme

We study the convergence of point and line stationary iterative methods for solving the linear system arising from a fourth-order 9-point compact finite difference discretization of the two-dimensional convection-diffusion equation with constant coefficients. We present new techniques to bound the spectral radii of iteration matrices in terms of the cell Reynolds numbers. We also derive analytic formulas for the spectral radii for special values of the cell Reynolds numbers and study asymptotic behaviors of the analytic bounds. The results provide rigorous justification for the numerical experiments conducted elsewhere, which show good stability for the fourth-order compact scheme. In addition, we compare the 9-point scheme with the traditional 5-point difference discretization schemes and conduct some numerical experiments to supplement our analyses.

[1]  Murli M. Gupta,et al.  A single cell high order scheme for the convection‐diffusion equation with variable coefficients , 1984 .

[2]  Jun Zhang.  Preconditioned iterative methods and finite difference schemes for convection-diffusion , 2000, Appl. Math. Comput..

[3]  S. Dennis,et al.  Compact h4 finite-difference approximations to operators of Navier-Stokes type , 1989 .

[4]  Thomas A. Manteuffel,et al.  An analysis of block successive overrelaxation for a class of matrices with complex spectra , 1988 .

[5]  Jun Zhang.  On convergence and performance of iterative methods with fourth‐order compact schemes , 1998 .

[6]  H. F. Jordan,et al.  Is SOR Color-Blind? , 1986 .

[7]  David J. Evans,et al.  Block iterative methods for the nine-point approximation to the convection-diffusion equation , 1996, Int. J. Comput. Math..

[8]  Murli M. Gupta,et al.  A Compact Multigrid Solver for Convection-Diffusion Equations , 1997 .

[9]  W. F. SPOTZ,et al.  Iterative and Parallel Performance of High-Order Compact Systems , 1998, SIAM J. Sci. Comput..

[10]  J. Gillis,et al.  Matrix Iterative Analysis , 1961 .

[11]  A. Segal.  Aspects of Numerical Methods for Elliptic Singular Perturbation Problems , 1982 .

[12]  G. Carey,et al.  High‐order compact scheme for the steady stream‐function vorticity equations , 1995 .

[13]  P. Wesseling.  An Introduction to Multigrid Methods , 1992 .

[14]  Jun Zhang,et al.  An explicit fourth‐order compact finite difference scheme for three‐dimensional convection–diffusion equation , 1998 .

[15]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[16]  T. A. Manteuffel.  Optimal parameters for linear second-degree stationary iterative methods , 1982 .

[17]  David M. Young,et al.  Applied Iterative Methods , 2004 .

[18]  G. Golub,et al.  ITERATIVE METHODS FOR CYCLICALLY REDUCED NON-SELF-ADJOINT LINEAR SYSTEMS , 1990 .

[19]  B. Fornberg,et al.  A compact fourth‐order finite difference scheme for the steady incompressible Navier‐Stokes equations , 1995 .

[20]  Chen Greif,et al.  Iterative Solution of Cyclically Reduced Systems Arising from Discretization of the Three-Dimensional Convection-Diffusion Equation , 1998, SIAM J. Sci. Comput..

[21]  Jun Zhang,et al.  Accuracy, robustness and efficiency comparison in iterative computation of convection diffusion equation with boundary layers , 2000 .

[22]  Irfan Altas,et al.  A high accuracy defect-correction multigrid method for the steady incompressible Navier-Stokes equations , 1994 .

[23]  Jun Zhang,et al.  On cyclic reduction and finite difference schemes , 2002 .

[24]  Jun Zhang.  Convergence of iterative methods for a fourth-order discretization scheme , 1997 .