A taxonomy and comparison of parallel block multi-level preconditioners for the incompressible Navier-Stokes equations.

Shadid, John Nicolas (Sandia National Laboratories, Albuquerque, NM) ; Elman, Howard (University of Maryland, College Park, MD) ; Shuttleworth, Robert R. (University of Maryland, College Park, MD) ; Howle, Victoria E. ; Tuminaro, Raymond Stephen

Sandia National Laboratories

关键词: Navier-Stokes Equations;    Finite Element Method Finite Element Method.;    99 General And Miscellaneous//Mathematics, Computing, And Information Science;    Mesh Generation.;    Approximations;   

DOI  :  10.2172/920807
RP-ID  :  SAND2007-2761
RP-ID  :  AC04-94AL85000
RP-ID  :  920807

美国|英语

来源: UNT Digital Library

【 摘 要 】

In recent years, considerable effort has been placed on developing efficient and robust solution algorithms for the incompressible Navier-Stokes equations based on preconditioned Krylov methods. These include physics-based methods, such as SIMPLE, and purely algebraic preconditioners based on the approximation of the Schur complement. All these techniques can be represented as approximate block factorization (ABF) type preconditioners. The goal is to decompose the application of the preconditioner into simplified sub-systems in which scalable multi-level type solvers can be applied. In this paper we develop a taxonomy of these ideas based on an adaptation of a generalized approximate factorization of the Navier-Stokes system first presented in [25]. This taxonomy illuminates the similarities and differences among these preconditioners and the central role played by efficient approximation of certain Schur complement operators. We then present a parallel computational study that examines the performance of these methods and compares them to an additive Schwarz domain decomposition (DD) algorithm. Results are presented for two and three-dimensional steady state problems for enclosed domains and inflow/outflow systems on both structured and unstructured meshes. The numerical experiments are performed using MPSalsa, a stabilized finite element code.