{"id":"2e1a7223-5bd2-498c-be6b-35ec9e8f5ba6","arxiv_id":"2605.21444","paper_version":2,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"CoDeS supplies a compression-aligned tensor stress from a modified-Helmholtz solve that activates only at shocks and remains inactive in expansions, contacts, and shear.","lead":"The paper introduces the CoDeS method, which applies a directional tensor stress derived from compressive flow directions to regularize shocks in compressible fluid simulations. A smart generalist might read it because the approach aims to preserve accuracy at contacts, shear layers, and vortices while handling shocks, potentially improving high-resolution CFD for engineering applications.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.3","headline":"No significant objection identified","rationale":"Reader's weakest assumption was flagged from the abstract alone. The complete text and listed test suite address exactly that assumption with concrete multi-dimensional evidence, including oblique-shock configurations. No internal inconsistency, hidden assumption in the tensor construction, or unsupported extrapolation is visible in the reported results or code availability.","tokens_in":1815,"tokens_out":258,"duration_ms":27232,"concrete_test":"Clone https://github.com/xubonan/code_for_CoDeS, rerun the Mach-3 slot jet at the paper's resolution, and compare the density and pressure fields along a line crossing the oblique shock against a reference high-order scheme; check for any new oscillations exceeding 1% of the jump.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The full manuscript supplies numerical results for the multidimensional Riemann problem, Mach-3 slot jet, and 3D Taylor-Green vortex that directly test the gating and modified-Helmholtz construction at oblique shocks and curved compressive fronts. These cases show the tensor stress remains localized to compressive principal directions, vanishes in shear and expansion, and produces results comparable to or better than high-order WENO/TENO references without reported non-physical oscillations or interface degradation.","agreement_with_reader":"disagree"},"referee_report":{"model":"grok-4.3","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.","tokens_in":1889,"tokens_out":476,"duration_ms":26538,"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.","major_comments":[],"minor_comments":[{"comment":"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.","section":null},{"comment":"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.","section":null}],"recommendation":"minor_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"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.","responses":[],"tokens_in":1458,"tokens_out":82,"duration_ms":21098,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The paper's core move is replacing the scalar IGR pressure with a tensor stress Pi_Sigma = sigma M, where M comes from the compressive eigenvectors of the velocity gradient and sigma solves a modified Helmholtz equation. The source is gated by volumetric and principal-strain compression so the term drops out in expansion, rotation, and contacts. That construction is new relative to the scalar versions cited.\n\nThe tests cover the expected ground: 1D tubes, multidimensional Riemann problems, oblique shocks in the slot jet, and the 3D Taylor-Green vortex. The results show the stress stays localized to compressive fronts, produces no visible oscillations at contacts or shear layers, and at matched resolution the 3D vortex is at least as energetic as the seventh-order WENO/TENO runs. Code and plotting scripts are on GitHub, which lets anyone check the implementation directly.\n\nThe main limitation is that the quantitative support is mostly side-by-side contour plots and statements of comparability rather than tabulated L1 or L2 errors against exact solutions or systematic grid-convergence studies in the multi-dimensional cases. That is not fatal for a methods paper, but it leaves the accuracy claim resting on visual inspection and the reference-scheme comparison.\n\nThis is aimed at CFD groups that already run high-order finite-volume codes and want a shock regularizer that does not degrade interfaces or vorticity. A reader working on alternative stabilization techniques would find the gating logic and the tensor form worth looking at. The work is clearly defined, the numerics are reproducible, and the central claim survives the multi-D tests that were run. I would send it to peer review.","headline":"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.","tokens_in":2371,"tokens_out":404,"would_cite":false,"duration_ms":22082,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"CoDeS supplies compression-selective shock regularization that preserves high-order resolution of contacts and vortical structures.","keywords":["CoDeS method","shock regularization","compressible flows","entropic stress tensor","principal compression directions","modified Helmholtz equation","high-order finite volume","vortical structures"],"falsifier":"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.","tokens_in":2717,"feed_emoji":"🌊","tokens_out":727,"duration_ms":44165,"temperature":0.7,"pith_summary":"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.","feed_headline":"Directional tensor regularizes shocks without harming contacts or vortices","feed_subtitle":"By aligning stress to compressive principal directions and gating by strain, CoDeS stays off in expansions and shear while activating at sho","key_machinery":"The tensor stress Π_Σ = σ M aligned to the compressive eigenspace of the velocity gradient, with σ from modified-Helmholtz and gated by compression measures.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"fun_headline_variants":["Tensor stress follows principal compression directions","CoDeS regularizes shocks with gated compression tensor","Entropic tensor stress targets compressive shocks only","Compression-directional stress for selective shock regularization"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Tensor stress follows principal compression directions","CoDeS regularizes shocks with gated compression tensor","Entropic tensor stress targets compressive shocks only","Compression-directional stress for selective shock regularization"]},"model":"grok-4.3","cost_usd":0.007307,"raw_usage":{"total_tokens":3429,"prompt_tokens":797,"num_sources_used":0,"completion_tokens":52,"cost_in_usd_ticks":73074500,"prompt_tokens_details":{"text_tokens":797,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":2580,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":797,"tokens_out":52,"duration_ms":28919,"temperature":1.0,"reasoning_tokens":2580,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-30T17:01:21.203729+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"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.","supporting_citations":[],"review_version":2}