{"id":"faf9148d-2ecd-4cb1-9ed7-3f3be8dff1eb","arxiv_id":"2506.13231","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":0,"one_line_summary":"A new Entropy-Stable/Double-Flux numerical flux for multi-component compressible Navier-Stokes equations is claimed to preserve kinetic energy, suppress interface oscillations, and satisfy the second law of thermodynamics.","lead":"This paper combines entropy-stable flux construction with the Double-Flux method to simulate multi-component compressible flows with variable specific heat ratios. It claims an oscillation-free, thermodynamically consistent scheme, backed by a semi-discrete entropy inequality proof and OpenFOAM benchmark cases.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 4.5's entropy-stability proof rests on Lemma 4.1, which assumes [[r(Y)]]=0 across species interfaces, and the final dissipation claim also relies on a false Lemma 4.2.","rationale":"The reader's weakest-assumption analysis correctly identifies Lemma 4.1 as a load-bearing failure: at the species interfaces used in the paper, r(Y) is discontinuous, so [[r]]=0 is false and the proof of Eq. (43) is invalid. My independent reading confirms that this is not a cosmetic gap; the proof of Theorem 4.5 passes through this identity at the point where the entropy flux jump is expanded. I additionally find that Lemma 4.2, which supplies the final dissipation inequality, is itself false for unequal species gas constants, giving a two-species counterexample where the asserted nonpositive quantity is positive. This means the central entropy-stability claim is not only unproved but is contradicted by direct evaluation of the final expression in a simple interface configuration. The numerical sections do not test the entropy inequality directly, so they provide no countervailing evidence. The method may still behave well in practice, but the paper's principal theoretical assertion is not supported. Since the reader already recommended REJECT and my analysis reinforces that conclusion, the verdict should remain unchanged.","tokens_in":19718,"tokens_out":3909,"duration_ms":42174,"concrete_test":"Re-derive Eq. (43)–(46) without Lemma 4.1, replacing [[rρ]] with r̄[[ρ]]+ρ̄[[r]], and track the resulting extra terms into the final entropy-difference expression. Then evaluate that expression for the two-species moving-interface configuration with H2/N2 (Y_H2 jumping between 0.1 and 0.9, r_H2=4124, r_N2=296.8, uniform ρ, v=100 m/s, p=1 atm). If the final expression is positive, the claimed inequality (46) is violated and the flux is not entropy-stable as stated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that flux (42) satisfies a semi-discrete entropy inequality—is unsupported because Theorem 4.5 uses two false ingredients. First, Eq. (43) invokes Lemma 4.1, [[r(Y)ρ]] = r[[ρ]], which is derived from [[r]]=0. For the paper's own motivating cases (H2/N2, He/air), r(Y)=ΣY_i r_i is discontinuous at the interface, so [[r]]≠0 and the jump identity must read [[rρ]] = r̄[[ρ]] + ρ̄[[r]]. The missing ρ̄[[r]] term propagates through the proof and changes the final balance. Second, even granting Lemma 4.1, the proof ends by invoking Lemma 4.2, whose inequality Σ r_i Ȳ_i [[ln Y_i]] ρ^ln ≤ 0 is false when the r_i differ. Counterexample with two species: r1=2, r2=1, Y1:0.1→0.9, Y2:1−Y1. The sum is 2·0.5·ln9 + 1·0.5·ln(1/9) = 0.5 ln9 ≈ 1.099 > 0. Thus Eq. (46) does not establish entropy stability; in fact the expression can be positive, meaning the flux is not provably entropy-stable for variable-r mixtures. The numerical benchmarks, which are qualitative and do not monitor the discrete entropy budget, cannot repair this. The theorem as stated is therefore not merely missing a minor assumption; its proof collapses at the load-bearing step.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a numerical flux for multi-component compressible Navier-Stokes equations that combines an entropy-stable formulation with the double-flux model. The central claim, stated in the abstract and in Section 4.2, is that the proposed flux satisfies a semi-discrete entropy inequality and is entropy-stable, while preserving kinetic energy and suppressing pressure oscillations at material interfaces. The method is implemented in OpenFOAM with adaptive mesh refinement and is tested on a moving interface, a two-dimensional shock-bubble interaction, and a three-dimensional under-expanded hydrogen jet. The paper also introduces a hybrid dissipation strategy intended to blend low and high dissipation mechanisms while maintaining entropy stability.","tokens_in":20042,"tokens_out":4875,"duration_ms":49931,"significance":"If the theoretical claim were valid, the scheme would be a valuable contribution to multi-component compressible flow simulation: it promises thermodynamic consistency, kinetic energy preservation, and oscillation-free behavior at species interfaces, which are typically conflicting goals. The numerical benchmark suite is reasonably broad, includes quantitative comparisons (Mach disk dimensions and jet penetration in Section 6.4), and reports an OpenFOAM implementation that could be useful to practitioners. However, the central proof of entropy stability rests on two lemmas that are false for the intended applications, so the paper's main theoretical conclusion is unsupported. The numerical experiments do not monitor a discrete entropy budget and therefore cannot compensate for the faulty proof.","major_comments":[{"comment":"Lemma 4.1 assumes [[r(Y)]] = 0, but r(Y) = Σ_i Y_i r_i is discontinuous at species interfaces, which are exactly the cases targeted by the paper (H2/N2 in Section 6.2, He/air in Section 6.3). The correct jump identity is [[rρ]] = \\bar r [[ρ]] + \\bar ρ [[r]]. Since Eq. (43) of Theorem 4.5 uses Lemma 4.1 directly, the proof of the semi-discrete entropy inequality is invalid for the intended applications.","section":"Section 4.1, Lemma 4.1 (Eq. 34)"},{"comment":"The claimed inequality Σ_i r_i \\bar Y_i [[ln Y_i]] ρ^ln ≤ 0 is false when the species gas constants differ. For two species with r1=2, r2=1 and Y1: 0.1→0.9, Y2=1−Y1, the sum equals 0.5 ln 9 ≈ 1.099, which is positive. The proof in Eq. (37) passes from Σ_i r_i Y_i [[ln Y_i]] to Σ_i [[ln Y_i]] r ρ^ln without a valid inequality, and the subsequent bound is therefore a non-sequitur. Since Eq. (46) in Theorem 4.5 uses this lemma as the final dissipation step, the entropy-stability conclusion is unsupported.","section":"Section 4.1, Lemma 4.2 (Eq. 35)"},{"comment":"Because Theorem 4.5 relies on Lemmas 4.1 and 4.2, both of which fail for variable-r mixtures, the proof that the flux (42) is entropy-stable collapses. The numerical experiments do not monitor the discrete entropy residual, so they cannot remedy this. The claim in the abstract and Section 4.2 that the flux 'satisfies a semi-discrete entropy inequality' is therefore not established.","section":"Section 4.2, Theorem 4.5"},{"comment":"The statement that the dissipation term −½ (vR−vL)^T R |Λ| T R^T (vR−vL) is a negative quadratic form requires proving that R |Λ| T R^T is symmetric positive semi-definite. The manuscript does not provide this proof for the modified scaling based on γ*, and the derivation in Eq. (53) introduces R^{-1} in a way that is not the standard entropy-variable dissipation structure, so the entropy-stability claim for the hybrid flux (54) is not justified.","section":"Section 4.3, Eq. (55)"}],"minor_comments":[{"comment":"The word 'challange' in the second paragraph should be 'challenge'.","section":"Introduction"},{"comment":"The definition of c*_{pi} is unclear: the integration variable appears as T in both the integrand and the denominator, and the reference state T0 is not specified. Please clarify the notation.","section":"Section 3.3, Eq. (32)"},{"comment":"The notation ρ^ln and Y_i^ln in Eq. (37) is not defined, and the inequality steps are not explained; this makes the already problematic proof even harder to follow.","section":"Section 4.1, proof of Lemma 4.2"},{"comment":"The entropy inequality is written for the inviscid entropy flux only, while the Navier-Stokes equations include viscous terms that contribute to entropy production; the statement should clarify the role of the viscous fluxes.","section":"Section 2.1, Eq. (11)"},{"comment":"The reported convergence order Δt^{2.91} is close but not equal to the formal third order of SSP-RK3; the text should explain how the measured order is obtained and why it differs from 3.","section":"Section 6.1, Figure 3"},{"comment":"The boundary condition at the inlet is described only as Dirichlet, but the velocity profile and turbulence specification are not given, which limits reproducibility of the three-dimensional jet simulation.","section":"Section 6.4"}],"recommendation":"reject","confidential_remarks":"The paper contains a substantial implementation effort and an interesting scheme, but the core proof of entropy stability is invalid because the key lemmas are false for the very cases the method is designed to handle. The numerical results, while impressive in scope, cannot serve as evidence for the entropy-stability claim without a discrete entropy diagnostic. I would consider a resubmission if the authors can correct the proof or reformulate the flux so that the entropy inequality genuinely holds for variable-r mixtures."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the punchline: the Entropy-Stable/Double-Flux combination is genuinely new, and the numerical campaign is serious; but the paper's central claim — the rigorous proof that flux (42) satisfies a semi-discrete entropy inequality — is wrong as written. I checked the proof chain, and it breaks at two load-bearing steps.\n\nWhat's new and good: wrapping the Double-Flux frozen γ* and e0* variables around an entropy-stable base flux, plus the pressure-switched hybrid LF/Roe dissipation, is a combination I don't find in the cited literature. The mixture-variable flux in Eq. 42 requiring only two logarithmic-mean evaluations instead of n+1 is a practical, credible efficiency gain. The benchmarks are extensive and mostly well executed: moving H2/N2 interface, shock-helium-bubble interaction with AMR against Terashima and Tryggvason's data, and a 3D under-expanded hydrogen jet with Mach-disk metrics close to Hamzehloo and Ruggles/Ekoto. The qualitative behavior — no pressure/velocity oscillations at the material interface, crisper oblique shocks than Double-Flux alone — matches what the scheme is supposed to deliver.\n\nThe soft spot is the theory, and it is not minor. Lemma 4.1 asserts [[r(Y)ρ]] = r̄[[ρ]] using [[r]] = 0, which is false at species interfaces: in the paper's own H2/N2 and He/air tests, r(Y) jumps, so the missing ρ̄[[r]] term rides through the proof of Theorem 4.5. Lemma 4.2, which produces the non-positive right-hand side, is also false when the r_i differ: with r1=2, r2=1 and Y1 going 0.1→0.9, the claimed sum is 0.5 ln 9 > 0. The stress-test counterexample holds up, so the theorem fails exactly where flagged. At best the flux is entropy-conservative in the equal-r limit, consistent with their n=1 test. The single-component entropy-convergence test in §6.1 cannot exercise these terms, and the multi-component cases never monitor the discrete entropy budget, so the numerics don't rescue the theory. Also, the paper concedes in §3.3 that Double-Flux is inherently non-conservative but never quantifies the conservation error in any benchmark; data and code are only 'available upon reasonable request.'\n\nWho this is for: anyone working on entropy-stable or multi-component interface schemes. The method may well behave well in practice — the numerics are suggestive, and the combination is worth knowing — but the thermodynamic guarantee is not established.\n\nRecommendation: send it to review. The paper deserves referee time; the trade-off it targets is real, and the defect, while load-bearing, is repairable in principle. A corrected proof carrying the [[r]] terms, or a fixed Lemma 4.2 under different dissipation, would patch it; failing that, the claim should be downgraded to numerical evidence plus entropy consistency only in the single-component limit. Don't accept the proof as stated.","headline":"The ES/DF combination is genuinely new and the numerics are serious, but Theorem 4.5's entropy-stability proof rests on two lemmas that are false for variable-r mixtures, so the paper's central claim is unsupported as written.","tokens_in":20571,"tokens_out":7099,"would_cite":true,"duration_ms":67983,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65M08","76N15","65M12"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims a new interface flux that joins Entropy-Stable and Double-Flux methods satisfies a semi-discrete entropy inequality, preserves kinetic energy, and suppresses pressure oscillations at species interfaces with variable…","keywords":["Entropy-Stable schemes","Double-Flux method","multi-component compressible flow","finite-volume methods","kinetic energy preservation","pressure oscillations","adaptive mesh refinement","supersonic flow simulation"],"falsifier":"Evaluate the identity $[\\![r(Y)\\rho]\\!] = \\bar{r}[\\![\\rho]\\!]$ at a face separating pure hydrogen from pure nitrogen at equal pressure and temperature: the jump in $r$ is large, and the equality omits the $\\bar{\\rho}[\\![r]\\!]$ term, so the identity is false exactly at the interfaces the method targets. A numerical corollary would be to run the one-dimensional moving-interface test and record the discrete entropy residual each stage; a scheme delivering the claimed inequality should never show entropy decreasing across the species jump.","tokens_in":19481,"feed_emoji":"💨","tokens_out":10211,"duration_ms":110391,"temperature":0.7,"pith_summary":"The paper's central claim is that a numerical flux can combine two previously incompatible ideas: the Entropy-Stable flux, which guarantees the second law of thermodynamics in discrete form, and the Double-Flux method, which freezes each cell's effective specific-heat ratio and reference energy during a timestep to stop spurious pressure oscillations at species interfaces. The proposed flux is said to satisfy a semi-discrete entropy inequality, to preserve kinetic energy, and to remain low-dissipation while keeping material interfaces oscillation-free. A sympathetic reader would care because uniting these properties in one flux would remove a known trade-off in multi-component compressible flow simulation: entropy stability without interface noise, or oscillation-free Double-Flux without thermodynamic consistency. The paper supports the claim with a theorem, a hybrid dissipation strategy, and benchmark tests ranging from moving interfaces through shock-bubble interaction to an under-expanded hydrogen jet.","feed_headline":"Entropy-stable flux kills interface pressure spikes","feed_subtitle":"Combining two flux families promises oscillation-free, thermodynamically consistent multi-component flow.","key_machinery":"The load-bearing object is the interface flux of Equation (42), assembled from logarithmic and arithmetic means of the reconstructed left and right states, with $\\gamma^*$ and $e_0^*$ frozen per cell as in Double-Flux. Entropy stability is certified through the shuffle condition of entropy-conservative flux theory, which relates the jump of entropy variables to the jump of the entropy potential $\\psi = r(Y)\\rho v$. The dissipation layer uses the eigenvector scaling relation $A_0 = R T R^T$ and the hybrid eigenvalue matrix $|\\Lambda| = (1-\\theta)|\\lambda| + \\theta|\\lambda_{\\max}|I$, blending Roe-type and Lax-Friedrichs-type dissipation by a local pressure-jump indicator. A key identity used in the proof is Lemma 4.1, $[\\![r(Y)\\rho]\\!] = \\bar{r}[\\![\\rho]\\!]$, which assumes the mixture gas constant does not jump across the face.","core_discovery":"On the paper's own terms, the central discovery is a new inviscid numerical flux, $F^{\\mathrm{es}^*}$ in Equation (42), that splices the Entropy-Stable flux for multi-component Euler equations together with the Double-Flux freezing of $\\gamma^*$ and $e_0^*$ per cell. Contracting the semi-discrete scheme with entropy variables reduces the entropy production to the jump expression $\\sum_i r_i Y_i [\\![\\ln Y_i]\\!] \\rho^{\\ln}$, which Lemma 4.2 shows is non-positive; hence the flux is claimed to satisfy the entropy inequality and be Entropy-Stable. The same flux is shown to meet the structural condition for kinetic-energy preservation, and a hybrid dissipation term built from eigenvector scaling is added in Equation (54) without breaking the entropy inequality. If these claims hold, the method delivers a thermodynamically consistent, oscillation-free, low-dissipation discretization for multi-component compressible flows with variable specific-heat ratios.","pith_inferences":["A testable repair of the proof would add dissipation or a correction proportional to $[\\![r]\\!]$ at species interfaces, restoring the jump identity and letting the discrete entropy inequality be checked directly.","Because Double-Flux freezes $\\gamma^*$ and $e_0^*$ per cell over each stage, the scheme is not fully conservative; measuring accumulated mass and energy errors on long interface advection would quantify the trade-off between oscillation-free capture and conservation.","The same flux structure could be ported to real-gas or transcritical mixtures by replacing the frozen perfect-gas closure with tabulated thermodynamics, an extension suggested by the paper's own real-gas references.","If the entropy-stability claim survives the proof repair, the scheme is a natural candidate for large-eddy simulation of reacting high-speed jets, where low dissipation and thermodynamic consistency both matter."],"forward_implications":["Moving species interfaces with different specific-heat ratios are captured without the pressure and velocity oscillations that plain Entropy-Stable fluxes produce.","The semi-discrete scheme obeys the second law in an integral sense, so entropy production is non-negative even near shocks.","Kinetic-energy preservation keeps numerical dissipation low, which should improve under-resolved turbulent and mixing simulations.","The hybrid dissipation adds Lax-Friedrichs-type damping near strong shocks and Roe-type low dissipation elsewhere, preserving sharp oblique shocks.","Benchmark results reproduce reference shock-bubble trajectories and Mach-disk positions for the under-expanded hydrogen jet, matching previous simulations and experiments."],"supporting_citations":[{"why":"supplies the shuffle condition and entropy-conservative flux framework on which the stability proof is built.","marker":"[12]"},{"why":"supplies the kinetic-energy-preserving entropy-conservative flux that Equation (27) extends to multi-component flow.","marker":"[14]"},{"why":"supplies the multi-component Entropy-Stable formulation and eigenvector-scaling construction that the new flux starts from.","marker":"[19]"},{"why":"introduces the Double-Flux idea for multifluid flows, the source of the oscillation-free interface behavior.","marker":"[20]"},{"why":"extends Double-Flux to reacting flows and justifies freezing $\\gamma^*$ and $e_0^*$ per cell during a timestep.","marker":"[29]"},{"why":"provides the stable logarithmic-mean evaluation used in the numerical flux.","marker":"[38]"},{"why":"supplies the entropy-variable mapping and multi-component entropy-stable discretization that defines the variables used in the proof.","marker":"[8]"},{"why":"provides the eigenvector scaling theorem used to build the entropy-stable dissipation term.","marker":"[36]"},{"why":"provides the reference shock-bubble trajectory data used to validate the scheme.","marker":"[55]"}],"fun_headline_variants":["Entropy-stable double-flux kills pressure spikes","Flux combo ensures oscillation-free multi-component flow","Energy-preserving flux tames shock-interface interactions","Double-flux and entropy stability merge for robust flow","Thermodynamic consistency meets low dissipation in new flux"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The proof of entropy stability assumes the mixture gas constant does not change across a cell face, but at the species interfaces the method is designed for — hydrogen against nitrogen, helium against air — that constant jumps sharply, so the key identity used in the theorem does not hold there.","fun_headline_variants_meta":{"raw":{"variants":["Entropy-stable double-flux kills pressure spikes","Flux combo ensures oscillation-free multi-component flow","Energy-preserving flux tames shock-interface interactions","Double-flux and entropy stability merge for robust flow","Thermodynamic consistency meets low dissipation in new flux"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000167,"raw_usage":{"total_tokens":1261,"prompt_tokens":954,"completion_tokens":307,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":570,"completion_tokens_details":{"reasoning_tokens":232}},"tokens_in":570,"tokens_out":307,"duration_ms":4330,"temperature":1.0,"reasoning_tokens":232,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T00:37:33.161896+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the identity $[\\![r(Y)\\rho]\\!] = \\bar{r}[\\![\\rho]\\!]$ at a face separating pure hydrogen from pure nitrogen at equal pressure and temperature: the jump in $r$ is large, and the equality omits the $\\bar{\\rho}[\\![r]\\!]$ term, so the identity is false exactly at the interfaces the method targets. A numerical corollary would be to run the one-dimensional moving-interface test and record the discrete entropy residual each stage; a scheme delivering the claimed inequality should never show entropy decreasing across the species jump.","supporting_citations":[{"cited_title":"Tadmor, The numerical viscosity of entropy stable schemes for systems of conservation laws","cited_arxiv_id":null,"evidence_quote":"supplies the shuffle condition and entropy-conservative flux framework on which the stability proof is built."},{"cited_title":"Gouasmi, K","cited_arxiv_id":null,"evidence_quote":"supplies the multi-component Entropy-Stable formulation and eigenvector-scaling construction that the new flux starts from."},{"cited_title":"Abgrall, S","cited_arxiv_id":null,"evidence_quote":"introduces the Double-Flux idea for multifluid flows, the source of the oscillation-free interface behavior."},{"cited_title":"Billet, R","cited_arxiv_id":null,"evidence_quote":"extends Double-Flux to reacting flows and justifies freezing $\\gamma^*$ and $e_0^*$ per cell during a timestep."},{"cited_title":"Ismail, P","cited_arxiv_id":null,"evidence_quote":"provides the stable logarithmic-mean evaluation used in the numerical flux."},{"cited_title":"Renac, Entropy stable, robust and high-order dgsem for the compressible multicomponent euler equations, Journal of Computational Physics 445 (2021) 110584","cited_arxiv_id":null,"evidence_quote":"supplies the entropy-variable mapping and multi-component entropy-stable discretization that defines the variables used in the proof."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the eigenvector scaling theorem used to build the entropy-stable dissipation term."},{"cited_title":"Terashima, G","cited_arxiv_id":null,"evidence_quote":"provides the reference shock-bubble trajectory data used to validate the scheme."}],"review_version":1}