pith. sign in

arxiv: 1907.01397 · v2 · pith:REDKJDIPnew · submitted 2019-07-01 · 🧮 math.NA · cs.NA

A conforming discontinuous Galerkin finite element method: Part II

Xiu Ye , Shangyou Zhang This is my paper

Pith reviewed 2026-05-25 12:24 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords conforming discontinuous Galerkinpolytopal meshesweak gradienterror estimatesfinite element methoddiscontinuous approximation
0
0 comments X

The pith

Conforming discontinuous Galerkin methods achieve optimal error estimates on general polytopal meshes through appropriate weak gradient design.

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

The paper extends a conforming discontinuous Galerkin finite element method from simplicial meshes to general polytopal meshes. The extension works by designing the weak gradient so the method stays conforming even though the approximation functions are discontinuous. Two formulations are given that treat boundary conditions differently. Optimal-order error estimates are proved for the resulting approximations in a discrete H1 norm and in the L2 norm, with numerical tests confirming the rates.

Core claim

By designing the weak gradient appropriately, the conforming DG method extends to general polytopal meshes while remaining conforming, and the approximations satisfy optimal-order error estimates in both the discrete H1 norm and the L2 norm.

What carries the argument

The weak gradient ∇_w, designed to enforce conformity on polytopal meshes while permitting discontinuous trial functions.

If this is right

  • The method applies directly to meshes made of general polytopes without requiring a simplicial triangulation.
  • Optimal convergence holds in the discrete H1 norm for the two boundary-condition variants.
  • Optimal convergence also holds in the L2 norm.
  • The formulations remain simple in structure despite using discontinuous spaces.

Where Pith is reading between the lines

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

  • The approach may allow more flexible mesh generation for domains with complex geometry.
  • Similar weak-gradient constructions could be tested on other classes of discontinuous methods or higher-order elements.
  • The two boundary-condition treatments may differ in implementation cost for large-scale problems.

Load-bearing premise

The weak gradient can be designed so the method remains conforming on arbitrary polytopal meshes.

What would settle it

A computation on a polytopal mesh in which the chosen weak gradient produces suboptimal convergence rates in the discrete H1 or L2 norm would falsify the error estimates.

Figures

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

Figure 3
Figure 3. Figure 3 [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 6
Figure 6. Figure 6 [PITH_FULL_IMAGE:figures/full_fig_p013_6.png] view at source ↗
Figure 6
Figure 6. Figure 6 [PITH_FULL_IMAGE:figures/full_fig_p014_6.png] view at source ↗
read the original abstract

A conforming discontinuous Galerkin (DG) finite element method has been introduced in [21] on simplicial meshes, which has the flexibility of using discontinuous approximation and the simplicity in formulation of the classic continuous finite element method. The goal of this paper is to extend the conforming DG finite element method in \cite{cdg1} so that it can work on general polytopal meshes by designing weak gradient $\nabla_w$ appropriately. Two different conforming DG formulations on polytopal meshes are introduced which handle boundary conditions differently. Error estimates of optimal order are established for the corresponding conforming DG approximation in both a discrete $H^1$ norm and the $L^2$ norm. Numerical results are presented to confirm the theory.

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 / 0 minor

Summary. The manuscript extends the conforming discontinuous Galerkin method of [21] from simplicial to general polytopal meshes by constructing a weak gradient operator ∇_w. Two formulations are introduced that differ in boundary-condition treatment. Optimal-order error estimates are claimed in a discrete H¹ seminorm and the L² norm, with numerical results presented for confirmation.

Significance. A conforming DG scheme that retains optimal convergence on arbitrary polytopal meshes while preserving the formulation simplicity of continuous FEM would be a useful addition to the literature on polytopal discretizations, particularly if the weak-gradient construction avoids hidden regularity assumptions beyond those stated for the mesh family.

major comments (2)
  1. [Abstract] Abstract (paragraph on extension from [21]): the central claim that an appropriate ∇_w renders both formulations conforming on general polytopal meshes is load-bearing for the optimal-order estimates in the discrete H¹ and L² norms; the manuscript must explicitly record the precise mesh hypotheses (star-shapedness, uniform shape-regularity, etc.) under which the local approximation properties of ∇_w are proved, because any implicit restriction would invalidate uniformity of the constants in the error bounds.
  2. The derivation of the error estimates (presumably in the sections following the formulation) relies on the consistency and stability properties induced by ∇_w; without the explicit definition of ∇_w and the accompanying approximation lemma, it is impossible to verify that the discrete H¹ seminorm controls the error without additional mesh assumptions.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thoughtful review and for highlighting the importance of mesh hypotheses and clarity in the weak-gradient construction. We address each major comment below and will revise the manuscript accordingly where needed.

read point-by-point responses

  1. Referee: [Abstract] Abstract (paragraph on extension from [21]): the central claim that an appropriate ∇_w renders both formulations conforming on general polytopal meshes is load-bearing for the optimal-order estimates in the discrete H¹ and L² norms; the manuscript must explicitly record the precise mesh hypotheses (star-shapedness, uniform shape-regularity, etc.) under which the local approximation properties of ∇_w are proved, because any implicit restriction would invalidate uniformity of the constants in the error bounds.

    Authors: We agree that the mesh hypotheses must be stated explicitly to guarantee uniformity of constants. The local approximation properties of ∇_w are proved under the standard assumptions that each polytopal element is star-shaped with respect to a ball whose radius is bounded below by a fixed fraction of the element diameter and that the mesh family is uniformly shape-regular. These hypotheses are already used in the proofs but were not restated in the abstract. In the revision we will add a concise sentence to the abstract and to the opening of Section 3 that records these precise conditions. revision: yes

  2. Referee: [—] The derivation of the error estimates (presumably in the sections following the formulation) relies on the consistency and stability properties induced by ∇_w; without the explicit definition of ∇_w and the accompanying approximation lemma, it is impossible to verify that the discrete H¹ seminorm controls the error without additional mesh assumptions.

    Authors: The weak gradient ∇_w is defined explicitly in Section 3 as the unique solution of a local variational problem on each element; the approximation lemma (Lemma 3.1) then quantifies its accuracy under the star-shapedness and shape-regularity assumptions already mentioned. Consistency of the method follows directly from the definition of ∇_w, while stability in the discrete H¹ seminorm is proved in Theorem 4.1. We will add forward references from the error-analysis sections (4 and 5) to Section 3 and Lemma 3.1 so that the logical chain is immediately visible. revision: partial

Circularity Check

0 steps flagged

No circularity; extension via new weak-gradient design and error analysis is independent

full rationale

The paper cites [21] only for the original simplicial-mesh idea and then introduces two new conforming DG formulations on polytopal meshes together with an explicit new construction of the weak gradient ∇_w. Optimal-order estimates in the discrete H¹ and L² norms are derived from this construction and standard approximation properties; no equation or claim reduces a stated prediction to a fitted quantity or to a self-citation chain by construction. The derivation therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the ability to define a weak gradient that restores conformity on polytopal meshes; this is presented as a design choice rather than derived from prior axioms.

axioms (2)
  • domain assumption A weak gradient operator can be constructed on polytopal meshes that makes the DG method conforming while preserving optimal approximation properties.
    This is the key technical step stated in the abstract for extending the method from simplicial to general polytopal meshes.
  • standard math Standard Sobolev-space assumptions and mesh-regularity conditions hold for the polytopal meshes under consideration.
    Implicit in any finite-element error analysis; required for the stated optimal-order estimates.

pith-pipeline@v0.9.0 · 5642 in / 1364 out tokens · 26289 ms · 2026-05-25T12:24:20.409927+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. A conforming DG method for the biharmonic equation on polytopal meshes

    math.NA 2019-07 unverdicted novelty 6.0

    A conforming DG finite element method is developed for the biharmonic equation on polytopal meshes, with optimal error estimates established in a discrete H2 norm and L2 estimates ranging from sub-optimal to optimal b...

Reference graph

Works this paper leans on

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

  1. [1]

    Arnold, F

    D. Arnold, F. Brezzi, B. Cockburn and D. Marini, Unified analysis of discontinuous Galerkin methods for elliptic problems, SIAM J. Numer. Anal., 39 (2002), 1749-1779

    work page 2002

  2. [2]

    Babuˇ ska, The finite element method with penalty, Math

    I. Babuˇ ska, The finite element method with penalty, Math. Comp., 27 (1973), 221-228

    work page 1973

  3. [3]

    Brenner, L

    S. Brenner, L. Owens and L. Sung, A weakly over-penalized symmetric interior penalty method, Ele. Trans. Numer. Anal., 30 (2008), 107-127

    work page 2008

  4. [4]

    Cockburn, J

    B. Cockburn, J. Gopalakrishnan, and R. Lazarov, Unified hybridization of discontinuous Galerkin, mixed, and conforming Galerkin methods for second order elliptic problems, SIAM J. Numer. Anal., 47 (2009), 1319-136

    work page 2009

  5. [5]

    Cockburn and C

    B. Cockburn and C. Shu, The local discontinuous Galerkin finite element method for convection- diffusion systems, SIAM J. Numer. Anal., 35 (1998), 2440-2463

    work page 1998

  6. [6]

    Douglas Jr

    J. Douglas Jr. and T. Dupont, Interior penalty procedures for elliptic and parabolic Galerkin methods, Computing Methods in Applied Sciences, (1976), 207-216. 15 Table 6.3 Error profiles and convergence rates for (6.1) on polygonal grids shown in Figure 6.2 level ∥uh−Q0u∥ rate |||uh−u||| rate dim byP1 elements with strongly enforced boundary condition 6 0....

    work page 1976

  7. [7]

    Lipnikov, G

    K. Lipnikov, G. Manzini, F. Brezzi and A. Buffa, The mimetic finite difference method for the 3D magnetostatic field problems on polyhedral meshes, J. Comput. Phys., 230 (2011), 305-328

    work page 2011

  8. [8]

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

    work page 2018

  9. [9]

    L. Mu, J. Wang, and X. Ye, weak Galerkin finite element method for second-order elliptic problems on polytopal meshes, International Journal of Numerical Analysis and Modeling, 12 (2015), 31-53

    work page 2015

  10. [10]

    L. Mu, X. Wang and X. Ye, A modified weak Galerkin finite element method for the Stokes equations, J. Comput. Appl. Math., 275 (2015), 79-90

    work page 2015

  11. [11]

    L. Mu, J. Wang, Y. Wang and X. Ye, A weak Galerkin mixed finite element method for bihar- monic equations, Numerical Solution of Partial Differential Equations: Theory, Algorithms, and Their Applications, 45 (2013), 247-277

    work page 2013

  12. [12]

    L. Mu, J. Wang, and X. Ye, A weak Galerkin finite element method for biharmonic equations on polytopal meshes, Numerical Methods for Partial Differential Equations, 30 (2014), 1003-1029

    work page 2014

  13. [13]

    Pietro and A

    D. Pietro and A. Ern, Hybrid high-order methods for variable-diffusion problems on general meshes, Comptes Rendus Mathmatique, 353 (2015), 31-34

    work page 2015

  14. [14]

    Reed and T

    W. Reed and T. Hill. Triangular mesh methods for the neutron transport equation. Technical Report LA-UR-73-0479, Los Alamos Scientific Laboratory, Los Alamos, NM, 1973

    work page 1973

  15. [15]

    Wang and X

    J. Wang and X. Ye, A weak Galerkin finite element method for second-order elliptic problems. J. Comput. Appl. Math. 241 (2013), 103-115. 16 Table 6.4 Error profiles and convergence rates for (6.1) on polygonal grids shown in Figure 6.2 level ∥uh−Q0u∥ rate |||uh−u||| rate dim byP4 elements with strongly enforced boundary condition 2 0.7295E-04 3.68 0.4484E-0...

    work page 2013

  16. [16]

    Wang and X

    J. Wang and X. Ye, A Weak Galerkin mixed finite element method for second-order elliptic problems, Math. Comp., 83 (2014), 2101-2126

    work page 2014

  17. [17]

    X. Wang, N. Malluwawadu, F Gao and T. McMillan, A modified weak Galerkin finite element method, J. Comput. Appl. Math., 217 (2014), 319-327

    work page 2014

  18. [18]

    Wheeler, An elliptic collocation-finite element method with interior penalties

    M. Wheeler, An elliptic collocation-finite element method with interior penalties. SIAM J. Numer. Anal., 15 (1978), 152-161

    work page 1978

  19. [19]

    Ye and S

    X. Ye and S. Zhang, A conforming discontinuous Galerkin finite element method, International Journal of Numerical Analysis and Modeling, accepted

    work page