pith. sign in

arxiv: 2604.20018 · v1 · submitted 2026-04-21 · 🧮 math.NA · cs.NA

Pressure-Robust H(div)-Conforming HDG Methods for the Steady Stokes Equations with an Application to Tangential Boundary Control

Gang Chen , Wenyi Liu , Yangwen Zhang This is my paper

Pith reviewed 2026-05-10 01:12 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords HDG methodspressure-robustStokes equationsdivergence-freeH(div) conformingboundary controlhybridizable discontinuous Galerkinfinite element methods
0
0 comments X

The pith

H(div)-conforming HDG methods for Stokes equations produce exactly divergence-free velocities and pressure robustness with low pressure regularity.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper constructs a family of hybridizable discontinuous Galerkin methods that conform to the H(div) space for the steady Stokes equations, using BDM or RT velocity spaces paired with hybrid traces. The discrete velocities satisfy the divergence-free condition exactly inside each element. This property delivers pressure robustness, so that velocity accuracy holds independently of pressure approximation quality or viscosity size. Optimal convergence rates in energy and L2 norms are proved for the BDM variants, with corresponding results and weaker pressure bounds for the RT variants. The same framework supplies error estimates for a tangential boundary control problem and is confirmed by two- and three-dimensional tests.

Core claim

We develop pressure-robust H(div)-conforming HDG methods for the steady Stokes equations based on BDM and RT velocity spaces with either discontinuous or continuous hybrid traces. The discrete velocities are exactly divergence-free, which yields pressure robustness. The consistency argument requires only low pressure regularity rather than H1. For the BDM variants we derive optimal energy-norm estimates and optimal L2-velocity convergence, while for the RT variants we obtain optimal velocity convergence and weaker pressure estimates. We also prove a uniform spectral equivalence for the pressure Schur complement and apply the BDM discontinuous-trace scheme to obtain error estimates for the控制,

What carries the argument

H(div)-conforming HDG discretization using BDM or RT velocity spaces together with hybrid trace spaces, which enforces that the discrete velocity is exactly divergence-free elementwise.

If this is right

  • Optimal energy-norm and L2-velocity error estimates hold for the BDM variants.
  • Optimal velocity convergence with weaker pressure estimates holds for the RT variants.
  • The hybridized system admits a uniformly spectrally equivalent pressure Schur complement suitable for iterative solvers.
  • Error estimates for control, state, and adjoint variables are available for the tangential boundary control problem using the BDM discontinuous-trace scheme.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The exact mass conservation may improve long-time accuracy in incompressible flow computations where cumulative divergence errors matter.
  • The low-regularity pressure assumption suggests the methods could handle discontinuous or singular pressures more reliably than schemes requiring H1 pressure.
  • The framework could be tested on related div-constrained problems such as Darcy flow or magnetohydrodynamics to check whether similar robustness carries over.

Load-bearing premise

The BDM or RT spaces together with the chosen hybrid trace spaces must produce a consistent and stable discretization of the Stokes problem under only low pressure regularity.

What would settle it

On a fixed mesh, compute the elementwise divergence of the discrete velocity solution; if it is nonzero in any element, or if the velocity error grows as viscosity is driven toward zero, the exact divergence-free and pressure-robustness claims are falsified.

read the original abstract

We develop a family of $H(\mathrm{div})$-conforming hybridizable discontinuous Galerkin methods for the steady Stokes equations based on BDM and RT velocity spaces with either discontinuous or continuous hybrid traces. In contrast to our earlier pressure-robust HDG method for tangential boundary control, the present analysis does not require the pressure to belong to $H^1$; instead, the consistency argument only assumes low pressure regularity. The discrete velocities are exactly divergence-free, which yields pressure robustness. For the BDM variants we derive optimal energy-norm estimates and optimal $L^2$-velocity convergence, while for the RT variants we obtain optimal velocity convergence and weaker pressure estimates. We also analyze the hybridized linear system and prove a uniform spectral equivalence for the pressure Schur complement relevant to iterative solvers. As an application, we revisit the Stokes tangential boundary control problem and derive error estimates for the control, state, and adjoint variables using the BDM discontinuous-trace scheme. Two- and three-dimensional numerical experiments confirm the predicted convergence rates, the exact divergence-free property, and the robustness of the method with respect to the viscosity parameter.

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

0 major / 3 minor

Summary. The manuscript develops a family of H(div)-conforming hybridizable discontinuous Galerkin (HDG) methods for the steady Stokes equations based on BDM and RT velocity spaces with either discontinuous or continuous hybrid traces. The discrete velocities are exactly divergence-free, yielding pressure robustness. The consistency argument requires only low pressure regularity (not H^1). Optimal energy-norm and L^2-velocity estimates are derived for BDM variants; RT variants obtain optimal velocity convergence with weaker pressure estimates. The hybridized system is analyzed with a uniform spectral equivalence result for the pressure Schur complement. Error estimates are derived for the tangential boundary control problem using the BDM discontinuous-trace scheme. Two- and three-dimensional numerical experiments confirm the predicted rates, exact divergence-free property, and viscosity robustness.

Significance. If the central claims hold, the work makes a solid contribution to pressure-robust discretizations for incompressible flows by relaxing the pressure regularity assumption relative to earlier HDG approaches while preserving the exact divergence-free property. The Schur-complement analysis supports practical iterative solvers, and the boundary-control application demonstrates utility beyond the model problem. Numerical confirmation of rates and robustness properties is a positive feature.

minor comments (3)
  1. The abstract and introduction could more explicitly contrast the low-regularity consistency argument with the H^1 assumption in the authors' prior tangential-control HDG paper to highlight the technical advance.
  2. In the numerical experiments, include a brief statement on how the stabilization parameters are chosen (e.g., scaling with mesh size or viscosity) to facilitate reproducibility of the reported rates.
  3. Figure captions would benefit from indicating the polynomial degree and the precise norm plotted (energy vs. L2) so that readers can immediately match the plots to the theorems.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive evaluation of our manuscript and the recommendation for minor revision. The referee summary correctly identifies the main contributions, including the relaxed pressure regularity assumption in the consistency analysis, the exact divergence-free property of the discrete velocities, the optimal convergence results for the BDM variants, the Schur complement analysis, and the application to the tangential boundary control problem.

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained via new consistency and approximation arguments

full rationale

The paper develops new H(div)-conforming HDG schemes for Stokes using BDM/RT spaces, with the exact divergence-free property following directly from the local divergence constraint and test-function choice in the hybridized formulation. Consistency holds under low pressure regularity (no H^1 assumption), and error estimates (energy-norm, L2-velocity) are obtained from standard approximation properties of the discrete spaces. The hybridized system analysis (Schur complement spectral equivalence) and control estimates are likewise derived independently. The single self-citation to prior work serves only for contrast and to revisit the boundary-control application; it is not load-bearing for the central claims or estimates. No self-definitional reductions, fitted inputs renamed as predictions, or ansatz smuggling occur.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Abstract-only review yields no explicit free parameters, invented entities, or ad-hoc axioms beyond standard finite-element assumptions; the method relies on established BDM/RT spaces and Sobolev-space theory for Stokes.

axioms (1)
  • standard math Standard Sobolev regularity and finite-element approximation properties for BDM and RT spaces on simplicial meshes
    Invoked implicitly for consistency, stability, and error estimates of the HDG discretization.

pith-pipeline@v0.9.0 · 5513 in / 1415 out tokens · 41696 ms · 2026-05-10T01:12:01.786900+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

38 extracted references · 38 canonical work pages

  1. [1]

    N. Arada. Optimal control of shear-thinning fluids.SIAM J. Control Optim.50(4):2515–2542, 2012

    work page 2012

  2. [2]

    G. R. Barrenechea, M. Bosy, V. Dolean, F. Nataf, and P.-H. Tournier. Hybrid discontinuous Galerkin discretisation and domain decomposition preconditioners for the Stokes problem. Comput. Methods Appl. Math.19(4):703–722, 2019

    work page 2019

  3. [3]

    Banas, M

    L. Banas, M. Klein, and A. Prohl. Control of interface evolution in multiphase fluid flows.SIAM J. Control Optim.52(4):2284–2318, 2014

    work page 2014

  4. [4]

    Boffi, F

    D. Boffi, F. Brezzi, L. F. Demkowicz, R. G. Durán, R. S. Falk, and M. Fortin.Mixed Finite Elements, Compatibility Conditions, and Applications. Lecture Notes in Mathematics 1939, Springer, Berlin, 2008

    work page 1939

  5. [5]

    Boffi, F

    D. Boffi, F. Brezzi, and M. Fortin.Mixed Finite Element Methods and Applications. Springer Series in Computational Mathematics 44, Springer, Heidelberg, 2013

    work page 2013

  6. [6]

    Z. Cai, X. Ye, and S. Zhang. Discontinuous Galerkin finite element methods for interface problems: a priori and a posteriori error estimations.SIAM J. Numer. Anal.49(5):1761–1787, 2011. 34 Pressure-RobustH(div)-Conforming HDG Methods

    work page 2011

  7. [7]

    G. Chen, M. Feng, and X. Xie. Robust globally divergence-free weak Galerkin methods for Stokes equations.J. Comput. Math.34(5):549–572, 2016

    work page 2016

  8. [8]

    G. Chen, W. Gong, M. Mateos, J. R. Singler, and Y. Zhang. A new global divergence-free and pressure-robust HDG method for tangential boundary control of Stokes equations.Comput. Methods Appl. Mech. Engrg.405:115837, 2023

    work page 2023

  9. [9]

    G. Chen, J. R. Singler, and Y. Zhang. An HDG method for Dirichlet boundary control of convection dominated diffusion PDEs.SIAM J. Numer. Anal.57(4):1919–1946, 2019

    work page 1919

  10. [10]

    Chen and X

    G. Chen and X. Xie. Analysis of a class of globally divergence-free HDG methods for stationary Navier–Stokes equations.Sci. China Math., 2023

    work page 2023

  11. [11]

    Cockburn and F.-J

    B. Cockburn and F.-J. Sayas. Divergence-conforming HDG methods for Stokes flows.Math. Comp.83(288):1571–1598, 2014

    work page 2014

  12. [12]

    M. Dauge. Stationary Stokes and Navier–Stokes systems on two- or three-dimensional domains with corners. I. Linearized equations.SIAM J. Math. Anal.20(1):74–97, 1989

    work page 1989

  13. [13]

    D. A. Di Pietro and A. Ern.Mathematical Aspects of Discontinuous Galerkin Methods. Mathé- matiques & Applications 69, Springer, Heidelberg, 2012

    work page 2012

  14. [14]

    Du and F.-J

    S. Du and F.-J. Sayas.An Invitation to the Theory of the Hybridizable Discontinuous Galerkin Method. SpringerBriefs in Mathematics, Springer, Cham, 2019

    work page 2019

  15. [15]

    Ern and J.-L

    A. Ern and J.-L. Guermond.Theory and Practice of Finite Elements. Applied Mathematical Sciences 159, Springer, New York, 2004

    work page 2004

  16. [16]

    R. D. Falgout and U. M. Yang. hypre: A library of high performance preconditioners. In P. M. A. Sloot, A. G. Hoekstra, C. J. K. Tan, and J. J. Dongarra, editors,Computational Science — ICCS 2002, pages 632–641. Springer, Berlin, Heidelberg, 2002

    work page 2002

  17. [17]

    G. Fu, Y. Jin, and W. Qiu. Parameter-free superconvergentH(div)-conforming HDG methods for the Brinkman equations.IMA J. Numer. Anal.39(2):957–982, 2019

    work page 2019

  18. [18]

    Girault and P.-A

    V. Girault and P.-A. Raviart.Finite Element Methods for Navier–Stokes Equations: Theory and Algorithms. Springer Series in Computational Mathematics 5, Springer, Berlin, 1986

    work page 1986

  19. [20]

    W. Gong, W. Hu, M. Mateos, J. Singler, X. Zhang, and Y. Zhang. A new HDG method for Dirichlet boundary control of convection diffusion PDEs II: low regularity.SIAM J. Numer. Anal.56(4):2262–2287, 2018

    work page 2018

  20. [21]

    W. Gong, W. Hu, M. Mateos, J. R. Singler, and Y. Zhang. Analysis of a hybridizable discon- tinuous Galerkin scheme for the tangential control of the Stokes system.ESAIM Math. Model. Numer. Anal.54(6):2229–2264, 2020

    work page 2020

  21. [22]

    W. Gong, M. Mateos, J. Singler, and Y. Zhang. Analysis and approximations of Dirichlet boundary control of Stokes flows in the energy space.SIAM J. Numer. Anal.60(1):450–474, 2022. 35 Gang Chen, Wenyi Liu, Yangwen Zhang

    work page 2022

  22. [23]

    Guennebaud, B

    G. Guennebaud, B. Jacob, et al. Eigen v3.http://eigen.tuxfamily.org, 2010

    work page 2010

  23. [24]

    Hintermüller, T

    M. Hintermüller, T. Keil, and D. Wegner. Optimal control of a semidiscrete Cahn–Hilliard– Navier–Stokes system with nonmatched fluid densities.SIAM J. Control Optim.55(3):1954–1989, 2017

    work page 1954

  24. [25]

    Hinze and K

    M. Hinze and K. Kunisch. Second order methods for optimal control of time-dependent fluid flow.SIAM J. Control Optim.40(3):925–946, 2001

    work page 2001

  25. [26]

    W. Hu, J. Shen, J. R. Singler, Y. Zhang, and X. Zheng. A superconvergent hybridizable discontinuous Galerkin method for Dirichlet boundary control of elliptic PDEs.Numer. Math. 144(2):375–411, 2020

    work page 2020

  26. [27]

    R. B. Kellogg and J. E. Osborn. A regularity result for the Stokes problem in a convex polygon. J. Functional Analysis21(4):397–431, 1976

    work page 1976

  27. [28]

    Kuberry and H

    P. Kuberry and H. Lee. Analysis of a fluid-structure interaction problem recast in an optimal control setting.SIAM J. Numer. Anal.53(3):1464–1487, 2015

    work page 2015

  28. [29]

    P. L. Lederer, C. Lehrenfeld, and J. Schöberl. Hybrid discontinuous Galerkin methods with relaxed H(div)-conformity for incompressible flows. Part I.SIAM J. Numer. Anal.56(4):2070– 2094, 2018

    work page 2070

  29. [30]

    Lehrenfeld and J

    C. Lehrenfeld and J. Schöberl. High order exactly divergence-free hybrid discontinuous Galerkin methods for unsteady incompressible flows.Comput. Methods Appl. Mech. Engrg.307:339–361, 2016

    work page 2016

  30. [31]

    Mathew, I

    G. Mathew, I. Mezic, S. Grivopoulos, U. Vaidya, and L. Petzold. Optimal control of mixing in Stokes fluid flows.J. Fluid Mech.580:261–281, 2007

    work page 2007

  31. [32]

    J.-C. Nédélec. A new family of mixed finite elements inR3.Numer. Math.50(1):57–81, 1986

    work page 1986

  32. [33]

    Raviart and J

    P.-A. Raviart and J. M. Thomas. A mixed finite element method for 2nd order elliptic problems. InMathematical Aspects of Finite Element Methods, Lecture Notes in Mathematics 606, pages 292–315. Springer, Berlin, 1977

    work page 1977

  33. [34]

    Rhebergen and G

    S. Rhebergen and G. N. Wells. A hybridizable discontinuous Galerkin method for the Navier– Stokes equations with pointwise divergence-free velocity field.J. Sci. Comput.76(3):1484–1501, 2018

    work page 2018

  34. [35]

    Tachim Medjo

    T. Tachim Medjo. Iterative methods for a class of robust control problems in fluid mechanics. SIAM J. Numer. Anal.39(5):1625–1647, 2001/02

    work page 2001

  35. [36]

    Tachim Medjo and L

    T. Tachim Medjo and L. R. Tcheugoue Tebou. Adjoint-based iterative method for robust control problems in fluid mechanics.SIAM J. Numer. Anal.42(1):302–325, 2004

    work page 2004

  36. [37]

    P. You, Z. Ding, and J. Zhou. Optimal boundary control of the Stokes fluids with point velocity observations.SIAM J. Control Optim.36(3):981–1004, 1998

    work page 1998

  37. [38]

    Zheng, G

    X. Zheng, G. Chen, and X. Xie. A divergence-free weak Galerkin method for quasi-Newtonian Stokes flows.Sci. China Math.60(8):1515–1528, 2017

    work page 2017

  38. [39]

    Cockburn and F.-J

    B. Cockburn and F.-J. Sayas. Divergence-conforming HDG methods for Stokes flows.Math. Comp.83(288):1571–1598, 2014. 36

    work page 2014