Pith. sign in

REVIEW 2 minor 2 cited by

A Compression-Directional Entropic Stress Method for Shock-Regularized Compressible Flow

T0 review · 0 major / 2 minor · reviewed 2026-06-30 · grok-4.3

Pith's one-line read CoDeS supplies compression-selective shock regularization that preserves high-order resolution of contacts and vortical structures.

desk verdict CoDeS turns scalar entropic regularization into a compression-aligned tensor that stays off in shear and expansion, and the tests plus released code make the claim hold up. read the letter →

arxiv 2605.21444 v2 pith:LKNPP7G7 submitted 2026-05-20 physics.flu-dyn

classification physics.flu-dyn
keywords CoDeSmethodshockregularizationcompressibleflowsentropicstresstensorprincipalcompressiondirectionsmodifiedHelmholtzequationhigh-orderfinitevolumevorticalstructures
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper introduces the Compression-Directional Entropic Stress method to regularize shocks in compressible flow calculations. Instead of a uniform scalar addition, it builds a stress tensor aligned with the principal directions of compression, obtained by solving a modified Helmholtz equation and gated so it activates only under compression. This setup makes the regularization disappear in expansions, rigid rotations, and across contacts. A reader would care if it allows high-resolution schemes to keep sharp interfaces and energetic vortices even when shocks are present. The reported tests across one to three dimensions support that the method stays localized to compressive regions and preserves more energy in vortical flows than standard high-order references.

What carries the argument

The tensor stress Π_Σ = σ M aligned to the compressive eigenspace of the velocity gradient, with σ from modified-Helmholtz and gated by compression measures.

What would settle it

Running a simulation of an oblique shock wave or a curved shock front using the method and observing either instability or unexpected artifacts at the resolutions reported in the tests would falsify the central claim.

Watch

Extended reading notes

Core claim

CoDeS replaces scalar multidimensional entropic pressure with a tensor stress aligned with the principal directions of compression. The stress has the form Π_Σ = σ M, where σ comes from a modified-Helmholtz equation and M is built from the compressive eigenspace of the symmetric velocity-gradient tensor. The source is gated by volumetric and principal-strain compression so the regularization vanishes in smooth expansion, rigid-body rotation, and ideal contacts while recovering the compressive one-dimensional mechanism at planar shocks. The same tensor is used in momentum and energy fluxes.

Load-bearing premise

Gating the source term by volumetric and principal-strain compression combined with the modified-Helmholtz solve for sigma will produce stable and accurate results across all tested multi-dimensional configurations without introducing new artifacts at oblique shocks or curved fronts.

Editorial extensions

If this is right

  • The regularization vanishes in smooth expansion, rigid-body rotation, and ideal contacts.
  • It supplies localized stress at shocks while remaining weak in shear- and vorticity-dominated regions.
  • The three-dimensional Taylor--Green results at matched resolutions are comparable to or more energetic than seventh-order WENO/TENO references.
  • CoDeS provides a compression-selective shock regularization compatible with high-order finite-volume resolution of contacts, interfaces, shear layers, and vortical structures.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The selective nature of the regularization may enable more accurate long-term evolution of structures in flows that combine shocks with turbulence.
  • Applying similar directional gating to other numerical artifacts could improve fidelity in multi-physics simulations involving both discontinuities and smooth features.
  • Further tests on problems with strong shock curvature would help confirm the absence of artifacts beyond the cases already examined.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

0 major / 2 minor

Summary. The paper introduces the Compression-Directional Entropic Stress (CoDeS) method, which replaces scalar multidimensional entropic pressure with a tensor stress Π_Σ = σ M aligned to the compressive eigenspace of the velocity-gradient tensor. The scalar σ is obtained from a modified-Helmholtz solve whose source is gated by volumetric and principal-strain compression indicators, ensuring the regularization vanishes in expansion, rigid rotation, and contacts while recovering the 1D IGR mechanism at planar shocks. The same tensor is inserted into both momentum and energy fluxes. Numerical tests are reported on the Sod tube, double rarefaction, multidimensional Riemann problem, viscous shock tube, two-fluid triple point, Mach-3 slot jet, and 3D Taylor–Green vortex; results indicate that the stress remains localized to compressive fronts, produces no visible artifacts at oblique shocks or curved fronts, and yields solutions comparable to or more energetic than seventh-order WENO/TENO references at matched resolution. All code, case files, and plotting scripts are released on GitHub.

Significance. If the reported behavior holds under quantitative scrutiny, CoDeS supplies a compression-selective regularization that is compatible with high-order finite-volume discretizations of contacts, interfaces, shear layers, and vortical structures. The public release of the complete implementation, case settings, and figure-generation code is a clear strength that supports reproducibility and community verification. The approach could reduce the need for ad-hoc limiters or artificial viscosity in multi-dimensional compressible-flow simulations while preserving the underlying high-order scheme in smooth regions.

minor comments (2)
  1. The abstract states that 3D Taylor–Green results are 'comparable to or more energetic' than WENO/TENO references but supplies no L2 or L∞ error norms, convergence rates, or direct comparison tables; adding a short quantitative summary table in the abstract or §4 would strengthen the claim without lengthening the manuscript.
  2. Notation for the compressive projector M and the modified-Helmholtz operator is introduced in the abstract but not cross-referenced to the first appearance of the governing equations in the main text; a single forward reference would improve readability.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the detailed and positive summary of our work on the Compression-Directional Entropic Stress (CoDeS) method, as well as for recognizing its potential utility and the value of the public code release. The recommendation for minor revision is noted. No specific major comments were provided in the report.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity

full rationale

The derivation introduces an explicit tensor construction Π_Σ = σ M with σ from a modified-Helmholtz equation and M from the compressive eigenspace of the velocity-gradient tensor, gated by volumetric and principal-strain compression. No equation reduces the output stress to a fitted parameter, renamed input, or self-citation chain; the construction is stated directly from the symmetric velocity-gradient tensor and recovers the 1D IGR mechanism only at planar shocks by design. Numerical tests on Sod tube, multidimensional Riemann, Mach-3 jet, and 3D Taylor-Green supply independent verification that the stress localizes correctly without artifacts, confirming the chain is self-contained.

Assumptions & free parameters 0 free parameters · 0 assumptions · 0 invented entities

Abstract provides insufficient detail to enumerate free parameters or axioms; the modified-Helmholtz equation for sigma and the eigenspace construction are introduced without stated constants or background lemmas.

how reviews work

0 comments
Cite this review

Pith. "Pith review of A Compression-Directional Entropic Stress Method for Shock-Regularized Compressible Flow." pith.science (2026). https://pith.science/paper/LKNPP7G7

@misc{pith2026260521444,
  author       = {Pith},
  title        = {Pith review of: A Compression-Directional Entropic Stress Method for Shock-Regularized Compressible Flow},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/LKNPP7G7}},
  note         = {Machine review of arXiv:2605.21444}
}
abstract

We introduce the Compression-Directional Entropic Stress (CoDeS) method inspired by information geometric regularization. CoDeS replaces scalar multidimensional entropic pressure with a tensor stress aligned with the principal directions of compression. The stress has the form $\boldsymbol{\Pi}_{\Sigma}=\sigma\boldsymbol{M}$, where $\sigma$ is obtained from a modified-Helmholtz equation and $\boldsymbol{M}$ is constructed from the compressive eigenspace of the symmetric velocity-gradient tensor. The source is gated by volumetric and principal-strain compression, so the regularization vanishes in smooth expansion, rigid-body rotation, and ideal contacts, while recovering the compressive one-dimensional IGR mechanism at planar shocks. The same tensor stress is used in the conservative momentum flux and the stress-work energy flux. CoDeS is tested on one-, two-, and three-dimensional problems including smooth expansion, double rarefaction, the Sod shock tube, multidimensional Riemann flow, a viscous shock tube, a two-fluid triple point, a Mach-3 slot jet, and a supersonic Taylor--Green vortex. The results show that CoDeS remains inactive in expansive and contact regions, supplies localized stress at shocks, and concentrates regularization along compressive wave structures while remaining weak in shear- and vorticity-dominated regions. At matched resolutions, the three-dimensional Taylor--Green results are comparable to or more energetic than seventh-order WENO/TENO references. These results indicate that CoDeS provides a compression-selective shock regularization compatible with high-order finite-volume resolution of contacts, interfaces, shear layers, and vortical structures. All the code, case settings, and code for plotting figures of this paper are available at https://github.com/xubonan/code\_for\_CoDeS.

Figures

Figures reproduced from arXiv: 2605.21444 by the authors.

Figure 1
Figure 1. Smooth isentropic simple-wave expansion at 𝑡 = 0.1 on the 𝑛 = 320 grid. The CoDeS solution remains visually indistinguishable from the exact Euler solution in 𝜌, 𝑢, and 𝑝, while the CoDeS entropic stress remains zero to roundoff. In contrast, scalar IGR generates a finite stress inside the smooth expansion fan. Bonan Xu and Chihyung Wen: Preprint submitted to Elsevier Page 22 of 21 [PITH_FULL_IMAGE:figures/full_fig… view at source ↗
Figure 1
Figure 1. Smooth isentropic simple-wave expansion at 𝑡 = 0.1 on the 𝑛 = 320 grid. The CoDeS solution remains visually indistinguishable from the exact Euler solution in 𝜌, 𝑢, and 𝑝, while the CoDeS entropic stress remains zero to roundoff. In contrast, scalar IGR generates a finite stress inside the smooth expansion fan. Bonan Xu et al.: Preprint submitted to Elsevier Page 23 of 22 [PITH_FULL_IMAGE:figures/full_fig_p023_1.png] view at source ↗
Figure 2
Figure 2. 𝐿1 density-error convergence for the smooth isentropic simple-wave expansion. CoDeS recovers the high-order behavior of the underlying seventh-order finite-volume discretization, whereas scalar IGR converges at a substantially lower rate because it activates a finite stress in the smooth expansion region. Bonan Xu and Chihyung Wen: Preprint submitted to Elsevier Page 23 of 21 [PITH_FULL_IMAGE:figures/full_fig_p023_2.png] view at source ↗
Figures from the paper (27 more)
Figure 2
Figure 2. Figure 2: 𝐿1 density-error convergence for the smooth isentropic simple-wave expansion. CoDeS recovers the high-order behavior of the underlying seventh-order finite-volume discretization, whereas scalar IGR converges at a substantially lower rate because it activates a finite s…
Figure 3
Figure 3. Figure 3: Double-rarefaction problem computed on the 𝑛 = 200 grid with seventh-order reconstruction and 𝐶𝛼 = 2. The CoDeS and scalar-IGR solutions are compared with the exact Euler solution for 𝜌, 𝑢, and 𝑝, and the corresponding entropic stress is shown in panel (d). CoDeS prese…
Figure 3
Figure 3. Figure 3: Double-rarefaction problem computed on the 𝑛 = 200 grid with seventh-order reconstruction and 𝐶𝛼 = 2. The CoDeS and scalar-IGR solutions are compared with the exact Euler solution for 𝜌, 𝑢, and 𝑝, and the corresponding entropic stress is shown in panel (d). CoDeS prese…
Figure 4
Figure 4. Figure 4: Sod shock tube at 𝑡 = 0.2 on the 𝑛 = 400 grid with 𝐶𝛼 = 2, seventh-order linear upwind reconstruction, and SSP-RK3 time integration. The exact Euler solution is compared with CoDeS and scalar IGR for 𝜌, 𝑢, and 𝑝, and the corresponding entropic stress is shown in panel …
Figure 4
Figure 4. Figure 4: Sod shock tube case at 𝑡 = 0.2 on the 𝑛 = 400 grid with 𝐶𝛼 = 2, seventh-order linear upwind reconstruction, and SSP-RK3 time integration. The exact Euler solution is compared with CoDeS and scalar IGR for 𝜌, 𝑢, and 𝑝, and the corresponding entropic stress is shown in p…
Figure 5
Figure 5. Figure 5: Pressure-error behavior for the two-dimensional isentropic vortex. Panel (a) shows the 𝐿∞ pressure-error convergence after four periods, 𝑡 = 400. Panel (b) shows the corresponding pressure-error history on the 𝑁 = 200 mesh over the first four periods. CoDeS maintains t…
Figure 5
Figure 5. Figure 5: Pressure-error behavior for the two-dimensional isentropic vortex. Panel (a) shows the 𝐿∞ pressure-error convergence after four periods, 𝑡 = 400. Panel (b) shows the corresponding pressure-error history on the 𝑁 = 200 mesh over the first four periods. CoDeS maintains t…
Figure 6
Figure 6. Figure 6: Density fields for the perturbed two-dimensional Riemann problem at 𝑡 = 0.8 on a 500 × 500 grid. All three calculations reproduce the large-scale interacting Riemann structure. CoDeS preserves a sharper and more coherent central roll-up while maintaining clean shock tr…
Figure 6
Figure 6. Figure 6: Density fields for the perturbed two-dimensional Riemann problem at 𝑡 = 0.8 on a 500 × 500 grid. All three calculations reproduce the large-scale interacting Riemann structure. CoDeS preserves a sharper and more coherent central roll-up while maintaining clean shock tr…
Figure 7
Figure 7. Figure 7: Activation diagnostics for the perturbed two-dimensional Riemann problem at 𝑡 = 0.8. The CoDeS stress magnitude is concentrated primarily along shocks and strong compressive wave fronts, while remaining small in the central shear-layer roll-up and over most smooth regi…
Figure 7
Figure 7. Figure 7: Activation diagnostics for the perturbed two-dimensional Riemann problem at 𝑡 = 0.8. The CoDeS stress magnitude is concentrated primarily along shocks and strong compressive wave fronts, while remaining small in the central shear-layer roll-up and over most smooth regi…
Figure 8
Figure 8. Figure 8: Viscous shock-tube problem at 𝑡 = 1.0 on the 1280 × 640 grid. Shown are density 𝜌, pressure 𝑝, spanwise vorticity 𝜔𝑧 , and the CoDeS stress magnitude ‖𝚷Σ‖𝐹 . The density and pressure fields show the reflected shock system, oblique compression waves, and post-shock stru…
Figure 8
Figure 8. Figure 8: Viscous shock-tube problem at 𝑡 = 1.0 on the 1280 × 640 grid. Shown are density 𝜌, pressure 𝑝, spanwise vorticity 𝜔𝑧 , and the CoDeS stress magnitude ‖𝚷Σ‖𝐹 . The density and pressure fields show the reflected shock system, oblique compression waves, and post-shock stru…
Figure 9
Figure 9. Figure 9: Density evolution for the two-fluid triple-point problem at 𝑇 = 1, 2, 3, and 4. The sequence shows the deformation of the material interface, the development of baroclinically generated roll-up, and the propagation of the associated compressive wave system. 0 3.5 7 x 0…
Figure 9
Figure 9. Figure 9: Density evolution for the two-fluid triple-point problem at 𝑡 = 1, 2, 3, and 4. The sequence shows the deformation of the material interface, the development of baroclinically generated roll-up, and the propagation of the associated compressive wave system. 0 3.5 7 x 0…
Figure 10
Figure 10. Figure 10: Final-time diagnostics for the two-fluid triple-point problem at 𝑇 = 4. Shown are pressure 𝑝, volume fraction 𝛼𝓁 , the density-gradient indicator log(1 + |∇𝜌|), and the CoDeS stress magnitude ‖𝚷Σ‖𝐹 . The pressure and density-gradient fields identify the compressive wa…
Figure 10
Figure 10. Figure 10: Final-time diagnostics for the two-fluid triple-point problem at 𝑡 = 4. Shown are pressure 𝑝, volume fraction 𝛼𝓁 , the density-gradient indicator log(1 + |∇𝜌|), and the CoDeS stress magnitude ‖𝚷Σ‖𝐹 . The pressure and density-gradient fields identify the compressive wa…
Figure 11
Figure 11. Figure 11: Density evolution for the Mach–3 slot jet at 𝑇 = 2, 4, 6, and 8. The sequence shows the formation of the jet core, shock-cell structure, barrel-shock system, and shear-layer roll-up while the plume remains stable through the final time. Bonan Xu and Chihyung Wen: Prep…
Figure 11
Figure 11. Figure 11: Density evolution for the Mach–3 slot jet at 𝑡 = 2, 4, 6, and 8. The sequence shows the formation of the jet core, shock-cell structure, barrel-shock system, and shear-layer roll-up while the plume remains stable through the time duration. Bonan Xu et al.: Preprint su…
Figure 12
Figure 12. Figure 12: Density-gradient indicator for the Mach–3 slot jet at 𝑇 = 2, 4, 6, and 8. The Schlieren-type diagnostic highlights the shock-cell pattern, compression fronts, slot-lip shear layers, and downstream wave interactions. Bonan Xu and Chihyung Wen: Preprint submitted to Els…
Figure 12
Figure 12. Figure 12: Density-gradient indicator for the Mach–3 slot jet at 𝑡 = 2, 4, 6, and 8. The Schlieren-type diagnostic highlights the shock-cell pattern, compression fronts, slot-lip shear layers, and downstream wave interactions. Bonan Xu et al.: Preprint submitted to Elsevier Page…
Figure 13
Figure 13. Figure 13: CoDeS stress magnitude ‖𝚷Σ‖𝐹 for the Mach–3 slot jet at 𝑇 = 2, 4, 6, and 8. The stress is concentrated near compressive structures, including the inlet compression region, barrel-shock system, and downstream shock cells, while remaining weak over much of the ambient f…
Figure 13
Figure 13. Figure 13: CoDeS stress magnitude ‖𝚷Σ‖𝐹 for the Mach–3 slot jet at 𝑡 = 2, 4, 6, and 8. The stress is concentrated near compressive structures, including the inlet compression region, barrel-shock system, and downstream shock cells, while remaining weak over much of the ambient f…
Figure 14
Figure 14. Figure 14: Boundary-localized flux-fallback activation near the left slot-inlet boundary at 𝑇 = 2, 4, 6, and 8. The plotted quantity is the cell-centered diagnostic of the maximum adjacent fallback weight 𝜔𝑓 . Activation is confined to the prescribed boundary layer near the phys…
Figure 14
Figure 14. Figure 14: Boundary-localized flux-fallback activation near the left slot-inlet boundary at 𝑡 = 2, 4, 6, and 8. The plotted quantity is the cell-centered diagnostic of the maximum adjacent fallback weight 𝜔𝑓 . Activation is confined to the prescribed boundary layer near the phys…
Figure 15
Figure 15. Figure 15: Nondimensional solenoidal dissipation 𝜖𝑠 for the supersonic Taylor–Green vortex. CoDeS, WENO-7/LF, and TENO-7/LF are compared at the available grid resolutions. The CoDeS sequence shows increasing peak solenoidal dissipation under refinement, indicating improved resol…
Figure 15
Figure 15. Figure 15: Nondimensional solenoidal dissipation 𝜖𝑠 for the supersonic Taylor–Green vortex. CoDeS, WENO-7/LF, and TENO-7/LF are compared at the available grid resolutions. The CoDeS sequence shows increasing peak solenoidal dissipation under refinement, indicating improved resol…

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 2 Pith papers

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

  1. Information geometric regularization for computing sensitivities of flows with shocks

    math.NA 2026-08 conditional novelty 6.0 of 10

    The paper derives forward and adjoint sensitivity equations for the information geometrically regularized Euler equations and numerically verifies them against finite differences and automatic differentiation.

  2. Discontinuous Galerkin Semidiscretization of the Information Geometric Regularized Compressible Euler Equations

    math.NA 2026-08 conditional novelty 6.0 of 10

    A DG discretization of the IGR-regularized Euler equations stabilizes shocks and preserves fine-scale flow features in 1D/2D benchmarks without limiters or artificial viscosity.

Reference graph

Works this paper leans on

35 extracted references · 35 canonical work pages · cited by 2 Pith papers

  1. [1]

    Godunov, A difference scheme for numerical computation of discontinuous solutions of fluid dynamics, Mat

    S. Godunov, A difference scheme for numerical computation of discontinuous solutions of fluid dynamics, Mat. Sb 47 (1959) 271–306

  2. [3]

    P. L. Roe, Approximate riemann solvers, parameter vectors, and difference schemes, Journal of computational physics 43 (1981) 357–372

  3. [5]

    E. F. Toro, M. Spruce, W. Speares, Restoration of the contact surface in the hll-riemann solver, Shock waves 4 (1994) 25–34

  4. [6]

    Harten, B

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

  5. [7]

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

  6. [8]

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

  7. [9]

    Jiang, C.-W

    G.-S. Jiang, C.-W. Shu, Efficient implementation of weighted eno schemes, Journal of computational physics 126 (1996) 202–228

  8. [10]

    A.K.Henrick,T.D.Aslam,J.M.Powers, Mappedweightedessentiallynon-oscillatoryschemes:achievingoptimalordernearcriticalpoints, Journal of Computational Physics 207 (2005) 542–567

Show all 35 references
  1. [11]

    Borges, M

    R. Borges, M. Carmona, B. Costa, W. S. Don, An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws, Journal of computational physics 227 (2008) 3191–3211

  2. [12]

    Suresh, H

    A. Suresh, H. T. Huynh, Accurate monotonicity-preserving schemes with runge–kutta time stepping, Journal of Computational Physics 136 (1997) 83–99. Bonan Xu et al.:Preprint submitted to ElsevierPage 21 of 22 A Compression-Directional Entropic Stress Method for Shock-Regularize...

  3. [13]

    L. Fu, X. Y. Hu, N. A. Adams, A family of high-order targeted eno schemes for compressible-fluid simulations, Journal of Computational Physics 305 (2016) 333–359

  4. [14]

    L. Fu, X. Y. Hu, N. A. Adams, Targeted eno schemes with tailored resolution property for hyperbolic conservation laws, Journal of Computational Physics 349 (2017) 97–121

  5. [15]

    L. Fu, X. Y. Hu, N. A. Adams, A new class of adaptive high-order targeted eno schemes for hyperbolic conservation laws, Journal of Computational Physics 374 (2018) 724–751

  6. [17]

    F.Ducros,V.Ferrand,F.Nicoud,C.Weber,D.Darracq,C.Gacherieu,T.Poinsot, Large-eddysimulationoftheshock/turbulenceinteraction, Journal of Computational Physics 152 (1999) 517–549

  7. [18]

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

  8. [19]

    B.Fiorina,S.K.Lele, Anartificialnonlineardiffusivitymethodforsupersonicreactingflowswithshocks, JournalofComputationalPhysics 222 (2007) 246–264

  9. [20]

    Kawai, S

    S. Kawai, S. K. Lele, Localized artificial diffusivity scheme for discontinuity capturing on curvilinear meshes, Journal of Computational Physics 227 (2008) 9498–9526

  10. [21]

    Kawai, S

    S. Kawai, S. K. Shankar, S. K. Lele, Assessment of localized artificial diffusivity scheme for large-eddy simulation of compressible turbulent flows, Journal of computational physics 229 (2010) 1739–1762

  11. [22]

    Guermond, R

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

  12. [23]

    P.-O.Persson,J.Peraire, Sub-cellshockcapturingfordiscontinuousgalerkinmethods, in:44thAIAAaerospacesciencesmeetingandexhibit, 2006, p. 112

  13. [24]

    Johnsen, J

    E. Johnsen, J. Larsson, A. V. Bhagatwala, W. H. Cabot, P. Moin, B. J. Olson, P. S. Rawat, S. K. Shankar, B. Sjögreen, H. C. Yee, et al., Assessment of high-resolution methods for numerical simulations of compressible turbulence with shock waves, Journal of Computational Physic...

  14. [25]

    Pirozzoli, Numerical methods for high-speed flows, Annual review of fluid mechanics 43 (2011) 163–194

    S. Pirozzoli, Numerical methods for high-speed flows, Annual review of fluid mechanics 43 (2011) 163–194

  15. [26]

    R. Cao, F. Schäfer, Information geometric regularization of the barotropic euler equation, arXiv preprint arXiv:2308.14127 (2023)

  16. [27]

    R.Cao, F.Schäfer, Informationgeometricregularization ofunidimensional pressurelesseulerequations yieldsglobalstrong solutions, arXiv preprint arXiv:2411.15121 (2024)

  17. [28]

    Barham, B

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

  18. [29]

    Radhakrishnan, B

    A. Radhakrishnan, B. Wilfong, S. H. Bryngelson, F. Schäfer, Shocks without shock capturing: Information geometric regularization of finite volume methods for navier–stokes-like problems, arXiv preprint arXiv:2604.06546 (2026)

  19. [30]

    Taylor, R

    S. Taylor, R. J. Spiteri, S. Gaudreault, Thermodynamically constrained information geometric regularization for compressible flows, 2026. doi:10.48550/arXiv.2604.27286.arXiv:2604.27286

  20. [31]

    VonNeumann, R

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

  21. [32]

    E.J.Caramana,M.J.Shashkov,P.P.Whalen, Formulationsofartificialviscosityformulti-dimensionalshockwavecomputations, Journalof Computational Physics 144 (1998) 70–97

  22. [33]

    Harten, P

    A. Harten, P. D. Lax, B. v. Leer, On upstream differencing and godunov-type schemes for hyperbolic conservation laws, SIAM review 25 (1983) 35–61

  23. [34]

    Einfeldt, On godunov-type methods for gas dynamics, SIAM Journal on numerical analysis 25 (1988) 294–318

    B. Einfeldt, On godunov-type methods for gas dynamics, SIAM Journal on numerical analysis 25 (1988) 294–318

  24. [35]

    V. V. Rusanov, The calculation of the interaction of non-stationary shock waves and obstacles, USSR Computational Mathematics and Mathematical Physics 1 (1962) 304–320

  25. [36]

    P. D. Lax, Weak solutions of nonlinear hyperbolic equations and their numerical computation, Communications on pure and applied mathematics 7 (1954) 159–193

  26. [37]

    S. H. Bryngelson, K. Schmidmayer, V. Coralic, J. C. Meng, K. Maeda, T. Colonius, Mfc: An open-source high-order multi-component, multi-phase, and multi-scale compressible flow solver, Computer Physics Communications 266 (2021) 107396

  27. [38]

    Gottlieb, C.-W

    S. Gottlieb, C.-W. Shu, Total variation diminishing runge-kutta schemes, Mathematics of computation 67 (1998) 73–85. Bonan Xu et al.:Preprint submitted to ElsevierPage 22 of 22 A Compression-Directional Entropic Stress Method for Shock-Regularized Compressible Flow 0.0 0.5 1.0...

Pith tools

Reviewed June 30, 2026 · model on record in the stance chip above.