p-Multilevel Preconditioners for HHO Discretizations of the Stokes Equations with Static Condensation

doi:10.1007/s42967-021-00142-5

Abstract

Abstract: We propose a p-multilevel preconditioner for hybrid high-order (HHO) discretizations of the Stokes equation, numerically assess its performance on two variants of the method, and compare with a classical discontinuous Galerkin scheme. An efficient implementation is proposed where coarse level operators are inherited using $ L^2 $-orthogonal projections defined over mesh faces and the restriction of the fine grid operators is performed recursively and matrix-free. Both h- and k-dependency are investigated tackling two- and three-dimensional problems on standard meshes and graded meshes. For the two HHO formulations, featuring discontinuous or hybrid pressure, we study how the combination of p-coarsening and static condensation influences the V-cycle iteration. In particular, two different static condensation procedures are considered for the discontinuous pressure HHO variant, resulting in global linear systems with a different number of unknowns and matrix non-zero entries. Interestingly, we show that the efficiency of the solution strategy might be impacted by static condensation options in the case of graded meshes.

Key words: Stokes equations, Divergence free constraint, Hybrid high-order, Discontinuous Galerkin, p-multigrid, Static condensation

CLC Number:

Lorenzo Botti, Daniele A. Di Pietro. p-Multilevel Preconditioners for HHO Discretizations of the Stokes Equations with Static Condensation[J]. Communications on Applied Mathematics and Computation, 2022, 4(3): 783-822.

TrendMD

References

[1] Adler, J.H., Benson, T.R., MacLachlan, S.P.: Preconditioning a mass-conserving discontinuous Galerkin discretization of the Stokes equations. Numer. Linear Algebra Appl. 24(3), e2047 (2017)
[2] Aghili, J., Boyaval, S., Di Pietro, D.A.: Hybridization of mixed high-order methods on general meshes and application to the Stokes equations. Comput. Methods Appl. Math. 15(2), 111–134 (2015)
[3] Antonietti, P.F., Giani, S., Houston, P.:

\begin{document}$ hp $\end{document}

-Version composite discontinuous Galerkin methods for elliptic problems on complicated domains. SIAM J. Sci. Comput. 35(3), A1417–A1439 (2013)
[4] Antonietti, P.F., Sarti, M., Verani, M.: Multigrid algorithms for

\begin{document}$ hp $\end{document}

-discontinuous Galerkin discretizations of elliptic problems. SIAM J. Numer. Anal. 53(1), 598–618 (2015)
[5] Antonietti, P.F., Cangiani, A., Collis, J., Dong, Z., Georgoulis, E.H., Giani, S., Houston, P.: Review of discontinuous Galerkin finite element methods for partial differential equations on complicated domains. In: Building Bridges: Connections and Challenges in Modern Approaches to Numerical Partial Differential Equations, volume 114 of Lect. Notes Comput. Sci. Eng., pp. 279–308. Springer, Cham (2016)
[6] Antonietti, P.F., Houston, P., Pennesi, G., Süli, E.: An agglomeration-based massively parallel non-overlapping additive Schwarz preconditioner for high-order discontinuous Galerkin methods on polytopic grids. Math. Comput. 89(325), 2047–2083 (2020)
[7] Ayuso de Dios, B., Brezzi, F., Marini, L., Xu, J., Zikatanov, L.: A simple preconditioner for a discontinuous Galerkin method for the Stokes problem. J. Sci. Comput. 58, 517–547 (2012)
[8] Balay, S., Abhyankar, S., Adams, M.F., Brown, J., Brune, P., Buschelman, K., Dalcin, L., Dener, A., Eijkhout, V., Gropp, W.D., Karpeyev, D., Kaushik, D., Knepley, M.G., May, D.A., McInnes, L.C., Mills, R.T., Munson, T., Rupp, K., Sanan, P., Smith, B.F., Zampini, S., Zhang, H., Zhang, H.: PETSc Web page. https://www.mcs.anl.gov/petsc (2019)
[9] Bassi, F., Rebay, S.: A high-order accurate discontinuous finite element method for the numerical solution of the compressible Navier-Stokes equations. J. Comput. Phys. 131(2), 267–279 (1997)
[10] Bassi, F., Rebay, S., Mariotti, G., Pedinotti, S., Savini, M.: A high-order accurate discontinuous finite element method for inviscid and viscous turbomachinery flows. In: Decuypere, R., Dibelius, G. (eds) Proceedings of the 2nd European Conference on Turbomachinery Fluid Dynamics and Thermodynamics, pp. 99–109 (1997)
[11] Bassi, F., Crivellini, A., Di Pietro, D.A., Rebay, S.: An artificial compressibility flux for the discontinuous Galerkin solution of the incompressible Navier-Stokes equations. J. Comput. Phys. 218(2), 794–815 (2006)
[12] Bassi, F., Ghidoni, A., Rebay, S., Tesini, P.: High-order accurate

\begin{document}$ p $\end{document}

-multigrid discontinuous Galerkin solution of the Euler equations. Int. J. Numer. Methods Fluids 60(8), 847–865 (2009)
[13] Bassi, F., Botti, L., Colombo, A., Rebay, S.: Agglomeration based discontinuous Galerkin discretization of the Euler and Navier-Stokes equations. Comput. Fluids 61, 77–85 (2012)
[14] Bassi, F., Botti, L., Colombo, A., Di Pietro, D.A., Tesini, P.: On the flexibility of agglomeration based physical space discontinuous Galerkin discretizations. J. Comput. Phys. 231(1), 45–65 (2012)
[15] Bassi, F., Botti, L., Colombo, A.: Agglomeration-based physical frame DG discretizations: an attempt to be mesh free. Math. Models Methods Appl. Sci. 24(8), 1495–1539 (2014)
[16] Botti, L., Di Pietro, D.A.: Numerical assessment of hybrid high-order methods on curved meshes and comparison with discontinuous Galerkin methods. J. Comput. Phys. 370, 58–84 (2018)
[17] Botti, L., Colombo, A., Bassi, F.:

\begin{document}$ h $\end{document}

-Multigrid agglomeration based solution strategies for discontinuous Galerkin discretizations of incompressible flow problems. J. Comput. Phys. 347, 382–415 (2017)
[18] Botti, L., Di Pietro, D.A., Droniou, J.: A hybrid high-order discretisation of the Brinkman problem robust in the Darcy and Stokes limits. Comput. Methods Appl. Mech. Eng. 341, 278–310 (2018)
[19] Botti, L., Colombo, A., Crivellini, A., Franciolini, M.:

\begin{document}$ h $\end{document}

\begin{document}$ p $\end{document}

\begin{document}$ hp $\end{document}

-Multilevel discontinuous Galerkin solution strategies for elliptic operators. Int. J. Comput. Fluid Dyn. 33(9), 362–370 (2019)
[20] Botti, L., Di Pietro, D.A., Droniou, J.: A hybrid high-order method for the incompressible Navier-Stokes equations based on Temam’s device. J. Comput. Phys. 376, 786–816 (2019)
[21] Cangiani, A., Dong, Z., Georgoulis, E.H., Houston, P.:

\begin{document}$ hp $\end{document}

-Version Discontinuous Galerkin Methods on Polygonal and Polyhedral Meshes. Springer Briefs in Mathematics. Springer, Cham (2017)
[22] Charrier, D.E., May, D.A., Schnepp, S.M.: Symmetric interior penalty discontinuous Galerkin discretizations and block preconditioning for heterogeneous stokes flow. SIAM J. Sci. Comput. 39(6), B1021–B1042 (2017)
[23] Cockburn, B., Shu, C.-W.: TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws II. General framework. Math. Comput. 52(186), 411–435 (1989)
[24] Cockburn, B., Shu, C.-W.: The Runge-Kutta local projection

\begin{document}$ P^1 $\end{document}

-discontinuous-Galerkin finite element method for scalar conservation laws. RAIRO Modél. Math. Anal. Numér. 25(3), 337–361 (1991)
[25] Cockburn, B., Shu, C.-W.: The Runge-Kutta discontinuous Galerkin method for conservation laws V. Multidimensional systems. J. Comput. Phys. 141(2), 199–224 (1998)
[26] Cockburn, B., Lin, S.Y., Shu, C.-W.: TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws III. One-dimensional systems. J. Comput. Phys. 84(1), 90–113 (1989)
[27] Cockburn, B., Hou, S., Shu, C.-W.: The Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws IV. The multidimensional case. Math. Comput. 54(190), 545–581 (1990)
[28] Cockburn, B., Kanschat, G., Schötzau, D., Schwab, C.: Local discontinuous Galerkin methods for the Stokes system. SIAM J. Numer. Anal. 40(1), 319–343 (2002)
[29] Cockburn, B., Dubois, O., Gopalakrishnan, J., Tan, S.: Multigrid for an HDG method. IMA J. Numer. Anal. 34(4), 1386–1425 (2013)
[30] Cockburn, B., Di Pietro, D.A., Ern, A.: Bridging the Hybrid High-Order and Hybridizable Discontinuous Galerkin methods. ESAIM Math. Model. Numer. Anal. 50(3), 635–650 (2016)
[31] Di Pietro, D.A.: Analysis of a discontinuous Galerkin approximation of the Stokes problem based on an artificial compressibility flux. Int. J. Numer. Methods Fluids 55(8), 793–813 (2007)
[32] Di Pietro, D.A., Droniou, J.: The Hybrid High-order Method for Polytopal Meshes. Design, Analysis, and Applications. Number 19 in Modeling, Simulation and Application. Springer International Publishing (2020)
[33] Di Pietro, D.A., Ern, A.: Discrete functional analysis tools for discontinuous Galerkin methods with application to the incompressible Navier-Stokes equations. Math. Comput. 79(271), 1303–1330 (2010)
[34] Di Pietro, D.A., Ern, A.: Mathematical Aspects of Discontinuous Galerkin Methods. Mathematics and Applications, vol. 69. Springer, Heidelberg (2012)
[35] Di Pietro, D.A., Krell, S.: Benchmark session: the 2D Hybrid High-Order method. In: Cancès, C., Omnes, P. (eds) Finite Volumes for Complex Applications VIII: Methods and Theoretical Aspects. pp. 91–106 (2017)
[36] Di Pietro, D.A., Krell, S.: A hybrid high-order method for the steady incompressible Navier-Stokes problem. J. Sci. Comput. 74(3), 1677–1705 (2018)
[37] Di Pietro, D.A., Ern, A., Lemaire, S.: An arbitrary-order and compact-stencil discretization of diffusion on general meshes based on local reconstruction operators. Comput. Methods Appl. Math. 14(4), 461–472 (2014)
[38] Di Pietro, D.A., Ern, A., Lemaire, S.: Building Bridges: Connections and Challenges in Modern Approaches to Numerical Partial Differential Equations, A Review of Hybrid High-order Methods: Formulations, Computational Aspects, Comparison with Other Methods, no 114 in Lecture Notes in Computational Science and Engineering, pp. 205–236. Springer (2016)
[39] Di Pietro, D.A., Hülsemann, F., Matalon, P., Mycek, P., Rüde, U., Ruiz, D.: An

\begin{document}$ h $\end{document}

-multigrid method for hybrid high-order discretizations. SIAM J. Sci. Comput. (2021). (To appear)
[40] Elman, H.C., Silvester, D.J., Wathen, A.J.: Finite Elements and Fast Iterative Solvers: with Applications in Incompressible Fluid Dynamics. Numerical Mathematics and Scientific Computation, 2nd edn. Oxford University Press, Oxford (2014)
[41] Fabien, M.S., Knepley, M.G., Mills, R.T., Rivière, B.M.: Manycore parallel computing for a hybridizable discontinuous Galerkin nested multigrid method. SIAM J. Sci. Comput. 41(2), C73–C96 (2019)
[42] Fidkowski, K.J., Oliver, T.A., Lu, J., Darmofal, D.L.:

\begin{document}$ p $\end{document}

-Multigrid solution of high-order discontinuous Galerkin discretizations of the compressible Navier-Stokes equations. J. Comput. Phys. 207(1), 92–113 (2005)
[43] Franciolini, M., Botti, L., Colombo, A., Crivellini, A.:

\begin{document}$ p $\end{document}

-Multigrid matrix-free discontinuous Galerkin solution strategies for the under-resolved simulation of incompressible turbulent flows. Comput. Fluids 206, 104558 (2020)
[44] Franciolini, M., Fidkowski, K.J., Crivellini, A.: Efficient discontinuous Galerkin implementations and preconditioners for implicit unsteady compressible flow simulations. Comput. Fluids 203, 104542 (2020)
[45] Kanschat, G., Mao, Y.: Multigrid methods for Hdiv-conforming discontinuous Galerkin methods for the Stokes equations. J. Numer. Math. 23(1), 51–66 (2015)
[46] Kronbichler, M., Wall, W.A.: A performance comparison of continuous and discontinuous Galerkin methods with fast multigrid solvers. SIAM J. Sci. Comput. 40(5), A3423–A3448 (2018)
[47] Labeur, R.J., Wells, G.N.: Energy stable and momentum conserving hybrid finite element method for the incompressible Navier-Stokes equations. SIAM J. Sci. Comput. 34(2), A889–A913 (2012)
[48] Nastase, C.R., Mavriplis, D.J.: High-order discontinuous Galerkin methods using an

\begin{document}$ hp $\end{document}

-multigrid approach. J. Comput. Phys. 213(1), 330–357 (2006)
[49] Nguyen, N.C., Peraire, J., Cockburn, B.: An implicit high-order hybridizable discontinuous Galerkin method for the incompressible Navier-Stokes equations. J. Comput. Phys. 230(4), 1147–1170 (2011)
[50] Quarteroni, A., Sacco, R., Saleri, F.: Numerical Mathematics, vol. 37. Springer, Berlin (2000)
[51] Rhebergen, S., Wells, G.: Preconditioning of a hybridized discontinuous Galerkin finite element method for the Stokes equations. J. Sci. Comput. 77, 1936–1952 (2018)
[52] Rhebergen, S., Wells, G.N.: A hybridizable discontinuous Galerkin method for the Navier-Stokes equations with pointwise divergence-free velocity field. J. Sci. Comput. 76(3), 1484–1501 (2018)
[53] Rønquist, E.M., Patera, A.T.: Spectral element multigrid I. Formulation and numerical results. J. Sci. Comput. 2(4), 389–406 (1987)
[54] Shahbazi, K., Mavriplis, D.J., Burgess, N.K.: Multigrid algorithms for high-order discontinuous Galerkin discretizations of the compressible Navier-Stokes equations. J. Comput. Phys. 228(21), 7917–7940 (2009)
[55] Smith, B.F., Bjørstad, P.E., Gropp, W.D.: Domain Decomposition: Parallel Multilevel Methods for Elliptic Partial Differential Equations. Cambridge University Press, Cambridge (1996)