Premium
Performance of fully coupled algebraic multilevel domain decomposition preconditioners for incompressible flow and transport
Author(s) -
Lin Paul T.,
Sala Marzio,
Shadid John N.,
Tuminaro Ray S.
Publication year - 2006
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.1624
Subject(s) - preconditioner , domain decomposition methods , multigrid method , mathematics , linear system , rate of convergence , robustness (evolution) , cholesky decomposition , finite element method , mathematical optimization , computer science , key (lock) , eigenvalues and eigenvectors , mathematical analysis , partial differential equation , biochemistry , chemistry , physics , computer security , quantum mechanics , gene , thermodynamics
Abstract This study investigates algebraic multilevel domain decomposition preconditioners of the Schwarz type for solving linear systems associated with Newton–Krylov methods. The key component of the preconditioner is a coarse approximation based on algebraic multigrid ideas to approximate the global behaviour of the linear system. The algebraic multilevel preconditioner is based on an aggressive coarsening graph partitioning of the non‐zero block structure of the Jacobian matrix. The scalability of the preconditioner is presented as well as comparisons with a two‐level Schwarz preconditioner using a geometric coarse grid operator. These comparisons are obtained on large‐scale distributed‐memory parallel machines for systems arising from incompressible flow and transport using a stabilized finite element formulation. The results demonstrate the influence of the smoothers and coarse level solvers for a set of 3D example problems. For preconditioners with more than one level, careful attention needs to be given to the balance of robustness and convergence rate for the smoothers and the cost of applying these methods. For properly chosen parameters, the two‐ and three‐level preconditioners are demonstrated to be scalable to 1024 processors. Copyright © 2006 John Wiley & Sons, Ltd.