An Augmented Lagrangian Preconditioner for the 3D Stationary Incompressible Navier--Stokes Equations at High Reynolds Number

From MaRDI portal
Publication:5240802

DOI10.1137/18M1219370zbMATH Open1448.65261arXiv1810.03315OpenAlexW3102568135WikidataQ127126513 ScholiaQ127126513MaRDI QIDQ5240802FDOQ5240802


Authors: P. E. Farrell, Lawrence Mitchell, Florian Wechsung Edit this on Wikidata


Publication date: 29 October 2019

Published in: SIAM Journal on Scientific Computing (Search for Journal in Brave)

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 mathrmRe=4.5 to 3 at mathrmRe=1000 and 5 at mathrmRe=5000.


Full work available at URL: https://arxiv.org/abs/1810.03315




Recommendations




Cites Work


Cited In (47)

Uses Software





This page was built for publication: An Augmented Lagrangian Preconditioner for the 3D Stationary Incompressible Navier--Stokes Equations at High Reynolds Number

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5240802)