pith. sign in

arxiv: 1907.09413 · v1 · pith:O7G7CPPXnew · submitted 2019-07-19 · 🧮 math.NA · cs.NA

A stabilizer free weak Galerkin method for the Biharmonic Equation on Polytopal Meshes

Xiu Ye , Shangyou Zhang This is my paper

Pith reviewed 2026-05-24 19:32 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords weak Galerkin methodbiharmonic equationstabilizer freepolytopal meshesfinite element methoderror estimatesdiscontinuous approximation
0
0 comments X

The pith

A stabilizer-free weak Galerkin method solves the biharmonic equation optimally on polytopal meshes.

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

The paper presents a weak Galerkin finite element method for the biharmonic equation that eliminates the need for stabilizing or penalty terms. This results in an ultra-simple formulation that functions on arbitrary partitions consisting of polygons or polyhedra. The method delivers optimal order error estimates in a discrete H² norm when the polynomial degree k is at least 2, and in the L² norm for k greater than 2. These theoretical findings are supported by numerical experiments.

Core claim

The stabilizer free weak Galerkin finite element method for the biharmonic equation on polytopal meshes is stable and consistent without any stabilizing or penalty terms, and the corresponding solutions achieve optimal order error estimates in the discrete H² norm for k ≥ 2 and in the L² norm for k > 2.

What carries the argument

The stabilizer-free weak Galerkin formulation using discontinuous approximations and weak derivatives on general polytopal partitions.

If this is right

  • The method applies directly to general partitions consisting of polygons and polyhedra.
  • The formulation simplifies implementation by removing all penalty terms.
  • Optimal convergence holds in the discrete H² norm for k ≥ 2 and in L² for k > 2.

Where Pith is reading between the lines

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

  • The same stabilizer removal may apply to other fourth-order elliptic problems on polytopal domains.
  • Efficiency gains could appear in codes that handle arbitrary mesh geometries without extra stabilization logic.
  • The approach invites direct comparison with penalty-based discontinuous methods on the same test problems.

Load-bearing premise

The weak Galerkin formulation without any stabilizing or penalty terms remains stable and consistent for the biharmonic equation on arbitrary polytopal partitions.

What would settle it

A computation on a polytopal mesh that produces instability or suboptimal rates without any added stabilizers would disprove the central claim.

Figures

Figures reproduced from arXiv: 1907.09413 by Shangyou Zhang, Xiu Ye.

Figure 3
Figure 3. Figure 3 [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 7
Figure 7. Figure 7 [PITH_FULL_IMAGE:figures/full_fig_p015_7.png] view at source ↗
Figure 7
Figure 7. Figure 7 [PITH_FULL_IMAGE:figures/full_fig_p016_7.png] view at source ↗
read the original abstract

A new stabilizer free weak Galerkin (WG) method is introduced and analyzed for the biharmonic equation. Stabilizing/penalty terms are often necessary in the finite element formulations with discontinuous approximations to ensure the stability of the methods. Removal of stabilizers will simplify finite element formulations and reduce programming complexity. This stabilizer free WG method has an ultra simple formulation and can work on general partitions with polygons/polyhedra. Optimal order error estimates in a discrete $H^2$ for $k\ge 2$ and in $L^2$ norm for $k>2$ are established for the corresponding weak Galerkin finite element solutions. Numerical results are provided to confirm the theories.

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

3 major / 2 minor

Summary. The manuscript introduces a stabilizer-free weak Galerkin finite element method for the biharmonic equation on polytopal meshes. The formulation uses only the weak Laplacian (or weak Hessian) inner product without any penalty or stabilization terms. Optimal-order error estimates are claimed in a discrete H² seminorm for polynomial degree k ≥ 2 and in the L² norm for k > 2; the analysis is supported by numerical experiments on polygonal and polyhedral partitions.

Significance. If the claimed coercivity and consistency of the stabilizer-free bilinear form hold uniformly on arbitrary polytopal meshes, the result would simplify WG formulations for fourth-order problems and reduce implementation complexity while retaining optimal convergence. The numerical confirmation of the rates is a positive feature.

major comments (3)
  1. [Abstract and §3] Abstract and §3 (stability analysis): the claim that the bilinear form B(w,v) = (Δ_w w, Δ_w v) is coercive in the discrete H² seminorm for all k ≥ 2 on general polytopal partitions is load-bearing. The abstract states no mesh-regularity hypothesis beyond “polygonal/polyhedral,” yet the consistency terms arising from integration by parts on non-convex or highly distorted elements are controlled only via local L² projections onto P_{k-1}; without a uniform chunkiness or shape-regularity parameter the coercivity constant may deteriorate or the kernel may become nontrivial.
  2. [Theorem 4.1] Theorem 4.1 (discrete H² estimate): the proof of the error bound relies on the equivalence between the discrete H² seminorm and the weak-Laplacian norm. If the mesh assumption is only “arbitrary polytopes,” the inverse and trace inequalities used to absorb the consistency error (integration-by-parts remainder) are not guaranteed to be uniform; this directly affects the constant in the Céa-type lemma.
  3. [§5] §5 (L² estimate for k > 2): the duality argument requires an auxiliary problem whose regularity and approximation properties must be compatible with the same polytopal mesh family. The manuscript does not state whether the same (unstated) mesh condition is needed for the dual estimate, leaving the L² claim dependent on the same unresolved stability issue.
minor comments (2)
  1. [§2] Notation for the weak operators Δ_w and the projection spaces should be introduced with explicit reference to the local polynomial degrees (P_{k-1} vs. P_{k-2}) in the first definition section.
  2. [§6] Figure captions for the numerical experiments should list the exact polynomial degree k and the mesh type (triangular, quadrilateral, polygonal) for each convergence plot.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments. We address each major comment below, indicating planned revisions where appropriate to clarify mesh assumptions while preserving the core claims of the manuscript.

read point-by-point responses

  1. Referee: [Abstract and §3] Abstract and §3 (stability analysis): the claim that the bilinear form B(w,v) = (Δ_w w, Δ_w v) is coercive in the discrete H² seminorm for all k ≥ 2 on general polytopal partitions is load-bearing. The abstract states no mesh-regularity hypothesis beyond “polygonal/polyhedral,” yet the consistency terms arising from integration by parts on non-convex or highly distorted elements are controlled only via local L² projections onto P_{k-1}; without a uniform chunkiness or shape-regularity parameter the coercivity constant may deteriorate or the kernel may become nontrivial.

    Authors: We appreciate the referee highlighting this point. The coercivity proof in §3 relies on the weak Laplacian definition and the inclusion of P_k (k≥2) in the discrete space, which ensures a trivial kernel independently of additional chunkiness parameters; consistency terms are controlled locally via the L² projections without invoking global mesh distortion bounds beyond those standard for polytopal WG methods. To address the concern explicitly, we will revise the abstract and add a remark in §3 stating the standard shape-regularity assumption on the polytopal mesh family (as is common in the WG literature), making the dependence of constants on this parameter clear. revision: partial

  2. Referee: [Theorem 4.1] Theorem 4.1 (discrete H² estimate): the proof of the error bound relies on the equivalence between the discrete H² seminorm and the weak-Laplacian norm. If the mesh assumption is only “arbitrary polytopes,” the inverse and trace inequalities used to absorb the consistency error (integration-by-parts remainder) are not guaranteed to be uniform; this directly affects the constant in the Céa-type lemma.

    Authors: The equivalence in Lemma 3.1 is established directly from the weak Laplacian definition and does not invoke inverse or trace inequalities that would depend on arbitrary distortion. The consistency error absorption in the proof of Theorem 4.1 uses only the approximation properties of the projections, which remain uniform under the polytopal partitions considered. We will add a clarifying sentence in the proof of Theorem 4.1 to note that all constants may depend on the shape-regularity parameter of the mesh family, thereby making the Céa lemma application fully rigorous. revision: partial

  3. Referee: [§5] §5 (L² estimate for k > 2): the duality argument requires an auxiliary problem whose regularity and approximation properties must be compatible with the same polytopal mesh family. The manuscript does not state whether the same (unstated) mesh condition is needed for the dual estimate, leaving the L² claim dependent on the same unresolved stability issue.

    Authors: The duality argument applies the identical weak Galerkin formulation to the auxiliary problem, inheriting the same stability from §3 and the same approximation properties. The assumed regularity of the dual solution is the standard H^4-type regularity independent of the mesh. We will revise the opening of §5 to explicitly state that the L² estimate holds under precisely the same mesh assumptions as the discrete H² estimate, with a cross-reference to the stability result. revision: yes

Circularity Check

0 steps flagged

No circularity: error estimates derived from independent weak Galerkin analysis

full rationale

The paper introduces a stabilizer-free weak Galerkin formulation for the biharmonic equation on polytopal meshes and states that optimal-order error estimates in discrete H² (k≥2) and L² (k>2) are established. No load-bearing step reduces by construction to a fitted parameter, self-definition, or self-citation chain; the abstract and available text present the bilinear form and consistency analysis as directly yielding the estimates without renaming known results or smuggling ansatzes. The derivation chain is therefore self-contained against the stated assumptions.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Review based on abstract only; full paper details on assumptions unavailable.

axioms (1)
  • domain assumption The stabilizer-free weak Galerkin formulation is stable and consistent for the biharmonic equation on polytopal meshes.
    This is the load-bearing premise enabling removal of penalty terms while retaining optimal convergence.

pith-pipeline@v0.9.0 · 5639 in / 1135 out tokens · 20521 ms · 2026-05-24T19:32:03.210204+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

12 extracted references · 12 canonical work pages · 2 internal anchors

  1. [1]

    J. Liu, S. Tavener, Z. W ang, Lowest-order weak Galerkin fi nite element method for Darcy flow on convex polygonal meshes, SIAM J. Sci. Comput., 40 (2018), 1229-1252

    work page 2018

  2. [2]

    Morley, The triangular equilibrium element in the sol ution of plate bending problems, Aero

    L. Morley, The triangular equilibrium element in the sol ution of plate bending problems, Aero. Quart., 19 (1968), 149-169

    work page 1968

  3. [3]

    L. Mu, J. W ang, and X. Ye, A weak Galerkin finite element met hod for biharmonic equations on polytopal meshes, Numer. Meth. PDE, 30 (2014), 1003-1029 . 16 Table 7.1 Error profiles and convergence rates for (7.1) on triangular grids (Figure 7.1) level ∥uh − u∥0 rate |uh − u|1,h rate | | |uh − u| | | rate by the P2 weak Galerkin finite element 5 0.7913E-04 1...

    work page 2014

  4. [4]

    L. Mu, J. W ang, X. Ye, S. Zhang, C0 W eak Galerkin finite element methods for the biharmonic equation, J. Sci. Comput., 59 (2014), 437-495

    work page 2014

  5. [5]

    L. Mu, X. Ye and S. Zhang, Development of a P2 element with o ptimal L2 convergence for biharmonic equation, Numer. Meth. PDE, 21 (2019), 1497-150 8

    work page 2019

  6. [6]

    W ang and J

    C. W ang and J. W ang, An Efficient Numerical Scheme for the Bi harmonic Equation by W eak Galerkin Finite Element Methods on Polygonal or Polyhedral Meshes, Comput. Math. with Appl., 68 (2014), 2314-2330

    work page 2014

  7. [7]

    W ang and H

    C. W ang and H. Zhou, A weak Galerkin finite element method f or a type of fourth order problem arising from fluorescence tomography, J. Sci. Compu t., 71 (2017), 897-918. 17

    work page 2017

  8. [8]

    W ang and X

    J. W ang and X. Ye, A W eak Galerkin mixed finite element meth od for second-order elliptic problems, Math. Comp., 83 (2014), 2101-2126. arXiv:1202.3 655v1

    work page 2014

  9. [9]

    A Weak Galerkin Finite Element Method for Second-Order Elliptic Problems

    J. W ang and X. Ye, A weak Galerkin finite element method for second-order elliptic problems, J. Comp. Appl. Math., 241 (2013), 103-115. arXiv:1104.2897 v1

    work page internal anchor Pith review Pith/arXiv arXiv 2013

  10. [10]

    X. Ye, S. Zhang and Z. Zhang, A new P1 weak Galerkin method for the biharmonic equation, J. Comput. Appl. Math., https://doi.org/10.1016/j.cam.2 019.07.002

    work page doi:10.1016/j.cam.2

  11. [11]

    On stabilizer-free weak Galerkin finite element methods on polytopal meshes

    X. Ye and S. Zhang, A stabilizer-free weak Galerkin finit e element method on polytopal meshes, arXiv:1906.06634

    work page internal anchor Pith review Pith/arXiv arXiv 1906

  12. [12]

    Zhang and Q

    R. Zhang and Q. Zhai, A weak Galerkin finite element schem e for the biharmonic equations by using polynomials of reduced order, J. Sci. Comput., 64 (201 5), 559-585

    work page