pith. sign in

arxiv: 2605.19454 · v1 · pith:WBFAHLAInew · submitted 2026-05-19 · 🧮 math.NA · cs.NA

A Unified Transmissibility-Based Interior Penalty DG Method for Heterogeneous and Anisotropic Diffusion

Gregory Etangsale , Vincent Fontaine , Anis Younes This is my paper

Pith reviewed 2026-05-20 03:01 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords discontinuous Galerkininterior penaltyhybridized methodsheterogeneous diffusionanisotropic diffusiontransmissibilitya priori error estimatesnumerical stability
0
0 comments X

The pith

Eliminating the skeletal unknown from a hybridized interior penalty method produces a unified primal DG scheme stable for any diffusion contrast and anisotropy.

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

This paper derives a primal discontinuous Galerkin formulation for heterogeneous and anisotropic diffusion by exact algebraic elimination of the skeletal unknown from a compact hybridized interior penalty method. The resulting Unified Interior Penalty DG scheme incorporates transmissibility-based weights and two stabilization terms that act on the primal jump and the jump of the normal diffusive flux, scaling respectively with the harmonic mean and the inverse arithmetic mean of the face-wise transmissibilities. It unifies several earlier interior penalty approaches while proving that the scheme remains consistent, coercive, and bounded with stability properties independent of material contrast and anisotropy. Quasi-optimal a priori error estimates in the energy norm are established for all variants. A sympathetic reader would care because the construction supplies a parameter-robust numerical method for flow problems in complex media.

Core claim

The UIP-DG scheme is obtained by exact algebraic elimination of the skeletal unknown in a compact hybridized interior penalty method. It involves transmissibility-based weights together with stabilization terms on the primal jump and on the jump of the normal diffusive flux, which scale with the harmonic mean and the inverse arithmetic mean of the face-wise transmissibilities. This construction unifies several interior penalty approaches previously introduced independently and yields a robust method whose stability properties are independent of the diffusion contrast and anisotropy. Consistency, coercivity, and boundedness of the formulation are proved, together with quasi-optimal energynorm

What carries the argument

The Unified Interior Penalty DG (UIP-DG) scheme with transmissibility-based weights and dual stabilization terms obtained by exact algebraic elimination of the skeletal unknown from the hybridized interior penalty formulation.

If this is right

  • The UIP-DG scheme unifies several interior penalty DG methods that had been introduced independently.
  • Stability properties and quasi-optimal error estimates hold independently of the diffusion contrast and anisotropy.
  • The two stabilization terms ensure coercivity and boundedness for all variants of the scheme.
  • Numerical experiments confirm that the theoretical convergence rates are attained in practice.

Where Pith is reading between the lines

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

  • The transmissibility weighting suggests a direct link to finite-volume schemes used in reservoir simulation, potentially allowing hybrid implementations.
  • Elimination of the skeletal unknown may simplify code structure in existing DG solvers by removing the need to handle additional interface degrees of freedom.
  • The same algebraic-elimination technique could be applied to other hybridized DG methods for different elliptic operators.
  • Extending the analysis to time-dependent or nonlinear diffusion would test whether the contrast-independent properties carry over.

Load-bearing premise

The hybridized interior penalty formulation is well-defined and exact algebraic elimination of the skeletal unknown preserves the stability and consistency properties of the original hybrid method.

What would settle it

Compute the energy-norm error of the UIP-DG solution on a sequence of successively refined meshes for a heterogeneous anisotropic diffusion problem with large contrast; if the observed convergence rate falls below the predicted quasi-optimal rate, the stability and error claims are falsified.

Figures

Figures reproduced from arXiv: 2605.19454 by Anis Younes, Gregory Etangsale, Vincent Fontaine.

Figure 1
Figure 1. Figure 1: Two heterogeneous anisotropic benchmark configurations on the four-subdomain [PITH_FULL_IMAGE:figures/full_fig_p021_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: Test 1 – Representation of the discrete piecewise linear solution [PITH_FULL_IMAGE:figures/full_fig_p022_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Test 2 – (Left) A non-uniform triangular mesh. (Right) The SUIP-DG solution for [PITH_FULL_IMAGE:figures/full_fig_p024_4.png] view at source ↗
read the original abstract

We derive a primal discontinuous Galerkin (DG) formulation for heterogeneous and anisotropic diffusion, obtained by exact algebraic elimination of the skeletal unknown in a compact hybridized interior penalty (H-IP) method. The resulting Unified Interior Penalty DG (UIP-DG) scheme involves transmissibility-based weights inherited from the hybrid formulation, together with two stabilization terms acting respectively on the primal jump and on the jump of the normal diffusive flux. These penalties scale, respectively, with the harmonic mean and with the inverse arithmetic mean of the face-wise transmissibilities. This construction provides a unified perspective on several interior penalty approaches previously introduced independently, while yielding a robust method with stability properties independent of the diffusion contrast and anisotropy. We prove consistency, coercivity, and boundedness of the formulation, and derive quasi-optimal energy-norm a priori error estimates for all variants. Numerical experiments confirm the theoretical claims.

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

0 major / 3 minor

Summary. The manuscript derives a primal discontinuous Galerkin formulation for heterogeneous and anisotropic diffusion by exact algebraic elimination of the skeletal unknown from a hybridized interior penalty (H-IP) method. The resulting Unified Interior Penalty DG (UIP-DG) scheme incorporates transmissibility-based weights together with two stabilization terms, one acting on the primal jump (scaled by the harmonic mean of face-wise transmissibilities) and one on the jump of the normal diffusive flux (scaled by the inverse arithmetic mean). The authors prove consistency, coercivity, and boundedness of the formulation with stability properties independent of diffusion contrast and anisotropy, derive quasi-optimal energy-norm a priori error estimates for all variants, and present numerical experiments confirming the theoretical claims. This construction unifies several previously introduced interior penalty approaches.

Significance. If the claims hold, the work provides a valuable unified perspective on interior penalty DG methods for diffusion problems, delivering a robust scheme whose stability and error estimates are independent of contrast and anisotropy. The explicit algebraic elimination that preserves consistency, coercivity, and boundedness, together with the detailed proofs of these properties and the quasi-optimal estimates, constitute clear strengths. The transmissibility-based weights and the two distinct penalty scalings offer a coherent framework that could be useful in applications involving high-contrast or anisotropic media.

minor comments (3)
  1. [Theorem 4.3] In the statement of the main theorem on quasi-optimal estimates, the dependence of the hidden constant on the polynomial degree and mesh regularity parameters should be made explicit to clarify the range of applicability.
  2. [Section 5] The numerical experiments section would benefit from an additional table summarizing the observed convergence rates for the two stabilization variants across all test cases, including the high-anisotropy and high-contrast configurations.
  3. [Section 2.2] Notation for the face-wise transmissibility tensor T_f should be introduced once in the preliminaries and then used consistently; the current alternation between T_f and the harmonic/arithmetic means occasionally obscures the scaling arguments.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the positive assessment and recommendation for minor revision. The provided summary accurately captures the derivation of the UIP-DG scheme via algebraic elimination from the H-IP method, the transmissibility-based weights, the two stabilization terms, and the contrast- and anisotropy-independent stability and quasi-optimal error estimates.

read point-by-point responses

  1. Referee: No specific major comments are listed in the report; the referee summary and significance section are entirely positive and correctly describe the manuscript.

    Authors: We appreciate the recognition that the explicit elimination preserves consistency, coercivity, and boundedness, and that the transmissibility-based framework unifies prior interior-penalty approaches while delivering robust stability. Since no concrete issues or requested changes were identified, we see no need to alter the manuscript at this stage. revision: no

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper derives the primal UIP-DG formulation explicitly by algebraic elimination of the skeletal unknown from a compact hybridized interior penalty (H-IP) method. It supplies the consistency identity, coercivity proof, boundedness argument, and quasi-optimal error estimates directly in the manuscript, with stability shown independent of contrast and anisotropy via the transmissibility-based weights and dual penalty scalings. No load-bearing step reduces to a fitted input renamed as prediction, self-definitional construction, or unverified self-citation chain; the hybridization is taken as a well-defined starting point whose elimination preserves the stated properties, making the central claims self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The paper introduces no new free parameters, axioms beyond standard assumptions on the diffusion tensor, or invented entities. The central construction relies on the existence of a well-posed hybridized IP method and the algebraic invertibility of the local skeletal coupling, both treated as background facts.

axioms (1)
  • domain assumption The diffusion tensor is symmetric positive definite and bounded away from zero and infinity on each element.
    Required for the transmissibility definitions and for coercivity of the bilinear form.

pith-pipeline@v0.9.0 · 5679 in / 1341 out tokens · 27261 ms · 2026-05-20T03:01:37.486991+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

20 extracted references · 20 canonical work pages

  1. [1]

    Riviere, Discontinuous Galerkin methods for solving elliptic and parabolic equations: theory and implementation, SIAM, 2008

    B. Riviere, Discontinuous Galerkin methods for solving elliptic and parabolic equations: theory and implementation, SIAM, 2008. 25

    work page 2008

  2. [2]

    D. A. Di Pietro, A. Ern, Mathematical Aspects of Discontinuous Galerkin Methods, Mathématiques et Applications, Springer, Berlin, 2011

    work page 2011

  3. [3]

    D. N. Arnold, F. Brezzi, B. Cockburn, L. D. Marini, Unified analysis of dis- continuous galerkin methods for elliptic problems, SIAM J. Numer. Anal. 39 (5) (2002) 1749–1779.doi:10.1137/S0036142901384162

    work page doi:10.1137/s0036142901384162 2002

  4. [4]

    Dryja, On discontinuous Galerkin methods for elliptic problems with discontinuous coefficients, Comput

    M. Dryja, On discontinuous Galerkin methods for elliptic problems with discontinuous coefficients, Comput. Methods Appl. Math. 3 (1) (2003) 76– 85.doi:10.2478/cmam-2003-0007

    work page doi:10.2478/cmam-2003-0007 2003

  5. [5]

    A. Ern, A. F. Stephansen, P. Zunino, A discontinuous Galerkin method with weighted averages for advection–diffusion equations with locally small and anisotropic diffusivity, IMA J. Numer. Anal. 29 (2) (2009) 235–256. doi:10.1093/imanum/drm050

    work page doi:10.1093/imanum/drm050 2009

  6. [6]

    Douglas, T

    J. Douglas, T. F. Dupont, Interior penalty procedures for elliptic and parabolic galerkin methods, in: Computing Methods in Applied Sciences, Vol. 58 of Lecture Notes in Physics, Springer, 1976, pp. 207–216

    work page 1976

  7. [7]

    Cockburn, J

    B. Cockburn, J. Gopalakrishnan, R. Lazarov, Unified hybridization of dis- continuous galerkin, mixed, and continuous galerkin methods for second order elliptic problems, SIAM J. Numer. Anal. 47 (2) (2009) 1319–1365. doi:10.1137/070706616

    work page doi:10.1137/070706616 2009

  8. [8]

    Lehrenfeld, Hybrid discontinuous galerkin methods for incompressible flow problems, Ph.D

    C. Lehrenfeld, Hybrid discontinuous galerkin methods for incompressible flow problems, Ph.D. thesis, RWTH Aachen University (2010)

    work page 2010

  9. [9]

    B. Cockburn, Static condensation, hybridization, and the devising of the HDG methods, in: Building Bridges: Connections and Challenges in Mod- ern Approaches to Numerical Partial Differential Equations, Springer In- ternational Publishing, Cham, 2016, Ch. 3, pp. 129–177. URLhttps://doi.org/10.1007/978-3-319-41640-3_5

    work page doi:10.1007/978-3-319-41640-3_5 2016

  10. [10]

    N. C. Nguyen, J. Peraire, B. Cockburn, An implicit high-order hybridizable discontinuous galerkin method for linear convection–diffusion equations, J. Comput. Phys. 228 (9) (2009) 3232–3254.doi:10.1016/j.jcp.2009.01. 030

    work page doi:10.1016/j.jcp.2009.01 2009

  11. [11]

    G. N. Wells, Analysis of an interface stabilized finite element method: the advection-diffusion-reaction equation, SIAM J. Numer. Anal. 49 (1) (2011) 87–109.doi:10.1137/090775464

    work page doi:10.1137/090775464 2011

  12. [12]

    M. S. Fabien, M. G. Knepley, B. M. Rivière, Families of interior penalty hybridizable discontinuous Galerkin methods for second order elliptic prob- lems, Journal of Numerical Mathematics 28 (3) (2020) 161–174.doi: 10.1515/jnma-2019-0027. 26

    work page doi:10.1515/jnma-2019-0027 2020

  13. [13]

    Burman, P

    E. Burman, P. Zunino, A domain decomposition method based on weighted interior penalties for advection-diffusion-reaction problems, SIAM J. Nu- mer. Anal. 44 (4) (2006) 1612–1638.doi:10.1137/050634736

    work page doi:10.1137/050634736 2006

  14. [14]

    Etangsale, M

    G. Etangsale, M. Fahs, V. Fontaine, N. Rajaonison, Improved error esti- mates of hybridizable interior penalty methods using a variable penalty for highly anisotropic diffusion problems, Computers & Mathematics with Applications 119 (2022) 89–99.doi:10.1016/j.camwa.2022.05.029

    work page doi:10.1016/j.camwa.2022.05.029 2022

  15. [15]

    Codina, J

    R. Codina, J. Principe, J. Baiges, Subscales on the element boundaries in the variational two-scale finite element method, Comput. Methods Appl. Mech. Engrg. 198 (2009) 838–852.doi:10.1016/j.cma.2008.10.020

    work page doi:10.1016/j.cma.2008.10.020 2009

  16. [16]

    S. C. Brenner, L. R. Scott, The Mathematical Theory of Finite Element Methods, 3rd Edition, Springer, 2008

    work page 2008

  17. [17]

    S. C. Brenner, Poincaré–Friedrichs inequalities for piecewiseH 1 func- tions, SIAM J. Numer. Anal. 41 (1) (2003) 306–324.doi:10.1137/ S0036142902401311

    work page 2003

  18. [18]

    Schöberl, C++11 implementation of finite elements in ngsolve, Institute for Analysis and Scientific Computing Report (2014)

    J. Schöberl, C++11 implementation of finite elements in ngsolve, Institute for Analysis and Scientific Computing Report (2014)

    work page 2014

  19. [19]

    URLhttps://ngsolve.org

    NGSolve development team, Ngsolve: High order finite element library (2024). URLhttps://ngsolve.org

    work page 2024

  20. [20]

    R. B. Kellogg, On the Poisson equation with intersecting interfaces, in: Applicable Analysis, Vol. 4, Taylor & Francis, 1975, pp. 101–129. 27

    work page 1975