{"id":"8253ee3b-2250-496d-98c9-1fa3d8b4b01c","arxiv_id":"2608.00746","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":3,"one_line_summary":"The paper derives a WENO-C scheme with relative-smoothness weights and coupled mass-volume fraction reconstruction, satisfying consistency, conservation, equilibrium, and bound preservation for compressible N-phase phase-field flows.","lead":"This paper introduces a new WENO reconstruction scheme for compressible multiphase flow simulations that stays consistent, conservative, and bound-preserving for any number of phases. It matters because standard high-order WENO schemes can create unphysical oscillations at material interfaces, and the new scheme is designed to remove those errors while keeping high accuracy.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Fully discrete bound-preservation guarantee is imported from [20], and the headline applications run at CFL=0.2, above the stated sufficient CFL=1/12, leaving the strongest BP demonstration outside the proven regime.","rationale":"I read the paper as a construction-and-verification paper for a finite-volume WENO reconstruction for compressible multiphase flows. The algebraic construction is carefully built: Eq. (15) defines relative-smoothness weights that are invariant under affine scaling, so the same weights serve mass and volume fraction; Eq. (16) couples the mass and volume-fraction reconstructions; Eq. (17) and Eq. (20) apply shared positivity limiters; Eq. (19) normalizes the volume fractions to unit sum; and Eq. (21) enforces the energy bound. These steps preserve consistency and conservation by using the same limiter or scaling values across components and by redefining U_C through Eq. (6). The numerical verification is extensive: two-phase and five-phase advection, smooth and discontinuous convergence tests, two large applications with and without the Phase-Field mechanism, and an AMR demonstration. I do not see an internal flaw in the WENO-C derivation. The load-bearing soft spot is the fully discrete bound-preservation guarantee. The reconstruction-level admissibility is established, but the claim that the time-stepped solution stays in U is imported from [20] and applies only under a CFL condition. Section 4.2 reports the sufficient CFL as 1/12, then runs the application problems at CFL=0.2. Because no admissibility diagnostics are reported, the central property is not demonstrated in the most impressive simulations. This is exactly the reader's weakest assumption, and I agree with that identification. The concern is addressable without changing the construction: re-run at CFL=1/12 to have a fully proven bound-preserving demonstration, or report a diagnostic showing membership in U at the larger CFL, or add a short proof that the necessary CFL for this flux and reconstruction is at least 0.2. Since the reader already assigned CONDITIONAL, my read does not change the verdict. The lack of released code and data and the absence of a direct comparison with the bound-preserving WENO in [64] are secondary verification gaps, not correctness objections.","tokens_in":24673,"tokens_out":16508,"duration_ms":153981,"concrete_test":"Add a runtime diagnostic that checks, at every Runge-Kutta stage, that the cell-averaged state and the reconstructed states U_L, U_R, U_C satisfy the admissible-set conditions in Eq. (8): mass positivity, volume fractions in [0,1], sum of volume fractions equal to 1, and energy bound psi >= 0. Re-run Section 4.2.1 (shock-induced cavity collapse) and Section 4.2.2 (shock-vessel-bubble interaction) at both CFL=0.2 and CFL=1/12. If violations occur at CFL=0.2 but not at CFL=1/12, the bound-preservation claim must be restricted to the lower CFL; if no violations occur at either CFL, the imported sufficient condition is conservative and the concern is resolved for these configurations.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim 'bound-preserving' is a fully discrete statement, not just a property of the reconstruction. The paper proves in Section 3.4 that WENO-C produces reconstructed states U_L, U_R, U_C in the admissible set U, and then imports the time-stepping guarantee from the author's earlier theorem [20]. That theorem requires bound-preserving numerical fluxes and a CFL condition; Section 4.2 states the sufficient CFL for the fifth-order scheme is 1/12. Nevertheless, both application runs (shock-induced cavity collapse and shock-vessel-bubble interaction) are performed at CFL=0.2, justified only by preliminary tests that the sufficient condition is more restrictive than necessary. No diagnostic is reported showing that all states remained in U during those runs. If the CFL condition is practically relevant, the advertised bound-preserving behavior of the headline simulations is not actually established. This is not an internal inconsistency in WENO-C: under CFL<=1/12 the construction appears sound, and the algebraic consistency, conservation, and equilibrium steps are carefully built. The gap is that the strongest empirical support for the scheme's central property lies outside the proven regime, so the paper's claims should be qualified unless admissibility at CFL=0.2 is verified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript proposes a finite-volume WENO reconstruction scheme, WENO-C, for compressible two-/N-phase flows with a Phase-Field mechanism. The paper first analyzes why the standard WENO-JS scheme produces errors in low-density phases, attributing the failure to the scaling of the smoothness indicators under affine transformations (Eq. 14). It then introduces modified WENO weights based on relative smoothness between stencil candidates (Eq. 15) and combines these with a coupled reconstruction of phasic masses and volume fractions, positivity limiters, volume-fraction normalization, and an energy-bound limiter (Section 3.4). The central claim is that WENO-C is the first WENO scheme satisfying consistency, conservation, equilibrium, and bound preservation for compressible multiphase flows, and, for N phases, consistency of reduction and volume-fraction summation to unity. Numerical tests verify these properties in advection, smooth and shock-interface problems, and applications to shock-induced cavity collapse and shock-vessel-bubble interaction, including an AMR variant. A variant, WENO-CL, is also presented for problems with non-zero non-unity constant volume fractions.","tokens_in":24818,"tokens_out":9996,"duration_ms":101582,"significance":"If the claims are correct, the paper is a valuable practical contribution: it provides a reconstruction procedure that algebraically enforces the coupled requirements that previous WENO implementations for multiphase flows satisfied only partially. The construction is transparent and largely self-contained: the consistency argument for the modified weights is explicit, the coupled mass/volume-fraction reconstruction is clearly motivated, and the positivity and summation steps are stated as algebraic operations rather than fitted parameters. The numerical tests are extensive and include round-off-level measurements of the targeted errors, which is a real strength. The main caveat is that the fully discrete bound-preservation property is not re-proved here but imported from the author's earlier theorem in [20]; the headline application runs are performed at CFL=0.2, above the sufficient CFL=1/12 quoted from [20], and no admissibility diagnostics are reported for those runs. Thus the strongest empirical demonstration of the central property is conditional on an unverified CFL assumption. This is a fixable gap rather than a flaw in the reconstruction itself.","major_comments":[{"comment":"The paragraph beginning \"Following the analysis in [20]\" states a sufficient CFL condition of CFL<=1/12 for bound preservation, but both application problems (shock-induced cavity collapse and shock-vessel-bubble interaction) are run at CFL=0.2, justified only by the remark that preliminary tests show the sufficient condition is more restrictive than necessary. Since the paper's central claim is a fully discrete bound-preservation property and the fully discrete theorem is imported from [20], the headline simulations are outside the proven regime. Please either rerun these problems at CFL<=1/12, or report time-resolved diagnostics (for example, the minimum of alpha_p rho_p, the minimum of the energy function psi(U) from Eq. (9), and the maximum overshoot in sum_p alpha_p) over the full runs at CFL=0.2, and qualify the bound-preservation claim accordingly. As written, the advertised bound-preserving behavior of the two application problems is not actually established.","section":"Section 4.2"},{"comment":"The admissible set U in Eq. (8) depends on user-prescribed phasic minimum pressures (P_min)_p, and the energy limiter in Eq. (21) is applied relative to those values, but the manuscript nowhere reports the values of (P_min)_p used in the numerical tests. Without this information the energy-constraint step is not precisely specified and the reported bound-preservation verifications cannot be reproduced. Please list the prescribed minimum pressures for each phase in the verification and application problems, or state clearly how they were chosen.","section":"Section 4"},{"comment":"The order-of-accuracy statement for WENO-C is argued informally: the text states that the positivity limiter does not influence the order in smooth regions, citing [66,67,68], and that Eq. (19) and Eq. (20) act only in interfacial regions, where the order drops to first order. This reasoning is plausible and the numerical convergence tests in Sections 4.1.4 and 4.1.5 support it, but the statement as written is not a proof. Please make the locality assumption explicit, or add a brief derivation showing that the operations in Eqs. (17)-(21) are inactive or reduce to identity in the smooth regions of the tested problems.","section":"Section 3.4"}],"minor_comments":[{"comment":"The phrase \"serves as the building blog for these variants\" appears to be a typo; it should likely read \"building block\".","section":"Section 1"},{"comment":"In the consistent limiter formula, the case where all gradients S^{(alpha_p)}_{i+1/2} vanish simultaneously leads to a 0/0 ratio in the expressions for theta^{(alpha_p)}. The text should specify that the limiter is skipped (or set to a defined value such as 0) when the denominator is zero, since this case occurs in constant-volume-fraction regions.","section":"Section 3.4.1, Eq. (25)"},{"comment":"The air-helium shock tube is performed at CFL=0.4, which is also above the quoted sufficient CFL=1/12. If this test is intended as a demonstration of fully discrete bound preservation, it should either be run at a CFL within the proven regime or accompanied by admissibility diagnostics; otherwise, the CFL statement in Section 4.2 should be coordinated with the CFL values used in all tests.","section":"Section 4.1.5"},{"comment":"Reference [35] contains a typo in the page range (\"1—-36\"), and the author name in reference [15] is formatted inconsistently. A careful proofreading pass over the reference list is recommended.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The construction in the manuscript is sound and the numerical experiments are extensive; the main weakness is the verification gap for the fully discrete bound-preservation claim at CFL=0.2, which is fixable by rerunning or adding diagnostics. The paper depends heavily on the author's prior theorem in [20], so I would ask the editor to confirm that [20] indeed proves the exact CFL condition quoted here and that its hypotheses are satisfied by the HLL/HLLC fluxes used in the present paper. The novelty claim of being the first WENO scheme satisfying all four requirements should also be checked carefully against recent work such as [64] and the DG-based bound-preserving literature, since the paper's own references suggest active competition in this area."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nQuick take: this is a solid numerical methods paper, and the central construction is genuinely new. The WENO-JS diagnosis in Eq. (14)—that g = a f + b reconstructs like WENO-JS on f but with epsilon scaled by a^2—explains both the consistency failure and the oscillation of low-density masses in one clean line. The proposed relative-smoothness weights in Eq. (15) remove the epsilon sensitivity, and the coupled reconstruction of (alpha_p rho_p) and alpha_p with shared weights and limiters is an elegant way to keep consistency, conservation, and the admissibility set intact. The N-phase extension and the WENO-CL variant are natural, and the numerical verification is unusually thorough: five-phase advection, order tests, and two substantial application problems with and without the phase-field mechanism.\n\nWhere it is soft: the fully discrete bound-preservation guarantee is imported from the author's earlier theorem in [20]. That is legitimate—it is a published result—but the two headline applications run at CFL = 0.2, while the sufficient CFL stated from [20] is 1/12. The paper says preliminary tests show the sufficient condition is too restrictive, and that may be true; however, no diagnostics are reported showing that the states stayed in the admissible set during those runs. The strongest empirical demonstration of the central property therefore sits outside the proven regime. I would not call this fatal—the construction under CFL <= 1/12 appears sound—but the claims should be qualified, or one page of admissibility diagnostics for the CFL = 0.2 runs should be added. A direct comparison with the bound-preserving WENO of Zhang and Cheng [64] would also make the 'first WENO scheme satisfying the four requirements' claim easier to accept; right now the claim rests on comparing with WENO-JS and MUSCL.\n\nThe order-of-accuracy statement for the modified weights is verified numerically, not proven. That is a minor point for a methods paper, but worth noting. No code or data are released, which also weakens reproducibility, though that is common for this literature.\n\nThe math, the data, and the citation pattern look solid. The paper is carefully argued, honestly acknowledges its own limitations, and the central contribution is clear. It deserves a serious referee. I would recommend peer review, with the CFL caveat as the main requested revision.","headline":"A genuinely useful WENO reconstruction for compressible multiphase flows, with a crisp explanation of why WENO-JS misbehaves; the main blemish is that the headline applications run above the proven bound-preserving CFL.","tokens_in":25461,"tokens_out":2306,"would_cite":true,"duration_ms":20666,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65M08","76N10","35L65"],"pacs":["47.55.-t","47.40.-x","02.70.-c"],"model":"deepseek-v4-flash","headline":"The paper claims to be the first WENO scheme for compressible multiphase flows that enforces consistency, conservation, equilibrium, and physical bound preservation simultaneously, and verifies this for two-phase and arbitrary-N-phase flow.","keywords":["Multiphase flows","Compressible flows","Phase-Field methods","WENO schemes","Bound-preserving schemes","Shock-interface interactions","Consistency of reduction","Volume fraction summation"],"falsifier":"Run the two-phase air-water advection problem of Section 4.1.1 at $t=1$ and measure the consistency error $((\\alpha_1 \\rho_1)/\\rho_1 - \\alpha_1)$ and temperature; the paper's claim is that these are at round-off, so any systematic growth above machine precision would falsify the consistency claim. For the bound-preservation claim, run the shock-induced cavity collapse with the energy-bound projection removed at the theorem's stated $CFL \\leq 1/12$ and at $CFL=0.2$: the paper predicts failure whenever the energy bound is skipped, so a completed simulation without the projection would contradict it.","tokens_in":24356,"feed_emoji":"🫧","tokens_out":8121,"duration_ms":65063,"temperature":0.7,"pith_summary":"The paper sets out to fix a specific failure mode: standard WENO reconstructions, when applied to compressible multiphase flows, produce errors in the phasic mass at low-density material interfaces and break thermal equilibrium, and in N-phase settings they create fictitious phases, local voids, or overfilling. The author proposes a WENO scheme (\"WENO-C\") that reconstructs masses and volume fractions in a coupled way using scale-invariant, relative-smoothness weights, then enforces positivity, volume-fraction summation to unity, and an energy bound by projection. The claim is that this is the first WENO scheme satisfying the four critical requirements — consistency, conservation, equilibrium, and bound preservation — for compressible multiphase flows, with the N-phase extra requirements of reduction consistency and unity summation. If the claim is right, high-order WENO methods become usable for shock-interface problems, including large-density-ratio and multi-phase applications, without the unphysical states that currently cause simulation failures.","feed_headline":"Relative-smoothness weights fix WENO for multiphase flows","feed_subtitle":"Coupled mass-volume fraction reconstruction cuts consistency errors to round-off and stops fictitious phases.","key_machinery":"The central object is the modified WENO weight $\\omega^{(s)} = d^{(s)} / (d^{(s)} + \\sum_{k \\neq s} d^{(k)} (IS^{(s,k)})^2)$, where $IS^{(s,k)}$ is the ratio of the smoothness indicators of stencils $s$ and $k$ (set to 1 if both vanish). Because ratios are scale-invariant, the weights sense discontinuities in the phasic mass even when the phasic density is tiny. Around this weight formula, the scheme couples the reconstruction of $\\alpha_p \\rho_p$ and $\\alpha_p$ using the same weights, applies a positivity limiter, rescales to make $\\sum_p \\alpha_p = 1$, re-limits interior states, and projects the full conservative state onto the energy bound; the conservation requirement is kept by reconstructing the cell polynomial so that its Gauss-Lobatto quadrature matches the cell average.","core_discovery":"On its own terms, the paper's central discovery is a diagnosis and a fix. It shows that the standard WENO-JS weights fail for quantities like the phasic mass $g = a f + b$ because the smoothness indicators scale as $a^2$ and the epsilon-paranoia parameter effectively becomes $\\epsilon/a^2$, which is not small when $a$ is a small density; hence non-smooth stencils are not deactivated and the mass oscillates. The fix replaces absolute smoothness with relative smoothness ratios between stencil pairs, making the weights scale-invariant, and then builds the reconstruction so that the volume fraction and the mass use the same weights, positivity and unity-sum steps, and an energy-bound projection, guaranteeing by construction that reconstructed face and interior states lie in the admissible set.","pith_inferences":["The relative-smoothness weight formula suggests a general design principle: any WENO variant whose smoothness indicators are homogeneous of degree two will inherit the $\\epsilon/a^2$ failure when applied to scaled quantities; replacement by ratios should be a drop-in fix for those variants too.","The fact that the demonstration runs use $CFL=0.2$ while the sufficient bound-preservation condition from the earlier theorem is $CFL \\leq 1/12$ hints that the theoretical CFL condition is conservative; if so, the practically stable range of WENO-C is wider than proven, but that remains open.","The energy-bound projection, not merely mass positivity, appears to be what lets high-order schemes survive shock-interface impacts; the failure table suggests that low-order dissipative schemes can mask the missing bound, which points to a testing protocol for future high-order multiphase methods.","WENO-CL, which preserves non-zero non-unity constant volume fractions, is positioned for miscible or reacting flows; one testable extension is to use the same limiter for mass-fraction transport in combustion models, where the sum-to-unity constraint is on mass fractions."],"forward_implications":["WENO-JS should not be applied directly to phasic masses in low-density phases; the paper's two-phase advection tests show oscillations in mass and temperature that WENO-C removes to round-off.","For N-phase flows with $N \\geq 3$, a consistent WENO reconstruction need not produce fictitious phases, local voids, or overfilling; the five-phase and four-phase advection tests verify this.","High-order WENO schemes can be made bound-preserving for compressible multiphase flows, which the paper argues is necessary because skipping the energy-bound step made the shock-induced collapse simulations fail.","The modified weights can be adopted by other WENO variants (e.g., WENO-M, WENO-Z, TENO), since the analysis operates on the weight construction rather than on a particular variant.","The scheme integrates with adaptive mesh refinement and with both the Allaire and Kapila five-equation models, giving it a wide operational range."],"supporting_citations":[{"why":"States the theorem that reconstructed states in the admissible set and bound-preserving fluxes give a bound-preserving scheme, which the new scheme verifies by construction.","marker":"[20]"},{"why":"Provides the base WENO-JS weights whose epsilon-scaling failure is analyzed and fixed.","marker":"[25]"},{"why":"Furnishes the consistent and conservative phase-field model and the consistency requirement on mass and volume fraction.","marker":"[18]"},{"why":"Supplies the consistent limiter used in the WENO-CL variant and the reduction-consistent formulation.","marker":"[19]"},{"why":"Establishes the primitive-variable reconstruction used to maintain pressure equilibrium.","marker":"[26]"},{"why":"Gives the positivity limiter formula used to enforce mass and volume-fraction bounds.","marker":"[66]"},{"why":"Shows how positivity limiters can be added to high-order schemes without destroying accuracy.","marker":"[67]"},{"why":"Defines the five-equation model with $K_p=0$ used in the application problems.","marker":"[2]"},{"why":"Defines the five-equation model with $K_p=\\rho c^2/(\\rho_p c_p^2)-1$ used in the application problems.","marker":"[28]"}],"fun_headline_variants":["Relative smoothness keeps WENO stable for multiphase","WENO fix stops spurious phases in compressible flows","Scale-invariant WENO weights for consistent multiphase","New WENO scheme guarantees bounds in multiphase flows","Relative smoothness cures WENO mass oscillations"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The scheme's guarantee that the fully discrete solution stays in-bounds inherits the bound-preservation theorem from earlier work; the paper checks the reconstructed states are admissible but does not re-prove that theorem, and some demonstration runs use a larger CFL number than the theorem's sufficient condition.","fun_headline_variants_meta":{"raw":{"variants":["Relative smoothness keeps WENO stable for multiphase","WENO fix stops spurious phases in compressible flows","Scale-invariant WENO weights for consistent multiphase","New WENO scheme guarantees bounds in multiphase flows","Relative smoothness cures WENO mass oscillations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00068,"raw_usage":{"total_tokens":3085,"prompt_tokens":935,"completion_tokens":2150,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":551,"completion_tokens_details":{"reasoning_tokens":2066}},"tokens_in":551,"tokens_out":2150,"duration_ms":14032,"temperature":1.0,"reasoning_tokens":2066,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T15:17:38.392547+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the two-phase air-water advection problem of Section 4.1.1 at $t=1$ and measure the consistency error $((\\alpha_1 \\rho_1)/\\rho_1 - \\alpha_1)$ and temperature; the paper's claim is that these are at round-off, so any systematic growth above machine precision would falsify the consistency claim. For the bound-preservation claim, run the shock-induced cavity collapse with the energy-bound projection removed at the theorem's stated $CFL \\leq 1/12$ and at $CFL=0.2$: the paper predicts failure whenever the energy bound is skipped, so a completed simulation without the projection would contradict it.","supporting_citations":[{"cited_title":"Huang and E","cited_arxiv_id":null,"evidence_quote":"States the theorem that reconstructed states in the admissible set and bound-preserving fluxes give a bound-preserving scheme, which the new scheme verifies by construction."},{"cited_title":"Efficient implementation of weighted eno schemes.J","cited_arxiv_id":null,"evidence_quote":"Provides the base WENO-JS weights whose epsilon-scaling failure is analyzed and fixed."},{"cited_title":"Huang and E","cited_arxiv_id":null,"evidence_quote":"Furnishes the consistent and conservative phase-field model and the consistency requirement on mass and volume fraction."},{"cited_title":"Huang and E","cited_arxiv_id":null,"evidence_quote":"Supplies the consistent limiter used in the WENO-CL variant and the reduction-consistent formulation."},{"cited_title":"Johnsen and T","cited_arxiv_id":null,"evidence_quote":"Establishes the primitive-variable reconstruction used to maintain pressure equilibrium."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the positivity limiter formula used to enforce mass and volume-fraction bounds."},{"cited_title":"Zhang and C.-W","cited_arxiv_id":null,"evidence_quote":"Shows how positivity limiters can be added to high-order schemes without destroying accuracy."},{"cited_title":"Allaire, S","cited_arxiv_id":null,"evidence_quote":"Defines the five-equation model with $K_p=0$ used in the application problems."},{"cited_title":"Kapila, R","cited_arxiv_id":null,"evidence_quote":"Defines the five-equation model with $K_p=\\rho c^2/(\\rho_p c_p^2)-1$ used in the application problems."}],"review_version":2}