pith. sign in

arxiv: 2604.25101 · v1 · submitted 2026-04-28 · 🧮 math.NA · cs.NA· math-ph· math.MP

A Posteriori Error Estimation for Parabolic Equations with Enriched Galerkin Finite Element Methods

Hyun-Geun Shin , Yi-Yung Yang , Sanghyun Lee This is my paper

Pith reviewed 2026-05-07 15:48 UTC · model grok-4.3

classification 🧮 math.NA cs.NAmath-phmath.MP MSC 65N3065N1565M15
keywords a posteriori error estimationenriched Galerkin methodparabolic equationsresidual-based estimatorsadaptive mesh refinementfinite element methodsreliability and efficiency
0
0 comments X

The pith

A residual-based estimator is reliable and efficient for enriched Galerkin discretizations of linear parabolic equations.

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

The paper establishes a residual-based a posteriori error estimation framework for the enriched Galerkin finite element method on linear parabolic equations. It proves that the estimator both bounds the true error from above and is bounded by the true error from below, under standard assumptions on meshes and solutions. The authors embed the estimator in an adaptive mesh refinement algorithm and verify through examples that it delivers efficient error reduction without prior knowledge of the exact solution. This matters because it supplies a practical tool for controlling approximation quality in time-dependent problems where uniform meshes waste computation on smooth regions.

Core claim

The paper proves that a residual-based a posteriori error estimator for the enriched Galerkin finite element discretization of linear parabolic equations is both reliable and efficient. It further shows that the estimator integrates directly into an adaptive mesh refinement strategy that produces meshes yielding reliable error control, as confirmed by numerical tests on several model problems.

What carries the argument

The residual-based a posteriori error estimator, which sums local residual contributions from the enriched Galerkin weak form over space-time elements to produce an upper and lower bound on the global error.

If this is right

  • The estimator can be inserted into existing adaptive codes to automatically refine meshes only where the local residual is large.
  • Efficiency of the estimator implies that the adaptive procedure stops when the true error falls below a prescribed tolerance.
  • The approach inherits the local conservation property of the underlying enriched Galerkin scheme, allowing the estimator to respect local mass balance.
  • Numerical examples demonstrate that the resulting adaptive meshes achieve comparable accuracy to uniform meshes at substantially lower degrees of freedom.

Where Pith is reading between the lines

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

  • The same residual construction could be tested on nonlinear parabolic equations to check whether the reliability proof extends without major changes.
  • Because enriched Galerkin methods are cheaper than discontinuous Galerkin methods, the estimator may reduce overall computational cost for long-time simulations compared with existing adaptive DG frameworks.
  • The framework supplies a template for deriving residual estimators for other locally conservative finite element schemes applied to evolution equations.

Load-bearing premise

The residual estimator remains reliable and efficient when the mesh satisfies standard shape-regularity conditions, the time-stepping scheme meets usual stability requirements, and the solution has sufficient regularity for the parabolic problem.

What would settle it

A concrete computation on a sequence of successively refined meshes where the ratio of the true error to the residual estimator either diverges to infinity or collapses to zero for a known exact solution of a linear parabolic equation would falsify reliability or efficiency.

Figures

Figures reproduced from arXiv: 2604.25101 by Hyun-Geun Shin, Sanghyun Lee, Yi-Yung Yang.

Figure 1
Figure 1. Figure 1: Example 1.1-1.2. Convergence tests with uniform and adaptive mesh. view at source ↗
Figure 2
Figure 2. Figure 2: Example 1.2. The h–adaptive mesh at t = 0.1, 0.25, and 0.5 in the fifth cycle. The convergence of the errors is illustrated in Figure 1a, and the convergence rate achieve the optimal order with 1.1 in ∥ · ∥l∞(0,T;H1 ) norm. We note the error with the adaptivity is less than the uniform mesh with the similar number of Dofs. The comparison is shown in Figure1a. Moreover, we emphasize that we obtain the bette… view at source ↗
Figure 3
Figure 3. Figure 3: Example 1.1-1.2 The Error estimations η N at t = 0.5 with uniform and adaptive mesh. 5.2.3. Example 1.3 Coarsening with tolerance and effectivity index In this experiment, we perform adaptive mesh refinement as described in Algorithm 1, incor￾porating coarsening and a tolerance of τ = 1 × 10−3 . The coarsening and refinement parameters 21 view at source ↗
Figure 4
Figure 4. Figure 4: , with a minimal mesh size of hmin = 2 −10 . (a) t = 0.1 (b) t = 0.25 (c) t = 0.5 view at source ↗
Figure 5
Figure 5. Figure 5: Example 1.3. The results of adaptive mesh with view at source ↗
Figure 6
Figure 6. Figure 6: Example 2. The numbers of degree of freedom at each time step that required to satisfy view at source ↗
Figure 7
Figure 7. Figure 7: Example 2. The h–adaptive mesh using EG-Q1 (the first row) and EG-Q2 (the second row) with τ = 5 × 10−4 , θcoarse = 40%, and θrefine = 20% at t = 0.1, 0.25, and 0.5 (columns). (a) Error (b) Effectivity index view at source ↗
Figure 8
Figure 8. Figure 8: Example 2. The error and effectivity index using EG-Q1 and EG-Q2 with adaptive mesh refinement with τ = 5 × 10−4 , θcoarse = 40%, and θrefine = 20%. 25 view at source ↗
read the original abstract

This paper introduces a novel a posteriori error estimation framework for the enriched Galerkin (EG) finite element method applied to linear parabolic equations. While the EG method has been recognized for its local conservation property and computational efficiency compared to discontinuous Galerkin methods, its mathematical analysis in the context of a posteriori error estimation for parabolic problems remains unexplored. In this work, we prove reliability and efficiency using the residual-based approach. Furthermore, we integrate these error estimators into an adaptive mesh refinement strategy, demonstrating their effectiveness in achieving efficient and reliable error control through several numerical examples. The proposed approach provides a significant advancement in the mathematical foundation and practical applicability of the EG method for time-dependent problems.

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

Summary. The manuscript develops a residual-based a posteriori error estimator for the enriched Galerkin finite element discretization of linear parabolic equations. It proves reliability and efficiency of the estimator under standard mesh regularity, time-stepping, and solution regularity assumptions in appropriate Bochner spaces, and demonstrates its use within an adaptive mesh refinement algorithm through several numerical examples that show bounded effectivity indices.

Significance. If the proofs hold, the work supplies a missing theoretical foundation for reliable and efficient adaptive computations with the enriched Galerkin method on time-dependent problems. The local conservation property of EG is thereby made practically usable via rigorous error control, and the combination of standard residual decomposition with Clément-type interpolants adapted to the enrichment space plus numerical confirmation of effectivity indices constitutes a solid contribution to the numerical analysis of parabolic problems.

minor comments (2)
  1. [Abstract] The abstract states that reliability and efficiency are proved but does not name the precise mesh-regularity or time-step restrictions under which the constants remain independent of h and Δt; a single sentence in the abstract or introduction clarifying these assumptions would improve accessibility.
  2. [Numerical experiments] In the numerical section, the effectivity indices are reported to remain bounded, yet the tables or figures do not explicitly tabulate the dependence on the enrichment degree or on the time-step size; adding one column or subplot would make the efficiency claim more transparent.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript and the recommendation for minor revision. The referee's summary accurately reflects the development of a residual-based a posteriori error estimator for enriched Galerkin discretizations of linear parabolic equations, along with the proofs of reliability and efficiency and the numerical demonstration of adaptive mesh refinement.

Circularity Check

0 steps flagged

No significant circularity in the derivation chain

full rationale

The paper derives reliability and efficiency bounds for a residual-based a posteriori estimator on the enriched Galerkin discretization of linear parabolic equations. The proof proceeds via standard residual decomposition, integration by parts, and Clément-type interpolants adapted to the enrichment space, under explicitly stated assumptions on mesh regularity, time-stepping, and solution regularity in Bochner spaces. These steps are independent of the target result and do not reduce to fitted parameters, self-definitions, or load-bearing self-citations. Numerical examples validate the effectivity indices but are not part of the proof. The approach is a direct application of established residual techniques to a new discretization context, remaining self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

No free parameters, axioms, or invented entities are described in the abstract; standard finite-element assumptions on meshes and function spaces are presumed but not itemized.

pith-pipeline@v0.9.0 · 5420 in / 1094 out tokens · 40535 ms · 2026-05-07T15:48:12.731292+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

45 extracted references

  1. [1]

    A posteriori error estimation in finite element analysis, Computer Methods in Applied Mechanics and Engineering 142 (1) (1997) 1–88

    1997

  2. [2]

    S. Sun, J. Liu, A Locally Conservative Finite Element Method Based on Piecewise Constant Enrichment of the Continuous Galerkin Method, SIAM Journal on Scientific Computing 31 (4) (2009) 2528–2548

    2009

  3. [3]

    Becker, E

    R. Becker, E. Burman, P. Hansbo, M. Larson, A reduced P1-discontinuous Galerkin method (2004)

    2004

  4. [4]

    Lee, Y .-J

    S. Lee, Y .-J. Lee, M. F. Wheeler, A Locally Conservative Enriched Galerkin Approximation and Efficient Solver for Elliptic and Parabolic Problems, SIAM Journal on Scientific Computing 38 (3) (2016) A1404–A1429

    2016

  5. [5]

    S. Lee, M. F. Wheeler, Adaptive enriched Galerkin methods for miscible displacement problems with entropy residual stabilization, Journal of Computational Physics 331 (2017) 19–37

    2017

  6. [6]

    S. Lee, M. F. Wheeler, Enriched galerkin methods for two-phase flow in porous media with capillary pressure, Journal of Computational Physics 367 (2018) 65–86

    2018

  7. [7]

    Kadeethum, H

    T. Kadeethum, H. M. Nick, S. Lee, F. Ballarin, Enriched galerkin discretization for modeling poroelasticity and permeability alteration in heterogeneous porous media, Journal of Computational Physics 427 (2021) 110030

    2021

  8. [8]

    Lee, S.-Y

    S. Lee, S.-Y . Yi, Locking-free and locally-conservative enriched galerkin method for poroelasticity, Journal of Scientific Computing 94 (1) (2023) 26

    2023

  9. [9]

    J. Choo, S. Lee, Enriched galerkin finite elements for coupled poromechanics with local mass conservation, Com- puter Methods in Applied Mechanics and Engineering 341 (2018) 311–332

    2018

  10. [10]

    Kuzmin, S

    D. Kuzmin, S. Lee, Y .-Y . Yang, Bound-preserving and entropy stable enriched galerkin methods for nonlinear hyperbolic equations, Journal of Computational Physics (2025) 114323

    2025

  11. [11]

    S.-Y . Yi, S. Lee, Physics-preserving enriched galerkin method for a fully-coupled thermo-poroelasticity model, Numerische Mathematik 156 (3) (2024) 949–978

    2024

  12. [12]

    Kadeethum, H

    T. Kadeethum, H. Nick, S. Lee, F. Ballarin, Flow in porous media with low dimensional fractures by employing enriched galerkin method, Advances in Water Resources 142 (2020) 103620

    2020

  13. [13]

    S. Lee, A. Mikeli ´c, M. F. Wheeler, T. Wick, Phase-field modeling of proppant-filled fractures in a poroelastic medium, Computer Methods in Applied Mechanics and Engineering 312 (2016) 509–541

    2016

  14. [14]

    Adjerid, J

    S. Adjerid, J. E. Flaherty, I. Babu ˇska, A posteriori error estimation for the finite element method-of-lines solution of parabolic problems, Mathematical Models and Methods in Applied Sciences 9 (02) (1999) 261–286

    1999

  15. [15]

    A. V . Gaevskaya, S. I. Repin, A Posteriori Error Estimates for Approximate Solutions of Linear Parabolic Prob- lems., Differential Equations 41 (7) (2005)

    2005

  16. [16]

    Picasso, Adaptive finite elements for a linear parabolic problem, Computer Methods in Applied Mechanics and Engineering 167 (3-4) (1998) 223–237

    M. Picasso, Adaptive finite elements for a linear parabolic problem, Computer Methods in Applied Mechanics and Engineering 167 (3-4) (1998) 223–237

    1998

  17. [17]

    Z. Chen, J. Feng, An adaptive finite element algorithm with reliable and efficient error control for linear parabolic problems, Mathematics of computation 73 (247) (2004) 1167–1193

    2004

  18. [18]

    Morin, R

    P. Morin, R. H. Nochetto, K. G. Siebert, Data oscillation and convergence of adaptive FEM, SIAM Journal on Numerical Analysis 38 (2) (2000) 466–488

    2000

  19. [19]

    Makridakis, R

    C. Makridakis, R. H. Nochetto, Elliptic reconstruction and a posteriori error estimates for parabolic problems, SIAM journal on numerical analysis 41 (4) (2003) 1585–1594

    2003

  20. [20]

    Lakkis, C

    O. Lakkis, C. Makridakis, Elliptic reconstruction and a posteriori error estimates for fully discrete linear parabolic problems, Mathematics of computation 75 (256) (2006) 1627–1658

    2006

  21. [21]

    Bansch, F

    E. Bansch, F. Karakatsani, C. Makridakis, A posteriori error control for fully discrete Crank–Nicolson schemes, SIAM Journal on numerical analysis 50 (6) (2012) 2845–2872

    2012

  22. [22]

    B ¨ansch, F

    E. B ¨ansch, F. Karakatsani, C. Makridakis, The effect of mesh modification in time on the error control of fully discrete approximations for parabolic equations, Applied numerical mathematics 67 (2013) 35–63

    2013

  23. [23]

    Karakatsani, A posteriori error estimates for fully discrete fractional-stepϑ-approximations for parabolic equa- tions, IMA Journal of Numerical Analysis 36 (3) (2016) 1334–1361

    F. Karakatsani, A posteriori error estimates for fully discrete fractional-stepϑ-approximations for parabolic equa- tions, IMA Journal of Numerical Analysis 36 (3) (2016) 1334–1361

    2016

  24. [24]

    O. J. Sutton, Long-time L ∞(L2) a posteriori error estimates for fully discrete parabolic problems, IMA Journal of Numerical analysis 40 (1) (2020) 498–529

    2020

  25. [25]

    W. Wang, L. Yi, Delay-dependent elliptic reconstruction and optimal L ∞(L2) a posteriori error estimates for fully discrete delay parabolic problems, Mathematics of Computation 91 (338) (2022) 2609–2643

    2022

  26. [26]

    E. H. Georgoulis, C. G. Makridakis, Lower bounds, elliptic reconstruction and a posteriori error control of parabolic problems, IMA Journal of Numerical Analysis 43 (6) (2023) 3212–3242

    2023

  27. [27]

    T. Linβ, M. Ossadnik, G. Radojev, A unified approach to maximum-norm a posteriori error estimation for second- order time discretizations of parabolic equations, IMA Journal of Numerical Analysis 44 (3) (2024) 1644–1659. 26

    2024

  28. [28]

    Akrivis, C

    G. Akrivis, C. Makridakis, R. Nochetto, A posteriori error estimates for the Crank–Nicolson method for parabolic equations, Mathematics of computation 75 (254) (2006) 511–531

    2006

  29. [29]

    Kim, E.-J

    D. Kim, E.-J. Park, Adaptive Crank-Nicolson methods with dynamic finite-element spaces for parabolic problems, Discrete and Continuous Dynamical Systems-B 10 (4) (2008) 873–886

    2008

  30. [30]

    J. S. Gupta, R. K. Sinha, G. M. M. Reddy, J. Jain, A posteriori error analysis of the Crank–Nicolson finite element method for linear parabolic interface problems: A reconstruction approach, Journal of Computational and Applied Mathematics 340 (2018) 173–190

    2018

  31. [31]

    Cangiani, E

    A. Cangiani, E. H. Georgoulis, O. J. Sutton, Adaptive non-hierarchical Galerkin methods for parabolic problems with application to moving mesh and virtual element methods, Mathematical Models and Methods in Applied Sciences 31 (04) (2021) 711–751

    2021

  32. [32]

    J. Dai, L. Chen, M. Yang, A posteriori error estimates of the weak Galerkin finite element methods for parabolic problems, Journal of Computational and Applied Mathematics 445 (2024) 115822

    2024

  33. [33]

    T. Ray, R. K. Sinha, An adaptive immersed finite element method for linear parabolic interface problems with nonzero flux jump, Calcolo 60 (2) (2023) 21

    2023

  34. [34]

    Schmidt, K

    A. Schmidt, K. Siebert, The finite element toolbox ALBERTA (2005)

    2005

  35. [35]

    E. H. Georgoulis, O. Lakkis, J. M. Virtanen, A posteriori error control for discontinuous Galerkin methods for parabolic problems, SIAM journal on numerical analysis 49 (2) (2011) 427–458

    2011

  36. [36]

    Verf ¨urth, A posteriori error estimation techniques for finite element methods, OUP Oxford, 2013

    R. Verf ¨urth, A posteriori error estimation techniques for finite element methods, OUP Oxford, 2013

    2013

  37. [37]

    Rivi `ere, Discontinuous Galerkin methods for solving elliptic and parabolic equations: theory and implementa- tion, SIAM, 2008

    B. Rivi `ere, Discontinuous Galerkin methods for solving elliptic and parabolic equations: theory and implementa- tion, SIAM, 2008

    2008

  38. [38]

    M. F. Wheeler, An Elliptic Collocation-Finite Element Method with Interior Penalties, SIAM Journal on Numerical Analysis 15 (1) (1978) 152–161

    1978

  39. [39]

    Dawson, S

    C. Dawson, S. Sun, M. F. Wheeler, Compatible algorithms for coupled flow and transport, Computer Methods in Applied Mechanics and Engineering 193 (23) (2004) 2565–2580

    2004

  40. [40]

    S. Sun, M. F. Wheeler, Symmetric and nonsymmetric discontinuous Galerkin methods for reactive transport in porous media, SIAM Journal on Numerical Analysis 43 (1) (2005) 195–219

    2005

  41. [41]

    Rivi `ere, M

    B. Rivi `ere, M. F. Wheeler, V . Girault, A Priori Error Estimates for Finite Element Methods Based on Discontinuous Approximation Spaces for Elliptic Problems, SIAM Journal on Numerical Analysis 39 (3) (2002) 902–931

    2002

  42. [42]

    Dolejˇs´ı, M

    V . Dolejˇs´ı, M. Feistauer, Discontinuous Galerkin method: Analysis and Applications to Compressible flow, V ol. 48, Springer, 2015

    2015

  43. [43]

    J. T. Oden, M. Ainsworth, A posteriori error estimation in finite element analysis, Pure and Applied Mathematics: A Wiley Series of Texts, Monographs and Tracts, Wiley-Interscience, 2011

    2011

  44. [44]

    Arndt, W

    D. Arndt, W. Bangerth, M. Feder, M. Fehling, R. Gassm ¨oller, T. Heister, L. Heltai, M. Kronbichler, M. Maier, P. Munch, J.-P. Pelteret, S. Sticko, B. Turcksin, D. Wells, Thedeal.IILibrary, version 9.4, Journal of Numerical MathematicsAccepted (2022)

    2022

  45. [45]

    D ¨orfler, A convergent adaptive algorithm for poisson’s equation, SIAM Journal on Numerical Analysis 33 (3) (1996) 1106–1124

    W. D ¨orfler, A convergent adaptive algorithm for poisson’s equation, SIAM Journal on Numerical Analysis 33 (3) (1996) 1106–1124. 27

    1996