An algebraic multilevel parallelizable preconditioner for large-scale CFD problems

Abstract An efficient parallelizable preconditioner for solving large-scale CFD problems is presented. It is adapted to coarse-grain parallelism and can be used for both shared and distributed-memory parallel computers. The proposed preconditioner consists of two independent approximations of the system matrix. The first one is a block-diagonal, fully parallelizable approximation of the given system. The second matrix is coarser than the original one and is built using algebraic multi-grid methods. The preconditioner is used to compute the steady solution of the compressible Navier-Stokes equations for subsonic laminar flows, on shared-memory computers, for a moderate number of processors. The coupled two-level preconditioner is robust and has a large potential for parallelism. Interesting savings in computational time for parallel computations are obtained when comparing the two-level preconditioner with the well-known incomplete Gaussian factorization.

[1]  E. Cuthill,et al.  Reducing the bandwidth of sparse symmetric matrices , 1969, ACM '69.

[2]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[3]  Michel Fortin,et al.  Parallelizable block diagonal preconditioners for the compressible Navier-Stokes equations , 1994 .

[4]  P. Sonneveld.  CGS, A Fast Lanczos-Type Solver for Nonsymmetric Linear systems , 1989 .

[5]  Alan George,et al.  The Evolution of the Minimum Degree Ordering Algorithm , 1989, SIAM Rev..

[6]  Bruce Hendrickson,et al.  A Multi-Level Algorithm For Partitioning Graphs , 1995, Proceedings of the IEEE/ACM SC95 Conference.

[7]  A method for finite element parallel viscous compressible flow calculations , 1994 .

[8]  J. A. George.  Computer implementation of the finite element method , 1971 .

[9]  A. H. Sherman,et al.  Comparative Analysis of the Cuthill–McKee and the Reverse Cuthill–McKee Ordering Algorithms for Sparse Matrices , 1976 .

[10]  Daniel B. Szyld,et al.  A Block Ordering Method for Sparse Matrices , 1990, SIAM J. Sci. Comput..

[11]  Roland W. Freund,et al.  A Transpose-Free Quasi-Minimal Residual Algorithm for Non-Hermitian Linear Systems , 1993, SIAM J. Sci. Comput..

[12]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[13]  L. Dutto.  The effect of ordering on preconditioned GMRES algorithm, for solving the compressible Navier-Stokes equations , 1993 .

[14]  Bruce Hendrickson,et al.  The Chaco user`s guide. Version 1.0 , 1993 .