Robust and Efficient Multilevel-ILU Preconditioned Newton-GMRES for Incompressible Navier-Stokes

11/14/2020
by   Qiao Chen, et al.
0

We introduce a new preconditioned Newton-GMRES method for solving the nonlinear systems arising from incompressible Navier-Stokes equations. When the Reynolds number is relatively high, these systems often involve millions of degrees of freedom (DOFs), and the nonlinear systems are difficult to converge, partially due to their saddle-point structure. In this work, we propose to alleviate these issues by introducing a new preconditioned nonlinear solver called HILUNG, or Hierarchical Incomplete-LU preconditioned Newton-Gmres method. A key novelty of HILUNG is to incorporate an accurate and stable multilevel preconditioner called HILUCSI, which is particularly effective for solving saddle-point problems. HILUCSI enables robust and rapid convergence of the inner iterations in Newton-GMRES. We also introduce physics-aware sparsifying operators, adaptive refactorization and thresholding, and iterative refinement, to improve efficiency without compromising robustness. We show that HILUNG can robustly solve the standard 2D driven-cavity problem with Re 5000, while other nonlinear solvers failed to converge at Re 1000 with a similar configuration. HILUNG also improved the efficiency over another state-of-the-art multilevel ILU preconditioner and a multi-threaded direct solver by more than an order of magnitude for the 3D flow-over-cylinder problem with one million DOFs and enabled the efficient solution with about ten million DOFs using only 60GB of RAM while others fail due to nonrobustness or memory limitation.

READ FULL TEXT

page 1

page 12

research
11/22/2019

HILUCSI: Simple, Robust, and Fast Multilevel ILU for Large-Scale Saddle-Point Problems from PDEs

Incomplete factorization is a widely used preconditioning technique for ...
research
03/24/2020

Advanced Newton Methods for Geodynamical Models of Stokes Flow with Viscoplastic Rheologies

Strain localization and resulting plasticity and failure play an importa...
research
04/22/2022

Robust and efficient primal-dual Newton-Krylov solvers for viscous-plastic sea-ice models

We present a Newton-Krylov solver for a viscous-plastic sea-ice model. T...
research
09/02/2020

Approximate Generalized Inverses with Iterative Refinement for ε-Accurate Preconditioning of Singular Systems

We introduce a new class of preconditioners to enable flexible GMRES to ...
research
06/07/2021

Smooth Implicit Hybrid Upwinding for Compositional Multiphase Flow in Porous Media

In subsurface multiphase flow simulations, poor nonlinear solver perform...
research
07/29/2020

Localized Nonlinear Solution Strategies for Efficient Simulation of Unconventional Reservoirs

Accurate and efficient numerical simulation of unconventional reservoirs...
research
04/16/2019

Short Overview of Early Developments of the Hardy Cross Type Methods for Computation of Flow Distribution in Pipe Networks

Hardy Cross originally proposed a method for analysis of flow in network...

Please sign up or login with your details

Forgot password? Click here to reset