{"id":"38c9a9e0-3fbf-49d4-a976-1cc6e31d9cd8","arxiv_id":"2607.27426","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"low","formal_verification":"none","parameter_count":4,"one_line_summary":"Constant-transfer bubble pressure and vapor mass can be advanced at CHyQMOM radius–velocity nodes to close polydisperse bubbly-flow equations without extra mixed moments, matching EL screens to 1.5% RMS.","lead":"The authors add heat and mass transfer to an Euler–Euler moment model of polydisperse bubbles by evolving pressure and vapor mass only at existing quadrature nodes. That keeps the transported moment set small while matching Monte Carlo and Euler–Lagrange bubble-screen pressures to about 1.5% RMS.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"No significant objection identified beyond the already-disclosed pb/mv collocation and dilute-regime limits.","rationale":"The reader’s strongest claim matches the abstract and Secs. 4–5 evidence, and the weakest assumption (delta pinning of pb/mv at CHyQMOM nodes, Sec. 2.5) is correctly identified as the main modeling soft spot. That choice is intentional cost control and is partially validated by MC node/mean comparisons and by the 1.5% EE–EL screen pressure error; the paper also discloses non-realizability under strong shocks/small R0. No internal inconsistency, missing baseline, or untested leap undermines the stated dilute-screen results. Novelty is incremental (constant-transfer on existing CHyQMOM) but the verification package and artifacts are solid for a methods paper. Stress-testing does not surface a stronger load-bearing failure mode, so the ACCEPT verdict and low correctness risk stand.","tokens_in":14673,"tokens_out":631,"duration_ms":47065,"concrete_test":"Extract the centerline liquid-pressure history from the nb=51 EE bubble-screen run; drive an offline ensemble of Ns≥1000 independent constant-transfer Keller–Miksis ODEs (same IC statistics as Sec. 5) with that pl(t); compare the MC mean of R^3 pbw to the QBMM four-node closure over t/t0∈[0,20]. If relative RMS on that closure exceeds ~2–3% (comparable to Table 1 pb error), the collocation under screen forcing is weaker than the 1.5% EE–EL figure suggests.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central engineering claim holds under the paper’s own evidence. The load-bearing modeling choice is still the one the reader flags: after Eq. 23 (Sec. 2.5), pb and mv are pinned as deltas at the four (R,Ṙ) CHyQMOM nodes per R0 bin, with no mixed-moment hierarchy, and are stored outside the conservative state (Sec. 3) so they are not fluxed—only corrected through ∂Rj,i/∂t from the moment equations (25)–(27). That approximation could in principle bias the ensemble closures R^3 pbw and R^2 Ṙ that enter mixture pressure (2) and void fraction (6) once conditional variance of pb|R becomes large. However, the manuscript already stress-tests the claim where it matters: monodisperse harmonic MC (Table 1, Figs. 2–3) bounds mean R/pb/mv error; nb-convergence and σ-sweeps quantify screen closure error; polytropic contrast shows the expected damping; and the 3D screen gives 1.5% RMS vs the mean of 40 volume-averaged EL runs. Those checks, plus open pinned code, make the collocation a disclosed applicability bound rather than an untested hole in the reported results.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The manuscript extends conditional hyperbolic quadrature-of-moments (CHyQMOM) modeling of polydisperse bubbly flows by advancing constant-transfer ODEs for bubble pressure pb and vapor mass mv at the four joint (R, Ṙ) nodes per equilibrium-radius bin, rather than closing pb with a polytropic relation. Node values supply the ensemble averages needed by the mixture pressure and void-fraction equations without transporting mixed pb/mv moments. The model is implemented in MFC, verified against monodisperse Monte Carlo under harmonic forcing, tested on 1D bubble screens for nb and initial-distribution sensitivity, contrasted with the polytropic closure, and validated in 3D against the mean of 40 volume-averaged Euler–Lagrange realizations (1.5% relative RMS center-pressure error).","tokens_in":14993,"tokens_out":1219,"duration_ms":24746,"significance":"If the reported checks hold, this is a useful and incremental advance for Euler–Euler sub-grid bubbly-flow modeling: it brings heat and mass transfer into an existing CHyQMOM framework at modest cost, documents open code and seeds, and supplies concrete verification (MC node evolution, Table 1 RMS errors), closure-error quantification, and a nontrivial 3D EE–EL comparison. The suppression of high-frequency polytropic ringing under the tested screen conditions is a clear, falsifiable modeling consequence. Strengths include reproducible archival code, explicit disclosure of the delta-collocation assumption, and multi-level evidence rather than a single demonstration case.","major_comments":[{"comment":"Sec. 2.5 (after Eq. 23) and Sec. 3: pb and mv are collocated as deltas at the (R, Ṙ) nodes and stored outside the conserved state, so they are not fluxed—only corrected via ∂Rj,i/∂t from Eqs. (25)–(27). The 1.5% EE–EL agreement and MC checks support the claim under the tested dilute, moderate-forcing screens, but the manuscript does not bound when conditional variance of pb|R would invalidate the collocation (e.g., stronger collapse or broader thermal states). A short applicability statement or one additional diagnostic (node-wise pb spread vs. an MC reference under a harder drive) would make the central closure claim more durable without changing the method.","section":"Sec. 2.5, Eqs. (23)–(27); Sec. 3"},{"comment":"Sec. 6 states instability at large pressure ratios / small R0 (non-realizable moments), shared with the prior polytropic implementation, and notes failures for strong shocks past immersed boundaries. This is disclosed, but it is load-bearing for the claimed sub-grid-scale scope. The paper should either (i) quantify the pressure-ratio / R0 envelope of the presented screen and harmonic cases, or (ii) show that the variance floor (Eq. 28) and RK treatment do not silently degrade the ensemble closures near that envelope. Without that, the 1.5% EE–EL result risks being read as general beyond the dilute acoustic-screen regime actually tested.","section":"Sec. 6; Eq. (28)"}],"minor_comments":[{"comment":"Table 1: the 1.4× error increase for constant-transfer vs. polytropic is attributed to an inexact initial condition; the pseudo-polytropic control is helpful, but a one-sentence note on how pb0, mv0 are set from the isothermal assumption (Eq. 23) would clarify reproducibility of the MC comparison.","section":"Table 1; Eq. (23)"},{"comment":"Figs. 4–6: RMS errors are relative to nb=91; stating the absolute center-pressure scale (or peak p/p0) in the caption would help readers judge whether 10^{-2}–level relative errors are dynamically important.","section":"Figs. 4–6"},{"comment":"Fig. 8: the EE curve exits the EL min–max envelope by ~4% of p0 in early phases; a brief remark in Sec. 5.5 on whether this is sampling, grid coarsening, or collocation bias would prevent over-reading the 1.5% RMS alone.","section":"Sec. 5.5; Fig. 8"},{"comment":"Notation: ∂Rj,i/∂t vs. Ṙj,i (Eqs. 24–25) is important but easy to miss; emphasizing once in the text that equality holds only in the ensemble average would aid readers implementing the method.","section":"Sec. 2.5, Eqs. (24)–(25)"},{"comment":"Minor prose: abstract and introduction repeat the four-node / no-mixed-moments claim several times; a single crisp statement would tighten the front matter. Also fix “arXiv:2607.27426v1 … 29 Jul 2026” date anomaly if it is a typesetting artifact in the source.","section":"Abstract; Introduction"}],"recommendation":"minor_revision","confidential_remarks":"The work is a natural and well-executed extension of the authors’ prior polytropic QBMM and EL papers; novelty is real but incremental. Fit for a multiphase/fluids computational journal is good. I would not block on the collocation assumption given the disclosed checks and open code; minor revision to bound applicability is enough. No integrity or citation-pattern concerns."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The one thing worth knowing is that this is a clean engineering fix, not a new theory. They take Preston-style constant-transfer ODEs for bubble pressure and vapor mass and advance them only at the four existing second-order CHyQMOM (R, Ṙ) nodes per R0 bin, so the mixture pressure and void-fraction closures pick up heat/mass transfer without enlarging the transported moment set. That is the actual new result relative to their own polytropic QBMM line and to Fox-style conditional inversion.\n\nWhat they do well is the evidence package. Monodisperse harmonic forcing is checked against Ns=1000 Monte Carlo with tabulated RMS on R, pb, mv and a pseudo-polytropic control that isolates the initial-condition error. Bubble-screen runs quantify nb convergence and σR/σṘ/σR0 sensitivity. The polytropic contrast shows the expected suppression of high-frequency ringing. The 3D screen lands at 1.5% relative RMS versus the mean of 40 volume-averaged Euler–Lagrange realizations, with open code at a pinned commit. Citations sit where they should (Preston, Fox/CHyQMOM, their prior MFC/QBMM/EL work). Math and numerics are standard and readable; the ∂Rj,i/∂t correction path for the non-conserved node states is spelled out.\n\nSoft spots are real but already on the table and proportionate. After Eq. 23 they pin pb and mv as deltas at the (R, Ṙ) nodes and store them outside the conservative state, so they are not fluxed—only corrected through the moment equations. That can bias R3 pbw once conditional variance of pb|R grows. They also inherit the usual dilute/spherical/no-slip/no-breakup limits and the known non-realizable-moment instability at strong shocks and small R0. None of that refutes the reported tests; it bounds the regime. Free parameters (variance floor, initial widths, nb, ensemble sizes) are ordinary for this class of method.\n\nThis is for people who actually run dilute polydisperse bubbly CFD in naval hydro or cavitation and need a usable heat/mass-transfer SGS inside an Euler–Euler solver. It deserves a serious referee. I would engage with it and expect to cite the coupling and the EL comparison if I am working in the same lane.","headline":"Solid incremental SGS upgrade: node-local constant-transfer inside CHyQMOM without mixed moments, backed by MC, screen sweeps, and 1.5% EE–EL agreement, with the delta-collocation of pb/mv as a disclosed bound rather than a hidden hole.","tokens_in":15689,"tokens_out":609,"would_cite":true,"duration_ms":11668,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Constant-transfer heat and mass transfer can be closed inside a quadrature bubble model by advancing pressure and vapor mass only at existing radius–velocity nodes, without new mixed moments.","keywords":["bubbly flow","population balance","quadrature-based moment methods","sub-grid-scale modeling","heat and mass transfer","compressible multiphase flow","CHyQMOM","bubble screen"],"falsifier":"Repeat the three-dimensional bubble-screen comparison with the same forcing and initial polydispersity; if the Euler–Euler center-pressure relative RMS error against a well-converged ensemble of volume-averaged Euler–Lagrange runs substantially exceeds the reported 1.5%, or if high-frequency oscillations reappear while the polytropic model remains quieter, the closure claim fails.","tokens_in":15467,"feed_emoji":"💧","tokens_out":970,"duration_ms":17500,"temperature":0.7,"pith_summary":"Ensemble-averaged models of many bubbles need statistics of an evolving population, but earlier quadrature moment methods closed bubble pressure with a polytropic relation that ignores heat and mass transfer at the wall. This paper puts the constant-transfer equations for internal bubble pressure and vapor mass into a conditional hyperbolic quadrature method. Four joint radius–radial-velocity nodes are obtained per equilibrium-radius bin; pressure and vapor mass are marched only at those nodes and then used to close the mixture pressure and void-fraction equations. No mixed pressure or vapor-mass moments are added to the transported set. Monte Carlo checks confirm the nodes and mean bubble variables under harmonic forcing. In bubble-screen tests the constant-transfer closure removes the high-frequency pressure oscillations seen with the polytropic model, and three-dimensional comparisons reach 1.5% relative RMS error against the mean of forty volume-averaged Euler–Lagrange runs. A sympathetic reader cares because this keeps heat and mass transfer inside a cheap Euler–Euler framework while staying close to more expensive Lagrangian truth under the conditions examined.","feed_headline":"Bubble heat and mass transfer closed at four quadrature nodes","feed_subtitle":"Constant-transfer model matches Lagrangian bubble screens to 1.5% RMS without extra mixed moments","key_machinery":"Node-collocated constant-transfer closure: second-order conditional hyperbolic inversion yields four (R, Ṙ) abscissae per R0 bin; pb and mv are integrated at those abscissae alone and their nodal values close the barred mixture terms (R3 pbw, R2 Ṙ, etc.).","core_discovery":"Advancing the constant-transfer equations for bubble pressure and vapor mass only at the four joint radius–radial-velocity CHyQMOM nodes per equilibrium-radius bin supplies the ensemble averages needed by the compressible mixture equations without enlarging the transported moment set. Under the tested conditions the resulting model matches Monte Carlo means, suppresses the high-frequency polytropic pressure oscillations, and agrees with the mean of forty volume-averaged Euler–Lagrange bubble-screen realizations to 1.5% relative RMS error.","pith_inferences":["If the delta-function pinning of pb and mv remains accurate only while the joint (R, Ṙ) distribution stays near-Gaussian, strong shock-driven collapses may still force a return to higher-dimensional moment inversions.","The reported stability limit at large pressure ratios and small R0 is likely the practical barrier to using the model for shock–bubble-cloud interaction until realizability limiters or sub-stepping are added.","The same collocation pattern could be tried for other internal bubble variables (e.g., dissolved-gas mass) without immediately exploding the moment vector."],"forward_implications":["Heat and mass transfer can be retained in Euler–Euler bubbly-flow simulations without transporting extra mixed pb/mv moments.","Under the tested screen conditions the constant-transfer closure removes the high-frequency natural-frequency ringing that the polytropic closure produces.","Closure error grows mainly with broader equilibrium-radius distributions and remains modest under the reported discretizations of R0.","The same nodal-collocation idea can be dropped into existing high-order compressible multiphase solvers that already carry CHyQMOM moments."],"fun_headline_variants":["Four CHyQMOM nodes close bubble heat and mass transfer","Constant-transfer equations advanced only at quadrature nodes","No mixed moments: pressure and vapor mass at four nodes","Quad-node constant-transfer beats polytropic oscillations","1.5% RMS match to EL bubble screens via four-node closure"],"cache_read_input_tokens":128,"weakest_assumption_plain":"Bubble pressure and vapor mass are treated as delta functions sitting exactly on the existing radius–velocity quadrature nodes, so they never receive their own independent moment hierarchy.","fun_headline_variants_meta":{"raw":{"variants":["Four CHyQMOM nodes close bubble heat and mass transfer","Constant-transfer equations advanced only at quadrature nodes","No mixed moments: pressure and vapor mass at four nodes","Quad-node constant-transfer beats polytropic oscillations","1.5% RMS match to EL bubble screens via four-node closure"]},"model":"grok-4.5","effort":"low","cost_usd":0.003764,"raw_usage":{"total_tokens":1189,"prompt_tokens":788,"num_sources_used":0,"completion_tokens":67,"cost_in_usd_ticks":37644000,"prompt_tokens_details":{"text_tokens":788,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":334,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":788,"tokens_out":67,"duration_ms":6179,"temperature":1.0,"reasoning_tokens":334,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-31T01:15:40.302874+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Repeat the three-dimensional bubble-screen comparison with the same forcing and initial polydispersity; if the Euler–Euler center-pressure relative RMS error against a well-converged ensemble of volume-averaged Euler–Lagrange runs substantially exceeds the reported 1.5%, or if high-frequency oscillations reappear while the polytropic model remains quieter, the closure claim fails.","supporting_citations":[],"review_version":1}