[petsc-users] Field Split PC for Fully-Coupled 3d stationary incompressible Navier-Stokes Solution Algorithm

[Klaij, Christiaan]

Hi Fabian,

After reading your thread, I think we have the same
discretization. Some thoughts:

- The stabilization by pressure-weighted interpolation should
  lead to a diagonally dominant system that can be solved by
  algebraic multigrid (done for example in CFX), see
  http://dx.doi.org/10.1080/10407790.2014.894448
  http://dx.doi.org/10.1016/j.jcp.2008.08.027

- If you go for fieldsplit with mild tolerances, the
  preconditioner will become variable and you might need FGMRES
  or GCR instead.

- As you point out, most preconditioners are designed for stable
  FEM. My experience is that those conclusions sometimes (but not
  always) hold for FVM. Due to the specific form of the
  stabilization, other choices might be feasible. More research
  is needed :-)

Chris


dr. ir. Christiaan Klaij
CFD Researcher
Research & Development
E mailto:C.Klaij at marin.nl
T +31 317 49 33 44


MARIN
2, Haagsteeg, P.O. Box 28, 6700 AA Wageningen, The Netherlands
T +31 317 49 39 11, F +31 317 49 32 45, I www.marin.nl