{"id":"8e4ca3de-8d52-4b5a-9e7d-936c265448c2","arxiv_id":"2607.29224","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A low-Mach-number flow solver for closed systems uses Newton iteration on the nonlinear equation of state so that global mass stays exactly constant as the thermodynamic pressure evolves.","lead":"This paper presents a numerical method for simulating slow-moving flows of real fluids, such as supercritical CO2, inside sealed containers, keeping the total mass exactly constant while the background pressure changes. It matters because it makes accurate heat-transfer simulations feasible in regimes where fully compressible solvers are too slow.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Mass-conservation update (Eq. 14d) is mathematically sound, but the uncharacterized constant-coefficient pressure splitting in §3.3 (Eqs. 17–18) is the weakest load-bearing link for the framework's claimed accuracy and robustness in real-gas regimes.","rationale":"The reader's weakest_assumption correctly identifies the pressure-splitting projection as the load-bearing approximation for the framework's overall accuracy and robustness, even though the global mass-conservation property itself is guaranteed by construction via the Newton iteration on Eq. (14d). The paper explicitly flags the stability constraint on ρ* in a footnote but gives no quantitative support—no reported ρ* values, no convergence or sensitivity tests, and no error estimate relative to a variable-coefficient solver. This is a genuine gap that affects the credibility of the validation results in the strongly nonlinear regimes claimed as the method's target. The mass-preservation claim is nonetheless sound and well-supported by the derivation, so the concern does not justify rejection; it warrants a conditional acceptance pending the additional characterization. The reader's verdict of CONDITIONAL remains appropriate, and my analysis does not move it.","tokens_in":46,"tokens_out":9673,"duration_ms":229236,"concrete_test":"Re-run the transcritical CO2 channel (Sec. 5.2.2) and the VdW differentially heated cavity (Sec. 5.1.2) using a variable-coefficient Poisson solver (or an iterative solver that resolves the full 1/ρ^k+1 pressure gradient) and compare mean velocity, mean temperature, Nusselt number, and TKE budgets. Additionally, report the chosen ρ* for each run, and perform a sensitivity test by repeating the CO2 case with ρ* = ρ_min and ρ* = 0.5 ρ_min; if the statistics shift by more than the typical DNS sampling error, the splitting error is significant and the accuracy claim is overstated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central mass-preservation claim rests on the Newton–Raphson update (Eq. 14d), which enforces M0 = ∫ρ(p0,T)dV at every RK stage for any EoS with tabulated isothermal compressibility. That part is constructionally sound and does not depend on the pressure splitting. However, the abstract also claims accuracy and robustness for strongly nonlinear real-gas flows. The segregated algorithm's efficiency and stability hinge on the pressure-splitting projection in §3.3: Eq. (17a) replaces the variable-density pressure gradient with a constant-coefficient term using a reference density ρ*, plus an extrapolated pressure p*1 from the previous two steps (Eq. 17b). The footnote states ρ* must be ≤ min(ρ) for stability, citing Dong & Shen (2012) and Demou et al. (2019), but the paper reports no actual ρ* values, no sensitivity study, and no error estimate for the splitting at the tested density ratios (≈7 in the VdW cavity, ≈3 in the CO2 channel). The correction step (Eq. 18c) enforces the divergence constraint exactly regardless of the splitting error, so global mass is preserved, but the momentum error can bias the velocity field, temperature advection, and thus all flow statistics and heat-transfer predictions. At larger density ratios or stronger pseudo-boiling peaks, where β and χ vary sharply, the linearly extrapolated p*1 may be far from the true hydrodynamic pressure, and the constant-coefficient Poisson solve may not be a reliable approximation. This is the primary load-bearing risk for the practical validity of the proposed framework beyond the tested cases.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a low-Mach-number numerical method for closed systems with a general, possibly tabulated, equation of state. The thermodynamic pressure p0 is advanced at each Runge–Kutta stage by a Newton–Raphson iteration that enforces global mass conservation, using the isothermal compressibility to form the Jacobian. The method is implemented in the CaNS finite-difference framework with a segregated temperature/momentum update, WENO for the temperature equation, and a constant-coefficient pressure-splitting projection so that FFT-based Poisson solvers can be retained. The method of manufactured solutions is used for verification, followed by validation against a differentially heated cavity (ideal gas and van der Waals fluid), ideal-gas turbulent channel flow, and transcritical CO2 turbulent channel flow.","tokens_in":22984,"tokens_out":7858,"duration_ms":80514,"significance":"If the results hold, the main contribution is a genuinely EoS-agnostic, mass-preserving update of the thermodynamic pressure for closed-system low-Mach-number flows. Equation (14d) is parameter-free and reduces exactly to the classical ideal-gas update in the linear limit, which is a clean and useful property. The use of isothermal compressibility allows direct coupling to tabulated thermodynamic data, and the validation against the transcritical CO2 DNS is encouraging. The paper therefore addresses a real gap in low-Mach-number real-gas solvers, where closed-system global mass conservation is often only approximate.","major_comments":[{"comment":"The abstract and §4.2 state that the manufactured-solution study verifies thermodynamic-pressure evolution, but Table 2 reports only L2(eT), L2(eu), and L2(ev) with convergence orders. No error norm, convergence order, or pointwise comparison for p0 is reported. Since the p0 update in Eq. (14d) is the central novelty, the verification is incomplete without at least L2(p0) or a comparison of p0(t) against the manufactured p0,M(t). Please add this result or explicitly state what aspect of the thermodynamic-pressure evolution is being verified.","section":"§4.2, Table 2"},{"comment":"The constant-coefficient pressure splitting replaces the variable-density pressure gradient with 1/ρ*∇p1 + (1/ρ − 1/ρ*)∇p*1. The paper gives only the stability constraint ρ* ≤ min(ρ) in footnote 2, citing Dong & Shen (2012) and Demou et al. (2019), but does not report the ρ* values used in the VdW differentially heated cavity (density ratio 7.34) or the CO2 channel (density ratio ≈3), nor a sensitivity study, nor an error estimate for the splitting. This is load-bearing for the claimed accuracy and robustness in real-gas regimes: an inaccurate hydrodynamic pressure due to the splitting can bias the velocity field, temperature advection, and turbulence statistics even though global mass conservation remains exact. Please report the actual ρ* values and either provide a sensitivity study with respect to ρ*/minρ or an error estimate for the tested conditions.","section":"§3.3, Eqs. (17a)–(18c)"},{"comment":"The transcritical channel is initialized at p0(0) = 8.4 MPa, which the authors state was chosen so that the steady-state thermodynamic pressure is close to the constant reference pressure of 8 MPa used by Wan et al. (2025). This tuning weakens the independence of the validation. Because the central claim is about closed-system mass preservation, please report the thermodynamic-pressure history p0(t) and the global mass residual over the transient and statistically stationary periods. If possible, show sensitivity to the initial pressure (e.g., starting at 8.0 MPa) or otherwise quantify how the closed-system evolution differs from the constant-pressure reference.","section":"§5.2.2, supercritical CO2 channel"},{"comment":"The manufactured velocity field is constructed from continuity, Eq. (23), while the projection method enforces the thermodynamic divergence constraint, Eq. (8). The text states that non-zero residuals are added to the temperature transport equation, but it does not state whether the corresponding manufactured source term also enters the divergence constraint in Eq. (8). If the source is included only in the temperature equation, the enforced divergence differs from the manufactured velocity divergence by an extra βS/(ρCp) term, and the verification of the projection step is incomplete. Please clarify the construction or adjust the manufactured fields so that Eq. (8) is satisfied consistently with the added source.","section":"§4.1, Eqs. (8) and (23)"}],"minor_comments":[{"comment":"The paper claims that the Newton–Raphson procedure converges in at most two to three iterations for all cases considered, but no iteration counts or convergence tolerance are reported. Adding a short table or statement of the typical iteration count would support this claim.","section":"§3.1, Eqs. (14a)–(14d)"},{"comment":"The turbulent transcritical CO2 case is presented at a single grid resolution. A coarse-grid comparison, even qualitative, would strengthen the confidence in the turbulence statistics.","section":"§5.2.2"},{"comment":"The caption states that the ordinate multiplier is shown on the axis, but the printed axis label does not display the numerical multiplier. Including the multiplier on the axis would improve readability.","section":"Figure 11"},{"comment":"The reference density ρ* is introduced with only a footnote on the stability condition. It would help to state explicitly how ρ* is selected in practice and whether it is fixed or updated during the simulation.","section":"§3.3, Eq. (17a)"}],"recommendation":"major_revision","confidential_remarks":"The central mass-conservation construction is sound and the paper has a clear contribution. The missing p0 convergence data and the uncharacterized pressure-splitting approximation are, in my view, fixable by additions rather than by reworking the method. The transcritical validation is useful but would be notably strengthened by reporting p0(t) and mass conservation, given the deliberate pressure initialization. I would not reject the paper; after a careful revision addressing the verification and sensitivity questions, it could be acceptable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The core contribution is real and clean: for closed-system low-Mach real-gas flows, the thermodynamic pressure is updated by a Newton–Raphson iteration on global mass, with the derivative expressed through isothermal compressibility (Eq. 14d). That is a neat, parameter-free idea that works for tabulated equations of state and reduces exactly to the ideal-gas formula. The thermodynamic derivations in Eqs. (7)–(9) and Appendices A–B checked out consistently. The validation is also credible: ideal-gas cavity against Le Quéré, channel flows against Nicoud, and the transcritical CO2 case against Wan et al. shows good agreement in mean profiles and TKE budgets. The VdW cavity at Mach ~1e-7 is a good demonstration of why this class of solver matters.\n\nWhere the paper is softer: the verification section claims to verify thermodynamic-pressure evolution, but Table 2 reports errors only for temperature and velocity, no p0 error. That is a simple omission and easy to fix, but it matters because the paper's selling point is precisely the p0 update. Second, the transcritical validation initializes p0 to 8.4 MPa specifically so the steady state lands near the reference 8 MPa. The authors discuss this honestly, and the agreement with Wan et al. suggests the comparison is fair, but it is a tuned setup rather than an independent prediction. Third, the pressure-splitting projection in §3.3 is the weakest load-bearing link. The stability condition on ρ* ≤ min(ρ) is cited rather than characterized, no ρ* values are reported, and there is no error estimate or sensitivity study for density ratios up to ~7. The reader's stress-test note is right about this being the main practical risk. For the tested cases it works, but robustness at larger density contrasts is unquantified. I would not call this a fatal flaw; it is an underdocumented part of an otherwise sound method. No code or data are shipped, which is a legitimate minor complaint for a method paper.\n\nNet assessment: this is a well-executed incremental advance, not a field re-organizer. The math is solid, the validation is genuinely external, and the central claim holds up. The soft spots are real but repairable with modest additions. I would send this to peer review; a competent referee can verify the derivation in an afternoon and the missing p0 error and ρ* details are straightforward to request.","headline":"Solid, useful method paper: the Newton–Raphson thermodynamic-pressure update is a genuine, well-derived contribution, but the pressure-splitting projection is undertested and the verification omits the very quantity it claims to verify.","tokens_in":23474,"tokens_out":1213,"would_cite":true,"duration_ms":11653,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65M06","76D05","76F65"],"pacs":[],"model":"deepseek-v4-flash","headline":"A low-Mach-number flow solver now enforces exact global mass conservation for real fluids with nonlinear equations of state.","keywords":["low-Mach-number flows","real-gas equation of state","thermodynamic pressure","global mass conservation","pressure-correction projection","transcritical CO2","method of manufactured solutions","variable-property turbulence"],"falsifier":"Run a closed-system differentially heated cavity with the same van der Waals fluid at progressively larger wall temperature differences, producing wall density ratios of roughly 10, 20, and 50, and compare mean Nusselt number and velocity profiles with a variable-coefficient Poisson solver as ground truth. If the mass-preserving solution's velocity statistics depart from the reference as the density ratio grows, the pressure-splitting approximation is the cause.","tokens_in":1477,"feed_emoji":"⚖️","tokens_out":2196,"duration_ms":65480,"temperature":0.7,"pith_summary":"This paper sets out a low-Mach-number numerical framework for closed systems filled with real fluids governed by general nonlinear equations of state. Its central aim is to keep the system's total mass exactly constant by evolving the spatially uniform thermodynamic pressure so that the equation of state integrates to the initial mass at every Runge-Kutta stage. The thermodynamic update is done with a Newton-Raphson iteration whose derivative uses isothermal compressibility, so it works with tabulated thermodynamic data rather than a specific analytic equation of state. Because the thermodynamic state is advanced before momentum and the pressure projection is recast as a constant-coefficient Poisson solve, the method stays computationally efficient. If it holds up, closed-system simulations of supercritical fluids, where thermodynamic-pressure evolution changes heat-capacity peaks, no longer need to accept formal mass drift.","feed_headline":"Low-Mach solver enforces exact total mass for real-gas flows","feed_subtitle":"A thermodynamic-pressure update built on isothermal compressibility makes mass conservation exact for any equation of state.","key_machinery":"The load-bearing identity is the Newton-Raphson update p0^(m+1) = p0^m + (M0 - ∫ρ dV) / (∫ρχ dV), where χ = (1/ρ) ∂ρ/∂p|T is the isothermal compressibility. This converts the global mass constraint into a root-finding problem whose derivative is a standard thermodynamic quantity, so the method is agnostic to the equation of state's analytical form. The second mechanism is the pressure-splitting projection: the reciprocal density in the pressure gradient is replaced by a constant reference reciprocal density, chosen at or below the minimum density, plus a correction built from a linearly extrapolated pressure. That turns a variable-coefficient Poisson equation into a constant-coefficient one","core_discovery":"The central claim is that global mass conservation in a closed low-Mach-number system can be enforced exactly and efficiently for any single-phase equation of state by solving the nonlinear constraint M0 = ∫ρ(p0,T)dV at each stage. The paper derives a Newton-Raphson update for p0 whose residual derivative is the volume integral of ρ times isothermal compressibility, a quantity available from analytic equations of state or property tables; in the ideal-gas limit the update reduces to the standard explicit formula. The algorithm separates thermodynamic from hydrodynamic updates: temperature is advanced first, p0 is corrected, density and properties come from the equation of state, and only the","pith_inferences":["The exact-mass property is independent of the projection's accuracy: p0 is forced to satisfy ∫ρ dV = M0 even if the pressure-splitting velocity error grows at large density ratios, so future work should separate these two errors in convergence tests.","Because the update uses only ρ and χ, the same iteration could be applied to mixtures or to tabulated single-fluid models beyond single-component single-phase cases, though the paper only presents the latter.","A direct test of the splitting assumption would be a closed cavity with wall density ratios above the tested range (roughly 7 in the van der Waals case and 3 in the transcritical channel), comparing velocity statistics against a variable-coefficient Poisson or fully compressible reference; divergence of statistics with unchanged mass conservation would isolate the splitting error.","The Newton iteration converges in two to three iterations because the previous stage's p0 is a good initial guess; cases with very rapid pressure transients or noisy table interpolation could need more iterations or a safeguarded root-finder."],"forward_implications":["In the ideal-gas limit the Newton-Raphson update collapses to the known explicit thermodynamic-pressure formula, so existing low-Mach solvers are recovered as a special case.","Any fluid whose tables supply isothermal compressibility can be simulated without coding a new analytic equation of state; the paper demonstrates this with lookup-table CO2 properties.","Closed-system real-fluid simulations no longer need to prescribe a constant thermodynamic pressure; the pressure evolves and mass is conserved during thermal transients.","The method reaches regimes where fully compressible solvers are impractical, such as differentially heated cavity flows at characteristic Mach numbers near 10^-7.","Second-order spatial accuracy in velocity and temperature holds for a nonlinear equation of state, as verified by manufactured solutions."],"fun_headline_variants":["Exact mass preservation for any EOS in low-Mach flows","Newton-Raphson pressure step guarantees global mass in closed systems","Exact mass for real gases: compressibility-based pressure update","Low-Mach closed systems: exact mass via Newton-updated p0","Mass-conserving low-Mach solver for any equation of state"],"cache_read_input_tokens":24704,"weakest_assumption_plain":"The scheme assumes that replacing the variable-density pressure gradient with a constant-reference-density gradient plus a linearly extrapolated pressure correction remains accurate enough at the density contrasts of interest; the paper reports no error estimate for that splitting, so a failure there would corrupt velocities while leaving the global mass property intact.","fun_headline_variants_meta":{"raw":{"variants":["Exact mass preservation for any EOS in low-Mach flows","Newton-Raphson pressure step guarantees global mass in closed systems","Exact mass for real gases: compressibility-based pressure update","Low-Mach closed systems: exact mass via Newton-updated p0","Mass-conserving low-Mach solver for any equation of state"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000805,"raw_usage":{"total_tokens":3355,"prompt_tokens":708,"completion_tokens":2647,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":452,"completion_tokens_details":{"reasoning_tokens":2567}},"tokens_in":452,"tokens_out":2647,"duration_ms":17772,"temperature":1.0,"reasoning_tokens":2567,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T11:15:10.939380+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a closed-system differentially heated cavity with the same van der Waals fluid at progressively larger wall temperature differences, producing wall density ratios of roughly 10, 20, and 50, and compare mean Nusselt number and velocity profiles with a variable-coefficient Poisson solver as ground truth. If the mass-preserving solution's velocity statistics depart from the reference as the density ratio grows, the pressure-splitting approximation is the cause.","supporting_citations":[],"review_version":1}