Laboratory for High Performance Scientific Computing and Numerical Simulation
Department of Computer Science
University of Kentucky
773 Anderson Hall
Lexington, KY 40506-0046, USA
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.
Technical Report No. 308-00, Department of Computer Science, University of Kentucky, Lexington, KY, 2000. The second author's research work was supported in part by the U.S. National Science Foundation under grants CCR-9902022, CCR-9988165 and CCR-0043861.