pith. machine review for the scientific record. sign in

arxiv: 2604.06546 · v1 · submitted 2026-04-08 · 💻 cs.CE · cs.NA· math.NA

Shocks without shock capturing: Information geometric regularization of finite volume methods for Navier--Stokes-like problems

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

classification 💻 cs.CE cs.NAmath.NA
keywords information geometric regularizationfinite volume methodsshock wavesNavier-Stokes equationscomputational fluid dynamicsregularizationshock capturingnumerical methods for conservation laws
0
0 comments X

The pith

Information geometric regularization embeds into finite volume methods to resolve shocks as smooth adjustable profiles without artificial viscosity or limiters.

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

The paper shows how information geometric regularization can be added directly to standard finite volume discretizations of Navier-Stokes-like equations. Shocks, which normally create numerical oscillations, are replaced by smooth profiles whose width is controlled by a single parameter while small-scale flow structures remain intact. The approach is tested on standard one- and two-dimensional benchmark problems that contain both smooth regions and discontinuities. In these tests the regularized schemes produce the expected physical solutions and reach accuracy levels comparable to established WENO and localized artificial diffusion methods. The regularization also reduces the number of memory accesses and arithmetic operations required per time step.

Core claim

Embedding information geometric regularization into finite-volume schemes replaces shock singularities with smooth profiles of adjustable width without dissipating fine-scale features such as turbulence or acoustics, yielding solutions that recover the expected behavior on all canonical benchmarks while remaining competitive in accuracy with WENO and LAD schemes and requiring fewer memory accesses and arithmetic operations per time step.

What carries the argument

Information geometric regularization (IGR), which replaces shock singularities with smooth profiles of adjustable width without dissipating fine-scale features.

If this is right

  • Existing finite volume codes can treat discontinuous flows without adding custom limiter or Riemann-solver logic.
  • Turbulent or acoustic structures near shocks remain undamped, allowing direct simulation of mixed smooth-discontinuous regimes.
  • Lower arithmetic and memory cost per step permits either larger grids or shorter run times for the same accuracy.
  • The regularization width parameter can be chosen independently of mesh size, separating physical smoothing from numerical resolution.

Where Pith is reading between the lines

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

  • The same regularization could be inserted into other conservation-law discretizations such as discontinuous Galerkin or spectral element methods.
  • Because fine-scale features are preserved, the method may improve accuracy in multi-physics problems that couple fluid shocks to acoustics or chemistry.
  • The underlying information-geometric construction might imply additional discrete invariance properties that help long-time conservation on unstructured meshes.

Load-bearing premise

That the IGR regularization term can be inserted into existing finite volume codes for Navier-Stokes-like equations and will maintain numerical stability and physical fidelity in every flow regime without further tuning or post-processing.

What would settle it

A high-Mach-number shock-tube or blast-wave simulation in which the IGR-regularized finite volume scheme develops growing oscillations or violates discrete conservation to a greater degree than a standard WENO scheme on the same mesh.

Figures

Figures reproduced from arXiv: 2604.06546 by Anand Radhakrishnan, Benjamin Wilfong, Florian Sch\"afer, Spencer H. Bryngelson.

Figure 1
Figure 1. Figure 1: Comparison of localized artificial diffusivity (LAD) and information geometric regularization (IGR) approaches for shock treatment (modified from [23], with author permission.) 2.4 A discrete alternative: Limiters Rather than regularizing the PDE, an alternative approach changes the numerical discretization to handle discontinuities. ENO schemes [28] adaptively select the smoothest stencil at each cell int… view at source ↗
Figure 2
Figure 2. Figure 2: Information geometric regularization (IGR). [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Comparison of high-order polynomial reconstruction strategies in FVM near a jump discontinuity (orange). Blue horizontal segments are cell averages; the orange arrow marks the target face. (a) Fifth-order reconstruction from smooth data yields a single consistent polynomial. (b) The same stencil near a jump produces oscillatory left/right reconstructions. (c) (W)ENO selects among candidate sub-stencils (do… view at source ↗
Figure 4
Figure 4. Figure 4: demonstrates convergence in these three regimes on a one-dimensional test case with initial condition u(x) = 1.5 sin(2πx), ρ(x) ≡ 1, e(x) = 4, distinguishing behavior before and after shock formation. For a fixed α, we observe second-order convergence before shock formation and almost second-order convergence after shock formation. Since we are using a second-order accurate discretization of the entropic p… view at source ↗
Figure 5
Figure 5. Figure 5: Sine wave propagation comparing velocity profiles at final time T = 0.4. (a) and (b) compare all methods at m = 150 and m = 300, respectively. (c) shows IGR and LAD at m = 300 with α = 10∆x2 and 100∆x2 . The reference solution (dashed blue) is computed with characteristic WENO and m = 104 . 0 0.2 0.4 0.6 0.8 1 0 2 4 ρ(x) (a) m = 200 Initial (Sharp) LAD IGR WENO (Char.) WENO (Comp.) Reference 0 0.2 0.4 0.6 … view at source ↗
Figure 6
Figure 6. Figure 6: Shu–Osher problem comparing density profiles at resolutions m = 200, 300, and 400. The top row (a–c) shows the full domain; the bottom row (d–f) shows the post-shock oscillatory region. All methods are compared against a reference solution computed with characteristic WENO at m = 5000. 13 [PITH_FULL_IMAGE:figures/full_fig_p013_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Sod shock tube convergence study. Rows 1–2 (a–f): density profiles, full and zoomed near the contact discontinuity. Rows 3–4 (g–l): specific internal energy profiles, full and zoomed. Columns correspond to (m, ε) = (200, 0.02), (400, 0.01), and (800, 0.005). WENO uses the sharp initial condition; LAD and IGR use smoothed initial conditions. All results are compared against the exact Riemann solver solution… view at source ↗
Figure 8
Figure 8. Figure 8: Leblanc shock tube convergence study showing specific internal energy profiles. Each row corresponds to a different smoothing parameter (ε = 0, 0.025, 0.05, 0.1) and each column shows a different resolution (m = 450, 900, 1800). Runs not shown amount to simulations that break down due to negative densities or energies. All results are compared against the exact Riemann solver solution for the sharp initial… view at source ↗
Figure 9
Figure 9. Figure 9: Riemann test problem with N = 5002 grid points using (a) IGR/LF, (b) WENO5/LF, and (c) WENO5/HLLC with a (d) reference vanishing viscosity solution using WENO5/LF at N = 50002 grid points. (a) T = 0.05 (b) T = 0.10 (c) T = 0.15 (d) T = 0.20 1 6 11 16 21 Density ( ρ) [PITH_FULL_IMAGE:figures/full_fig_p018_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Double Mach problem with mx × my = 800 × 200 grid points at (a) T = 0.05, (b) T = 0.10, (c) T = 0.15, and (d) T = 0.20 using IGR and tanh smoothing. 18 [PITH_FULL_IMAGE:figures/full_fig_p018_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Convergence rate for error ∥e∥∞ in pressure with grid size for IGR with varying α and WENO5. 0 5 10 15 20 25 30 10−4 10−3 10−2 10−1 100 T [periods] ∥e∥∞ (a) N = 2002 IGR (α = 10∆x2 ) IGR (α = 2∆x2 ) IGR (α = 0) WENO + HLLC WENO + LF 0 5 10 15 20 25 30 T [periods] (b) N = 4002 [PITH_FULL_IMAGE:figures/full_fig_p020_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Growth of error ∥e∥∞ in pressure for an advecting isentropic vortex with grid sizes (a) N = 2002 and (b) N = 4002 versus the number of periods across the domain. 20 [PITH_FULL_IMAGE:figures/full_fig_p020_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Turbulent kinetic energy cascade for IGR/LF, WENO5/HLLC, and WENO5/LF for N = 1283 , 2563 , and 5123 grid points, normalized by a WENO5/HLLC solution that used 20483 grid points. 6.2.9 3D multi-jet configuration Many-engine rockets, such as the SpaceX Super Heavy, are powered by many smaller engines than the industry standard at the time of writing. This choice comes with a key engineering challenge: high… view at source ↗
Figure 14
Figure 14. Figure 14: Application of IGR to Mach 10 rocket simulations with strong shocks. Panel (a) shows an artistic rendering of a 33-jet configuration. Panels (b) and (c) show a comparison between IGR and WENO simulations of a simpler three-jet configuration. The IGR result in (b) avoids the grid-aligned artifacts present in the (c) WENO case. Images adapted from Wilfong et al. [23] with author permission. 7 Comparison to … view at source ↗
read the original abstract

Shock waves in high-speed fluid dynamics produce near-discontinuities in the fluid momentum, density, and energy. Most contemporary works use artificial viscosity or limiters as numerical mitigation of the Gibbs--Runge oscillations that result from traditional numerics. These approaches face a delicate balance in achieving sufficiently regular solutions without dissipating fine-scale features, such as turbulence or acoustics. Recent work by Cao and Sch\"afer introduces information geometric regularization (IGR), the first inviscid regularization method for fluid dynamics. IGR replaces shock singularities with smooth profiles of adjustable width, without dissipating fine-scale features. This work provides a strategy for the practical use of IGR in finite-volume-based numerical methods. We illustrate its performance on canonical test problems and compare it against established approaches based on limiters and Riemann solvers. Results show that the finite volume IGR approach recovers the expected solutions in all cases. Across canonical benchmarks, IGR achieves accuracy competitive with WENO and LAD shock-capturing schemes in both smooth and discontinuous flow regimes. The IGR approach is computationally light, with meaningfully fewer memory accesses and arithmetic operations per time step.

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

Summary. The manuscript proposes a strategy for embedding information geometric regularization (IGR) into standard finite-volume discretizations of Navier-Stokes-like systems. It replaces shock singularities with smooth profiles whose width is controlled by an information-geometric length scale, claims that this recovers expected solutions on canonical benchmarks without traditional shock-capturing, achieves accuracy competitive with WENO and LAD schemes in both smooth and discontinuous regimes, and incurs meaningfully lower memory-access and arithmetic cost per time step.

Significance. If the central claims hold, the work would be significant for high-speed flow simulation: it offers an inviscid regularization that avoids dissipating fine-scale features while simplifying the numerical treatment of discontinuities and reducing per-step cost relative to limiters or artificial-viscosity methods. The approach extends the prior IGR framework of Cao and Schäfer into a practical FV setting.

major comments (1)
  1. [Abstract] Abstract: the central claim that IGR can be directly embedded 'without additional tuning' and remains 'computationally light' is load-bearing, yet the regularization length scale is described as adjustable. Without explicit demonstration (e.g., a single fixed width that works across all reported benchmarks without per-problem adjustment or post-processing), the method risks reintroducing a tunable dissipation scale comparable to existing shock-capturing techniques, directly weakening the advertised advantage.
minor comments (1)
  1. [Abstract] Abstract: the performance assertions (recovery of expected solutions, competitive accuracy, reduced operation count) are stated without reference to specific test cases, error norms, or implementation details; these should be summarized concisely so readers can gauge the scope of the validation.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful review and constructive feedback. We address the single major comment below.

read point-by-point responses
  1. Referee: [Abstract] Abstract: the central claim that IGR can be directly embedded 'without additional tuning' and remains 'computationally light' is load-bearing, yet the regularization length scale is described as adjustable. Without explicit demonstration (e.g., a single fixed width that works across all reported benchmarks without per-problem adjustment or post-processing), the method risks reintroducing a tunable dissipation scale comparable to existing shock-capturing techniques, directly weakening the advertised advantage.

    Authors: We agree that the regularization length scale is an adjustable parameter, as stated in the abstract. The manuscript's reference to embedding 'without additional tuning' is intended to contrast IGR with traditional approaches that require per-problem calibration of limiters, Riemann solvers, or artificial viscosity coefficients; the IGR length scale is instead chosen once relative to mesh size to set a target shock width and then held fixed for the entire benchmark suite. However, the referee correctly identifies that the current text does not explicitly demonstrate this consistency. We will therefore revise the abstract to remove any ambiguity around parameter selection and add a concise statement (with supporting detail in the results section) confirming that the same relative length scale was used without per-benchmark retuning or post-processing. revision: yes

Circularity Check

0 steps flagged

No significant circularity; central FV embedding and benchmark results are independent of prior IGR definition

full rationale

The paper's core contribution is a practical embedding strategy for IGR into standard finite-volume discretizations, with performance claims (accuracy competitive with WENO/LAD, fewer operations, recovery of expected solutions) resting on numerical experiments across canonical problems rather than any self-referential derivation. The abstract attributes the IGR concept itself to prior work by Cao and Schäfer, but this citation supplies the regularization mechanism as an external input; the present results do not reduce to it by construction, nor do any quoted equations or claims equate a 'prediction' to a fitted parameter or rename a known result. No load-bearing step collapses to self-definition, fitted-input-as-prediction, or an unverified self-citation chain. The derivation chain is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review supplies no explicit free parameters, axioms, or invented entities; IGR itself is referenced as prior work.

pith-pipeline@v0.9.0 · 5520 in / 1070 out tokens · 70347 ms · 2026-05-10T18:27:31.163289+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Thermodynamically Constrained Information Geometric Regularization for Compressible Flows

    math.NA 2026-04 unverdicted novelty 6.0

    A thermodynamic extension of information geometric regularization for compressible flows introduces an anisotropic stress tensor and an elliptic equation that mitigates cusp singularities in simulations while preservi...

Reference graph

Works this paper leans on

62 extracted references · 62 canonical work pages · cited by 1 Pith paper

  1. [1]

    Shu, Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws, in: B

    C.-W. Shu, Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws, in: B. Cockburn, C. Johnson, C.-W. Shu, E. Tadmor (Eds.), Advanced Numerical Approximation of Nonlinear Hyperbolic Equations, volume 1697 of Lecture Notes in Mathematics, Springer, 1998, pp. 325–432

  2. [2]

    Gottlieb, C.-W

    D. Gottlieb, C.-W. Shu, On the Gibbs phenomenon and its resolution, SIAM Review39(1997) 644–668

  3. [3]

    Canuto, M

    C. Canuto, M. Hussaini, A. Quarteroni, T. A. Zang, Spectral Methods: Fundamentals in Single Domains, Springer, 2006. 24

  4. [4]

    VonNeumann, R

    J. VonNeumann, R. D. Richtmyer, A method for the numerical calculation of hydrodynamic shocks, Journal of Applied Physics21(1950) 232–237

  5. [5]

    Puppo, Numerical entropy production for central schemes, SIAM Journal on Scientific Computing25(2004) 1382–1415

    G. Puppo, Numerical entropy production for central schemes, SIAM Journal on Scientific Computing25(2004) 1382–1415

  6. [6]

    A. W. Cook, W. H. Cabot, Hyperviscosity for shock-turbulence interactions, Journal of Computational Physics203(2005) 379–385

  7. [7]

    Fiorina, S

    B. Fiorina, S. K. Lele, An artificial nonlinear diffusivity method for supersonic reacting flows with shocks, Journal of Computational Physics222(2007) 246–264

  8. [8]

    A. Mani, J. Larsson, P. Moin, Suitability of artificial bulk viscosity for large-eddy simulation of turbulent flows with shocks, Journal of Computational Physics228(2009) 7368–7374

  9. [9]

    G. E. Barter, D. L. Darmofal, Shock capturing with PDE-based artificial viscosity for DGFEM: Part I. formulation, Journal of Computational Physics229(2010) 1810–1827

  10. [10]

    Guermond, R

    J.-L. Guermond, R. Pasquetti, B. Popov, Entropy viscosity method for nonlinear conservation laws, Journal of Computational Physics230(2011) 4248–4267

  11. [11]

    O. P. Bruno, J. S. Hesthaven, D. V. Leibovici, FC-based shock-dynamics solver with neural- network localized artificial-viscosity assignment, Journal of Computational Physics: X15 (2022) 100110

  12. [12]

    Dolejˇ s´ ı, M

    V. Dolejˇ s´ ı, M. Feistauer, C. Schwab, On some aspects of the discontinuous Galerkin finite element method for conservation laws, Mathematics and Computers in Simulation61(2003) 333–346

  13. [13]

    Glaubitz, P

    J. Glaubitz, P. ¨Offner, H. Ranocha, T. Sonar, Smooth and compactly supported viscous sub-cell shock capturing for discontinuous Galerkin methods, Journal of Scientific Computing79(2019) 249–272

  14. [14]

    Van Leer, Towards the ultimate conservative difference scheme

    B. Van Leer, Towards the ultimate conservative difference scheme. V. a second-order sequel to Godunov’s method, Journal of computational Physics32(1979) 101–136

  15. [15]

    D. Ray, J. S. Hesthaven, An artificial neural network as a troubled-cell indicator, Journal of computational physics367(2018) 166–191

  16. [16]

    X.-D. Liu, S. Osher, T. Chan, Weighted essentially non-oscillatory schemes, Journal of Computational Physics115(1994) 200–212

  17. [17]

    Harten, B

    A. Harten, B. Engquist, S. Osher, S. R. Chakravarthy, Uniformly high order accurate essentially non-oscillatory schemes, III, Journal of computational physics131(1997) 3–47

  18. [18]

    Guermond, M

    J.-L. Guermond, M. Nazarov, B. Popov, I. Tomas, Second-order invariant domain preserving approximation of the euler equations using convex limiting, SIAM Journal on Scientific Computing40(2018) A3211–A3239

  19. [19]

    D. Kuzmin, Monolithic convex limiting for continuous finite element discretizations of hyperbolic conservation laws, Computer Methods in Applied Mechanics and Engineering361(2020) 112804. 25

  20. [20]

    Lozano, Watch your adjoints! Lack of mesh convergence in inviscid adjoint solutions, AIAA Journal57(2019) 3991–4006

    C. Lozano, Watch your adjoints! Lack of mesh convergence in inviscid adjoint solutions, AIAA Journal57(2019) 3991–4006

  21. [21]

    Bodony, A

    D. Bodony, A. Fikl, Adjoint-based sensitivity of shock-laden flows, Proceedings of the 2022 CTR summer program (2022)

  22. [22]

    R. Cao, F. Sch¨ afer, Information geometric regularization of the barotropic Euler equation, arXiv preprint arXiv:2308.14127 (2023)

  23. [23]

    Wilfong, A

    B. Wilfong, A. Radhakrishnan, H. Le Berre, D. J. Vickers, T. Prathi, N. Tselepidis, B. Dorschner, R. Budiardja, B. Cornille, S. Abbott, F. Sch¨ afer, S. H. Bryngelson, Simulating many-engine spacecraft: Exceeding 1 quadrillion degrees of freedom via information geometric regularization, in: Proceedings of SC ’25: The International Conference for High Perf...

  24. [24]

    R. Cao, F. Sch¨ afer, Information geometric regularization of unidimensional pressureless Euler equations yields global strong solutions, arXiv preprint arXiv:2411.15121 (2024)

  25. [25]

    arXiv preprint arXiv:2512.13948 (2025)

    W. Barham, B. K. Tran, B. S. Southworth, F. Sch¨ afer, Hamiltonian information geometric regularization of the compressible euler equations, arXiv preprint arXiv:2512.13948 (2025)

  26. [26]

    Guelmame, D

    B. Guelmame, D. Clamond, S. Junca, Hamiltonian regularisation of the unidimensional barotropic euler equations, Nonlinear Analysis: Real World Applications64(2022) 103455

  27. [27]

    Koreeda, H

    J. Koreeda, H. Yanagisawa, K. Maeno, H. Honma, S. Bystrov, V. Ivanov, F. Shugaev, Front structures of strong shock waves in air, in: Shock Waves @ Marseille II: Physico-Chemical Processes and Nonequilibrium Flow, Springer, 1995, pp. 263–268

  28. [28]

    Harten, B

    A. Harten, B. Engquist, S. Osher, S. R. Chakravarthy, Uniformly high order accurate essentially non-oscillatory schemes, III, Journal of Computational Physics71(1987) 231–303

  29. [29]

    van Leer, Towards the ultimate conservative difference scheme

    B. van Leer, Towards the ultimate conservative difference scheme. V. A second-order sequel to Godunov’s method, Journal of Computational Physics32(1979) 101–136

  30. [30]

    van Leer, Towards the ultimate conservative difference scheme

    B. van Leer, Towards the ultimate conservative difference scheme. I. The quest of monotonicity, Lecture Notes in Physics18(1973) 163–168

  31. [31]

    V. Arnold, Sur la g´ eom´ etrie diff´ erentielle des groupes de Lie de dimension infinie et ses applications ` a l’hydrodynamique des fluides parfaits, Annales de l’institut Fourier16(1966) 319–361

  32. [32]

    Khesin, G

    B. Khesin, G. Misio lek, K. Modin, Geometric hydrodynamics and infinite-dimensional newton’s equations, Bulletin of the American Mathematical Society58(2021) 377–442

  33. [33]

    Harten, P

    A. Harten, P. D. Lax, B. van Leer, On upstream differencing and Godunov-type schemes for hyperbolic conservation laws, SIAM Review25(1983) 35–61

  34. [34]

    E. F. Toro, M. Spruce, W. Speares, Restoration of the contact surface in the HLL-Riemann solver, Shock Waves4(1994) 25–34

  35. [35]

    E. F. Toro, Riemann Solvers and Numerical Methods for Fluid Dynamics: A Practical Intro- duction, 3 ed., Springer, 2009. 26

  36. [36]

    Chauvat, J.-M

    Y. Chauvat, J.-M. Moschetta, J. Gressier, Shock wave numerical structure and the carbuncle phenomenon, International Journal for Numerical Methods in Fluids47(2005) 903–909

  37. [37]

    Gottlieb, C.-W

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

  38. [38]

    P. S. Rawat, X. Zhong, On high-order shock-fitting and front-tracking schemes for numerical simulation of shock–disturbance interactions, Journal of Computational Physics229(2010) 6744–6780

  39. [39]

    C.-W. Shu, S. Osher, Efficient implementation of essentially non-oscillatory shock-capturing schemes, Journal of Computational Physics77(1988) 439–471

  40. [40]

    G. A. Sod, A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws, Journal of Computational Physics27(1978) 1–31

  41. [41]

    Loub` ere, M

    R. Loub` ere, M. J. Shashkov, A subcell remapping method on staggered polygonal grids for arbitrary Lagrangian–Eulerian methods, Journal of Computational Physics209(2005) 105–138

  42. [42]

    Liang, L

    T. Liang, L. Fu, A new high-order shock-capturing TENO scheme combined with skew- symmetric-splitting method for compressible gas dynamics and turbulence simulation, Computer Physics Communications302(2024) 109236

  43. [43]

    Woodward, P

    P. Woodward, P. Colella, The numerical simulation of two-dimensional fluid flow with strong shocks, Journal of Computational Physics54(1984) 115–173

  44. [44]

    C.-W. Shu, Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws, in: Advanced Numerical Approximation of Nonlinear Hyperbolic Equations: Lectures given at the 2nd Session of the Centro Internazionale Matematico Estivo (CIME) held in Cetraro, Italy, June 23–28, 1997, Springer, 2006, pp. 325–432

  45. [45]

    S. C. Spiegel, H. Huynh, J. R. DeBonis, A survey of the isentropic Euler vortex problem using high-order methods, in: 22nd AIAA computational fluid dynamics conference, 2015, p. 2444

  46. [46]

    M. E. Brachet, D. I. Meiron, S. A. Orszag, B. G. Nickel, R. H. Morf, U. Frisch, Small-scale structure of the Taylor–Green vortex, Journal of Fluid Mechanics130(1983) 411–452

  47. [47]

    D. J. Lusher, N. D. Sandham, Assessment of low-dissipative shock-capturing schemes for the compressible Taylor–Green vortex, AIAA Journal59(2021) 533–545

  48. [48]

    S. B. Pope, Turbulent Flows, Cambridge University Press, Cambridge, UK, 2000

  49. [49]

    Mehta, F

    M. Mehta, F. Canabal, S. B. Tashakkor, S. D. Smith, Numerical base heating sensitivity study for a four-rocket engine core configuration, Journal of Spacecraft and Rockets50(2013) 509–526

  50. [50]

    S. K. Godunov, A difference method for numerical calculation of discontinuous solutions of the equations of hydrodynamics, Matematicheskii Sbornik47(1959) 271–306

  51. [51]

    Colella, P

    P. Colella, P. R. Woodward, The piecewise parabolic method (PPM) for gas-dynamical simulations, Journal of Computational Physics54(1984) 174–201

  52. [52]

    P. L. Roe, Approximate Riemann solvers, parameter vectors, and difference schemes, Journal of Computational Physics43(1981) 357–372. 27

  53. [53]

    Liou, A sequel to AUSM, Part II: AUSM+-up for all speeds, Journal of Computational Physics214(2006) 137–170

    M.-S. Liou, A sequel to AUSM, Part II: AUSM+-up for all speeds, Journal of Computational Physics214(2006) 137–170

  54. [54]

    J. P. Boris, D. L. Book, Flux-corrected transport. I. SHASTA, a fluid transport algorithm that works, Journal of Computational Physics11(1973) 38–69

  55. [55]

    S. T. Zalesak, Fully multidimensional flux-corrected transport algorithms for fluids, Journal of Computational Physics31(1979) 335–362

  56. [56]

    Clain, S

    S. Clain, S. Diot, R. Loub` ere, A high-order finite volume method for systems of conservation laws—Multi-dimensional Optimal Order Detection (MOOD), Journal of Computational Physics 230(2011) 4028–4050

  57. [57]

    Tadmor, The numerical viscosity of entropy stable schemes for systems of conservation laws

    E. Tadmor, The numerical viscosity of entropy stable schemes for systems of conservation laws. I, Mathematics of Computation49(1987) 91–103

  58. [58]

    Kurganov, E

    A. Kurganov, E. Tadmor, New high-resolution central schemes for nonlinear conservation laws and convection–diffusion equations, Journal of Computational Physics160(2000) 241–282

  59. [59]

    Cockburn, C.-W

    B. Cockburn, C.-W. Shu, The Runge–Kutta discontinuous Galerkin method for conservation laws V: Multidimensional systems, Journal of Computational Physics141(1998) 199–224

  60. [60]

    J. J. Quirk, A contribution to the great Riemann solver debate, International Journal for Numerical Methods in Fluids18(1994) 555–574

  61. [61]

    Glimm, J

    J. Glimm, J. W. Grove, X. Li, K.-M. Shyue, Y. Zeng, Q. Zhang, Three-dimensional front tracking, SIAM Journal on Scientific Computing19(1998) 703–727

  62. [62]

    S. K. Lele, Compact finite difference schemes with spectral-like resolution, Journal of computa- tional physics103(1992) 16–42. 28