pith. sign in

arxiv: 2604.07202 · v1 · submitted 2026-04-08 · 🧮 math.NA · cs.NA

Robust preconditioning for an HDG discretization of the time-dependent Stokes equations

Pith reviewed 2026-05-10 17:50 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords discretizationcondensationpreconditionersrobuststaticstokestheoreticaltime-dependent
0
0 comments X

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.

The paper aims to construct preconditioners for the linear systems that result from statically condensing an HDG discretization of the time-dependent Stokes problem. These preconditioners are designed to remain effective no matter the values of physical parameters like viscosity or time step size, and discretization parameters like mesh size or polynomial degree. A sympathetic reader would care because efficient and reliable solvers are essential for simulating viscous flows in engineering and science, where changing parameters often cause standard methods to slow down dramatically. The work builds on prior analysis to prove uniform well-posedness and then gives conditions on a norm for the face unknowns that ensure the robustness carries over after condensation.

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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

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)
  1. [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.
  2. [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)
  1. [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.
  2. [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

2 responses · 0 unresolved

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
  1. 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

  2. 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

0 steps flagged

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

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on standard finite-element well-posedness theory and the authors' prior HDG analysis; no free parameters or invented entities are introduced in the abstract.

axioms (1)
  • domain assumption Uniform well-posedness of the HDG formulation before static condensation holds in appropriately defined norms
    Invoked to identify sufficient conditions on the face-space norm for preconditioner robustness.

pith-pipeline@v0.9.0 · 5445 in / 1238 out tokens · 23023 ms · 2026-05-10T17:50:12.660465+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Preconditioning of a hybridizable discontinuous Galerkin method for the coupled Stokes--Darcy system

    math.NA 2026-04 unverdicted novelty 5.0

    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

45 extracted references · 45 canonical work pages · cited by 1 Pith paper

  1. [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

  2. [2]

    An observation on the uniform preconditioners for the mixed Darcy problem.Numerical Methods for Partial Differential Equations, 36(6):1718–1734, 2020

    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

  3. [3]

    Springer, 2013

    Daniele Boffi, Franco Brezzi, and Michel Fortin.Mixed finite element methods and applications, vol- ume 44. Springer, 2013

  4. [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

  5. [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

  6. [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

  7. [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

  8. [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

  9. [9]

    Springer Nature, 2019

    Shukai Du and Francisco-Javier Sayas.An Invitation to the Theory of the Hybridizable Discontinuous Galerkin Method: Projections, Estimates, Tools. Springer Nature, 2019

  10. [10]

    Equivalence of local- and global- best approximations, a simple stable local commuting projector, and optimalhpapproximation esti- mates inH(div).IMA J

    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

  11. [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

  12. [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

  13. [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

  14. [14]

    Holst, and Martin W

    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

  15. [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

  16. [16]

    Local Fourier analysis of multigrid for hybridized and embedded discontinuous Galerkin methods.SIAM Journal on Scientific Computing, 43(5):S612– S636, 2021

    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

  17. [17]

    Lee, and Sander Rhebergen

    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

  18. [18]

    Preconditioning of a hybridizable discon- tinuous Galerkin method for Biot’s consolidation model.arXiv preprint arXiv:2508.12991, 2025

    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. [19]

    SIAM, 2002

    Nicholas J Higham.Accuracy and stability of numerical algorithms. SIAM, 2002

  20. [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

  21. [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

  22. [22]

    Preconditioning a hybridizable discon- tinuous Galerkin method for Navier-Stokes at high Reynolds number.arXiv preprint arXiv:2512.02971, 2025

    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. [23]

    Analysis of injection operators in geometric multigrid solvers for HDG methods.SIAM Journal on Numerical Analysis, 60(4):2293–2317, 2022

    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

  24. [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

  25. [25]

    Two-level Schwarz methods for hybridizable discon- tinuous Galerkin methods.Journal of Scientific Computing, 95(1):9, 2023

    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

  26. [26]

    Homogeneous multigrid for HDG applied to the Stokes equation.IMA Journal of Numerical Analysis, 44(5):3124–3152, 2024

    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

  27. [27]

    Uniform preconditioners for the time dependent Stokes problem.Numerische Mathematik, 98(2):305–327, 2004

    Kent-Andre Mardal and Ragnar Winther. Uniform preconditioners for the time dependent Stokes problem.Numerische Mathematik, 98(2):305–327, 2004

  28. [28]

    Preconditioning discretizations of systems of partial differ- ential equations.Numerical Linear Algebra with Applications, 18(1):1–40, 2011

    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

  29. [29]

    Uniform preconditioners for a parameter dependent saddle point problem with application to generalized Stokes interface equations.Numerische Mathematik, 105(1):159–191, 2006

    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

  30. [30]

    Analysis of a hybridized/interface stabilized finite element method for the Stokes equations.SIAM Journal on Numerical Analysis, 55(4):1982–2003, 2017

    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

  31. [31]

    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

    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

  32. [32]

    Preconditioning of a hybridized discontinuous Galerkin finite element method for the Stokes equations.Journal of Scientific Computing, 77(3):1936–1952, 2018

    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

  33. [33]

    An embedded–hybridized discontinuous Galerkin finite element method for the Stokes equations.Computer Methods in Applied Mechanics and Engineering, 358:112619, 2020

    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

  34. [34]

    Preconditioning for a pressure-robust HDG discretization of the Stokes equations.SIAM Journal on Scientific Computing, 44(1):A583–A604, 2022

    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

  35. [35]

    NETGEN an advancing front 2D/3D-mesh generator based on abstract rules.Com- puting and visualization in science, 1(1):41–52, 1997

    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

  36. [36]

    C++ 11 implementation of finite elements in NGSolve.Institute for analysis and scientific computing, Vienna University of Technology, 30, 2014

    Joachim Sch¨ oberl. C++ 11 implementation of finite elements in NGSolve.Institute for analysis and scientific computing, Vienna University of Technology, 30, 2014

  37. [37]

    Domain decomposition preconditioning for high order hybrid discontinuous Galerkin methods on tetrahedral meshes

    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

  38. [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

  39. [39]

    AIR algebraic multigrid for a space-time hybridizable discontinuous Galerkin discretization of advection (-diffusion).SIAM Journal on Scientific Computing, 43(5):A3393–A3416, 2021

    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

  40. [40]

    On fixed-point, krylov, and 2×2 block preconditioners for nonsymmetric problems.SIAM Journal on Matrix Analysis and Applications, 41(2):871–900, 2020

    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

  41. [41]

    Analysis of BDDC algorithms for Stokes problems with hybridizable discontinuous Galerkin discretizations.Electron

    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

  42. [42]

    BDDC algorithms for advection-diffusion problems with HDG discretiza- tions.Computers & Mathematics with Applications, 101:74–106, 2021

    Xuemin Tu and Jinjin Zhang. BDDC algorithms for advection-diffusion problems with HDG discretiza- tions.Computers & Mathematics with Applications, 101:74–106, 2021

  43. [43]

    Analysis of an interface stabilized finite element method: the advection-diffusion- reaction equation.SIAM journal on numerical analysis, 49(1):87–109, 2011

    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

  44. [44]

    Nonoverlapping spectral additive Schwarz methods for hybrid discontinuous Galerkin discretizations.IMA Journal of Numerical Analysis, 44(1):192–224, 2024

    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

  45. [45]

    Robust BDDC algorithms for the Brinkman problem with HDG dis- cretizations.Computer Methods in Applied Mechanics and Engineering, 400:115548, 2022

    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...