pith. sign in

arxiv: 2605.21444 · v2 · pith:LKNPP7G7new · submitted 2026-05-20 · ⚛️ physics.flu-dyn

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

Bonan Xu , Chihyung Wen , Peixu Guo This is my paper

Pith reviewed 2026-05-21 02:51 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords shock regularizationcompressible flowentropic stressfinite volume methodprincipal compression directionsTaylor-Green vortexmultidimensional shocks
0
0 comments X

The pith

A directional tensor stress regularizes shocks selectively while vanishing in expansions, contacts and shear.

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

The paper introduces the Compression-Directional Entropic Stress method (CoDeS) for finite-volume treatment of shock-dominated compressible flows. It replaces a scalar entropic pressure with a tensor stress aligned to the principal compression directions extracted from the symmetric velocity-gradient tensor. The stress is sourced from a modified Helmholtz equation and gated so that it activates only under compression and recovers the one-dimensional mechanism at planar shocks. Tests on one-, two-, and three-dimensional problems show the regularization concentrates at shocks, stays inactive in expansions and ideal contacts, and remains weak in shear and vorticity regions. At matched grid resolutions the three-dimensional Taylor-Green vortex yields energy levels comparable to or higher than those obtained with seventh-order WENO and TENO schemes.

Core claim

CoDeS constructs the stress tensor Π_Σ = σ M, where M is the compressive eigenspace matrix of the symmetric velocity-gradient tensor, σ is obtained from a modified-Helmholtz equation, and the entire term is gated by volumetric and principal-strain compression; the same tensor enters both the momentum and energy fluxes, recovering the one-dimensional IGR mechanism at shocks while remaining zero in smooth expansion, rigid-body rotation, and ideal contacts.

What carries the argument

The compressive eigenspace matrix M derived from the symmetric velocity-gradient tensor, which orients the entropic stress exactly along the principal axes of compression.

If this is right

  • CoDeS supplies localized stress only at compressive shocks.
  • The regularization concentrates along compressive wave structures while staying 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.
  • CoDeS remains compatible with high-order finite-volume resolution of contacts, interfaces, and vortical structures.

Where Pith is reading between the lines

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

  • The same directional gating idea could be applied to other stabilization techniques that currently rely on isotropic sensors.
  • Coarser grids might suffice for shock-capturing problems if the method preserves accuracy in smooth rotational and shear regions.
  • Further tests on strong multidimensional shock interactions would directly check whether the gating prevents the new instabilities the authors assume will not appear.

Load-bearing premise

The compressive eigenspace matrix M together with volumetric and principal-strain gating will not create new artifacts or instabilities in genuinely multidimensional shock interactions beyond those already present in the one-dimensional limit.

What would settle it

A simulation of a multidimensional shock interaction (such as the two-fluid triple point or Mach-3 slot jet) that exhibits oscillations, instabilities, or excess dissipation absent from the corresponding one-dimensional planar-shock test.

Figures

Figures reproduced from arXiv: 2605.21444 by Bonan Xu, Chihyung Wen, Peixu Guo.

Figure 1
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:figure… view at source ↗
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 stress in the smooth expansion region. Bonan Xu and Chihyung Wen: Preprint submitted to Elsevier Page 23 of 21 [PITH_FULL_IMAGE:figures/full… view at source ↗
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 preserves the rarefaction structure without generating appreciable stress in the expansion regions, whereas scalar IGR produces nonzero stress ne… view at source ↗
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 (d). Both CoDeS and scalar IGR generate comparable stress at the right-going shock, confirming that the CoDeS compression gate does not supp… view at source ↗
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 the lowest error over the tested resolutions and throughout the long-time advection, whereas scalar IGR exhibits substantially larger error g… view at source ↗
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 transitions, whereas scalar IGR visibly distorts parts of the vortical interaction region and WENO-5/LF gives a more diffuse roll-up with osci… view at source ↗
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 regions. The scalar IGR entropic pressure activates more broadly, including in vortical and shear-dominated regions, and exhibits sign-indefinit… view at source ↗
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 structures generated by shock-wall and shock–boundary-layer interaction. The vorticity field highlights the wall-generated shear layers and roll… view at source ↗
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 1.5 3 y 0.15 0.30 0.45 0.60 p (a) Pressure 𝑝. 0 3.5 7 x 0 1.5 3 y 0.0 0.2 0.4 0.6 0.8 1.0 αl (b) Volume fraction 𝛼𝓁 . 0 3.5 7 x 0 1.5 3 y 0… view at source ↗
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 waves and material interfaces, while the stress diagnostic shows that CoDeS is concentrated primarily along compressive wave structures and re… view at source ↗
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: Preprint submitted to Elsevier Page 29 of 21 [PITH_FULL_IMAGE:figures/full_fig_p029_11.png] view at source ↗
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 Elsevier Page 30 of 21 [PITH_FULL_IMAGE:figures/full_fig_p030_12.png] view at source ↗
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 flow and vortical shear-layer roll-up. Bonan Xu and Chihyung Wen: Preprint submitted to Elsevier Page 31 of 21 [PITH_FULL_IMAGE:figures/full… view at source ↗
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 physical inlet and does not replace the high-order flux in the interior jet, shock cells, or shear-layer roll-up. Bonan Xu and Chihyung Wen: Pre… view at source ↗
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 resolution of the vortical cascade. Bonan Xu and Chihyung Wen: Preprint submitted to Elsevier Page 33 of 21 [PITH_FULL_IMAGE:figures/full_fig_p0… view at source ↗
read the original 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.

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

3 major / 3 minor

Summary. The paper introduces the Compression-Directional Entropic Stress (CoDeS) method for shock regularization in compressible flows. It replaces scalar entropic pressure with a tensor stress Π_Σ = σ M, where M is constructed from the compressive eigenspace of the symmetric velocity-gradient tensor and σ is obtained from a modified-Helmholtz equation. The source is gated by volumetric and principal-strain compression indicators. The method is claimed to recover the 1-D IGR mechanism at planar shocks, remain inactive in expansive, contact, and shear regions, and provide localized regularization along compressive structures. Numerical tests on 1D Sod shock tube, 2D Riemann problems, viscous shock tube, triple point, Mach-3 jet, and 3D supersonic Taylor-Green vortex demonstrate the expected selective behavior and competitive performance against high-order WENO/TENO schemes.

Significance. If the selectivity claims hold under rigorous validation, CoDeS could represent a meaningful advance in shock-capturing techniques for high-resolution simulations of compressible flows involving shocks, interfaces, and turbulence. By aligning regularization with principal compression directions and gating it appropriately, the method aims to reduce numerical dissipation in non-compressive regions compared to traditional scalar approaches. The 3D Taylor-Green results suggesting maintained or higher energy levels are potentially important for under-resolved turbulent flows, but require quantitative substantiation to confirm significance.

major comments (3)
  1. [§3.2, definition of M] The construction of the compressive eigenspace matrix M from the symmetric velocity-gradient tensor, combined with volumetric and principal-strain gating, is asserted to keep the stress strictly compressive and inactive outside shocks. However, no analytic demonstration is provided that this holds for oblique, curved, or interacting shocks (e.g., in the multidimensional Riemann problem), which is load-bearing for the claim that the method remains weak in shear- and vorticity-dominated regions without introducing artifacts.
  2. [§5, Taylor-Green vortex test] The claim that at matched resolutions the three-dimensional Taylor-Green results are comparable to or more energetic than seventh-order WENO/TENO references lacks supporting quantitative data such as kinetic energy decay curves, L2 error norms, or direct residual comparisons. This is critical for evaluating the method's performance and the assertion of reduced interference with vortical structures.
  3. [Abstract and §5] The central claims of expected behavior and improved selectivity are supported only by qualitative descriptions of the test suite results. The manuscript supplies no quantitative error tables, convergence rates, or side-by-side comparisons against established regularization methods, which undermines the strength of the evidence for the method's advantages.
minor comments (3)
  1. [§2] Additional references to related work on tensor-based artificial viscosity or directional shock sensors would strengthen the background section.
  2. [Figure captions] The figure captions for the 3D results should include more details on the resolution and reference scheme parameters for easier comparison.
  3. [Notation] The modified-Helmholtz equation for σ should be presented with an explicit equation number to improve clarity.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed review of our manuscript on the CoDeS method. We address each major comment point by point below, providing our responses and indicating planned revisions to the manuscript.

read point-by-point responses

  1. Referee: [§3.2, definition of M] The construction of the compressive eigenspace matrix M from the symmetric velocity-gradient tensor, combined with volumetric and principal-strain gating, is asserted to keep the stress strictly compressive and inactive outside shocks. However, no analytic demonstration is provided that this holds for oblique, curved, or interacting shocks (e.g., in the multidimensional Riemann problem), which is load-bearing for the claim that the method remains weak in shear- and vorticity-dominated regions without introducing artifacts.

    Authors: We agree that a complete analytic demonstration for arbitrary oblique, curved, and interacting shocks is not supplied in the current manuscript. By construction, M is formed exclusively from the compressive eigenvectors of the symmetric velocity-gradient tensor, and the source term is multiplied by indicators that are identically zero for pure shear, rotation, and expansion. This ensures the stress vanishes outside compressive regions by design. The multidimensional Riemann problem test in §5 numerically confirms localization to shock structures without visible artifacts in adjacent shear layers. We will add a clarifying paragraph in §3.2 describing this design rationale and its implications for multidimensional flows. revision: partial

  2. Referee: [§5, Taylor-Green vortex test] The claim that at matched resolutions the three-dimensional Taylor-Green results are comparable to or more energetic than seventh-order WENO/TENO references lacks supporting quantitative data such as kinetic energy decay curves, L2 error norms, or direct residual comparisons. This is critical for evaluating the method's performance and the assertion of reduced interference with vortical structures.

    Authors: The referee correctly identifies the absence of quantitative support for the energy comparison. In the revised manuscript we will include time histories of total kinetic energy for the 3D supersonic Taylor-Green vortex at the reported resolutions, plotted against the seventh-order WENO and TENO reference solutions. We will also report L2 velocity norms at selected times to quantify the differences and substantiate the claim of maintained or higher energy levels due to reduced dissipation outside compressive regions. revision: yes

  3. Referee: [Abstract and §5] The central claims of expected behavior and improved selectivity are supported only by qualitative descriptions of the test suite results. The manuscript supplies no quantitative error tables, convergence rates, or side-by-side comparisons against established regularization methods, which undermines the strength of the evidence for the method's advantages.

    Authors: We acknowledge that the current presentation emphasizes qualitative field visualizations to demonstrate selective activation. For the Sod shock tube the numerical solution matches the exact Riemann solution to visual accuracy, providing implicit quantitative validation. We will revise §5 to include a summary table of L1 and L2 error norms for the 1D Sod and selected 2D Riemann cases, together with a short discussion of observed accuracy in smooth regions. Direct comparisons against other artificial-viscosity or entropic regularizations are outside the present scope and will be noted as future work. The abstract will be updated to reflect these additions and to moderate the language on selectivity advantages. revision: yes

Circularity Check

0 steps flagged

No significant circularity; construction is explicit modeling choice, not reduction to inputs

full rationale

The paper explicitly constructs the tensor M from the eigendecomposition of the symmetric velocity-gradient tensor and gates the source term using local volumetric and principal-strain compression indicators. These choices ensure by definition that regularization vanishes outside compressive regions and recovers the 1D IGR limit at normal shocks. However, this is an independent ansatz for the method rather than a derivation that reduces tautologically to fitted data or prior results within the paper. Numerical experiments on multiple test cases then demonstrate the intended behavior, but the core selectivity follows directly from the stated definitions without statistical forcing or self-referential loops. No load-bearing self-citations, uniqueness theorems, or fitted parameters renamed as predictions appear in the derivation chain.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 1 invented entities

The central construction rests on the kinematic definition of compression directions and the assumption that a Helmholtz solve can supply a suitable scalar multiplier without introducing new length scales that conflict with the underlying finite-volume scheme.

free parameters (1)
  • Helmholtz diffusion or relaxation coefficient
    The modified-Helmholtz equation for σ requires at least one coefficient whose value is not fixed by the abstract and may be chosen or calibrated.
axioms (1)
  • domain assumption The principal directions of the symmetric part of the velocity gradient correctly identify the local compression axes for regularization purposes.
    Invoked when constructing the matrix M from the compressive eigenspace.
invented entities (1)
  • Compression-directional entropic stress tensor Π_Σ no independent evidence
    purpose: To provide anisotropic regularization that vanishes outside compressive regions.
    New tensor field introduced by the paper; no independent experimental signature is claimed.

pith-pipeline@v0.9.0 · 5820 in / 1454 out tokens · 39722 ms · 2026-05-21T02:51:40.403187+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.