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 -