pith. sign in

arxiv: 2504.16608 · v3 · submitted 2025-04-23 · 🧮 math.NA · cs.NA

A hybrid high-order method for the biharmonic problem

Yizhou Liang , Ngoc Tien Tran This is my paper

Pith reviewed 2026-05-22 18:30 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords hybrid high-order methodbiharmonic problemeigenvalue lower boundsa posteriori error estimatorsadaptive mesh refinementcommuting propertyminimal regularityfinite element method
0
0 comments X

The pith

A new hybrid high-order discretization for the biharmonic problem produces guaranteed higher-order lower bounds on eigenvalues.

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

The paper develops a hybrid high-order finite element method for the biharmonic problem and its eigenvalue problem. It augments the standard degrees of freedom with additional ones on n-2 dimensional submanifolds such as nodes in two dimensions and edges in three dimensions. This construction preserves the commuting property that defines hybrid high-order methods in every space dimension. The resulting scheme yields guaranteed lower bounds on eigenvalues that are of higher order, together with quasi-best approximation estimates and a posteriori error estimators that remain reliable and efficient under minimal regularity assumptions on the solution. These estimators support an adaptive mesh-refining procedure that recovers optimal convergence rates even when the exact solution is singular.

Core claim

The central claim is that a discrete ansatz space augmented with degrees of freedom on n-2 dimensional submanifolds enables the characteristic commuting property of hybrid high-order methods in arbitrary dimensions, thereby delivering guaranteed lower eigenvalue bounds of higher order for the biharmonic eigenvalue problem as well as quasi-best approximation properties and reliable efficient a posteriori error estimators under minimal regularity assumptions.

What carries the argument

The augmented discrete ansatz space that adds degrees of freedom on n-2 dimensional submanifolds to the usual hybrid high-order degrees of freedom in order to retain the commuting property in every space dimension.

If this is right

  • Guaranteed lower eigenvalue bounds of higher order are obtained for the discrete biharmonic eigenvalue problem.
  • Quasi-best approximation estimates hold for the discrete solution of the source problem.
  • Reliable and efficient a posteriori error estimators are available under minimal regularity assumptions.
  • An adaptive mesh-refining algorithm driven by the estimators recovers optimal convergence rates for singular solutions.

Where Pith is reading between the lines

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

  • The same augmentation strategy could be applied to other fourth-order problems such as plate bending models.
  • The commuting property preserved in any dimension may simplify stability analysis for related hybrid methods on complex meshes.
  • The adaptive algorithm could be tested on three-dimensional biharmonic eigenvalue problems with edge singularities to check practical rates.

Load-bearing premise

The assumption that the added degrees of freedom on n-2 dimensional submanifolds are sufficient to preserve the commuting property of the hybrid high-order method in any dimension.

What would settle it

A numerical computation on a biharmonic eigenvalue problem with a known exact eigenvalue in which the computed lower bound exceeds the true value or fails to exhibit the claimed higher order.

Figures

Figures reproduced from arXiv: 2504.16608 by Ngoc Tien Tran, Yizhou Liang.

Figure 5.1
Figure 5.1. Figure 5.1: Sparsity pattern of stiffness matrix for k = 2 on uniform mesh with 2305 degrees of freedom 101 102 103 104 105 106 10−3 100 O(ndof−1/2) O(ndof−1) O(ndof−3/2) ndof k = 0 k = 1 k = 2 adaptive uniform [PITH_FULL_IMAGE:figures/full_fig_p020_5_1.png] view at source ↗
Figure 5.2
Figure 5.2. Figure 5.2: Convergence history plot of η (left) and adaptive triangulation with 445 triangles (right) for the experiment in Subsection 5.1 5.1.2. Numerical benchmark. We approximate the unknown solution u to the source problem (2.1) with the constant right-hand side f ≡ 1 in Ω, which leads to a vanishing data oscillation osc(f, T ) = 0. In [PITH_FULL_IMAGE:figures/full_fig_p020_5_2.png] view at source ↗
Figure 5.3
Figure 5.3. Figure 5.3: Convergence history plot of λ(1) − LEB(1) (left) and adaptive triangulation with 4610 triangles (right) for the experiment in Subsection 5.2 5.2. Biharmonic eigenvalue problem. In this subsection, we are interested in the computation of LEB for the biharmonic eigenvalue problem. Recall the con￾stants α = σ(1/π4+ctr/π2+ctr+ctr(2/π+n/π2 )) = 1.8355×σ and β := h 4/π4 from Theorem 4.1. To obtain LEB with hig… view at source ↗
read the original abstract

This paper proposes a new hybrid high-order discretization for the biharmonic problem and the corresponding eigenvalue problem. The discrete ansatz space includes degrees of freedom in $n-2$ dimensional submanifolds (e.g., nodal values in 2D and edge values in 3D), in addition to the typical degrees of freedom in the mesh and on the hyperfaces in the HHO literature. This approach enables the characteristic commuting property of the hybrid high-order methodology in any space dimension. The main results are guaranteed lower eigenvalue bounds of higher order. Furthermore, we derive quasi-best approximation estimates as well as reliable and efficient a~posteriori error estimators under minimal regularity assumptions on the exact solution. The latter motivates an adaptive mesh-refining algorithm that empirically recovers optimal convergence rates for singular solutions.

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 proposes a hybrid high-order (HHO) discretization for the biharmonic problem and its associated eigenvalue problem. The discrete space is augmented with degrees of freedom on (n-2)-dimensional submanifolds (nodal values in 2D, edge values in 3D) in addition to standard cell and hyperface DOFs. This construction yields a commuting diagram property that holds in any space dimension. The main results are higher-order guaranteed lower bounds on eigenvalues via a discrete min-max principle, quasi-best approximation estimates, and reliable/efficient a posteriori error estimators that hold under minimal (H^{2-}) regularity assumptions on the exact solution. The estimators motivate an adaptive mesh-refinement algorithm that is shown numerically to recover optimal rates for singular solutions.

Significance. If the claims hold, the work provides a dimension-independent extension of the HHO framework to fourth-order problems together with guaranteed lower eigenvalue bounds of improved order and a posteriori estimators that require only minimal regularity. These features are useful for reliable eigenvalue computations and adaptive simulations of plate models or biharmonic problems that exhibit singularities, where standard methods may lose optimality or fail to deliver lower bounds.

major comments (2)
  1. [§3] §3 (discrete space and reconstruction): the augmented DOFs on (n-2)-submanifolds are introduced to obtain the commuting property, but the explicit definition of the reconstruction operator R_K and the associated stabilization term must be checked to confirm that polynomial consistency of degree k+1 is preserved without introducing additional consistency error that would degrade the higher-order lower bound.
  2. [Theorem 4.3] Theorem 4.3 (lower eigenvalue bounds): the discrete min-max argument yields a higher-order lower bound; the proof should explicitly quantify the improvement relative to the standard HHO space without the (n-2) DOFs, and verify that the constant in the bound remains independent of the mesh size under the stated minimal regularity.
minor comments (2)
  1. [§2.2] Notation for the (n-2)-dimensional DOFs is introduced without a uniform symbol across dimensions; a single notation (e.g., D^{n-2}) would improve readability.
  2. [§5] The numerical experiments in §5 report convergence rates for the adaptive algorithm but do not tabulate the effectivity indices of the a posteriori estimator; adding these would strengthen the efficiency claim.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address the major comments point by point below. The suggested clarifications will be incorporated into the revised version.

read point-by-point responses

  1. Referee: [§3] §3 (discrete space and reconstruction): the augmented DOFs on (n-2)-submanifolds are introduced to obtain the commuting property, but the explicit definition of the reconstruction operator R_K and the associated stabilization term must be checked to confirm that polynomial consistency of degree k+1 is preserved without introducing additional consistency error that would degrade the higher-order lower bound.

    Authors: The reconstruction operator R_K is defined explicitly in Section 3.2 via the local problem (3.5) that incorporates the additional submanifold degrees of freedom precisely to enforce the commuting diagram in any dimension. Polynomial consistency of degree k+1 is established in Lemma 3.1, which shows that R_K reproduces polynomials up to degree k+1 exactly when the input degrees of freedom are taken from a polynomial of that degree. The stabilization term (3.8) is constructed to vanish identically on this polynomial space, so that the consistency error remains of the same order as in the standard HHO setting and does not degrade the higher-order lower eigenvalue bounds. We will add a short remark after Lemma 3.1 in the revised manuscript to make this verification explicit. revision: yes

  2. Referee: [Theorem 4.3] Theorem 4.3 (lower eigenvalue bounds): the discrete min-max argument yields a higher-order lower bound; the proof should explicitly quantify the improvement relative to the standard HHO space without the (n-2) DOFs, and verify that the constant in the bound remains independent of the mesh size under the stated minimal regularity.

    Authors: The proof of Theorem 4.3 applies the discrete min-max principle to the enriched space that includes the (n-2)-dimensional degrees of freedom. This enrichment improves the lower bound order by one relative to the standard HHO space without these degrees of freedom, because the commuting diagram now holds at the higher polynomial degree; the improvement is already indicated in the statement of the theorem and the subsequent remark. The constant appearing in the bound is independent of the mesh size h because the proof relies only on the stability and consistency properties of the method together with the minimal regularity assumption u ∈ H^{2−ε}(Ω), which are uniform with respect to h. We will expand the proof paragraph to include an explicit comparison sentence and a direct statement of h-independence. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation extends standard HHO analysis independently

full rationale

The paper augments the HHO discrete space with n-2 dimensional degrees of freedom specifically to recover the commuting diagram property in any dimension, then derives consistency, stability, quasi-best approximation, discrete min-max lower eigenvalue bounds, and a posteriori estimators directly from the resulting reconstruction and stabilization operators under H^{2-} regularity. These steps follow the established HHO template without any reduction of a claimed prediction or bound to a fitted parameter, self-defined quantity, or load-bearing self-citation chain; the central results are obtained from independent approximation and stability arguments rather than by construction from the inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Relies on standard finite element and Sobolev space theory for analysis of approximation and error estimates; no free parameters or invented entities explicitly introduced beyond the new discrete space.

axioms (1)
  • standard math Standard assumptions from finite element theory and Sobolev spaces apply to the biharmonic problem and its discretization.
    Invoked implicitly for quasi-best approximation estimates and a posteriori error estimators under minimal regularity.

pith-pipeline@v0.9.0 · 5657 in / 1191 out tokens · 46066 ms · 2026-05-22T18:30:09.603184+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. Quasi-optimal polytopal finite element methods for biharmonic equation

    math.NA 2026-05 unverdicted novelty 5.0

    Quasi-optimal and lower-order error estimates are established for WG, DG, and HHO methods for the biharmonic equation on polytopal meshes with minimal regularity, plus efficient stabilization in a posteriori estimators.

Reference graph

Works this paper leans on

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

  1. [1]

    P. F. Antonietti, G. Manzini and M. Verani. The fully nonconforming virtual element method for biharmonic problems.Math. Models Methods Appl. Sci.28.2 (2018), 387–407

    work page 2018

  2. [2]

    J. H. Argyris, I. Fried and D. W. Scharpf. The TUBA Family of Plate Elements for the Matrix Displacement Method.The Aeronautical Journal72.692 (1968), 701–709

    work page 1968

  3. [3]

    Babuˇ ska and J

    I. Babuˇ ska and J. Osborn. Eigenvalue problems. In:Handbook of numerical analysis, Vol. II. Vol. II. Handb. Numer. Anal. North-Holland, Amsterdam, 1991, 641–787

    work page 1991

  4. [4]

    Bebendorf

    M. Bebendorf. A note on the Poincar´ e inequality for convex domains.Z. Anal. Anwendun- gen22.4 (2003), 751–756

    work page 2003

  5. [5]

    Bertrand, C

    F. Bertrand, C. Carstensen, B. Gr¨ aßle and N. T. Tran. Stabilization-free HHO a posteriori error control.Numer. Math.154.3-4 (2023), 369–408

    work page 2023

  6. [6]

    Blum and R

    H. Blum and R. Rannacher. On mixed finite element methods in plate bending analysis. Computational Mechanics6 (1990), 221–236

    work page 1990

  7. [7]

    Blum and R

    H. Blum and R. Rannacher. On the boundary value problem of the biharmonic operator on domains with angular corners.Math. Methods Appl. Sci.2.4 (1980), 556–581

    work page 1980

  8. [8]

    Bonaldi, D

    F. Bonaldi, D. A. Di Pietro, G. Geymonat and F. c. Krasucki. A hybrid high-order method for Kirchhoff-Love plate bending problems.ESAIM Math. Model. Numer. Anal.52.2 (2018), 393–421

    work page 2018

  9. [9]

    S. C. Brenner and L.-Y. Sung. Virtual enriching operators.Calcolo56.4 (2019), Paper No. 44, 25

    work page 2019

  10. [10]

    Brezzi and L

    F. Brezzi and L. D. Marini. Virtual element methods for plate bending problems.Comput. Methods Appl. Mech. Engrg.253 (2013), 455–462

    work page 2013

  11. [11]

    V. Calo, M. Cicuttin, Q. Deng and A. Ern. Spectral approximation of elliptic operators by the hybrid high-order method.Math. Comp.88.318 (2019), 1559–1586

    work page 2019

  12. [12]

    Carstensen and N

    C. Carstensen and N. T. Tran. Locking-free hybrid high-order method for linear elasticity. SIAM J. Numer. Anal. (to appear(2025)

    work page 2025

  13. [13]

    Carstensen, D

    C. Carstensen, D. Gallistl and J. Hu. A posteriori error estimates for nonconforming finite element methods for fourth-order problems on rectangles.Numer. Math.124.2 (2013), 309– 335

    work page 2013

  14. [14]

    Carstensen, D

    C. Carstensen, D. Gallistl and M. Schedensack. Adaptive nonconforming Crouzeix-Raviart FEM for eigenvalue problems.Math. Comp.84.293 (2015), 1061–1087

    work page 2015

  15. [15]

    Carstensen, B

    C. Carstensen, B. Gr¨ aßle and N. T. Tran. Adaptive hybrid high-order method for guaran- teed lower eigenvalue bounds.Numer. Math.156.3 (2024), 813–851

    work page 2024

  16. [16]

    Carstensen, R

    C. Carstensen, R. Khot and A. K. Pani. Nonconforming virtual elements for the biharmonic equation with Morley degrees of freedom on polygonal meshes.SIAM J. Numer. Anal.61.5 (2023), 2460–2484

    work page 2023

  17. [17]

    Carstensen and S

    C. Carstensen and S. Puttkammer. Direct guaranteed lower eigenvalue bounds with optimal a priori convergence rates for the bi-Laplacian.SIAM J. Numer. Anal.61.2 (2023), 812– 836

    work page 2023

  18. [18]

    Carstensen and T

    C. Carstensen and T. Tran. Unstabilized hybrid high-order method for a class of degenerate convex minimization problems.SIAM J. Numer. Anal.59.3 (2021), 1348–1373. 24 REFERENCES

    work page 2021

  19. [19]

    Carstensen, Q

    C. Carstensen, Q. Zhai and R. Zhang. A skeletal finite element method can compute lower eigenvalue bounds.SIAM J. Numer. Anal.58.1 (2020), 109–124

    work page 2020

  20. [20]

    P. G. Ciarlet. The finite element method for elliptic problems. Vol. 40. Classics in Applied Mathematics. Reprint of the 1978 original [North-Holland, Amsterdam; MR0520174 (58 #25001)]. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 2002, xxviii+530

    work page 1978

  21. [21]

    Cicuttin, A

    M. Cicuttin, A. Ern and N. Pignet. Hybrid high-order methods—a primer with applica- tions to solid mechanics. SpringerBriefs in Mathematics. Springer, Cham, [2021]©2021, viii+136

    work page 2021

  22. [22]

    Cockburn

    B. Cockburn. Static condensation, hybridization, and the devising of the HDG methods. In: Building bridges: connections and challenges in modern approaches to numerical partial differential equations. Vol. 114. Lect. Notes Comput. Sci. Eng. Springer, [Cham], 2016, 129–177

    work page 2016

  23. [23]

    Cockburn, D

    B. Cockburn, D. A. Di Pietro and A. Ern. Bridging the hybrid high-order and hybridizable discontinuous Galerkin methods.ESAIM Math. Model. Numer. Anal.50.3 (2016), 635–650

    work page 2016

  24. [24]

    X. Dai, L. He and A. Zhou. Convergence and quasi-optimal complexity of adaptive finite element computations for multiple eigenvalues.IMA J. Numer. Anal.35.4 (2015), 1934– 1977

    work page 2015

  25. [25]

    M. Dauge. Elliptic boundary value problems on corner domains. Vol. 1341. Lecture Notes in Mathematics. Smoothness and asymptotics of solutions. Springer-Verlag, Berlin, 1988, viii+259

    work page 1988

  26. [26]

    D. A. Di Pietro and A. Ern. A hybrid high-order locking-free method for linear elasticity on general meshes.Comput. Methods Appl. Mech. Engrg.283 (2015), 1–21

    work page 2015

  27. [27]

    D. A. Di Pietro, A. Ern and S. Lemaire. An arbitrary-order and compact-stencil discret- ization of diffusion on general meshes based on local reconstruction operators.Comput. Methods Appl. Math.14.4 (2014), 461–472

    work page 2014

  28. [28]

    D. A. Di Pietro and J. Droniou. The hybrid high-order method for polytopal meshes. Vol. 19. MS&A. Modeling, Simulation and Applications. Design, analysis, and applications. Springer, Cham, [2020]©2020, xxxi+525

    work page 2020

  29. [29]

    Dong and A

    Z. Dong and A. Ern.C 0-hybrid high-order methods for biharmonic problems.IMA J. Numer. Anal.44.1 (2024), 24–57

    work page 2024

  30. [30]

    Dong and A

    Z. Dong and A. Ern. Hybrid high-order and weak Galerkin methods for the biharmonic problem.SIAM J. Numer. Anal.60.5 (2022), 2626–2656

    work page 2022

  31. [31]

    Z. Dong, L. Mascotto and O. J. Sutton. Residual-based a posteriori error estimates forhp- discontinuous Galerkin discretizations of the biharmonic problem.SIAM J. Numer. Anal. 59.3 (2021), 1273–1298

    work page 2021

  32. [32]

    D¨ orfler

    W. D¨ orfler. A convergent adaptive algorithm for Poisson’s equation.SIAM J. Numer. Anal. 33.3 (1996), 1106–1124

    work page 1996

  33. [33]

    Ern and P

    A. Ern and P. Zanotti. A quasi-optimal variant of the hybrid high-order method for elliptic partial differential equations withH −1 loads.IMA J. Numer. Anal.40.4 (2020), 2163– 2188

    work page 2020

  34. [34]

    Gallistl

    D. Gallistl. Morley finite element method for the eigenvalues of the biharmonic operator. IMA J. Numer. Anal.35.4 (2015), 1779–1811

    work page 2015

  35. [35]

    Gallistl and S

    D. Gallistl and S. Tian. A posteriori error estimates for nonconforming discretizations of singularly perturbed biharmonic operators.SMAI J. Comput. Math.10 (2024), 355–372

    work page 2024

  36. [36]

    E. H. Georgoulis and P. Houston. Discontinuous Galerkin methods for the biharmonic problem.IMA J. Numer. Anal.29.3 (2009), 573–594

    work page 2009

  37. [37]

    Guzm´ an, A

    J. Guzm´ an, A. Lischke and M. Neilan. Exact sequences on Worsey-Farin splits.Math. Comp.91.338 (2022), 2571–2608

    work page 2022

  38. [38]

    R. Khot, D. Mora and R. Ruiz-Baier. Virtual element methods for Biot-Kirchhoff poroelasti- city.Math. Comp.94.353 (2025), 1101–1146

    work page 2025

  39. [39]

    D. Li, C. Wang, J. Wang and S. Zhang. High order Morley elements for biharmonic equa- tions on polytopal partitions.J. Comput. Appl. Math.443 (2024), Paper No. 115757, 15

    work page 2024

  40. [40]

    Melzer and R

    H. Melzer and R. Rannacher. Spannungskonzentrationen in Eckpunkten der Kirchhoffschen Platte.Bauingenieur55 (1980), 181–184

    work page 1980

  41. [41]

    L. S. D. Morley. The Triangular Equilibrium Element in the Solution of Plate Bending Problems.Aeronautical Quarterly19.2 (1968), 149–169

    work page 1968

  42. [42]

    Mozolevski and E

    I. Mozolevski and E. S¨ uli. A priori error analysis for thehp-version of the discontinuous Galerkin finite element method for the biharmonic equation.Comput. Methods Appl. Math. 3.4 (2003), 596–607

    work page 2003

  43. [43]

    L. Mu, J. Wang and X. Ye. Weak Galerkin finite element methods for the biharmonic equation on polytopal meshes.Numer. Methods Partial Differential Equations30.3 (2014), 1003–1029. REFERENCES 25

    work page 2014

  44. [44]

    N. T. Tran. Discrete weak duality of hybrid high-order methods for convex minimization problems.SIAM J. Numer. Anal.62.4 (2024), 1492–1514

    work page 2024

  45. [45]

    N. T. Tran. Lower eigenvalue bounds with hybrid high-order methods.IMA Journal of Numerical Analysis(Mar. 2026), 1–21. eprint:https : / / academic . oup . com / imajna / advance-article-pdf/doi/10.1093/imanum/draf148/67341236/draf148.pdf

    work page doi:10.1093/imanum/draf148/67341236/draf148.pdf 2026

  46. [46]

    Veeser and P

    A. Veeser and P. Zanotti. Quasi-optimal nonconforming methods for symmetric elliptic problems. II—Overconsistency and classical nonconforming elements.SIAM J. Numer. Anal.57.1 (2019), 266–292

    work page 2019

  47. [47]

    ˇZen´ ıˇ sek

    A. ˇZen´ ıˇ sek. Polynomial approximation on tetrahedrons in the finite element method.J. Approximation Theory7 (1973), 334–351

    work page 1973

  48. [48]

    A. J. Worsey and G. Farin. Ann-dimensional Clough-Tocher interpolant.Constr. Approx. 3.2 (1987), 99–110

    work page 1987

  49. [49]

    Zhang and Q

    R. Zhang and Q. Zhai. A weak Galerkin finite element scheme for the biharmonic equations by using polynomials of reduced order.J. Sci. Comput.64.2 (2015), 559–585

    work page 2015

  50. [50]

    S. Zhang. A family of 3D continuously differentiable finite elements on tetrahedral grids. Appl. Numer. Math.59.1 (2009), 219–233

    work page 2009

  51. [51]

    J. Zhao, S. Chen and B. Zhang. The nonconforming virtual element method for plate bending problems.Math. Models Methods Appl. Sci.26.9 (2016), 1671–1687. Mathematical Institute, University of Oxford, Woodstock Rd, OX2 6GG Oxford, UK Email address:yizhou.liang@maths.ox.ac.uk Institute of Mathematics, University of Augsburg, Universit ¨atsstr. 12a, 86159 Au...

    work page 2016