pith. sign in

arxiv: 2604.09498 · v1 · submitted 2026-04-10 · 🧮 math.NA · cs.NA

New Scheme Adaption Strategy for Hyperbolic Conservation Laws

Shaoshuai Chu , Michael Herty , Alexander Kurganov This is my paper

Pith reviewed 2026-05-10 16:16 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords scheme adaptionhyperbolic conservation lawsSBM limiterssmoothness indicatorEuler equationsnumerical dissipationadaptive limiters
0
0 comments X

The pith

A new scheme adaptation varies SBM limiter parameters continuously using a smoothness indicator to transition smoothly between compressive and dissipative regimes.

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

The paper replaces a sharp threshold for separating rough and smooth regions with continuous variation of one limiting parameter in SBM-type limiters. This activates compressive and overcompressive limiters only near shocks and contacts while gradually switching to dissipative limiters in smooth areas. One- and two-dimensional tests on the Euler equations of gas dynamics show the result is higher wave resolution and lower overall numerical dissipation than the prior threshold-based adaptive approach. The method still relies on the same smoothness indicator to guide the adaptation.

Core claim

By varying one limiting parameter continuously in SBM-type limiters according to the smoothness indicator, the scheme activates compressive and overcompressive limiters only near discontinuities while transitioning smoothly to dissipative limiters elsewhere, resulting in higher resolution and reduced numerical dissipation in one- and two-dimensional tests for the Euler equations.

What carries the argument

SBM-type limiters with one continuously varied limiting parameter, driven by the Lohner smoothness indicator to detect rough and smooth solution regions and produce a gradual switch.

If this is right

  • Sharper resolution of shocks and contacts in gas-dynamics simulations
  • Lower numerical dissipation away from discontinuities compared with fixed or threshold-based limiters
  • Stable performance in both one- and two-dimensional settings without abrupt limiter switches
  • Direct applicability to other hyperbolic systems that already use SBM-type limiters

Where Pith is reading between the lines

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

  • The continuous-adaptation idea could be tested on limiters other than SBM types to check generality
  • Long-time integration tests might show whether the reduced dissipation improves accuracy in problems with many wave interactions
  • Implementation cost remains comparable to the threshold method since only one parameter is adjusted

Load-bearing premise

That continuously varying one limiting parameter produces a stable smooth transition between compressive and dissipative regimes without introducing new oscillations or instabilities, relying on the smoothness indicator to correctly identify regions.

What would settle it

A standard Euler test such as the double Mach reflection or shock-tube problem in which the continuously adapted scheme produces visible oscillations or lower resolution than the sharp-threshold version.

read the original abstract

We introduce a new scheme adaption strategy for one- and two-dimensional hyperbolic systems of conservation laws. The proposed approach builds upon the adaptive framework introduced in [S. Chu, A. Kurganov, and I. Menshov, Appl. Numer. Math., 209 (2025), pp.155--170], where we first employed the smoothness indicator from [R. Lohner, Comput. Methods. Appl. Mech. Eng., 61 (1987), pp.323--338] to automatically detect ``rough'' and smooth parts of the computed solution, and then used different limiters in the detected regions. This adaptive strategy was based on a threshold needed to sharply separate ``rough'' and smooth regions. In this paper, we propose a different adaption strategy. We use SBM-type limiters and vary one of the limiting parameters continuously to allow a smooth transition between the ``rough'' and smooth areas. This way, compressive and overcompressive limiters are activated in the shock and contact wave vicinities only, while we gradually switch to dissipative limiters in the smooth regions. A series of one- and two-dimensional numerical tests for the Euler equations of gas dynamics demonstrates that the new scheme adaption strategy leads to a higher resolution and reduced numerical dissipation.

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

4 major / 0 minor

Summary. The manuscript introduces a new scheme adaptation strategy for one- and two-dimensional hyperbolic systems of conservation laws. It modifies the threshold-based adaptive framework of Chu, Kurganov, and Menshov (2025) by using SBM-type limiters and continuously varying one limiting parameter according to the Lohner smoothness indicator, thereby producing a smooth transition from compressive/overcompressive behavior near shocks and contacts to dissipative behavior in smooth regions. A series of numerical tests on the Euler equations is presented to demonstrate higher resolution and reduced numerical dissipation.

Significance. If the continuous-parameter transition can be shown to remain stable, the method would provide a threshold-free alternative to existing adaptive limiting strategies for hyperbolic conservation laws, potentially improving shock-capturing accuracy in gas-dynamics simulations without introducing additional oscillations. The work correctly identifies the role of the Lohner indicator but its impact remains modest until quantitative validation and stability arguments are supplied.

major comments (4)
  1. [Abstract] Abstract: the central claim that the new strategy 'leads to a higher resolution and reduced numerical dissipation' is unsupported by any quantitative error measures (L1 or L2 norms), baseline comparisons with fixed-limiter or threshold-based schemes, specified grid resolutions, or convergence rates.
  2. [Method description] Method description: no explicit functional form, equation, or algorithm is supplied for the continuous dependence of the SBM limiting parameter on the Lohner smoothness indicator, preventing reproduction and analysis of the claimed smooth transition.
  3. [Theoretical analysis] Theoretical analysis: the manuscript contains no proof or even discussion that continuous variation of the single limiting parameter preserves the TVD property, a maximum principle, or monotonicity, nor does it address how the transition avoids generating new extrema; this is load-bearing for the stability of the scheme on hyperbolic systems.
  4. [Numerical tests] Numerical tests section: the Euler-equation experiments provide no stability analysis, no discussion of how the smoothness indicator triggers the parameter variation without artifacts, and no comparison data that would substantiate the claimed improvement over the prior threshold-based approach.

Simulated Author's Rebuttal

4 responses · 1 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major comment point by point below, indicating where revisions will be made to improve clarity, reproducibility, and support for the claims.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the central claim that the new strategy 'leads to a higher resolution and reduced numerical dissipation' is unsupported by any quantitative error measures (L1 or L2 norms), baseline comparisons with fixed-limiter or threshold-based schemes, specified grid resolutions, or convergence rates.

    Authors: We agree that the abstract claim would be strengthened by quantitative support. In the revised manuscript we will add L1 and L2 error norms (where exact solutions are available), direct comparisons against both fixed-limiter schemes and the original threshold-based adaptive method of Chu et al. (2025), and results at multiple specified grid resolutions together with observed convergence rates. revision: yes

  2. Referee: [Method description] Method description: no explicit functional form, equation, or algorithm is supplied for the continuous dependence of the SBM limiting parameter on the Lohner smoothness indicator, preventing reproduction and analysis of the claimed smooth transition.

    Authors: The manuscript describes the continuous variation conceptually but does not supply an explicit formula. We will insert the precise functional dependence (the mapping from the Lohner indicator value to the SBM limiting parameter) together with the corresponding algorithmic steps in the revised method section to enable reproduction. revision: yes

  3. Referee: [Theoretical analysis] Theoretical analysis: the manuscript contains no proof or even discussion that continuous variation of the single limiting parameter preserves the TVD property, a maximum principle, or monotonicity, nor does it address how the transition avoids generating new extrema; this is load-bearing for the stability of the scheme on hyperbolic systems.

    Authors: The work is primarily numerical and does not contain a rigorous proof that the continuous parameter variation preserves TVD or the maximum principle. We will add a brief discussion explaining why the smooth transition is expected to avoid new extrema on the basis of the known properties of the underlying SBM limiters, while acknowledging that a full theoretical guarantee remains open. revision: partial

  4. Referee: [Numerical tests] Numerical tests section: the Euler-equation experiments provide no stability analysis, no discussion of how the smoothness indicator triggers the parameter variation without artifacts, and no comparison data that would substantiate the claimed improvement over the prior threshold-based approach.

    Authors: The presented tests demonstrate stable behavior across the chosen Euler problems, but we accept that explicit discussion and comparisons are missing. In revision we will add a short stability subsection, describe the indicator-triggered parameter variation and the absence of visible artifacts, and include side-by-side quantitative and visual comparisons with the threshold-based scheme of Chu et al. (2025). revision: yes

standing simulated objections not resolved
  • A rigorous proof that continuous variation of the limiting parameter preserves the TVD property or maximum principle for general hyperbolic systems.

Circularity Check

0 steps flagged

No significant circularity: new continuous-variation strategy independently proposed and tested numerically

full rationale

The paper cites its authors' prior threshold-based adaptive framework but explicitly modifies it by introducing continuous variation of one SBM-limiter parameter driven by the external Lohner smoothness indicator. The central claim of higher resolution and lower dissipation is established solely through one- and two-dimensional numerical experiments on the Euler equations, with no derivation that reduces the outcome to a fitted parameter, self-defined quantity, or load-bearing self-citation. No uniqueness theorems, ansatzes, or renamings of known results are invoked in a circular manner; the method and its validation remain self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The approach rests on the reliability of the Lohner smoothness indicator for region detection and on the ability of SBM-type limiters to be tuned continuously from compressive to dissipative behavior; no new entities are postulated.

free parameters (1)
  • limiting parameter
    One parameter inside the SBM-type limiters that is varied continuously to control the transition between regimes.
axioms (2)
  • domain assumption The smoothness indicator from Lohner (1987) reliably separates rough and smooth solution regions.
    The new strategy inherits this detector from the cited 2025 framework.
  • domain assumption SBM-type limiters admit a continuous parameterization that transitions from compressive to dissipative without loss of stability.
    Assumed to enable the smooth switch described in the abstract.

pith-pipeline@v0.9.0 · 5531 in / 1445 out tokens · 83757 ms · 2026-05-10T16:16:43.438960+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

20 extracted references · 20 canonical work pages

  1. [1]

    Chertock, A., Chu, S., Herty, M., Kurganov, A., Luk ´aˇcov´a-Medviˇdov´a, M.: Local characteristic decomposition based central-upwind scheme. J. Comput. Phys.473, Paper No. 111718 (2023)

    work page 2023

  2. [2]

    Chu, S., Kurganov, A., Menshov, I.: New adaptive low-dissipation central-upwind schemes. Appl. Numer. Math.209, 155–170 (2025)

    work page 2025

  3. [3]

    Chu, S., Kurganov, A., Xin, R.: New low-dissipation central-upwind schemes. Part II. J. Sci. Comput.103(1), Paper No. 33 (2025)

    work page 2025

  4. [4]

    World Scientific Publishing Co

    Gottlieb, S., Ketcheson, D., Shu, C.-W.: Strong stability preserving Runge-Kutta and multistep time discretizations. World Scientific Publishing Co. Pte. Ltd., Hackensack, NJ (2011)

    work page 2011

  5. [5]

    SIAM Rev.43(1), 89–112 (2001)

    Gottlieb, S., Shu, C.-W., Tadmor, E.: Strong stability-preserving high-order time discretization methods. SIAM Rev.43(1), 89–112 (2001)

    work page 2001

  6. [6]

    Johnsen, E.: On the treatment of contact discontinuities using WENO schemes. J. Comput. Phys. 230, 8665–8668 (2011)

    work page 2011

  7. [7]

    Kurganov, A., Tadmor, E.: Solution of two-dimensional Riemann problems for gas dynamics without Riemann problem solvers. Numer. Methods Partial Differential Equations18(5), 584– 608 (2002)

    work page 2002

  8. [8]

    Lie, K.-A., Noelle, S.: On the artificial compression method for second-order nonoscillatory central difference schemes for systems of conservation laws. SIAM J. Sci. Comput.24(4), 1157– 1174 (2003)

    work page 2003

  9. [9]

    L ¨ohner, R.: An adaptive finite element scheme for transient problems in CFD. Comput. Methods Appl. Mech. Eng.61(2), 323–338 (1987)

    work page 1987

  10. [10]

    Qiu, J., Shu, C.-W.: On the construction, comparison, and local characteristic decomposition for high-order central WENO schemes. J. Comput. Phys.183(1), 187–209 (2002)

    work page 2002

  11. [11]

    W.: Classification of the Riemann problem for two-dimensional gas dynamics

    Schulz-Rinne, C. W.: Classification of the Riemann problem for two-dimensional gas dynamics. SIAM J. Math. Anal.24, 76–88 (1993)

    work page 1993

  12. [12]

    W., Collins, J

    Schulz-Rinne, C. W., Collins, J. P., Glaz, H. M.: Numerical solution of the Riemann problem for two-dimensional gas dynamics. SIAM J. Sci. Comput.14, 1394–1414 (1993)

    work page 1993

  13. [13]

    Shi, J., Zhang, Y.-T., Shu, C.-W.: Resolution of high order WENO schemes for complicated flow structures. J. Comput. Phys.186, 690–696 (2003)

    work page 2003

  14. [14]

    Acta Numer.29, 701–762 (2020)

    Shu, C.-W.: Essentially non-oscillatory and weighted essentially non-oscillatory schemes. Acta Numer.29, 701–762 (2020)

    work page 2020

  15. [15]

    Shu, C.-W., Osher, S.: Efficient implementation of essentially non-oscillatory shock-capturing schemes. J. Comput. Phys.77, 439–471 (1988)

    work page 1988

  16. [16]

    Shu, C.-W., Osher, S.: Efficient implementation of essentially nonoscillatory shock-capturing schemes. II. J. Comput. Phys.83(1), 32–78 (1989)

    work page 1989

  17. [17]

    A., Toro, E

    Titarev, V. A., Toro, E. F.: WENO schemes based on upwind and centred TVD fluxes. Comput. & Fluids.34, 705–720 (2005)

    work page 2005

  18. [18]

    F., Titarev, V

    Toro, E. F., Titarev, V. A.: TVD fluxes for the high-order ADER schemes. J. Sci. Comput.24, 285–309 (2005)

    work page 2005

  19. [19]

    Woodward, P., Colella, P.: The numerical solution of two-dimensional fluid flow with strong shocks. J. Comput. Phys.54(1), 115–173 (1984)

    work page 1984

  20. [20]

    Two-dimensional Riemann problems

    Zheng, Y.: Systems of conservation laws. Two-dimensional Riemann problems. Birkh ¨auser Boston, Inc., Boston, MA (2001)

    work page 2001