Robust preconditioning for an HDG discretization of the time-dependent Stokes equations
Pith reviewed 2026-05-10 17:50 UTC · model grok-4.3
The pith
Parameter-robust preconditioners are derived for the statically condensed hybridizable discontinuous Galerkin discretization of the time-dependent Stokes equations.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We extend the analysis to derive new preconditioners that remain robust with respect to all physical and discretization parameters. The construction relies on first establishing uniform well-posedness of the HDG formulation (before static condensation) through appropriately defined norms. Based on this result, we identify sufficient conditions that a norm on the face space must satisfy to guarantee parameter-robustness of the resulting preconditioner for the statically condensed HDG system. Numerical experiments in two and three dimensions verify our theoretical results.
What carries the argument
The sufficient conditions on the norm defined on the face space, which ensure parameter-robustness of the preconditioner for the statically condensed HDG system, derived from the uniform well-posedness of the pre-condensation HDG formulation.
Load-bearing premise
The HDG formulation before static condensation must be uniformly well-posed in appropriately defined norms, and the chosen norm on the face space must satisfy the sufficient conditions needed for the preconditioner to inherit that robustness.
What would settle it
If numerical experiments show that the number of iterations required by the preconditioned iterative solver grows without bound as the viscosity parameter tends to zero or as the time step size is varied independently, while holding other parameters fixed, the claim of parameter-robustness would be falsified.
read the original abstract
We present parameter-robust preconditioners for linear systems that arise after applying static condensation to a hybridizable discontinuous Galerkin (HDG) discretization of the time-dependent Stokes problem. Building upon the theoretical framework introduced in our previous work [SIAM Journal on Scientific Computing, 47(6):A3212-A3238, 2025], we extend the analysis to derive new preconditioners that remain robust with respect to all physical and discretization parameters. The construction relies on first establishing uniform well-posedness of the HDG formulation (before static condensation) through appropriately defined norms. Based on this result, we identify sufficient conditions that a norm on the face space must satisfy to guarantee parameter-robustness of the resulting preconditioner for the statically condensed HDG system. Numerical experiments in two and three dimensions verify our theoretical results.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops parameter-robust preconditioners for the statically condensed linear systems arising from an HDG discretization of the time-dependent Stokes equations. It extends the authors' prior framework by first proving uniform well-posedness of the pre-condensation HDG formulation via suitably chosen norms, then deriving sufficient conditions on a norm for the face space that transfer this robustness to the condensed system. The claims are supported by numerical experiments in two and three dimensions across relevant parameter regimes.
Significance. If the central claims hold, the work advances robust solver technology for HDG methods on time-dependent incompressible flows, addressing a practical need for preconditioners that remain effective under variation of viscosity, time-step size, polynomial degree, and mesh size. The explicit reliance on uniform well-posedness before condensation and the numerical verification across 2D/3D regimes constitute clear strengths; the extension of the cited SIAM J. Sci. Comput. framework without introducing new parameter dependencies is a positive technical feature.
major comments (2)
- [Section 3 (analysis of well-posedness and face-norm conditions)] The abstract and outline indicate that uniform well-posedness of the HDG formulation (before condensation) is established via appropriately defined norms, followed by sufficient conditions on the face-space norm. However, the manuscript does not appear to contain an explicit verification that the chosen face norm satisfies those sufficient conditions once the time-dependent term is included; this step is load-bearing for the parameter-robustness claim.
- [Section 5 (numerical results)] The numerical experiments are stated to confirm robustness, yet the presentation does not include a direct comparison of iteration counts or condition-number bounds against the theoretical predictions for the full range of physical and discretization parameters (e.g., small viscosity combined with large time steps). This weakens the link between theory and verification.
minor comments (2)
- [Section 2] Notation for the face-space norm and the condensed Schur complement should be introduced with explicit reference to the corresponding definitions in the authors' prior SIAM J. Sci. Comput. paper to improve readability.
- [Section 2.1] A few typographical inconsistencies appear in the statement of the time-dependent Stokes problem (e.g., placement of the time derivative relative to the viscous term).
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments on our manuscript. We appreciate the positive assessment of the work's significance and address each major comment below, outlining the revisions we will make to strengthen the presentation.
read point-by-point responses
-
Referee: [Section 3 (analysis of well-posedness and face-norm conditions)] The abstract and outline indicate that uniform well-posedness of the HDG formulation (before condensation) is established via appropriately defined norms, followed by sufficient conditions on the face-space norm. However, the manuscript does not appear to contain an explicit verification that the chosen face norm satisfies those sufficient conditions once the time-dependent term is included; this step is load-bearing for the parameter-robustness claim.
Authors: We thank the referee for highlighting this point. Section 3 establishes uniform well-posedness of the pre-condensation HDG formulation using suitably chosen norms and then derives sufficient conditions on the face-space norm that transfer robustness to the condensed system. The specific face norm employed is constructed by extending the norm from our prior SIAM J. Sci. Comput. framework to incorporate the time-dependent term, ensuring by design that the conditions hold. To make this verification fully explicit and address the load-bearing nature of the step, we will add a dedicated paragraph or short subsection in the revised Section 3 that directly confirms the chosen norm satisfies all stated sufficient conditions, including the contributions from the time-dependent term. revision: yes
-
Referee: [Section 5 (numerical results)] The numerical experiments are stated to confirm robustness, yet the presentation does not include a direct comparison of iteration counts or condition-number bounds against the theoretical predictions for the full range of physical and discretization parameters (e.g., small viscosity combined with large time steps). This weakens the link between theory and verification.
Authors: We agree that a more explicit link between the theoretical predictions and the numerical results would improve the manuscript. Section 5 already reports iteration counts and timings for the preconditioned systems across 2D and 3D test cases, covering wide ranges of viscosity, time-step size, polynomial degree, and mesh size. In the revised version, we will augment Section 5 with additional tables or figures that directly compare the observed iteration counts and estimated condition numbers to the theoretical bounds for the complete parameter space, with particular emphasis on the regime of small viscosity combined with large time steps. This will provide a clearer quantitative verification of the robustness claims. revision: yes
Circularity Check
Minor self-citation; central derivation remains independent
full rationale
The paper explicitly builds on the authors' prior SIAM J. Sci. Comput. framework for HDG preconditioning but adds new analysis for the time-dependent Stokes case, including uniform well-posedness via chosen norms and sufficient conditions on the face-space norm to transfer robustness after static condensation. No step reduces a prediction or uniqueness claim to a fitted input, self-referential definition, or unverified self-citation chain. The extension carries over without introducing new parameter dependencies that would force circularity, and numerical verification is external to the derivation.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Uniform well-posedness of the HDG formulation before static condensation holds in appropriately defined norms
Forward citations
Cited by 1 Pith paper
-
Preconditioning of a hybridizable discontinuous Galerkin method for the coupled Stokes--Darcy system
Parameter-robust preconditioners are designed and analyzed for the statically condensed HDG discretization of the Stokes-Darcy system by first proving uniform well-posedness and then applying operator-preconditioning ...
Reference graph
Works this paper leans on
-
[1]
LocalL 2-bounded commuting projections in FEEC.ESAIM Math
Douglas Arnold and Johnny Guzm´ an. LocalL 2-bounded commuting projections in FEEC.ESAIM Math. Model. Numer. Anal., 55(5):2169–2184, 2021
work page 2021
-
[2]
Trygve Bærland, Miroslav Kuchta, Kent-Andre Mardal, and Travis Thompson. An observation on the uniform preconditioners for the mixed Darcy problem.Numerical Methods for Partial Differential Equations, 36(6):1718–1734, 2020
work page 2020
-
[3]
Daniele Boffi, Franco Brezzi, and Michel Fortin.Mixed finite element methods and applications, vol- ume 44. Springer, 2013
work page 2013
-
[4]
Iterative techniques for time dependent Stokes problems
James H Bramble and Joseph E Pasciak. Iterative techniques for time dependent Stokes problems. Computers & Mathematics with Applications, 33(1-2):13–30, 1997
work page 1997
-
[5]
Some fast 3D finite element solvers for the generalized Stokes problem
J Cahouet and J-P Chabard. Some fast 3D finite element solvers for the generalized Stokes problem. International Journal for Numerical Methods in Fluids, 8(8):869–895, 1988
work page 1988
-
[6]
Multigrid for an HDG method.IMA Journal of Numerical Analysis, 34(4):1386–1425, 2014
Bernardo Cockburn, Olivier Dubois, Jay Gopalakrishnan, and Shuguang Tan. Multigrid for an HDG method.IMA Journal of Numerical Analysis, 34(4):1386–1425, 2014
work page 2014
-
[7]
Bernardo Cockburn, Jayadeep Gopalakrishnan, and Raytcho Lazarov. Unified hybridization of discon- tinuous Galerkin, mixed, and continuous Galerkin methods for second order elliptic problems.SIAM Journal on Numerical Analysis, 47(2):1319–1365, 2009
work page 2009
-
[8]
Springer Science & Business Media, 2011
Daniele Antonio Di Pietro and Alexandre Ern.Mathematical aspects of discontinuous Galerkin methods, volume 69. Springer Science & Business Media, 2011. 25
work page 2011
-
[9]
Shukai Du and Francisco-Javier Sayas.An Invitation to the Theory of the Hybridizable Discontinuous Galerkin Method: Projections, Estimates, Tools. Springer Nature, 2019
work page 2019
-
[10]
Alexandre Ern, Thirupathi Gudi, Iain Smears, and Martin Vohral´ ık. Equivalence of local- and global- best approximations, a simple stable local commuting projector, and optimalhpapproximation esti- mates inH(div).IMA J. Numer. Anal., 42(2):1023–1049, 2022
work page 2022
-
[11]
Guosheng Fu. Uniform auxiliary space preconditioning for HDG methods for elliptic operators with a parameter dependent low order term.SIAM Journal on Scientific Computing, 43(6):A3912–A3937, 2021
work page 2021
-
[12]
Guosheng Fu and Wenzheng Kuang. Uniform block-diagonal preconditioners for divergence-conforming HDG methods for the generalized Stokes equations and the linear elasticity equations.IMA Journal of Numerical Analysis, 43(3):1718–1741, 2023
work page 2023
-
[13]
G. P. Galdi.An introduction to the mathematical theory of the Navier-Stokes equations. Springer Monographs in Mathematics. Springer, New York, second edition, 2011. Steady-state problems
work page 2011
-
[14]
Evan Gawlik, Michael J. Holst, and Martin W. Licht. Local finite element approximation of Sobolev differential forms.ESAIM Math. Model. Numer. Anal., 55(5):2075–2099, 2021
work page 2075
-
[15]
On the equation divu=gand Bogovski ˘i’s operator in Sobolev spaces of negative order
Matthias Geiß ert, Horst Heck, and Matthias Hieber. On the equation divu=gand Bogovski ˘i’s operator in Sobolev spaces of negative order. InPartial differential equations and functional analysis, volume 168 ofOper. Theory Adv. Appl., pages 113–121. Birkh¨ auser, Basel, 2006
work page 2006
-
[16]
Yunhui He, Sander Rhebergen, and Hans De Sterck. Local Fourier analysis of multigrid for hybridized and embedded discontinuous Galerkin methods.SIAM Journal on Scientific Computing, 43(5):S612– S636, 2021
work page 2021
-
[17]
Esteban Henr´ ıquez, Jeonghun J. Lee, and Sander Rhebergen. Parameter-robust preconditioning for hybridizable symmetric discretizations.SIAM Journal on Scientific Computing, 47(6):A3212–A3238, 2025
work page 2025
-
[18]
Esteban Henr´ ıquez, Jeonghun J Lee, and Sander Rhebergen. Preconditioning of a hybridizable discon- tinuous Galerkin method for Biot’s consolidation model.arXiv preprint arXiv:2508.12991, 2025
-
[19]
Nicholas J Higham.Accuracy and stability of numerical algorithms. SIAM, 2002
work page 2002
-
[20]
Johannes Kraus, Philip L Lederer, Maria Lymbery, Kevin Osthues, and Joachim Sch¨ oberl. Hybridized discontinuous galerkin/hybrid mixed methods for a multiple network poroelasticity model with appli- cation in biomechanics.SIAM Journal on Scientific Computing, 45(6):B802–B827, 2023
work page 2023
-
[21]
Johannes Kraus, Philip L Lederer, Maria Lymbery, and Joachim Sch¨ oberl. Uniformly well-posed hy- bridized discontinuous Galerkin/hybrid mixed discretizations for Biot’s consolidation model.Computer Methods in Applied Mechanics and Engineering, 384:113991, 2021
work page 2021
-
[22]
Alexander D Lindsay, Sander Rhebergen, and Ben S Southworth. Preconditioning a hybridizable discon- tinuous Galerkin method for Navier-Stokes at high Reynolds number.arXiv preprint arXiv:2512.02971, 2025
-
[23]
Peipei Lu, Andreas Rupp, and Guido Kanschat. Analysis of injection operators in geometric multigrid solvers for HDG methods.SIAM Journal on Numerical Analysis, 60(4):2293–2317, 2022
work page 2022
-
[24]
Homogeneous multigrid for HDG.IMA Journal of Numerical Analysis, 42(4):3135–3153, 2022
Peipei Lu, Andreas Rupp, and Guido Kanschat. Homogeneous multigrid for HDG.IMA Journal of Numerical Analysis, 42(4):3135–3153, 2022. 26
work page 2022
-
[25]
Peipei Lu, Andreas Rupp, and Guido Kanschat. Two-level Schwarz methods for hybridizable discon- tinuous Galerkin methods.Journal of Scientific Computing, 95(1):9, 2023
work page 2023
-
[26]
Peipei Lu, Wei Wang, Guido Kanschat, and Andreas Rupp. Homogeneous multigrid for HDG applied to the Stokes equation.IMA Journal of Numerical Analysis, 44(5):3124–3152, 2024
work page 2024
-
[27]
Kent-Andre Mardal and Ragnar Winther. Uniform preconditioners for the time dependent Stokes problem.Numerische Mathematik, 98(2):305–327, 2004
work page 2004
-
[28]
Kent-Andre Mardal and Ragnar Winther. Preconditioning discretizations of systems of partial differ- ential equations.Numerical Linear Algebra with Applications, 18(1):1–40, 2011
work page 2011
-
[29]
Maxim A Olshanskii, J¨ org Peters, and Arnold Reusken. Uniform preconditioners for a parameter dependent saddle point problem with application to generalized Stokes interface equations.Numerische Mathematik, 105(1):159–191, 2006
work page 2006
-
[30]
Sander Rhebergen and Garth N Wells. Analysis of a hybridized/interface stabilized finite element method for the Stokes equations.SIAM Journal on Numerical Analysis, 55(4):1982–2003, 2017
work page 1982
-
[31]
Sander Rhebergen and Garth N Wells. A hybridizable discontinuous Galerkin method for the Navier–Stokes equations with pointwise divergence-free velocity field.Journal of Scientific Comput- ing, 76(3):1484–1501, 2018
work page 2018
-
[32]
Sander Rhebergen and Garth N Wells. Preconditioning of a hybridized discontinuous Galerkin finite element method for the Stokes equations.Journal of Scientific Computing, 77(3):1936–1952, 2018
work page 1936
-
[33]
Sander Rhebergen and Garth N Wells. An embedded–hybridized discontinuous Galerkin finite element method for the Stokes equations.Computer Methods in Applied Mechanics and Engineering, 358:112619, 2020
work page 2020
-
[34]
Sander Rhebergen and Garth N Wells. Preconditioning for a pressure-robust HDG discretization of the Stokes equations.SIAM Journal on Scientific Computing, 44(1):A583–A604, 2022
work page 2022
-
[35]
Joachim Sch¨ oberl. NETGEN an advancing front 2D/3D-mesh generator based on abstract rules.Com- puting and visualization in science, 1(1):41–52, 1997
work page 1997
-
[36]
Joachim Sch¨ oberl. C++ 11 implementation of finite elements in NGSolve.Institute for analysis and scientific computing, Vienna University of Technology, 30, 2014
work page 2014
-
[37]
Joachim Sch¨ oberl and Christoph Lehrenfeld. Domain decomposition preconditioning for high order hybrid discontinuous Galerkin methods on tetrahedral meshes. InAdvanced finite element methods and applications, pages 27–56. Springer, 2013
work page 2013
-
[38]
PhD thesis, University of Waterloo, 2021
Abdullah Ali Sivas.Preconditioning of hybridizable discontinuous Galerkin discretizations of the Navier- Stokes equations. PhD thesis, University of Waterloo, 2021
work page 2021
-
[39]
Abdullah Ali Sivas, Ben S Southworth, and Sander Rhebergen. AIR algebraic multigrid for a space-time hybridizable discontinuous Galerkin discretization of advection (-diffusion).SIAM Journal on Scientific Computing, 43(5):A3393–A3416, 2021
work page 2021
-
[40]
Ben S Southworth, Abdullah A Sivas, and Sander Rhebergen. On fixed-point, krylov, and 2×2 block preconditioners for nonsymmetric problems.SIAM Journal on Matrix Analysis and Applications, 41(2):871–900, 2020
work page 2020
-
[41]
Xuemin Tu, Bin Wang, and Jinjin Zhang. Analysis of BDDC algorithms for Stokes problems with hybridizable discontinuous Galerkin discretizations.Electron. Trans. Numer. Anal., 52:553–570, 2020. 27
work page 2020
-
[42]
Xuemin Tu and Jinjin Zhang. BDDC algorithms for advection-diffusion problems with HDG discretiza- tions.Computers & Mathematics with Applications, 101:74–106, 2021
work page 2021
-
[43]
Garth N Wells. Analysis of an interface stabilized finite element method: the advection-diffusion- reaction equation.SIAM journal on numerical analysis, 49(1):87–109, 2011
work page 2011
-
[44]
Yi Yu, Maksymilian Dryja, and Marcus Sarkis. Nonoverlapping spectral additive Schwarz methods for hybrid discontinuous Galerkin discretizations.IMA Journal of Numerical Analysis, 44(1):192–224, 2024
work page 2024
-
[45]
Jinjin Zhang and Xuemin Tu. Robust BDDC algorithms for the Brinkman problem with HDG dis- cretizations.Computer Methods in Applied Mechanics and Engineering, 400:115548, 2022. A Useful results Lemma 14.There exists a positive uniform constant¯csuch that |||¯vh|||v,h ≤¯c|||vh|||v,1 ∀vh ∈V h.(77) Proof.The proof follows the same steps as the proof of [32, L...
work page 2022
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.