REVIEW 3 major objections 2 minor 42 references
Simplified Weak Galerkin Finite Element Methods for Biharmonic Equations on Non-Convex Polytopal Meshes
T0 review · 3 major / 2 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read This paper claims that a stabilizer-free weak Galerkin method for the biharmonic equation achieves optimal H^2 and L^2 convergence on both convex and non-convex polytopal meshes, using bubble functions instead of stabilization terms.
desk verdict Off-by-one degree error in Lemma 4.4 invalidates the central norm equivalence as written, but the fix is straightforward and the paper deserves a refereed revision. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The key machinery is the discrete weak second-order partial derivative \(\$partial^{2}$_{ij,w,r,T}\), defined on each element by the integration-by-parts identity \((\$partial^{2}$_{ij,w} v, \phi)_T = (v_0, \$partial^{2}$_{ji}\phi)_T - \langle v_b n_i, \partial_j \phi\rangle_{\partial T} + \langle v_{g i}, \phi n_j\rangle_{\partial T}\) for all \(\phi \in P_r(T)\). On a non-convex polytopal element with \(N\) edges, the paper sets \(r = 2N + k - 2\) and builds a full bubble function \(\Phi_B = \prod_{i=1}^N $l_i^{2}$\), which vanishes on \(\partial T\), together with edge bubbles \(\varphi_{e_k} = \prod_{i \ne k} $l_i^{2}$\). These bubbles convert boundary jump terms into interior norms and yield the two-sided bound \(C_1 \|v\|_{2,h} \le \left|\left|\left|v\right|\right|\right| \le C_2 \|v\|_{2,h}\). That norm equivalence is what carries both the uniqueness proof and the optimal \($H^{2}$\) and \($L^{2}$\) error estimates.
What would settle it
Check the degree of \((v_b - v_0) l_k \varphi_{e_k}\) on a pentagon with \(k=2\): if it exceeds \(r = 2N + k - 2\), identity (4.10) does not follow from (2.2). One could then verify (4.10) numerically on a single non-convex element, or run the scheme on an L-shaped domain and compare the observed discrete \($H^{2}$\) error rate with the claimed \($h^{{k-1}}$\).
Extended reading notes
Core claim
The central claim is that, for the biharmonic equation with Dirichlet and Neumann boundary conditions, the WG scheme defined by \((\$partial^{2}$_w u_h, \$partial^{2}$_w v) = (f, v_0)\) with no stabilization term has a unique solution and satisfies \(\left|\left|\left|u - u_h\right|\right|\right| \le C $h^{{k-1}}$ \|u\|_{k+1}\) and \(\|e_0\| \le C $h^{{k+1}}$ \|u\|_{k+1}\), provided the discrete weak second-order partial derivative is computed in a polynomial space of degree \(r = 2N + k - 2\) on an \(N\)-edge element. The proof obtains these rates by establishing a norm equivalence between the WG energy norm and a discrete \($H^{2}$\) norm, with bubble functions providing the needed lower bound on the boundary jump terms.
Load-bearing premise
In Lemma 4.4, the proof assumes that the edge bubble test function \((v_b - v_0) l_k \varphi_{e_k}\) is a polynomial of degree at most \(r = 2N + k - 2\) in each \(N\)-edge element, so that the defining identity for the discrete weak derivative applies; a direct degree count gives one degree higher, so the key boundary identity (4.10) is not justified as written.
Editorial extensions
If this is right
- Stabilizer-free WG schemes become available on general polytopal meshes, so non-convex cells need no special treatment at the discretization stage.
- The method is symmetric and positive definite, preserving the size and global sparsity of the stiffness matrix while removing the usual stabilization parameters.
- The claimed convergence rates \(O(h^{k-1})\) in the discrete \(H^2\) norm and \(O(h^{k+1})\) in \(L^2\) are optimal for polynomial degree \(k\), with the analysis stated to extend to \(d \ge 3\) dimensions.
- The bubble-function technique gives a template for proving norm equivalence in other stabilizer-free weak Galerkin schemes for fourth-order problems.
Reading between the lines
- If the degree-count issue in the bubble test function were repaired by raising \(r\) or choosing another bubble, the same error equation would still give the claimed rates, so the numerical convergence is likely robust even though the written proof has a gap.
- The edge bubble construction could be reused to build stabilizer-free WG methods for other fourth-order operators such as the plate equation with variable coefficients, where boundary terms require the same kind of control.
- A direct numerical test on a single non-convex element, measuring the discrete \(H^2\) error against \(h^{k-1}\) for \(k=2,3\), would settle whether the claimed rates hold in practice independent of the proof details.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes a stabilizer-free weak Galerkin (WG) finite element method for the biharmonic equation on polytopal meshes without a convexity restriction. The discrete weak second-order derivative is computed with polynomial degree r = 2N + k - 2 on an N-edge element, and bubble-function arguments are used to prove a norm equivalence between the WG energy and a discrete H^2 norm. On this basis the paper claims existence, uniqueness, and optimal-order error estimates in the discrete H^2 norm and in L^2, thereby extending the stabilizer-free WG framework of Xiu and Zhang [42] to non-convex polytopal meshes.
Significance. If the results are correct, the paper would provide a useful extension of stabilizer-free WG methods: it removes the convexity restriction, keeps the scheme symmetric and positive definite, and avoids stabilizers. The constructive bubble-function strategy is transparent, and the claimed optimal rates are concrete and falsifiable. The main caveat is that the key norm-equivalence lemma contains a genuine degree mismatch as written, so the central estimates are not yet supported; however, the defect appears to be a locally repairable off-by-one error rather than a demonstrated counterexample. The paper would also be stronger if its proof were self-contained with respect to the extension lemmas currently delegated to the author's arXiv preprint [23].
major comments (3)
- [Section 4, Lemma 4.4 (Eq. (4.10))] The admissible test-function space for identity (2.2) is P_r(T) with r = 2N + k - 2, but the function phi = (v_b - v_0) l_k phi_{e_k} has degree at most 2N - 1 + k = 2N + k - 1, because l_k is linear, phi_{e_k} has degree 2N - 2, and v_b - v_0 is generically of degree k. Thus (2.2) cannot be applied to this phi, and the displayed identity (4.10) is unjustified. Consequently the bounds (4.11) and the lower half of the norm equivalence (4.9) are not established; since (4.9) underpins the uniqueness proof in Theorem 4.5, the H^2 error estimate in Theorem 6.3, and the L^2 estimate in Theorem 7.1, the main results are unsupported as written. The defect appears repairable: setting r = 2N + k - 1 makes phi admissible, keeps Lemma 4.1 valid (the test function Phi_B partial^2_{ij} v_0 has degree 2N + k - 2), and only improves the Q_r projection estimates.
- [Section 4, Lemma 4.4 and Lemma 4.2] The proof relies on extending v_b and the trace of v_0 from e_k to the element T and on Lemma 4.2, but Lemma 4.2 and the extension construction are attributed to the author's arXiv preprint [23], while Lemma 4.3 in this paper only describes the extension of v_g. Since the inequality (4.7) and the extension properties are load-bearing for (4.11), the paper is not self-contained. The revision should include these proofs or cite a peer-reviewed version containing them.
- [Section 6, Lemma 6.2] In bounding the terms ||(Q_0 u - Q_b u) n_i||_{partial T} and ||partial_i(Q_0 u) - Q_n(partial_i u)||_{partial T}, the proof passes directly to estimates involving u - Q_0 u on T. This omits the edge projection errors u - Q_b u and partial_i u - Q_n(partial_i u); a trace-plus-approximation argument for Q_b and Q_n is needed to obtain the factor h^{k-1} in (6.3). The missing step is standard and fixable, but as written the proof of (6.3), which is used in Theorems 6.3 and 7.1, is incomplete.
minor comments (2)
- [Section 4, Lemma 4.4] The same symbol v_0 is used for the original interior polynomial and for its edge-trace extension, which makes the substitution into (2.2) ambiguous; please use a separate notation such as tilde{v}_0 for the extension.
- [Equation (4.10)] The boundary contribution is displayed as +C h_T^{-1} int |v_b - v_0|^2 phi_{e_k} ds, whereas a direct calculation gives the signed term -h_T^{-1} n_i n_j int |v_b - v_0|^2 phi_{e_k} ds; the proof should take absolute values or specify the sum over i,j.
Circularity Check
Load-bearing norm-equivalence lemmas are deferred to the author's own preprint [23], but the central error estimates are otherwise derived, not assumed.
-
self citation load bearing
[Section 4, Lemmas 4.2 and 4.3; used in Lemma 4.4, Theorem 4.5, Theorem 6.3, and Theorem 7.1]
"Lemma 4.2. [23] For {v0, vb, vg} ∈ Vh, let ϕ = (vb − v0)lkϕek. The following inequality holds: (4.7) ‖ϕ‖2 T ≤ ChT ∫ ek (vb − v0)2ds. ... The extension preserves the polynomial nature of vg as demonstrated in [23]."
The norm equivalence (4.9) is the backbone of the paper: Theorem 4.5 uses it for uniqueness, Theorem 6.3 uses it to convert ||v||_{2,h} into |||v||| in (6.5)-(6.6), and Theorem 7.1 uses it in (7.9). To prove the lower bound in (4.9), Lemma 4.4 needs Lemmas 4.2 and 4.3, which bound volume integrals of (v_b - v_0) and (v_gi - ∂_i v_0) by edge integrals. These lemmas carry no proof in the present paper; they are attributed to [23], an arXiv preprint by the same author. The cited preprint is therefore load-bearing and is not independent, machine-checked, or reproduced evidence. This is a self-citation gap, not a reduction of the theorem to its own definition; the subsequent a priori estimates are still derived rather than assumed.
full rationale
The paper contains no fitted parameters, no data-driven predictions, and no renaming of a known empirical pattern; the WG scheme and its error analysis are conventional a priori derivations. The only circularity-adjacent issue is the reliance of Section 4 on the author's own earlier preprint [23]. Lemmas 4.2 and 4.3 are cited, not proved, and they are used to establish the norm equivalence (4.9) that later supports uniqueness and the H2 and L2 error estimates. Because [23] is a single-author preprint by the same author, this step does not have the independent evidentiary status of a machine-checked or externally reproduced proof. Nevertheless, the main estimates do not reduce to the theorem statement by construction: the discrete weak derivative, the bubble-function arguments, the error equation, and the duality argument are all carried out in this manuscript. A separate off-by-one concern appears in Lemma 4.4, where the test function φ = (v_b - v_0)l_k φ_{e_k} has degree 2N + k - 1 while r = 2N + k - 2; that would invalidate the displayed equality (4.10) as written. This is a correctness/consistency defect, not a circularity, because the proof does not assume the target error estimate. Overall, the derivation chain has a load-bearing self-citation but no definitional circularity, so score 4.
Assumptions & free parameters
free parameters (1)
- r (discrete weak derivative degree) =
2N + k - 2
assumptions (6)
- domain assumption Shape regularity of the polytopal partition T_h
- ad hoc to paper Lemma 4.2 from [23]
- ad hoc to paper Lemma 4.3 from [23] (edge-to-element polynomial extension)
- domain assumption H4-regularity of the dual problem (7.2)
- standard math Domain inverse inequality with bubble weight (4.6)
- domain assumption Regularity of the exact solution u in H^{k+1}
Cite this review
Pith. "Pith review of Simplified Weak Galerkin Finite Element Methods for Biharmonic Equations on Non-Convex Polytopal Meshes." pith.science (2026). https://pith.science/paper/YU3QHLD6
@misc{pith2026241211315,
author = {Pith},
title = {Pith review of: Simplified Weak Galerkin Finite Element Methods for Biharmonic Equations on Non-Convex Polytopal Meshes},
year = {2026},
howpublished = {\url{https://pith.science/paper/YU3QHLD6}},
note = {Machine review of arXiv:2412.11315}
}
abstract
This paper presents a simplified weak Galerkin (WG) finite element method for solving biharmonic equations avoiding the use of traditional stabilizers. The proposed WG method supports both convex and non-convex polytopal elements in finite element partitions, utilizing bubble functions as a critical analytical tool. The simplified WG method is symmetric and positive definite. Optimal-order error estimates are established for WG approximations in both the discrete $H^2$ norm and the $L^2$ norm.
Reference graph
Works this paper leans on
-
[42]
Y. Xiu and S. Zhang , A stabilizer free weak Galerkin finite element method for the bihar- monic equation on polytopal meshes , SIAM J. Numer. Anal., vol 58, No. 5, pp. 2572-2588, 2020. Department of Mathematics, University of Florida, Gainesvi lle, FL 32611, USA. Email address : chunmei.wang@ufl.edu
work page 2020
-
[23]
C. W ang, Auto-Stabilized Weak Galerkin Finite Element Methods for B iharmonic Equations on Polytopal Meshes without Convexity Assumptions , arXiv:2409.05887
-
[1]
W. Cao, C. W ang and J. W ang , An Lp-Primal-Dual Weak Galerkin Method for div-curl Systems, Journal of Computational and Applied Mathematics, vol. 42 2, 114881, 2023
work page 2023
-
[2]
W. Cao, C. W ang and J. W ang, An Lp-Primal-Dual Weak Galerkin Method for Convection- Diffusion Equations , Journal of Computational and Applied Mathematics, vol. 41 9, 114698, 2023
work page 2023
-
[3]
W. Cao, C. W ang and J. W ang , A New Primal-Dual Weak Galerkin Method for Elliptic Interface Problems with Low Regularity Assumptions , Journal of Computational Physics, vol. 470, 111538, 2022
work page 2022
-
[4]
S. Cao, C. W ang and J. W ang , A new numerical method for div-curl Systems with Low Regularity Assumptions, Computers and Mathematics with Applications, vol. 144, pp . 47-59, 2022
work page 2022
- [5]
-
[6]
D. Li, Y. Nie, and C. W ang , Superconvergence of Numerical Gradient for Weak Galerkin Finite Element Methods on Nonuniform Cartesian Partitions in Three Dimensions , Com- puters and Mathematics with Applications, vol 78(3), pp. 90 5-928, 2019
work page 2019
Show all 42 references
-
[7]
D. Li, C. W ang and J. W ang , An Extension of the Morley Element on General Polytopal Partitions Using Weak Galerkin Methods , Journal of Scientific Computing, 100, vol 27, 2024
2024
-
[8]
D. Li, C. W ang and S. Zhang , Weak Galerkin methods for elliptic interface problems on curved polygonal partitions, Journal of Computational and Applied Mathematics, pp. 115 995, 2024
2024
-
[9]
D. Li, C. W ang, J. W ang and X. Ye , Generalized weak Galerkin finite element methods for second order elliptic problems , Journal of Computational and Applied Mathematics, vol. 445, pp. 115833, 2024
2024
-
[10]
D. Li, C. W ang, J. W ang and S. Zhang , High Order Morley Elements for Biharmonic Equations on Polytopal Partitions , Journal of Computational and Applied Mathematics, Vol. 443, pp. 115757, 2024
2024
-
[11]
D. Li, C. W ang and J. W ang , Curved Elements in Weak Galerkin Finite Element Methods , Computers and Mathematics with Applications, Vol. 153, pp. 20-32, 2024
2024
-
[12]
D. Li, C. W ang and J. W ang , Generalized Weak Galerkin Finite Element Methods for Bi- harmonic Equations , Journal of Computational and Applied Mathematics, vol. 43 4, 115353, 2023
2023
-
[13]
D. Li, C. W ang and J. W ang , An Lp-primal-dual finite element method for first-order transport problems, Journal of Computational and Applied Mathematics, vol. 43 4, 115345, 2023
2023
-
[14]
Li and C
D. Li and C. W ang , A simplified primal-dual weak Galerkin finite element method for Fokker-Planck type equations , Journal of Numerical Methods for Partial Differential Equa - tions, vol 39, pp. 3942-3963, 2023
2023
-
[15]
D. Li, C. W ang and J. W ang , Primal-Dual Weak Galerkin Finite Element Methods for Transport Equations in Non-Divergence Form , Journal of Computational and Applied Math- ematics, vol. 412, 114313, 2022
2022
-
[16]
D. Li, C. W ang, and J. W ang , Superconvergence of the Gradient Approximation for Weak Galerkin Finite Element Methods on Rectangular Partitions , Applied Numerical Mathemat- ics, vol. 150, pp. 396-417, 2020
2020
-
[17]
B. Li, X. Xie and S. Zhang , BPS preconditioners for a weak Galerkin finite element metho d for 2D diffusion problems with strongly discontinuous coeffic ients, Computers & Mathematics with Applications, 76(4), pp.701-724, 2018. 18 CHUNMEI W ANG
2018
-
[18]
W ang, A Preconditioner for the FETI-DP Method for Mortar-Type Cro uzeix-Raviart Element Discretization , Applications of Mathematics, Vol
C. W ang, A Preconditioner for the FETI-DP Method for Mortar-Type Cro uzeix-Raviart Element Discretization , Applications of Mathematics, Vol. 59, 6, pp. 653-672, 2014
2014
-
[19]
W ang, New Discretization Schemes for Time-Harmonic Maxwell Equa tions by Weak Galerkin Finite Element Methods , Journal of Computational and Applied Mathematics, Vol
C. W ang, New Discretization Schemes for Time-Harmonic Maxwell Equa tions by Weak Galerkin Finite Element Methods , Journal of Computational and Applied Mathematics, Vol. 341, pp. 127-143, 2018
2018
-
[20]
W ang, Low Regularity Primal-Dual Weak Galerkin Finite Element Me thods for Ill-Posed Elliptic Cauchy Problems , Int
C. W ang, Low Regularity Primal-Dual Weak Galerkin Finite Element Me thods for Ill-Posed Elliptic Cauchy Problems , Int. J. Numer. Anal. Mod., vol. 19(1), pp. 33-51, 2022
2022
-
[21]
W ang, A Modified Primal-Dual Weak Galerkin Finite Element Method f or Second Order Elliptic Equations in Non-Divergence Form , Int
C. W ang, A Modified Primal-Dual Weak Galerkin Finite Element Method f or Second Order Elliptic Equations in Non-Divergence Form , Int. J. Numer. Anal. Mod., vol. 18(4), pp. 500- 523, 2021
2021
-
[22]
W ang, A New Primal-Dual Weak Galerkin Finite Element Method for Il l-posed Elliptic Cauchy Problems , Journal of Computational and Applied Mathematics, vol 371 , 112629, 2020
C. W ang, A New Primal-Dual Weak Galerkin Finite Element Method for Il l-posed Elliptic Cauchy Problems , Journal of Computational and Applied Mathematics, vol 371 , 112629, 2020
2020
-
[24]
W ang and J
C. W ang and J. W ang, A Primal-Dual Weak Galerkin Finite Element Method for Fokke r- Planck Type Equations , SIAM Numerical Analysis, vol. 58(5), pp. 2632-2661, 2020
2020
-
[25]
W ang and J
C. W ang and J. W ang , A Primal-Dual Finite Element Method for First-Order Transp ort Problems, Journal of Computational Physics, Vol. 417, 109571, 2020
2020
-
[26]
W ang and J
C. W ang and J. W ang , Primal-Dual Weak Galerkin Finite Element Methods for Ellip tic Cauchy Problems , Computers and Mathematics with Applications, vol 79(3), p p. 746-763, 2020
2020
-
[27]
W ang and J
C. W ang and J. W ang , A Primal-Dual Weak Galerkin Finite Element Method for Secon d Order Elliptic Equations in Non-Divergence form , Mathematics of Computation, Vol. 87, pp. 515-545, 2018
2018
-
[28]
W ang and J
C. W ang and J. W ang, Discretization of Div-Curl Systems by Weak Galerkin Finite Element Methods on Polyhedral Partitions , Journal of Scientific Computing, Vol. 68, pp. 1144-1171, 2016
2016
-
[29]
W ang and J
C. W ang and J. W ang, A Hybridized Formulation for Weak Galerkin Finite Element M eth- ods for Biharmonic Equation on Polygonal or Polyhedral Mesh es, International Journal of Numerical Analysis and Modeling, Vol. 12, pp. 302-317, 2015
2015
-
[30]
W ang and C
J. W ang and C. W ang, Weak Galerkin Finite Element Methods for Elliptic PDEs , Science China, Vol. 45, pp. 1061-1092, 2015
2015
-
[31]
W ang and J
C. W ang and J. W ang , An Efficient Numerical Scheme for the Biharmonic Equation by Weak Galerkin Finite Element Methods on Polygonal or Polyhe dral Meshes, Journal of Com- puters and Mathematics with Applications, Vol. 68, 12, pp. 2 314-2330, 2014
2014
-
[32]
W ang, J
C. W ang, J. W ang, R. W ang and R. Zhang, A Locking-Free Weak Galerkin Finite Element Method for Elasticity Problems in the Primal Formulation , Journal of Computational and Applied Mathematics, Vol. 307, pp. 346-366, 2016
2016
-
[33]
W ang, J
C. W ang, J. W ang, X. Ye and S. Zhang , De Rham Complexes for Weak Galerkin Finite Element Spaces, Journal of Computational and Applied Mathematics, vol. 39 7, pp. 113645, 2021
2021
-
[34]
W ang, J
C. W ang, J. W ang and S. Zhang , Weak Galerkin Finite Element Methods for Optimal Control Problems Governed by Second Order Elliptic Partial Differential Equations , Journal of Computational and Applied Mathematics, in press, 2024
2024
-
[35]
W ang, J
C. W ang, J. W ang and S. Zhang , A parallel iterative procedure for weak Galerkin methods for second order elliptic problems , International Journal of Numerical Analysis and Modeling , vol. 21(1), pp. 1-19, 2023
2023
-
[36]
W ang, J
C. W ang, J. W ang and S. Zhang , Weak Galerkin Finite Element Methods for Quad-Curl Problems, Journal of Computational and Applied Mathematics, vol. 42 8, pp. 115186, 2023
2023
-
[37]
W ang, and X
J. W ang, and X. Ye , A weak Galerkin mixed finite element method for second-order elliptic problems, Math. Comp., vol. 83, pp. 2101-2126, 2014
2014
-
[38]
W ang, X
C. W ang, X. Ye and S. Zhang , A Modified weak Galerkin finite element method for the Maxwell equations on polyhedral meshes , Journal of Computational and Applied Mathemat- ics, vol. 448, pp. 115918, 2024
2024
-
[39]
W ang and S
C. W ang and S. Zhang , A Weak Galerkin Method for Elasticity Interface Problems , Journal of Computational and Applied Mathematics, vol. 419, 114726 , 2023. 19
2023
-
[40]
W ang and L
C. W ang and L. Zikatanov , Low Regularity Primal-Dual Weak Galerkin Finite Element Methods for Convection-Diffusion Equations , Journal of Computational and Applied Math- ematics, vol 394, 113543, 2021
2021
-
[41]
W ang and H
C. W ang and H. Zhou, A Weak Galerkin Finite Element Method for a Type of Fourth Ord er Problem arising from Fluorescence Tomography, Journal of Scientific Computing, Vol. 71(3), pp. 897-918, 2017
2017
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.