Weighted BFBT Preconditioner for Stokes Flow Problems with Highly Heterogeneous Viscosity

Abstract

We present a weighted BFBT approximation (w-BFBT) to the inverse Schur complement of a Stokes system with highly heterogeneous viscosity. When used as part of a Schur complement-based Stokes preconditioner, we observe robust fast convergence for Stokes problems with smooth but highly varying (up to 10 orders of magnitude) viscosities, optimal algorithmic scalability with respect to mesh refinement, and only a mild dependence on the polynomial order of high-order finite element discretizations ($Q_k\times P_{k-1}^{disc}$, order $k\ge 2$). For certain difficult problems, we demonstrate numerically that w-BFBT significantly improves Stokes solver convergence over the widely used inverse viscosity-weighted pressure mass matrix approximation of the Schur complement. In addition, we derive theoretical eigenvalue bounds to prove spectral equivalence of w-BFBT. Using detailed numerical experiments, we discuss modifications to w-BFBT at Dirichlet boundaries that decrease the number of iterations. The overall algorithmic performance of the Stokes solver is governed by the efficacy of w-BFBT as a Schur complement approximation and, in addition, by our parallel hybrid spectral-geometric-algebraic multigrid (HMG) method, which we use to approximate the inverses of the viscous block and variable-coefficient pressure Poisson operators within w-BFBT. Building on the scalability of HMG, our Stokes solver achieves a parallel efficiency of 90% while weak scaling over a more than 600-fold increase from 48 to all 30,000 cores of TACC's Lonestar 5 supercomputer.

arXiv commentTo appear in SIAM Journal on Scientific Computing

References in corpus (2)

• Scalable and efficient algorithms for the propagation of uncertainty from data through inference to prediction for large-scale problems, with application to flow of the Antarctic ice sheet
• Solution of nonlinear Stokes equations discretized by high-order finite elements on nonconforming and anisotropic meshes, with application to ice sheet dynamics

Cited by in corpus (8)

• High-order Methods for a Pressure Poisson Equation Reformulation of the Navier-Stokes Equations with Electric Boundary Conditions
• Advanced Newton Methods for Geodynamical Models of Stokes Flow with Viscoplastic Rheologies
• On augmented finite element formulation for the Navier--Stokes equations with vorticity and variable viscosity
• Divergence conforming DG method for the optimal control of the Oseen equation with variable viscosity
• An AMG saddle point preconditioner with application to mixed Poisson problems on adaptive quad/cube meshes
• Velocity-vorticity-pressure formulation for the Oseen problem with variable viscosity
• Incorporating variable viscosity in vorticity-based formulations for Brinkman equations
• Higher-Order Discontinuous Galerkin Splitting Schemes for Fluids with Variable Viscosity