pith. sign in

arxiv: 2606.24186 · v2 · pith:YZ4U4CUJnew · submitted 2026-06-23 · 🧮 math.NA · cs.NA

E Scheme and Flux-Limiter Scheme, Revisited

Huazhong Tang This is my paper

Pith reviewed 2026-06-30 10:22 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords E schemeflux-limiter schemescalar conservation lawdiscrete entropy conditionmonotone fluxnumerical entropy fluxsecond-order accuracyquasi-linear hyperbolic equation
0
0 comments X

The pith

For scalar conservation laws the E scheme satisfies the discrete entropy condition for any convex entropy, yet the numerical entropy flux is not unique and two-point monotone fluxes are E fluxes but not conversely.

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

The paper re-examines the E scheme introduced by Osher and the flux-limiter scheme of Sweby for quasi-linear hyperbolic conservation laws. It establishes that, apart from conservative monotone schemes, E schemes meet the discrete entropy condition for every convex entropy while leaving the numerical entropy flux non-unique. Two-point monotone fluxes always qualify as E fluxes, but the converse fails, and three-or-more-point E fluxes need not be monotone while multi-point monotone fluxes need not be E fluxes. The flux-limiter construction, which relies on an E-flux splitting combined with the Lax-Wendroff scheme, is shown not to guarantee second-order accuracy in both space and time.

Core claim

For a scalar conservation law, except for the conservative monotone schemes, the E scheme is a type of numerical methods that satisfy the discrete entropy condition for any convex entropy, but numerical entropy flux is not unique. Two-point monotone flux is E flux, but conversely it may not necessarily be correct. Moreover, multi-point (three or more points) E flux may not necessarily be monotone flux, and multi-point monotone flux may not necessarily be E flux. Sweby's flux-limiter scheme for the quasi-linear conservation laws was built on the E flux-based splitting and the LW scheme. It may not be second-order accurate in both space and time.

What carries the argument

The E flux, a numerical flux satisfying the discrete entropy inequality for every convex entropy pair, together with its relationship to monotone fluxes and the E-flux splitting used to construct flux-limiter schemes.

If this is right

  • Two-point monotone fluxes are always E fluxes.
  • Multi-point E fluxes are not required to be monotone.
  • Multi-point monotone fluxes are not required to be E fluxes.
  • The numerical entropy flux for an E scheme is not unique.
  • Sweby's flux-limiter scheme built on E-flux splitting need not achieve second-order accuracy in space and time.

Where Pith is reading between the lines

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

  • The non-uniqueness of the entropy flux may permit additional design freedom when extending E schemes to systems of conservation laws.
  • Loss of second-order accuracy in the flux-limiter scheme could degrade shock resolution in long-time simulations even when the underlying E flux is used.
  • The decoupling between E and monotone properties at three or more points suggests that entropy satisfaction alone does not guarantee TVD behavior for wider stencils.

Load-bearing premise

The analysis assumes that the original definitions and constructions of the E scheme and flux-limiter scheme remain the correct starting points and that the quasi-linear hyperbolic setting introduces no further constraints altering the stated relationships.

What would settle it

A concrete three-point numerical flux that satisfies the discrete entropy condition for every convex entropy yet violates monotonicity, or a grid-refinement study on a linear advection problem showing that the flux-limiter scheme converges only at first order near discontinuities.

read the original abstract

This paper revisits the {\em E scheme} of Osher \cite{Osher-SINUM1984} and the {\em flux-limiter scheme} of Sweby for quasi-linear hyperbolic conservation laws \cite{Sweby-SINUM1984}. Part of existing results will be re-understood and some new results will be presented. For a scalar conservation law, except for the conservative monotone schemes, the E scheme is a type of numerical methods that satisfy the discrete entropy condition for any convex entropy, but numerical entropy flux is not unique. Two-point monotone flux is E flux, but conversely it may not necessarily be correct. Moreover, multi-point (three or more points) E flux may not necessarily be monotone flux, and multi-point monotone flux may not necessarily be E flux. Sweby's flux-limiter scheme for the quasi-linear conservation laws was built on the E flux-based splitting $f_{j+1}-f_j=f_{j+1} { -\hat{f}^{\text{\tiny E}}_{j+\frac12}+\hat{f}^{\text{\tiny E}}_{j+\frac12}}-f_j$ and the LW scheme. It may not be second-order accurate in both space and time.

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

Summary. The paper revisits the E scheme of Osher (1984) and the flux-limiter scheme of Sweby (1984) for quasi-linear hyperbolic conservation laws. It re-interprets prior results to claim that, for scalar conservation laws, E schemes (except conservative monotone schemes) satisfy the discrete entropy condition for any convex entropy with non-unique numerical entropy flux; two-point monotone fluxes are E fluxes but not conversely; multi-point (three or more points) E fluxes need not be monotone and multi-point monotone fluxes need not be E fluxes; and Sweby's flux-limiter scheme, built on E-flux splitting and the Lax-Wendroff scheme, may not be second-order accurate in both space and time.

Significance. If the claimed distinctions between E and monotone fluxes hold for multi-point stencils and the accuracy limitation of the flux-limiter scheme is demonstrated, the work would clarify the scope of entropy-satisfying schemes and high-resolution methods for conservation laws, aiding the design of numerical fluxes beyond two-point stencils.

major comments (2)
  1. [Sections discussing multi-point E fluxes and monotone fluxes] The central claim that multi-point E fluxes need not be monotone (and conversely) requires an explicit multi-point generalization of the Osher E-flux inequality from the 1984 two-point definition; without quoting or deriving this generalized inequality in the relevant section, the separation between the two properties for stencils of width three or more is not secured.
  2. [Section on flux-limiter scheme] The assertion that Sweby's flux-limiter scheme may not be second-order accurate in both space and time for quasi-linear conservation laws needs a concrete derivation or counter-example showing order reduction, as the construction via E-flux splitting and LW scheme is stated without the supporting error analysis or numerical verification.
minor comments (1)
  1. [Abstract] The abstract mentions analysis for quasi-linear laws but the entropy discussion is restricted to the scalar case; clarify whether variable-coefficient terms alter the stated relationships.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thorough review and valuable comments on our manuscript. We address the two major comments point by point below. We believe these clarifications will improve the paper and plan to incorporate the necessary revisions.

read point-by-point responses

  1. Referee: [Sections discussing multi-point E fluxes and monotone fluxes] The central claim that multi-point E fluxes need not be monotone (and conversely) requires an explicit multi-point generalization of the Osher E-flux inequality from the 1984 two-point definition; without quoting or deriving this generalized inequality in the relevant section, the separation between the two properties for stencils of width three or more is not secured.

    Authors: We agree that an explicit multi-point generalization of the E-flux inequality is essential to rigorously establish the distinction for wider stencils. In the original Osher paper, the E-flux is defined for two-point interactions. Our manuscript extends this conceptually to multi-point by considering the flux that satisfies the entropy condition across multiple cells. To address this, we will add a derivation in the revised manuscript that generalizes the inequality to multi-point stencils, showing how the entropy dissipation holds without requiring monotonicity. This will secure the claim that multi-point E fluxes need not be monotone and vice versa. revision: yes

  2. Referee: [Section on flux-limiter scheme] The assertion that Sweby's flux-limiter scheme may not be second-order accurate in both space and time for quasi-linear conservation laws needs a concrete derivation or counter-example showing order reduction, as the construction via E-flux splitting and LW scheme is stated without the supporting error analysis or numerical verification.

    Authors: The manuscript states that the scheme 'may not be' second-order accurate due to the E-flux splitting combined with the Lax-Wendroff scheme for quasi-linear laws. To strengthen this, we will include either a theoretical error analysis demonstrating the order reduction or a specific numerical counter-example in the revised version. For instance, we can consider a specific quasi-linear equation where the limiter and splitting lead to first-order behavior in certain regimes. revision: yes

Circularity Check

0 steps flagged

No circularity; re-analysis of external 1984 definitions

full rationale

The paper's claims rest on re-examination of Osher (1984) and Sweby (1984) definitions for E-schemes and flux-limiter schemes, which are independent external sources. No self-citations appear, no parameters are fitted to data then renamed as predictions, and no equations reduce by construction to their own inputs. Multi-point extensions are presented as new interpretive results based on the original entropy condition rather than forced by redefinition or self-referential steps. The derivation chain is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The paper rests on the prior definitions of E scheme and flux-limiter scheme from the two 1984 references plus the standard domain assumption that the target equations are quasi-linear hyperbolic conservation laws.

axioms (1)
  • domain assumption The schemes operate on quasi-linear hyperbolic conservation laws.
    Explicitly stated in the abstract as the setting for all claims.

pith-pipeline@v0.9.1-grok · 5737 in / 1361 out tokens · 42779 ms · 2026-06-30T10:22:36.711342+00:00 · methodology

Review history (2 revisions) →

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

8 extracted references

  1. [1]

    Engquist and S

    B. Engquist and S. Osher, Stable and entropy condition satisfying approximations for transonic flow calculations,Math. Comp., 34(1980), pp.45-75

    1980

  2. [2]

    Harten, J.M

    A. Harten, J.M. Hyman, and P.D. Lax, On finite difference approximations and entropy conditions for shocks,Comm. Pure Appl. Math., 29(1976), pp.297-322

    1976

  3. [3]

    LeVeque,Numerical Methods for Conservation Laws, Birkhäuser, 1992, pp.145

    R.J. LeVeque,Numerical Methods for Conservation Laws, Birkhäuser, 1992, pp.145

    1992

  4. [4]

    E. M. Murman, Analysis of embedded shock waves calculated by relaxation methods,AIAA J., 12(1974), pp.626-633

    1974

  5. [5]

    Osher, Riemann solvers, the entropy condition, and difference approximations,SIAM J

    S. Osher, Riemann solvers, the entropy condition, and difference approximations,SIAM J. Numer. Anal., 21(1984), pp.217-235

    1984

  6. [6]

    Sweby, High resolution schemes using flux limiters for hyperbolic conservation laws,SIAM J

    P.K. Sweby, High resolution schemes using flux limiters for hyperbolic conservation laws,SIAM J. Numer. Anal., 21(1984), pp.995-1011

    1984

  7. [7]

    Wu and H.Z

    K.L. Wu and H.Z. Tang, High-order accurate physical-constraints-preserving finite difference WENO schemes for special relativistic hydrodynamics,J. Comput. Phys., 298(2015), pp.539-564

    2015

  8. [8]

    Xu and X.X

    Z.F. Xu and X.X. Zhang, Bound-Preserving High-Order Schemes, inHandbook of Numerical Analysis, Vol.18: HandbookofNumericalMethodsforHyperbolicProblemsBasicandFundamentalIssues, edited by Rémi Abgrall & Chi-Wang Shu, 2017, pp.81-102. 14

    2017