TY - GEN

T1 - Design of a parallel hybrid direct/iterative solver for CFD problems

AU - Thies, Jonas

AU - Wubs, Fred

PY - 2011

Y1 - 2011

N2 - We discuss the parallel implementation of a hybrid direct/iterative solver for a special class of saddle point matrices arising from the discretization of the steady Navier-Stokes equations on an Arakawa C-grid, the F-matrices. The two-level method described here has the following properties: (i) it is very robust, even at comparatively high Reynolds Numbers, (ii) a single parameter controls fill and convergence, making the method straightforward to use, (iii) the convergence rate is independent of the number of unknowns, (iv) it can be implemented on distributed memory machines in a natural way, (v) the matrix on the second level has the same structure and numerical properties as the original problem, so the method can be applied recursively. The implementation focusses on generality, modularity, code reuse and recursiveness. The solver is implemented using building blocks of the Trilinos libraries. We show its performance on a parallel computer for the Navier-Stokes equations.

AB - We discuss the parallel implementation of a hybrid direct/iterative solver for a special class of saddle point matrices arising from the discretization of the steady Navier-Stokes equations on an Arakawa C-grid, the F-matrices. The two-level method described here has the following properties: (i) it is very robust, even at comparatively high Reynolds Numbers, (ii) a single parameter controls fill and convergence, making the method straightforward to use, (iii) the convergence rate is independent of the number of unknowns, (iv) it can be implemented on distributed memory machines in a natural way, (v) the matrix on the second level has the same structure and numerical properties as the original problem, so the method can be applied recursively. The implementation focusses on generality, modularity, code reuse and recursiveness. The solver is implemented using building blocks of the Trilinos libraries. We show its performance on a parallel computer for the Navier-Stokes equations.

UR - http://www.scopus.com/inward/record.url?scp=84856322097&partnerID=8YFLogxK

U2 - 10.1109/eScience.2011.60

DO - 10.1109/eScience.2011.60

M3 - Conference contribution

AN - SCOPUS:84856322097

SN - 9780769545974

T3 - Proceedings - 2011 7th IEEE International Conference on eScience, eScience 2011

SP - 387

EP - 394

BT - Proceedings - 2011 7th IEEE International Conference on eScience, eScience 2011

T2 - 7th IEEE International Conference on eScience, eScience 2011

Y2 - 5 December 2011 through 8 December 2011

ER -