pith. sign in

arxiv: 2509.18908 · v3 · submitted 2025-09-23 · 🧮 math.NA · cs.NA

Novel Adaptive Schemes for Hyperbolic Conservation Laws

Shaoshuai Chu , Pingyao Feng , Vadim A. Kolotilov , Alexander Kurganov , Vladimir V. Ostapenko This is my paper

Pith reviewed 2026-05-18 14:44 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords adaptive schemeshyperbolic conservation lawssmoothness indicatorovercompressive limiterMinmod limitercentral-upwind schemesfinite-difference methodsshock capturing
0
0 comments X

The pith

Adaptive schemes for hyperbolic conservation laws use a smoothness indicator to apply overcompressive limiters only at contacts, Minmod2 elsewhere in rough zones, and a fifth-order scheme in smooth regions.

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

The paper presents new adaptive numerical methods for one- and two-dimensional hyperbolic systems of conservation laws. It builds on a smoothness indicator to locate rough solution areas, then selectively applies an overcompressive limiter to capture contact discontinuities sharply while using a dissipative Minmod2 limiter in other rough regions to suppress oscillations. In smooth areas the method switches to a quasi-linear fifth-order finite-difference scheme. This combination is designed to resolve nonlinear shocks and linearly degenerate contacts without generating spurious oscillations or staircase artifacts. The resulting schemes are tested on challenging examples to show the practical gains in resolution quality.

Core claim

The central claim is that a smoothness indicator can be used both to detect rough regions and to distinguish contact discontinuities within them, allowing the overcompressive limiter to be applied only where it sharpens contacts, the Minmod2 limiter to control oscillations in the remaining rough zones, and the fifth-order scheme to be used safely in smooth zones, thereby producing solutions that resolve shocks and contacts sharply while avoiding nonphysical structures.

What carries the argument

The smoothness indicator that both flags rough regions and classifies contact discontinuities, driving the switch among the overcompressive limiter, the Minmod2 limiter, and the quasi-linear fifth-order finite-difference scheme.

If this is right

  • The adaptive schemes extend the earlier adaption strategy by adding a targeted distinction between contact discontinuities and other rough features.
  • In smooth regions the fifth-order scheme can be used without risking oscillations near discontinuities.
  • Numerical examples confirm sharper resolution of both nonlinear shocks and linearly degenerate contacts compared with non-adaptive approaches.
  • The same smoothness indicator serves dual purposes: region detection and discontinuity-type classification.

Where Pith is reading between the lines

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

  • The approach could be tested on systems with multiple interacting waves to check whether the indicator still separates contacts reliably under wave collisions.
  • Similar indicator-driven switching might be explored for other high-order reconstructions or for adaptive mesh refinement strategies.
  • The method may reduce the need for manual tuning of limiter parameters in practical simulations of gas dynamics or shallow-water flows.

Load-bearing premise

The smoothness indicator must correctly identify contact discontinuities and separate them from other rough areas so that the appropriate limiter is chosen in each case.

What would settle it

A numerical test case in which the smoothness indicator mislabels a contact discontinuity, causing either visible oscillations or staircase structures to appear in the computed solution.

Figures

Figures reproduced from arXiv: 2509.18908 by Alexander Kurganov, Pingyao Feng, Shaoshuai Chu, Vadim A. Kolotilov, Vladimir V. Ostapenko.

Figure 2.1
Figure 2.1. Figure 2.1: Density ρ computed with ∆x = 1/40 by the OLD and MODIFIED schemes (left) and zoom at x ∈ [9, 10.2] (right). We then refine the mesh to ∆x = 1/200 and compute the numerical solutions until the same final time t = 5 by both the OLD and MODIFIED schemes. The obtained densities are presented in [PITH_FULL_IMAGE:figures/full_fig_p007_2_1.png] view at source ↗
Figure 2.2
Figure 2.2. Figure 2.2: Density ρ computed with ∆x = 1/200 by the OLD and MODIFIED schemes (left) and zoom at x ∈ [9.2, 10.4] (right). 3 Novel 1-D Adaptive Schemes In this section, we propose novel 1-D adaptive schemes, which achieve higher resolution than the adaptive schemes from [8] and at the same time do not suffer from the occurrence of either “stair-like” structures or large spurious oscillations. To this end, we first c… view at source ↗
Figure 2.3
Figure 2.3. Figure 2.3: Density ρ computed with ∆x = 1/40 (top row) and ∆x = 1/200 (bottom row) by the OLD and MODIFIED schemes (left) and zoom at x ∈ [4, 8] (right). using the following quantities, which are similar to (2.5)–(2.6): E p j = |pj+1 − 2pj + pj−1| |pj+1 − pj | + |pj − pj−1| + ε [PITH_FULL_IMAGE:figures/full_fig_p008_2_3.png] view at source ↗
Figure 5.1
Figure 5.1. Figure 5.1: Example 1: Density ρ computed on a uniform mesh with ∆x = 1/40 by the NEW and OLD schemes (left) and zoom at x ∈ [9, 9.8] (right) [PITH_FULL_IMAGE:figures/full_fig_p013_5_1.png] view at source ↗
Figure 5.2
Figure 5.2. Figure 5.2: Example 1: Density ρ computed on a uniform mesh with ∆x = 1/200 by the NEW and OLD schemes—zoom at x ∈ [4, 8] (left) and x ∈ [9.2, 10.4] (right). We then measure the CPU times consumed by both of the studied schemes on a uniform mesh with ∆x = 1/4000 and the CPU time consumed by the OLD scheme is about 356% larger than the CPU time consumed by the NEW scheme. This demonstrates that the NEW scheme is not … view at source ↗
Figure 5.3
Figure 5.3. Figure 5.3: Example 2: Density ρ computed by the NEW and OLD schemes (left) and zoom at x ∈ [−1, 0] (right). Example 3—Blast Wave Problem. In the last 1-D example, we consider a strong shocks interaction problem from [46] with the initial data, (ρ, u, p) [PITH_FULL_IMAGE:figures/full_fig_p014_5_3.png] view at source ↗
Figure 5.4
Figure 5.4. Figure 5.4: Example 3: Density ρ computed by the NEW and OLD schemes (left) and zoom at x ∈ [0.55, 0.85]. 5.2 2-D Examples We now turn to the 2-D numerical examples. In Examples 4–6, we will consider three different 2-D Riemann problems from [25]; see also [35, 36, 47]. In all of the examples below, we will show the computed densities and the corresponding “rough” areas detected by the SIs. In the plots obtained by … view at source ↗
Figure 5.5
Figure 5.5. Figure 5.5: Example 4: Density ρ computed by the OLD (left) and NEW (right) schemes [PITH_FULL_IMAGE:figures/full_fig_p015_5_5.png] view at source ↗
Figure 5.6
Figure 5.6. Figure 5.6: Example 4: The “rough” areas detected by the OLD (left) and NEW (right) schemes at the final time-step [PITH_FULL_IMAGE:figures/full_fig_p015_5_6.png] view at source ↗
Figure 5.7
Figure 5.7. Figure 5.7: Example 5: Density ρ computed by the LDCU (left), OLD (middle), and NEW (right) schemes [PITH_FULL_IMAGE:figures/full_fig_p016_5_7.png] view at source ↗
Figure 5.8
Figure 5.8. Figure 5.8: Example 5: The “rough” areas detected by the OLD (left) and NEW (right) schemes at the final time-step [PITH_FULL_IMAGE:figures/full_fig_p016_5_8.png] view at source ↗
Figure 5.10
Figure 5.10. Figure 5.10: Example 6: The “rough” areas detected by the OLD (left) and NEW (right) schemes at the final time-step [PITH_FULL_IMAGE:figures/full_fig_p017_5_10.png] view at source ↗
Figure 5.9
Figure 5.9. Figure 5.9: Example 6: Density ρ computed by the LDCU (left), OLD (middle), and NEW (right) schemes [PITH_FULL_IMAGE:figures/full_fig_p017_5_9.png] view at source ↗
Figure 5.11
Figure 5.11. Figure 5.11: Example 7: Density ρ computed by the OLD (left) and NEW (right) schemes. Next, we measure the CPU times consumed by the NEW and OLD schemes. The obtained results show that the CPU time consumed by the OLD scheme is about 122% larger than the CPU time consumed by the NEW scheme. This demostrates that like in the 1-D case, the 2-D NEW scheme is substantially more efficient than the OLD one though the effi… view at source ↗
Figure 5.12
Figure 5.12. Figure 5.12: Example 7: The “rough” areas detected by the OLD (left) and NEW (right) schemes at the final time-step. system: ρt + (ρu)x + (ρv)y = 0, (ρu)t + (ρu2 + p)x + (ρuv)y = 0, (ρv)t + (ρuv)x + (ρv2 + p)y = ρ, Et + [u(E + p)]x + [v(E + p)]y = ρv, and then use the following initial conditions: (ρ, u, v, p) [PITH_FULL_IMAGE:figures/full_fig_p019_5_12.png] view at source ↗
Figure 5.13
Figure 5.13. Figure 5.13: Example 8: Density ρ computed by the OLD and NEW adaptive schemes [PITH_FULL_IMAGE:figures/full_fig_p020_5_13.png] view at source ↗
Figure 5.14
Figure 5.14. Figure 5.14: Example 8: The “rough” areas detected by the OLD (left) and NEW (right) schemes at the final time-step. After splitting the computational domain into the corresponding three areas—“contact”, “rough”, and “smooth”—we apply different schemes inside each of them. In the “contact” areas, we use the low-dissipation central-upwind (LDCU) scheme with the overcompressive SBM limiter, which helps to extremely sh… view at source ↗
read the original abstract

We introduce new adaptive schemes for the one- and two-dimensional hyperbolic systems of conservation laws. Our schemes are based on an adaption strategy recently introduced in [{\sc S. Chu, A. Kurganov, and I. Menshov}, Appl. Numer. Math., 209 (2025)]. As there, we use a smoothness indicator (SI) to automatically detect ``rough'' parts of the solution and employ in those areas the second-order finite-volume low-dissipation central-upwind scheme with an overcompressive limiter, which helps to sharply resolve nonlinear shock waves and linearly degenerate contact discontinuities. In smooth parts, we replace the limited second-order scheme with a quasi-linear fifth-order (in space and third-order in time) finite-difference scheme, recently proposed in [{\sc V. A. Kolotilov, V. V. Ostapenko, and N. A. Khandeeva}, Comput. Math. Math. Phys., 65 (2025)]. However, direct application of this scheme may generate spurious oscillations near ``rough'' parts, while excessive use of the overcompressive limiter may cause staircase-like nonphysical structures in smooth areas. To address these issues, we employ the same SI to distinguish contact discontinuities, treated with the overcompressive limiter, from other ``rough'' regions, where we switch to the dissipative Minmod2 limiter. Advantage of the resulting adaptive schemes are clearly demonstrated on a number of challenging numerical examples.

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

1 major / 2 minor

Summary. The manuscript introduces adaptive schemes for one- and two-dimensional hyperbolic conservation laws. A smoothness indicator detects rough regions where a second-order central-upwind scheme with overcompressive limiter is applied to resolve shocks and contacts; a fifth-order finite-difference scheme is used in smooth regions. The same indicator further distinguishes contacts (overcompressive limiter) from other rough areas (Minmod2 limiter) to avoid oscillations and staircasing. Advantages are illustrated on several numerical test cases.

Significance. If the smoothness indicator reliably separates contacts, the adaptive combination could usefully merge high-order accuracy away from discontinuities with sharp, stable resolution at shocks and contacts. The work logically extends the cited prior schemes by the overlapping author groups and supplies visual evidence from challenging examples. Quantitative error norms, convergence rates, and baseline comparisons would substantially increase the significance.

major comments (1)
  1. [Section 3] The central adaptation logic relies on the smoothness indicator distinguishing contact discontinuities (for the overcompressive limiter) from other rough regions (for the Minmod2 limiter). No explicit formula, threshold, or additional sensor for this contact-specific classification is supplied, making it impossible to verify that misclassification will not reintroduce oscillations or staircasing and thereby undermine the claimed advantages.
minor comments (2)
  1. [Abstract] Abstract: the sentence beginning 'Advantage of the resulting adaptive schemes are' is grammatically incorrect and should be rephrased.
  2. [Numerical Experiments] The numerical examples would benefit from at least one table of L1 or L2 errors and observed convergence rates, even if only for the smooth test problems.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the positive overall assessment of our work and for the constructive major comment, which helps us improve the clarity of the adaptation strategy. We address the point below and will revise the manuscript accordingly.

read point-by-point responses

  1. Referee: [Section 3] The central adaptation logic relies on the smoothness indicator distinguishing contact discontinuities (for the overcompressive limiter) from other rough regions (for the Minmod2 limiter). No explicit formula, threshold, or additional sensor for this contact-specific classification is supplied, making it impossible to verify that misclassification will not reintroduce oscillations or staircasing and thereby undermine the claimed advantages.

    Authors: We thank the referee for this observation. The manuscript states that the same smoothness indicator (SI) is used both to detect rough regions and to further classify contacts (treated with the overcompressive limiter) versus other discontinuities (treated with the Minmod2 limiter). However, we agree that Section 3 does not supply the explicit formula for the SI, the numerical thresholds that trigger the contact-specific switch, or any auxiliary sensor that prevents misclassification. In the revised manuscript we will add a dedicated subsection in Section 3 that presents (i) the precise mathematical definition of the SI, (ii) the threshold values and decision logic used to label a rough cell as a contact, and (iii) a brief discussion of how the chosen thresholds were calibrated on the test problems to avoid re-introduction of oscillations or staircasing. These additions will make the classification procedure fully reproducible and will directly address the referee’s concern. revision: yes

Circularity Check

0 steps flagged

Minor self-citations to overlapping-author prior schemes, but new adaptive distinction and external numerical tests remain independent

full rationale

The paper explicitly builds on an adaptation strategy and smoothness indicator from Chu, Kurganov, Menshov (2025) and a fifth-order scheme from Kolotilov, Ostapenko, Khandeeva (2025), both with author overlap. However, the central contribution is the novel routing logic that applies the same SI to send contact discontinuities to the overcompressive limiter while routing other rough regions to the Minmod2 limiter, then switches to the high-order scheme in smooth areas. This combination is presented as new and its advantages are shown via direct numerical experiments on external test problems rather than any fitted parameter, self-referential definition, or uniqueness theorem imported from the cited works. No equation or claim reduces by construction to the inputs; the self-citations supply reusable components without bearing the load of the reported resolution and oscillation-avoidance properties.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The approach rests on the assumption that the smoothness indicator provides a reliable classification of solution features and that the cited base schemes behave as described near discontinuities; no new free parameters or invented entities are introduced.

axioms (1)
  • domain assumption The smoothness indicator accurately detects rough parts and distinguishes contact discontinuities from other discontinuities or oscillations.
    Invoked to decide which limiter to apply in rough regions.

pith-pipeline@v0.9.0 · 5810 in / 1304 out tokens · 41777 ms · 2026-05-18T14:44:56.536533+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

48 extracted references · 48 canonical work pages

  1. [1]

    Ar \`a ndiga, A

    F. Ar \`a ndiga, A. Baeza, and R. Donat , Vector cell-average multiresolution based on H ermite interpolation , Adv. Comput. Math., 28 (2008), pp. 1--22

    work page 2008

  2. [2]

    Ben-Artzi and J

    M. Ben-Artzi and J. Falcovitz , Generalized R iemann problems in computational fluid dynamics , vol. 11 of Cambridge Monographs on Applied and Computational Mathematics, Cambridge University Press, Cambridge, 2003

    work page 2003

  3. [3]

    Boscheri and M

    W. Boscheri and M. Dumbser , Arbitrary- L agrangian- E ulerian one-step WENO finite volume schemes on unstructured triangular meshes , Commun. Comput. Phys., 14 (2013), pp. 1174--1206

    work page 2013

  4. [4]

    S. Z. Burstein and A. A. Mirin , Third order difference methods for hyperbolic equations , J. Comput. Phys., 5 (1970), pp. 547--571

    work page 1970

  5. [5]

    Chertock, S

    A. Chertock, S. Chu, M. Herty, A. Kurganov, and M. Luk\'a c ov\'a-Medvi d ov\'a , Local characteristic decomposition based central-upwind scheme , J. Comput. Phys., 473 (2023). Paper No. 111718

    work page 2023

  6. [6]

    Chertock, S

    A. Chertock, S. Chu, and A. Kurganov , Adaptive high-order A-WENO schemes based on a new local smoothness indicator , E. Asian. J. Appl. Math., 13 (2023), pp. 576--609

    work page 2023

  7. [7]

    S. Chu, M. Herty, and A. Kurganov , Novel local characteristic decomposition based path-conservative central-upwind schemes , J. Comput. Phys., 524 (2025). Paper No. 113692

    work page 2025

  8. [8]

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

    work page 2025

  9. [9]

    S. Chu, A. Kurganov, and R. Xin , Low-dissipation central-upwind schemes for compressible multifluids , J. Comput. Phys., 518 (2024). Paper No. 113311

    work page 2024

  10. [10]

    S. Chu, A. Kurganov, and R. Xin , New low-dissipation central-upwind schemes. P art II , J. Sci. Comput., 193 (2025). Paper No. 33

    work page 2025

  11. [11]

    Dewar, A

    J. Dewar, A. Kurganov, and M. Leopold , Pressure-based adaption indicator for compressible E uler equations , Numer. Methods Partial Differential Equations, 31 (2015), pp. 1844--1874

    work page 2015

  12. [12]

    Dumbser, O

    M. Dumbser, O. Zanotti, R. Loub\`ere, and S. Diot , A posteriori subcell limiting of the discontinuous G alerkin finite element method for hyperbolic conservation laws , J. Comput. Phys., 278 (2014), pp. 47--75

    work page 2014

  13. [13]

    Fu and C.-W

    G. Fu and C.-W. Shu , A new troubled-cell indicator for discontinuous G alerkin methods for hyperbolic conservation laws , J. Comput. Phys., 347 (2017), pp. 305--327

    work page 2017

  14. [14]

    Z. Gao, Q. Liu, J. S. Hesthaven, B.-S. Wang, W. S. Don, and X. Wen , Non-intrusive reduced order modeling of convection dominated flows using artificial neural networks with application to R ayleigh- T aylor instability , Commun. Comput. Phys., 30 (2021), pp. 97--123

    work page 2021

  15. [15]

    Gelb and E

    A. Gelb and E. Tadmor , Spectral reconstruction of piecewise smooth functions from their discrete data , M2AN Math. Model. Numer. Anal., 36 (2002), pp. 155--175

    work page 2002

  16. [16]

    Gelb and E

    A. Gelb and E. Tadmor , Adaptive edge detectors for piecewise smooth data based on the minmod limiter , J. Sci. Comput., 28 (2006), pp. 279--306

    work page 2006

  17. [17]

    Godlewski and P.-A

    E. Godlewski and P.-A. Raviart , Numerical approximation of hyperbolic systems of conservation laws , Springer-Verlag, New York, second ed., 2021

    work page 2021

  18. [18]

    Gottlieb, D

    S. Gottlieb, D. Ketcheson, and C.-W. Shu , Strong stability preserving R unge- K utta and multistep time discretizations , World Scientific Publishing Co. Pte. Ltd., Hackensack, NJ, 2011

    work page 2011

  19. [19]

    Gottlieb, C.-W

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

    work page 2001

  20. [20]

    Guermond, R

    J.-L. Guermond, R. Pasquetti, and B. Popov , Entropy viscosity method for nonlinear conservation laws , J. Comput. Phys., 230 (2011), pp. 4248--4267

    work page 2011

  21. [21]

    J. S. Hesthaven , Numerical methods for conservation laws: From analysis to algorithms , Comput. Sci. Eng. 18, SIAM, Philadelphia, 2018

    work page 2018

  22. [22]

    Karni, A

    S. Karni, A. Kurganov, and G. Petrova , A smoothness indicator for adaptive algorithms for hyperbolic systems , J. Comput. Phys., 178 (2002), pp. 323--341

    work page 2002

  23. [23]

    V. A. Kolotilov, A. A. Kurganov, V. V. Ostapenko, N. A. Khandeeva, and S. Chu , On the accuracy of shock-capturing schemes calculating gas-dynamic shock waves , Comput. Math. Math. Phys., 63 (2023), pp. 1341--1349

    work page 2023

  24. [24]

    V. A. Kolotilov, V. V. Ostapenko, and N. A. Khandeeva , Fifth-order finite-difference scheme in space with increased accuracy in shock influence areas , Comput. Math. Math. Phys., 65 (2025), pp. 901--916

    work page 2025

  25. [25]

    Kurganov and E

    A. Kurganov and E. Tadmor , Solution of two-dimensional R iemann problems for gas dynamics without R iemann problem solvers , Numer. Methods Partial Differential Equations, 18 (2002), pp. 584--608

    work page 2002

  26. [26]

    M. E. Ladonkina, O. A. Nekliudova, V. V. Ostapenko, and V. F. Tishkin , On the accuracy of discontinuous Galerkin method calculating gas-dynamic shock waves , Dokl. Math., 107 (2023), pp. 120--125

    work page 2023

  27. [27]

    R. J. LeVeque , Finite Volume Methods for Hyperbolic Problems , Cambridge Texts in Appl. Math., Cambridge University Press, Cambridge, UK, 2002

    work page 2002

  28. [28]

    Lie and S

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

    work page 2003

  29. [29]

    Liska and B

    R. Liska and B. Wendrof , Comparison of several diference schemes on 1 D and 2 D test problems for the E uler equations , SIAM J. Sci. Comput., 25 (2003), pp. 995--1017

    work page 2003

  30. [30]

    L\"ohner , An adaptive finite element scheme for transient problems in CFD , Comput

    R. L\"ohner , An adaptive finite element scheme for transient problems in CFD , Comput. Methods Appl. Mech. Eng., 61 (1987), pp. 323--338

    work page 1987

  31. [31]

    Puppo and M

    G. Puppo and M. Semplice , Numerical entropy and adaptivity for finite volume schemes , Commun. Comput. Phys., 10 (2011), pp. 1132--1160

    work page 2011

  32. [32]

    Qiu and C.-W

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

    work page 2002

  33. [33]

    Qiu and C.-W

    J. Qiu and C.-W. Shu , A comparison of troubled-cell indicators for R unge- K utta discontinuous G alerkin methods using weighted essentially nonoscillatory limiters , SIAM J. Sci. Comput., 27 (2005), pp. 995--1013

    work page 2005

  34. [34]

    V. V. Rusanov , Third-order accurate shock-capturing schemes for computing discontinuous solutions , Dokl. Akad. Nauk SSSR., 18 (1968), pp. 1303--1305

    work page 1968

  35. [35]

    C. W. Schulz-Rinne , Classifcation of the R iemann problem for two-dimensional gas dynamics , SIAM J. Math. Anal., 24 (1993), pp. 76--88

    work page 1993

  36. [36]

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

    work page 1993

  37. [37]

    Shi, Y.-T

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

    work page 2003

  38. [38]

    C.-W. Shu , Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws , in Advanced numerical approximation of nonlinear hyperbolic equations ( C etraro, 1997), vol. 1697 of Lecture Notes in Math., Springer, Berlin, 1998, pp. 325--432

    work page 1997

  39. [39]

    Shu , Essentially non-oscillatory and weighted essentially non-oscillatory schemes , Acta Numer., 5 (2020), pp

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

    work page 2020

  40. [40]

    Shu and S

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

    work page 1988

  41. [41]

    Osher , Efficient implementation of essentially nonoscillatory shock-capturing schemes

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

    work page 1989

  42. [42]

    E. F. Toro , Riemann solvers and numerical methods for fluid dynamics: A practical introduction , Springer-Verlag, Berlin, Heidelberg, third ed., 2009

    work page 2009

  43. [43]

    M. J. Vuik and J. K. Ryan , Automated parameters for troubled-cell indicators using outlier detection , SIAM J. Sci. Comput., 38 (2016), pp. A84--A104

    work page 2016

  44. [44]

    Wang, C.-W

    W. Wang, C.-W. Shu, H. C. Yee, D. V. Kotov, and B. Sj\" o green , High order finite difference methods with subcell resolution for stiff multispecies discontinuity capturing , Commun. Comput. Phys., 17 (2015), pp. 317--336

    work page 2015

  45. [45]

    X. Wen, W. S. Don, Z. Gao, and J. S. Hesthaven , An edge detector based on artificial neural network with application to hybrid compact- WENO finite difference scheme , J. Sci. Comput., 83 (2020). Paper No. 49

    work page 2020

  46. [46]

    Woodward and P

    P. Woodward and P. Colella , The numerical solution of two-dimensional fluid flow with strong shocks , J. Comput. Phys., 54 (1988), pp. 115--173

    work page 1988

  47. [47]

    Zheng , Systems of conservation laws

    Y. Zheng , Systems of conservation laws. Two-dimensional Riemann problems , Progress in Nonlinear Differential Equations and their Applications, Birkh\" a user Boston, Inc., Boston, MA, 2001

    work page 2001

  48. [48]

    write newline

    " write newline "" before.all 'output.state := FUNCTION fin.entry add.period write newline FUNCTION new.block output.state before.all = 'skip after.block 'output.state := if FUNCTION not #0 #1 if FUNCTION and 'skip pop #0 if FUNCTION or pop #1 'skip if FUNCTION new.block.checka empty 'skip 'new.block if FUNCTION field.or.null duplicate empty pop "" 'skip ...

    work page