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

Abstract

In Benzi & Olshanskii (SIAM J.~Sci.~Comput., 28(6) (2006)) 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 is tightly coupled to 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 $7$ at $\mathrm{Re}=4.5$ to $3$ at $\mathrm{Re}=1000$ and $5$ at $\mathrm{Re}=5000$.

References in corpus (4)

• Firedrake: automating the finite element method by composing abstractions
• Composing Scalable Nonlinear Algebraic Solvers
• Solver composition across the PDE/linear algebra barrier
• TSFC: a structure-preserving form compiler

Cited by in corpus (18)

• A Reynolds-robust preconditioner for the Scott-Vogelius discretization of the stationary incompressible Navier-Stokes equations
• PCPATCH: software for the topological construction of multigrid relaxation methods
• A constrained transport divergence-free finite element method for Incompressible MHD equations
• An augmented Lagrangian preconditioner for the magnetohydrodynamics equations at high Reynolds and coupling numbers
• A novel block non-symmetric preconditioner for mixed-hybrid finite-element-based flow simulations
• Fully Parallel Mesh I/O using PETSc DMPlex with an Application to Waveform Modeling
• Monolithic multigrid for implicit Runge-Kutta discretizations of incompressible fluid flow
• Recycling augmented Lagrangian preconditioner in an incompressible fluid solver
• Robust and Efficient Multilevel-ILU Preconditioning of Hybrid Newton-GMRES for Incompressible Navier-Stokes Equations
• Uniform block-diagonal preconditioners for divergence-conforming HDG Methods for the generalized Stokes equations and the linear elasticity equations
• Two conjectures on the Stokes complex in three dimensions on Freudenthal meshes
• Analysis and computation of a pressure-robust method for the rotation form of the stationary incompressible Navier-Stokes equations by using high-order finite elements
• Scalable augmented Lagrangian preconditioners for fictitious domain problems
• On the stability and convergence of Discontinuous Galerkin schemes for incompressible flow
• Bifurcation analysis of two-dimensional Rayleigh--Bénard convection using deflation
• Velocity-vorticity-pressure formulation for the Oseen problem with variable viscosity
• Parallel subspace correction methods for semicoercive and nearly semicoercive convex optimization with applications to nonlinear PDEs
• A thermodynamically consistent Johnson-Segalman-Giesekus model: numerical simulation of the rod climbing effect