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

Jonas Thies*, Fred Wubs

*Corresponding author for this work

Research output: Chapter in Book/Conference proceedings/Edited volumeConference contributionScientificpeer-review

5 Citations (Scopus)

Abstract

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.

Original languageEnglish
Title of host publicationProceedings - 2011 7th IEEE International Conference on eScience, eScience 2011
Pages387-394
Number of pages8
DOIs
Publication statusPublished - 2011
Externally publishedYes
Event7th IEEE International Conference on eScience, eScience 2011 - Stockholm, Sweden
Duration: 5 Dec 20118 Dec 2011

Publication series

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

Conference

Conference7th IEEE International Conference on eScience, eScience 2011
CountrySweden
CityStockholm
Period5/12/118/12/11

Fingerprint

Dive into the research topics of 'Design of a parallel hybrid direct/iterative solver for CFD problems'. Together they form a unique fingerprint.

Cite this