pith. sign in

arxiv: 2605.21764 · v1 · pith:PSGJJOQKnew · submitted 2026-05-20 · 🧮 math.NA · cs.NA

Quasi-optimal polytopal finite element methods for biharmonic equation

Ngoc Tien Tran This is my paper

Pith reviewed 2026-05-22 08:19 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords biharmonic equationweak Galerkindiscontinuous Galerkinhybrid high-orderpolytopal mesheserror estimatesa posteriori estimators
0
0 comments X

The pith

Weak Galerkin, discontinuous Galerkin and hybrid-high order methods deliver quasi-optimal error estimates for the biharmonic equation on general polytopal meshes under minimal regularity.

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

The paper proves quasi-optimal and lower-order error estimates for weak Galerkin, discontinuous Galerkin, and hybrid-high order finite element methods applied to the biharmonic equation. These estimates hold under only minimal regularity assumptions on the solution and on completely general polytopal meshes without extra smoothness or mesh constraints. A reader would care because the results support reliable numerical approximation of fourth-order problems on complex domains where solutions lack high smoothness. The work also establishes that the stabilization terms contribute efficiently to a posteriori error estimators.

Core claim

This paper establishes quasi-optimal and lower-order error estimates for weak Galerkin, discontinuous Galerkin, and hybrid-high order finite element methods for the biharmonic equation under minimal regularity assumptions on general polytopal meshes. Furthermore, it is shown that the stabilization is an efficient contribution in a posteriori error estimators.

What carries the argument

Quasi-optimal error estimates under minimal regularity on polytopal meshes, which bound the discretization error in terms of best approximation without requiring additional solution smoothness or mesh regularity.

If this is right

  • The methods can be applied directly to biharmonic problems with low-regularity data on complex polytopal domains.
  • Stabilization terms can be incorporated reliably into a posteriori error estimators for adaptive mesh refinement.
  • Lower-order estimates remain available as practical bounds when full quasi-optimality cannot be attained due to limited smoothness.

Where Pith is reading between the lines

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

  • Similar quasi-optimal analysis may carry over to other fourth-order elliptic equations such as the Kirchhoff plate model.
  • The approach could support development of adaptive algorithms that exploit the efficient stabilization contribution for automatic mesh adaptation.
  • Extensions to three-dimensional or time-dependent biharmonic-type problems on polytopal meshes become more feasible under the same minimal assumptions.

Load-bearing premise

The stabilization terms in the methods must remain efficient contributors to the error bounds even on completely general polytopal meshes without added regularity requirements.

What would settle it

Numerical computation of observed convergence rates for a biharmonic test problem with a solution of minimal regularity solved on a sequence of highly irregular polytopal meshes, to check whether the rates match the quasi-optimal theoretical predictions.

read the original abstract

This paper establishes quasi-optimal and lower-order error estimates for weak Galerkin, discontinuous Galerkin, and hybrid-high order finite element methods for the biharmonic equation under minimal regularity assumptions on general polytopal meshes. Furthermore, it is shown that the stabilization is an efficient contribution in a~posteriori error estimators.

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 manuscript develops quasi-optimal a priori error estimates together with lower-order estimates for weak Galerkin, discontinuous Galerkin, and hybrid-high-order discretizations of the biharmonic equation. The analysis is performed under minimal regularity assumptions on the solution and is claimed to hold on completely general polytopal meshes; the paper further asserts that the stabilization terms are efficient contributors to a posteriori error estimators.

Significance. If the central claims hold, the work would provide a unified treatment of three popular polytopal methods for a fourth-order elliptic problem, extending their applicability to domains that benefit from general polygonal/polyhedral meshes while retaining quasi-optimality and a posteriori efficiency. Such results are of interest for both theoretical analysis and practical computation on complex geometries.

major comments (2)
  1. [Abstract and §1] Abstract and §1: the headline claim that the quasi-optimal estimates hold 'on general polytopal meshes' without additional mesh regularity is not yet secured. Standard stability and approximation arguments for WG, DG and HHO schemes on polytopes rely on each element being star-shaped with respect to a ball whose radius is comparable to the diameter (the chunkiness parameter) together with a uniform bound on the number of faces. If §§3–5 invoke these quantities without declaring them part of the 'minimal' hypothesis, the Céa-type quasi-optimality result is conditional on hidden mesh assumptions that contradict the stated generality.
  2. [§4] §4 (or the corresponding a-priori analysis section): the proof of the consistency and stability of the discrete bilinear form for the biharmonic operator must be checked against the mesh hypotheses. If the estimates reduce to a shape-regularity or star-shapedness condition that is not listed among the 'minimal regularity assumptions,' the central quasi-optimality statement requires revision or an explicit additional hypothesis.
minor comments (2)
  1. [§2] Notation for the stabilization terms and the broken Sobolev spaces should be introduced once and used consistently across the three methods to improve readability.
  2. [§5] The lower-order estimates are stated without an explicit comparison to the quasi-optimal rates; a short remark clarifying when the lower-order result is the best that can be expected under the given regularity would be helpful.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and for highlighting the need to clarify the mesh regularity assumptions underlying our claims. We address the major comments point by point below and will revise the manuscript to make the hypotheses explicit.

read point-by-point responses

  1. Referee: [Abstract and §1] Abstract and §1: the headline claim that the quasi-optimal estimates hold 'on general polytopal meshes' without additional mesh regularity is not yet secured. Standard stability and approximation arguments for WG, DG and HHO schemes on polytopes rely on each element being star-shaped with respect to a ball whose radius is comparable to the diameter (the chunkiness parameter) together with a uniform bound on the number of faces. If §§3–5 invoke these quantities without declaring them part of the 'minimal' hypothesis, the Céa-type quasi-optimality result is conditional on hidden mesh assumptions that contradict the stated generality.

    Authors: We agree that the quasi-optimal error estimates rely on the standard assumptions for polytopal meshes, namely that each element is star-shaped with respect to a ball whose radius is comparable to the diameter (bounded chunkiness parameter) and that the number of faces per element is uniformly bounded. These conditions are the minimal ones required for the approximation properties of the discrete spaces and for the stability arguments in WG, DG, and HHO methods; they are routinely adopted in the polytopal finite-element literature. In the manuscript the phrase 'general polytopal meshes' with 'minimal regularity assumptions' primarily concerns the solution regularity, while the mesh hypotheses are the conventional ones implicit in the analysis. To eliminate any ambiguity we will revise the abstract, Section 1, and the preliminaries section to state these mesh conditions explicitly as part of the overall hypotheses. revision: yes

  2. Referee: [§4] §4 (or the corresponding a-priori analysis section): the proof of the consistency and stability of the discrete bilinear form for the biharmonic operator must be checked against the mesh hypotheses. If the estimates reduce to a shape-regularity or star-shapedness condition that is not listed among the 'minimal regularity assumptions,' the central quasi-optimality statement requires revision or an explicit additional hypothesis.

    Authors: The consistency and stability proofs in the a-priori analysis section do invoke the star-shapedness and bounded-face assumptions when controlling the stabilization terms and establishing the discrete inf-sup condition for the biharmonic operator. These steps follow the standard arguments for polytopal methods and do not require stronger regularity than the chunkiness and face-count bounds already mentioned. We will add explicit cross-references to the mesh hypotheses in this section and ensure they appear in the statement of assumptions at the beginning of the paper. With these clarifications the quasi-optimality result will be correctly qualified; no alteration of the estimates themselves is required. revision: yes

Circularity Check

0 steps flagged

No circularity detected; derivation relies on standard a priori/a posteriori analysis

full rationale

The paper claims quasi-optimal error estimates for WG/DG/HHO methods on polytopal meshes under minimal regularity, with stabilization efficiency in a posteriori estimators. No quoted steps reduce by construction to fitted inputs, self-definitions, or load-bearing self-citations that presuppose the target result. The abstract and described analysis chain use conventional consistency/stability arguments and Céa-type lemmas without renaming known patterns or smuggling ansatzes via prior self-work as the sole justification. The derivation remains self-contained against external benchmarks such as standard finite-element theory.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Review based solely on abstract; full technical assumptions and any hidden parameters are not visible.

axioms (1)
  • domain assumption Minimal regularity assumptions suffice for the error estimates on general polytopal meshes
    The central claims rest on this premise for both a priori rates and a posteriori efficiency.

pith-pipeline@v0.9.0 · 5558 in / 1256 out tokens · 37409 ms · 2026-05-22T08:19:35.806128+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

29 extracted references · 29 canonical work pages · 1 internal anchor

  1. [1]

    D. N. Arnold, F. Brezzi, B. Cockburn and L. D. Marini. Unified analysis of discontinuous Galerkin methods for elliptic problems. In:SIAM J. Numer. Anal.39.5 (2001/02), 1749– 1779

    work page 2001

  2. [2]

    Berger, R

    A. Berger, R. Scott and G. Strang. Approximate boundary conditions in the finite element method. In:Symposia Mathematica, Vol. X (Convegno di Geometria Differenziale, IN- DAM, Rome, 1971 & Convegno di Analisi Numerica, INDAM, Rome, 1972). Academic Press, London-New York, 1972, 295–313

    work page 1971

  3. [3]

    Bertrand, C

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

    work page 2023

  4. [4]

    Blum and R

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

    work page 1980

  5. [5]

    S. C. Brenner, T. Gudi and L.-y. Sung. An a posteriori error estimator for a quadratic C0-interior penalty method for the biharmonic problem. In:IMA J. Numer. Anal.30.3 (2010), 777–798

    work page 2010

  6. [6]

    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. In:SIAM J. Numer. Anal. 61.5 (2023), 2460–2484

    work page 2023

  7. [7]

    Carstensen and N

    C. Carstensen and N. Nataraj. A priori and a posteriori error analysis of the Crouzeix- Raviart and Morley FEM with original and modified right-hand sides. In:Comput. Methods Appl. Math.21.2 (2021), 289–315

    work page 2021

  8. [8]

    Carstensen and N

    C. Carstensen and N. T. Tran. Locking-free hybrid high-order method for linear elasticity. In:SIAM J. Numer. Anal.63.2 (2025), 827–853

    work page 2025

  9. [9]

    J. C´ ea. Approximation variationnelle des probl` emes aux limites. In:Ann. Inst. Fourier (Grenoble)14 (1964), 345–444

    work page 1964

  10. [10]

    Chaumont-Frelet

    T. Chaumont-Frelet. A new family of a posteriori error estimates for non-conforming finite element methods leading to stabilization-free error bounds. In:arXiv:2506.23381(2026)

    work page arXiv 2026

  11. [11]

    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

  12. [12]

    D. A. Di Pietro and A. Ern. Mathematical aspects of discontinuous Galerkin methods. Vol. 69. Math´ ematiques & Applications (Berlin) [Mathematics & Applications]. Springer, Heidelberg, 2012, xviii+384

    work page 2012

  13. [13]

    Dong and A

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

    work page 2022

  14. [14]

    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. In:IMA J. Numer. Anal.40.4 (2020), 2163– 2188

    work page 2020

  15. [15]

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

    work page 2009

  16. [16]

    E. H. Georgoulis, P. Houston and J. Virtanen. Ana posteriorierror indicator for discontinu- ous Galerkin approximations of fourth-order elliptic problems. In:IMA J. Numer. Anal. 31.1 (2011), 281–298

    work page 2011

  17. [17]

    Grisvard

    P. Grisvard. Elliptic problems in nonsmooth domains. Vol. 69. Classics in Applied Mathem- atics. Reprint of the 1985 original [MR0775683], With a foreword by Susanne C. Brenner. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 2011, xx+410

    work page 1985

  18. [18]

    T. Gudi. A new error analysis for discontinuous finite element methods for linear elliptic problems. In:Math. Comp.79.272 (2010), 2169–2189

    work page 2010

  19. [19]

    Guzm´ an, A

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

    work page 2022

  20. [20]

    O. A. Karakashian and F. Pascal. Convergence of adaptive discontinuous Galerkin approx- imations of second-order elliptic problems. In:SIAM J. Numer. Anal.45.2 (2007), 641– 665

    work page 2007

  21. [21]

    A hybrid high-order method for the biharmonic problem

    Y. Liang and N. T. Tran. A hybrid high-order method for the biharmonic problem. In: arXiv:2504.16608(2026)

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  22. [22]

    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. In:Comput. Methods Appl. Math.3.4 (2003), 596–607

    work page 2003

  23. [23]

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

    work page 2014

  24. [24]

    Veeser and P

    A. Veeser and P. Zanotti. Quasi-optimal nonconforming methods for symmetric elliptic problems. I—Abstract theory. In:SIAM J. Numer. Anal.56.3 (2018), 1621–1642

    work page 2018

  25. [25]

    Veeser and P

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

    work page 2019

  26. [26]

    Veeser and P

    A. Veeser and P. Zanotti. Quasi-optimal nonconforming methods for symmetric elliptic problems. III—Discontinuous Galerkin and other interior penalty methods. In:SIAM J. Numer. Anal.56.5 (2018), 2871–2894

    work page 2018

  27. [27]

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

    work page 1987

  28. [28]

    Zhang and Q

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

    work page 2015

  29. [29]

    S. Zhang. A family of 3D continuously differentiable finite elements on tetrahedral grids. In:Appl. Numer. Math.59.1 (2009), 219–233. REFERENCES 17 (N. T. Tran)Universit ¨at Augsburg, 86159 Augsburg, Germany Email address:ngoc1.tran@uni-a.de

    work page 2009