An augmented Lagrangian preconditioner for the 3D stationary incompressible Navier-Stokes equations at high Reynolds number

Author: 

Farrell, P
Mitchell, L
Wechsung, F

Publication Date: 

8 October 2019

Journal: 

SIAM Journal on Scientific Computing

Last Updated: 

2020-08-08T09:14:43.493+01:00

Issue: 

5

Volume: 

41

DOI: 

10.1137/18M1219370

page: 

A3073–A3096-

abstract: 

In [M. Benzi and M. A. Olshanskii, SIAM J. Sci. Comput., 28 (2006), pp. 2095--2113] a preconditioner of augmented Lagrangian type was presented for the two-dimensional stationary incompressible Navier--Stokes equations that exhibits convergence almost independent of Reynolds number. The algorithm relies on a highly specialized multigrid method involving a custom prolongation operator and for robustness requires the use of piecewise constant finite elements for the pressure. However, the prolongation operator and velocity element used do not directly extend to three dimensions: the local solves necessary in the prolongation operator do not satisfy the inf-sup condition. In this work we generalize the preconditioner to three dimensions, proposing alternative finite elements for the velocity and prolongation operators for which the preconditioner works robustly. The solver is effective at high Reynolds number: on a three-dimensional lid-driven cavity problem with approximately one billion degrees of freedom, the average number of Krylov iterations per Newton step varies from 4.5 at Re = 10 to 3 at Re = 1000 and 5 at Re = 5000.

Symplectic id: 

1037297

Submitted to ORA: 

Submitted

Publication Type: 

Journal Article