Local error estimates for a finite element method combining linear and nonlinear stabilization for the linear hyperbolic transport equation

Erik Burman; Fabian Heimann

arxiv: 2604.21775 · v1 · submitted 2026-04-23 · 🧮 math.NA · cs.NA

Local error estimates for a finite element method combining linear and nonlinear stabilization for the linear hyperbolic transport equation

Erik Burman , Fabian Heimann This is my paper

Pith reviewed 2026-05-09 20:58 UTC · model grok-4.3

classification 🧮 math.NA cs.NA

keywords finite element methodhyperbolic transport equationcontinuous interior penaltyartificial diffusionlocal error estimateslocalization principlestabilizationGalerkin method

0 comments

The pith

Combined linear and nonlinear stabilization localizes errors in finite element solutions of the hyperbolic transport equation.

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

The paper proves that a continuous Galerkin finite element method for the linear hyperbolic transport equation, using continuous interior penalty stabilization on smooth regions paired with nonlinear artificial diffusion on rough regions, satisfies a localization principle. This principle yields weighted stability and a priori error estimates of order O(h^{k+1/2}) in the L2 norm at final time, with the estimates depending only on local regularity. The analysis is semi-discrete and shows that errors generated in low-regularity zones remain confined relative to the convection direction and do not degrade accuracy in distant smooth zones when the stabilization is applied locally. A sympathetic reader cares because the result removes the need for globally smooth solutions to obtain optimal local convergence rates.

Core claim

The central claim is that the combined stabilization satisfies a localization principle in terms of weighted stability, from which follow a priori error bounds that scale as O(h^{k+1/2}) in the L2(Ω) norm at time T. The proof assumes the solution has sufficient local regularity only on the smooth subdomains where the linear continuous interior penalty term is active, while the nonlinear artificial diffusion term controls the remaining rough subdomains; under this setup, typical numerical errors generated in rough parts stay localized and do not pollute the smooth parts.

What carries the argument

The combination of linear continuous interior penalty stabilization (active on smooth parts) and nonlinear artificial diffusion (active on rough parts), which together produce the weighted stability estimates that enforce localization.

If this is right

Error bounds in smooth subdomains depend only on the local regularity there and are independent of the size or strength of rough features elsewhere.
Numerical errors generated near discontinuities or low-regularity zones remain confined and do not propagate to contaminate accuracy in smooth regions.
The standard O(h^{k+1/2}) convergence rate is recovered locally in any subdomain where the solution is sufficiently regular.
Stabilization parameters can be chosen region-wise according to local smoothness without losing the global localization property.

Where Pith is reading between the lines

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

The localization property suggests that adaptive mesh refinement can be restricted to rough features without needing global refinement to protect smooth accuracy.
Similar weighted-stability arguments may extend to other first-order hyperbolic systems or to time-dependent problems with moving interfaces, provided local regularity assumptions can be maintained.
The approach offers a practical route to shock-capturing in convection-dominated flows while preserving high-order accuracy away from shocks.

Load-bearing premise

The continuous solution must possess sufficient local regularity on the smooth subdomains where the continuous interior penalty stabilization is applied, with artificial diffusion correctly activated on the rough subdomains.

What would settle it

A numerical experiment on a domain containing one localized rough patch (such as a discontinuity) and one distant smooth region would falsify the claim if the L2 error measured only in the smooth region at final time fails to converge at rate O(h^{k+1/2}) or exhibits visible pollution from the rough patch.

read the original abstract

In this paper, we investigate the combination of a linear continuous interior penalty type and a non-linear artificial diffusion stabilisation applied to the transport problem, based on continuous Galerkin finite elements in space. This method was recently introduced and analysed for globally smooth solutions in [Burman 2023, SIAM J. Sci. Comput., 45, 1, A96-A122]. We provide a rigorous proof of a localisation principle in terms of weighted stability and a priori error bound results, which follow the widely known $\mathcal{O}(h^{k+1/2})$ scaling in the $L^2(\Omega; t=T)$ norm, where $k$ denotes the polynomial order of the finite element space and $h$ the mesh size. The analysis is semi-discrete in space and assumes sufficient local regularity of the continuous solution on the smooth part of the domain, where the continuous interior penalty stabilisation is active, whilst artificial diffusion operates on the remaining rough parts of the domain. Thereby, the analysis demonstrates that typical numerical errors in the rough part stay localised relative to the convection velocity and do not negatively affect the smooth parts of the solution, if the stabilisation combination is set up accordingly.

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.

Desk Editor's Note private letter to a colleague

This paper localizes the error analysis for the mixed CIP-plus-artificial-diffusion method, but the localization holds only under an a priori smooth/rough partition with assumed local regularity.

read the letter

The main thing to know is that Burman and Heimann have extended their 2023 combined stabilization method with a localization principle. They prove weighted stability and recover the usual O(h^{k+1/2}) a priori bound in the final-time L2 norm, showing that errors generated where artificial diffusion is active stay localized relative to the convection velocity and do not pollute the CIP-treated smooth regions. The proof follows standard energy methods with weights to control the cross terms between the two subdomains. This is a legitimate incremental advance over the global-smoothness case in the earlier paper. The analysis is semi-discrete and states its assumptions clearly: the continuous solution must have enough local regularity exactly where CIP is applied, and the two stabilizations must be placed on a known non-overlapping partition aligned with the flow. The paper does this part cleanly and without circularity. The soft spots are the assumptions themselves. The partition into smooth and rough regions has to be known in advance from the exact solution, and if regularity breaks near the interfaces the weighted estimates may not close. That is a real practical limitation even though the paper does not hide it. No time discretization is treated, and the stabilization parameters remain free choices that still need tuning. This is for specialists working on stabilized finite elements for hyperbolic transport, especially those thinking about local error control or adaptivity. A reader already familiar with CIP and artificial diffusion will get the most from the localization step. It deserves a serious referee because the technical extension is concrete and the proof strategy is proportionate to the claims.

Referee Report

1 major / 2 minor

Summary. The paper analyzes a hybrid stabilization approach for the linear hyperbolic transport equation discretized with continuous Galerkin finite elements in space. It combines linear continuous interior penalty (CIP) stabilization on smooth regions with nonlinear artificial diffusion on rough regions, extending prior global analysis for smooth solutions. Under local regularity assumptions on the smooth parts, the authors prove weighted stability and a priori error bounds of order O(h^{k+1/2}) in the L^2(Ω) norm at final time, claiming that errors from the rough region remain localized relative to the convection direction and do not pollute the smooth region when the stabilization regions are chosen accordingly.

Significance. If the localization principle is rigorously established, the result is significant for the design of adaptive or hybrid stabilization techniques in hyperbolic transport problems. It offers a pathway to handle localized discontinuities or low-regularity features without global error pollution, which is relevant for applications involving shocks or rough data. The work builds directly on the 2023 global analysis and recovers the standard suboptimal rate for CIP-type methods.

major comments (1)

[Abstract and weighted stability analysis] The localization result is derived under the assumption of a known, non-overlapping a priori partition into smooth/rough regions aligned with the convection velocity, with CIP applied only where sufficient local regularity holds. The abstract explicitly states that the analysis demonstrates localization 'if the stabilisation combination is set up accordingly,' but does not detail how interface cross terms in the weighted energy estimate are controlled when regularity fails near the partition boundary or when the partition is only approximate. This assumption appears load-bearing for the central claim that rough-part errors stay localized.

minor comments (2)

[Section on stability analysis] The semi-discrete setting is clear, but the transition from the continuous problem to the discrete weighted norms could be expanded for readability, particularly the precise definition of the weights used to localize the estimates.
[Introduction and conclusions] A brief remark on how the partition might be identified in practice (e.g., via a posteriori indicators) would help bridge the theoretical assumption to computational use, even if not required for the proof.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the detailed reading and for highlighting the importance of the localization principle. We address the major comment below, clarifying the role of the a priori partition and the control of interface terms in the weighted estimates.

read point-by-point responses

Referee: [Abstract and weighted stability analysis] The localization result is derived under the assumption of a known, non-overlapping a priori partition into smooth/rough regions aligned with the convection velocity, with CIP applied only where sufficient local regularity holds. The abstract explicitly states that the analysis demonstrates localization 'if the stabilisation combination is set up accordingly,' but does not detail how interface cross terms in the weighted energy estimate are controlled when regularity fails near the partition boundary or when the partition is only approximate. This assumption appears load-bearing for the central claim that rough-part errors stay localized.

Authors: The manuscript assumes an exact, non-overlapping a priori partition aligned with the convection direction, with the CIP stabilization active only on the smooth subdomain where the solution satisfies the required local regularity. The weighted energy estimate is constructed with weights that are constant along streamlines in the smooth region and increase in the upstream direction from the rough region; this choice ensures that any cross terms generated at the interface (from integration by parts in the transport term or from the stabilization forms) are absorbed by the local regularity assumption and the non-negativity of the weights. Because the partition boundaries are flow-aligned, information from the rough region cannot propagate upstream into the smooth region. The paper does not claim the result for approximate or misaligned partitions; the phrase 'if the stabilisation combination is set up accordingly' in the abstract refers precisely to this exact-partition hypothesis. We agree that the abstract is concise and will expand it in the revised version to explicitly mention the control of interface cross terms via the weighted estimates and flow alignment. revision: partial

Circularity Check

0 steps flagged

No significant circularity; localization proof is self-contained via energy estimates under explicit assumptions

full rationale

The paper derives weighted stability and O(h^{k+1/2}) a priori error bounds for the combined CIP/artificial-diffusion stabilization using standard energy methods on the semi-discrete formulation. The localization principle follows directly from the method definition, the non-overlapping partition of smooth/rough regions, and the assumed local regularity on the CIP-active subdomain; no step reduces by construction to a fitted parameter, self-referential definition, or load-bearing self-citation chain. The reference to Burman 2023 introduces the global-smooth case but is not invoked to justify the new local-regularity result. The derivation remains independent of the target error bound and is externally falsifiable via the stated assumptions and energy estimates.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard finite-element assumptions for hyperbolic problems plus paper-specific conditions on where each stabilizer is applied; no new physical entities are introduced.

free parameters (1)

stabilization parameters
Coefficients controlling the strength of the linear penalty and nonlinear diffusion are part of the method definition and must be chosen consistently with the localization requirement, though explicit values are not given in the abstract.

axioms (2)

domain assumption Sufficient local regularity of the continuous solution on the smooth part of the domain
Explicitly stated as required for the CIP-active region in the abstract.
ad hoc to paper The stabilization combination is set up accordingly so that rough-part errors stay localized
Condition stated in the abstract as necessary for the localization principle to hold.

pith-pipeline@v0.9.0 · 5513 in / 1401 out tokens · 42333 ms · 2026-05-09T20:58:56.378011+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

14 extracted references · 14 canonical work pages

[1]

S. C. Brenner and L. R. Scott , The mathematical theory of finite element methods , Springer, 2008, https://doi.org/10.1007/978-0-387-75934-0

work page doi:10.1007/978-0-387-75934-0 2008
[2]

A. N. Brooks and T. J. Hughes , Streamline upwind/petrov-galerkin formulations for convection dominated flows with particular emphasis on the incompressible navier-stokes equations , Computer Methods in Applied Mechanics and Engineering, 32 (1982), pp. 199--259, https://doi.org/https://doi.org/10.1016/0045-7825(82)90071-8, https://www.sciencedirect.com/sc...

work page doi:10.1016/0045-7825(82)90071-8 1982
[3]

E. Burman , Weighted error estimates for transient transport problems discretized using continuous finite elements with interior penalty stabilization on the gradient jumps , Vietnam Journal of Mathematics, 50 (2022), pp. 833--866, https://doi.org/10.1007/s10013-022-00550-x

work page doi:10.1007/s10013-022-00550-x 2022
[4]

E. Burman , Some observations on the interaction between linear and nonlinear stabilization for continuous finite element methods applied to hyperbolic conservation laws , SIAM Journal on Scientific Computing, 45 (2023), pp. A96--A122, https://doi.org/10.1137/21M1464154

work page doi:10.1137/21m1464154 2023
[5]

Burman and A

E. Burman and A. Ern , A continuous finite element method with face penalty to approximate Friedrichs' systems , ESAIM: Mod\'elisation math\'ematique et analyse num\'erique, 41 (2007), pp. 55--76, https://doi.org/10.1051/m2an:2007007, https://www.numdam.org/articles/10.1051/m2an:2007007/

work page doi:10.1051/m2an:2007007 2007
[6]

Burman and A

E. Burman and A. Ern , Continuous interior penalty hp-finite element methods for advection and advection-diffusion equations , Mathematics of Computation, 76 (2007), pp. 1119--1140

work page 2007
[7]

Burman, J

E. Burman, J. Guzmán, and D. Leykekhman , Weighted error estimates of the continuous interior penalty method for singularly perturbed problems , IMA Journal of Numerical Analysis, 29 (2008), pp. 284--314, https://doi.org/10.1093/imanum/drn001, https://doi.org/10.1093/imanum/drn001, https://arxiv.org/abs/https://academic.oup.com/imajna/article-pdf/29/2/284...

work page doi:10.1093/imanum/drn001 2008
[8]

Computer Methods in Applied Mechanics and Engineering , volume =

E. Burman and P. Hansbo , Edge stabilization for galerkin approximations of convection–diffusion–reaction problems , Computer Methods in Applied Mechanics and Engineering, 193 (2004), pp. 1437--1453, https://doi.org/https://doi.org/10.1016/j.cma.2003.12.032, https://www.sciencedirect.com/science/article/pii/S004578250400043X. Recent Advances in Stabilized...

work page doi:10.1016/j.cma.2003.12.032 2004
[9]

D. A. Di Pietro and A. Ern , Mathematical aspects of discontinuous Galerkin methods , vol. 69, Springer Science & Business Media, 2011, https://doi.org/10.1007/978-3-642-22980-0

work page doi:10.1007/978-3-642-22980-0 2011
[10]

Douglas and T

J. Douglas and T. Dupont , Interior penalty procedures for elliptic and parabolic galerkin methods , in Computing Methods in Applied Sciences, R. Glowinski and J. L. Lions, eds., Berlin, Heidelberg, 1976, Springer Berlin Heidelberg, pp. 207--216

work page 1976
[11]

Ern and J.-L

A. Ern and J.-L. Guermond , Weighting the edge stabilization , SIAM Journal on Numerical Analysis, 51 (2013), pp. 1655--1677, https://doi.org/10.1137/120867482, https://doi.org/10.1137/120867482, https://arxiv.org/abs/https://doi.org/10.1137/120867482

work page doi:10.1137/120867482 2013
[12]

Springer, Cham, [2021]©2021.https://doi.org/10.1007/978-3-030-56341-7

A. Ern and J.-L. Guermond , Finite Elements I , Springer Cham, 2021, https://doi.org/10.1007/978-3-030-56341-7

work page doi:10.1007/978-3-030-56341-7 2021
[13]

Finite Elements III

A. Ern and J.-L. Guermond , Finite Elements III , Springer Cham, 2021, https://doi.org/10.1007/978-3-030-57348-5

work page doi:10.1007/978-3-030-57348-5 2021
[14]

T. J. Hughes , The finite element method: linear static and dynamic finite element analysis , Courier Corporation, 2012

work page 2012

[1] [1]

S. C. Brenner and L. R. Scott , The mathematical theory of finite element methods , Springer, 2008, https://doi.org/10.1007/978-0-387-75934-0

work page doi:10.1007/978-0-387-75934-0 2008

[2] [2]

A. N. Brooks and T. J. Hughes , Streamline upwind/petrov-galerkin formulations for convection dominated flows with particular emphasis on the incompressible navier-stokes equations , Computer Methods in Applied Mechanics and Engineering, 32 (1982), pp. 199--259, https://doi.org/https://doi.org/10.1016/0045-7825(82)90071-8, https://www.sciencedirect.com/sc...

work page doi:10.1016/0045-7825(82)90071-8 1982

[3] [3]

E. Burman , Weighted error estimates for transient transport problems discretized using continuous finite elements with interior penalty stabilization on the gradient jumps , Vietnam Journal of Mathematics, 50 (2022), pp. 833--866, https://doi.org/10.1007/s10013-022-00550-x

work page doi:10.1007/s10013-022-00550-x 2022

[4] [4]

E. Burman , Some observations on the interaction between linear and nonlinear stabilization for continuous finite element methods applied to hyperbolic conservation laws , SIAM Journal on Scientific Computing, 45 (2023), pp. A96--A122, https://doi.org/10.1137/21M1464154

work page doi:10.1137/21m1464154 2023

[5] [5]

Burman and A

E. Burman and A. Ern , A continuous finite element method with face penalty to approximate Friedrichs' systems , ESAIM: Mod\'elisation math\'ematique et analyse num\'erique, 41 (2007), pp. 55--76, https://doi.org/10.1051/m2an:2007007, https://www.numdam.org/articles/10.1051/m2an:2007007/

work page doi:10.1051/m2an:2007007 2007

[6] [6]

Burman and A

E. Burman and A. Ern , Continuous interior penalty hp-finite element methods for advection and advection-diffusion equations , Mathematics of Computation, 76 (2007), pp. 1119--1140

work page 2007

[7] [7]

Burman, J

E. Burman, J. Guzmán, and D. Leykekhman , Weighted error estimates of the continuous interior penalty method for singularly perturbed problems , IMA Journal of Numerical Analysis, 29 (2008), pp. 284--314, https://doi.org/10.1093/imanum/drn001, https://doi.org/10.1093/imanum/drn001, https://arxiv.org/abs/https://academic.oup.com/imajna/article-pdf/29/2/284...

work page doi:10.1093/imanum/drn001 2008

[8] [8]

Computer Methods in Applied Mechanics and Engineering , volume =

E. Burman and P. Hansbo , Edge stabilization for galerkin approximations of convection–diffusion–reaction problems , Computer Methods in Applied Mechanics and Engineering, 193 (2004), pp. 1437--1453, https://doi.org/https://doi.org/10.1016/j.cma.2003.12.032, https://www.sciencedirect.com/science/article/pii/S004578250400043X. Recent Advances in Stabilized...

work page doi:10.1016/j.cma.2003.12.032 2004

[9] [9]

D. A. Di Pietro and A. Ern , Mathematical aspects of discontinuous Galerkin methods , vol. 69, Springer Science & Business Media, 2011, https://doi.org/10.1007/978-3-642-22980-0

work page doi:10.1007/978-3-642-22980-0 2011

[10] [10]

Douglas and T

J. Douglas and T. Dupont , Interior penalty procedures for elliptic and parabolic galerkin methods , in Computing Methods in Applied Sciences, R. Glowinski and J. L. Lions, eds., Berlin, Heidelberg, 1976, Springer Berlin Heidelberg, pp. 207--216

work page 1976

[11] [11]

Ern and J.-L

A. Ern and J.-L. Guermond , Weighting the edge stabilization , SIAM Journal on Numerical Analysis, 51 (2013), pp. 1655--1677, https://doi.org/10.1137/120867482, https://doi.org/10.1137/120867482, https://arxiv.org/abs/https://doi.org/10.1137/120867482

work page doi:10.1137/120867482 2013

[12] [12]

Springer, Cham, [2021]©2021.https://doi.org/10.1007/978-3-030-56341-7

A. Ern and J.-L. Guermond , Finite Elements I , Springer Cham, 2021, https://doi.org/10.1007/978-3-030-56341-7

work page doi:10.1007/978-3-030-56341-7 2021

[13] [13]

Finite Elements III

A. Ern and J.-L. Guermond , Finite Elements III , Springer Cham, 2021, https://doi.org/10.1007/978-3-030-57348-5

work page doi:10.1007/978-3-030-57348-5 2021

[14] [14]

T. J. Hughes , The finite element method: linear static and dynamic finite element analysis , Courier Corporation, 2012

work page 2012