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 →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
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.
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
- 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.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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)
- 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.
- 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
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
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
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 from the paper (27 more)
Forward citations
Cited by 2 Pith papers
-
Information geometric regularization for computing sensitivities of flows with shocks
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.
-
Discontinuous Galerkin Semidiscretization of the Information Geometric Regularized Compressible Euler Equations
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
-
[1]
S. Godunov, A difference scheme for numerical computation of discontinuous solutions of fluid dynamics, Mat. Sb 47 (1959) 271–306
work page 1959
-
[3]
P. L. Roe, Approximate riemann solvers, parameter vectors, and difference schemes, Journal of computational physics 43 (1981) 357–372
work page 1981
-
[5]
E. F. Toro, M. Spruce, W. Speares, Restoration of the contact surface in the hll-riemann solver, Shock waves 4 (1994) 25–34
work page 1994
- [6]
-
[7]
C.-W. Shu, S. Osher, Efficient implementation of essentially non-oscillatory shock-capturing schemes, Journal of computational physics 77 (1988) 439–471
work page 1988
-
[8]
X.-D. Liu, S. Osher, T. Chan, Weighted essentially non-oscillatory schemes, Journal of computational physics 115 (1994) 200–212
work page 1994
-
[9]
G.-S. Jiang, C.-W. Shu, Efficient implementation of weighted eno schemes, Journal of computational physics 126 (1996) 202–228
work page 1996
-
[10]
A.K.Henrick,T.D.Aslam,J.M.Powers, Mappedweightedessentiallynon-oscillatoryschemes:achievingoptimalordernearcriticalpoints, Journal of Computational Physics 207 (2005) 542–567
work page 2005
Show all 35 references
-
[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
2008
-
[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...
1997
-
[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
2016
-
[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
2017
-
[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
2018
-
[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
1999
-
[18]
A. W. Cook, W. H. Cabot, Hyperviscosity for shock-turbulence interactions, Journal of Computational Physics 203 (2005) 379–385
2005
-
[19]
B.Fiorina,S.K.Lele, Anartificialnonlineardiffusivitymethodforsupersonicreactingflowswithshocks, JournalofComputationalPhysics 222 (2007) 246–264
2007
-
[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
2008
-
[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
2010
-
[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
2011
-
[23]
P.-O.Persson,J.Peraire, Sub-cellshockcapturingfordiscontinuousgalerkinmethods, in:44thAIAAaerospacesciencesmeetingandexhibit, 2006, p. 112
2006
-
[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...
2010
-
[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
2011
-
[26]
R. Cao, F. Schäfer, Information geometric regularization of the barotropic euler equation, arXiv preprint arXiv:2308.14127 (2023)
2023
-
[27]
R.Cao, F.Schäfer, Informationgeometricregularization ofunidimensional pressurelesseulerequations yieldsglobalstrong solutions, arXiv preprint arXiv:2411.15121 (2024)
2024 arXiv
-
[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)
2025
-
[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)
2026 arXiv
- [30]
-
[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
1950
-
[32]
E.J.Caramana,M.J.Shashkov,P.P.Whalen, Formulationsofartificialviscosityformulti-dimensionalshockwavecomputations, Journalof Computational Physics 144 (1998) 70–97
1998
-
[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
1983
-
[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
1988
-
[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
1962
-
[36]
P. D. Lax, Weak solutions of nonlinear hyperbolic equations and their numerical computation, Communications on pure and applied mathematics 7 (1954) 159–193
1954
-
[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
2021
-
[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...
1998
Reviewed June 30, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.