Symmetric Interior Penalty Discontinuous Galerkin Discretisations and Block Preconditioning for Heterogeneous Stokes Flow
For computational fluid dynamics practitioners, this work provides a robust solver for heterogeneous Stokes flow, though it is an incremental improvement over existing multigrid methods.
This paper develops stable SIP discretizations for variable viscosity Stokes flow and proposes a block preconditioned iterative method. The hp-multigrid method using a Q2_1 coarse space achieves robust convergence with respect to viscosity jumps and mild dependence on polynomial order.
Provable stable arbitrary order symmetric interior penalty discontinuous Galerkin (SIP) discretisations of variable viscosity, incompressible Stokes flow utilising $Q^2_k$--$Q_{k-1}$ elements and hierarchical Legendre basis polynomials are developed and investigated.For solving the resulting linear system, a block preconditioned iterative method is proposed. The nested viscous problem is solved by a $hp$-multilevel preconditioned Krylov subspace method. For the $p$-coarsening, a twolevel method utilising element-block Jacobi preconditioned iterations as a smoother is employed. Piecewise bilinear ($Q^2_1$) and piecewise constant ($Q^2_0$) $p$-coarse spaces are considered. Finally, Galerkin $h$-coarsening is proposed and investigated for the two $p$-coarse spaces considered. Through a number of numerical experiments, we demonstrate that utilising the $Q^2_1$ coarse space results in the most robust $hp$-multigrid method for variable viscosity Stokes flow. Using this $Q^2_1$ coarse space we observe that the convergence of the overall Stokes solver appears to be robust with respect to the jump in the viscosity and only mildly depending on the polynomial order $k$. It is demonstrated and supported by theoretical results that the convergence of the SIP discretisations and the iterative methods rely on a sharp choice of the penalty parameter based on local values of the viscosity.