{"id":"a5e55eed-862c-4ff2-a245-85793f5ff910","arxiv_id":"2608.01828","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A unified simulation of Leidenfrost water droplets reveals that ambient humidity and internal circulation control evaporation, and that axisymmetry is invalid for large droplets.","lead":"Using direct numerical simulations, this paper shows that ambient humidity and the flow circulating inside a Leidenfrost droplet substantially change droplet shape and evaporation rate, and that the common axisymmetric modeling assumption fails for large droplets. The findings point toward more accurate predictions for spray cooling and other heat transfer applications, and warn against simplified vapor-only models.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Axisymmetry claim rests on simplified 3D model that omits humidity; full-3D support is missing.","rationale":"The paper's most novel and surprising conclusion is that the unrealistic constraint of axisymmetry explains the prolate shapes seen in axisymmetric simulations of large Leidenfrost droplets. This claim goes beyond the already-notable humidity+circulation results. The reader's weakest assumption correctly identifies the simplified 3D model (Section 8.2) as the load-bearing support. I agree: that model omits the mixed gas-vapour phase and humidity, which the paper itself demonstrates are first-order effects on evaporation and thermal gradients (Sections 6–7). Its axisymmetric form does not match the full MGV model, so the 3D result is not clearly transferable. The linear stability analysis shows instability but not the resulting shape; the 3D simulation is the only evidence linking instability to puddle-like shapes. This is a genuine gap, but it does not undermine the humidity/circulation results, which are supported by comparisons to experiments (Bleiker & Specht 2007, Burton et al. 2012) and by scaling-law recovery. Thus conditional acceptance remains appropriate: the humidity/circulation findings can be accepted with high confidence, but the axisymmetry attribution requires full 3D confirmation or a quantitative validation of the simplified model. The lack of code/data also hampers reproduction but is secondary to the physics gap.","tokens_in":25353,"tokens_out":5517,"duration_ms":62014,"concrete_test":"","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract's claim that the large-droplet discrepancy is due to the constraint of axisymmetry rests entirely on Section 8.2's simplified coupled NS+lubrication model (Eqs. 8.2–8.3). That model (i) uses a patching angle of π/4 and surface-tangent gradients, both acknowledged approximations; (ii) omits the mixed gas-vapour phase and ambient humidity, which Sections 6–7 show strongly affect evaporation and thermal gradients; and (iii) in axisymmetric form yields shapes even more prolate than the full MGV model (Fig. 14), so it does not quantitatively reproduce the full model. The 3D result, while suggestive, could therefore be an artifact of the omitted physics rather than a demonstration that axisymmetry is the cause for the full model. The linear stability analysis (Section 5) establishes instability but not the post-critical shape; only the simplified 3D simulation links instability to puddle-like shapes. Without either a full 3D simulation of the MGV model or at least a quantitative axisymmetric validation of the simplified model against the full model, the axisymmetry attribution is unsupported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops an axisymmetric finite-element model of Leidenfrost water droplets in a mixed gas-vapour environment, with internal droplet flow, evaporation, and a non-isothermal droplet. Simulations are run from droplet radii of about 0.05 mm to several millimetres, spanning the puddle, spherical, and take-off regimes. The authors compare four models: mixed gas-vapour with and without internal circulation, and pure vapour with and without circulation. They report that ambient humidity and internal circulation together change the evaporation flux distribution and droplet height, increase the global evaporation rate, and improve agreement with measured droplet lifetimes. For large droplets, the axisymmetric model with circulation produces prolate, 'wimple'-shaped droplets that disagree with experiments; an azimuthal linear stability analysis shows instability for R ≳ 0.2 mm, and 3D simulations of a simplified coupled Navier-Stokes/lubrication model yield puddle-like shapes. The abstract concludes that the axisymmetric constraint is the cause of the discrepancy.","tokens_in":25626,"tokens_out":5162,"duration_ms":65790,"significance":"If the conclusions hold, the paper makes a useful contribution by coupling ambient humidity, internal circulation, and evaporation in one model and by providing a quantitative bridge between the lubrication-theory regime (Sobac et al. 2014) and the small-droplet take-off regime (Sobac et al. 2025). The paper is not curve-fitted to the target experiments: thermophysical properties are taken from prior literature, and the model recovers the known scalings h~R^{-1/2} and J~R^{9/5}. The azimuthal stability analysis (Section 5) and the systematic viscosity- and humidity-sweep methodology (Table 1) are valuable even if the final attribution to axisymmetry is not yet fully proven. The comparison with experimental evaporation times (Fig. 12d) is a concrete falsifiable prediction. However, the paper's strongest advertised claim, namely that the large-droplet shape discrepancy is due to the axisymmetry constraint, is supported only by a reduced 3D model that omits the very humidity effects that Sections 6-7 show to be important.","major_comments":[{"comment":"The central claim in the abstract, that the large-droplet discrepancy 'is due to the unrealistic constraint of axisymmetry', is not yet supported quantitatively. The only 3D evidence comes from the reduced coupled NS-lubrication model, which (i) uses a patching angle of π/4 and surface-tangent gradients, (ii) omits the mixed gas-vapour phase and ambient humidity, and (iii) in axisymmetric form gives shapes that are even more prolate than the full MGV model (Fig. 14). The 3D result is suggestive, but the axisymmetric reduced model does not reproduce the full model, so the change to puddle-like shapes could be affected by the omitted physics or by the patching approximation rather than solely by relaxing axisymmetry. I recommend either a full 3D MGV simulation at least at one representative large-droplet condition, or an explicit validation of the reduced model against the full MGV axisymm","section":"§8.2, Eqs. (8.2)-(8.3), Fig. 14"},{"comment":"The azimuthal stability analysis establishes linear instability of the axisymmetric base state, but not that the nonlinear saturated state is the puddle-like shape seen in the reduced 3D model. The growth-rate calculation also uses material properties evaluated at the film temperature and a fourth-order Taylor expansion of the saturation pressure (Appendix A.1), so the reported critical radius of about 0.2 mm is an approximation. The paper should state explicitly that the linear analysis gives only the onset, not the post-critical shape, and should quantify the sensitivity of the critical radius to the film-temperature and Taylor-expansion choices. This is not a fatal flaw, but it is load-bearing for the abstract's attribution claim.","section":"§5, Appendix A.1"},{"comment":"The model overpredicts internal droplet velocities by about an order of magnitude relative to the PIV measurements of Bouillant et al. (2018). The paper acknowledges this, but the overprediction is nevertheless used to explain the prolate shapes and to argue that internal circulation is necessary for accurate drying kinetics. Since the magnitude of the internal circulation is a key ingredient in the proposed mechanism, the conclusions would be more robust if the authors showed that the main results are preserved when the circulation is artificially limited to velocities consistent with experiments (e.g., by a modest viscosity increase or by including surfactant-induced Marangoni stresses). Without such a sensitivity test, the quantitative claim that the MGV-with-circulation model 'best aligns' with experimental lifetimes may be partly fortuitous.","section":"§5, Fig. 7; §7, Fig. 12(d)"},{"comment":"The statement that ambient humidity has a 'significant impact' on geometry and drying kinetics is clearly demonstrated in the stable small-to-moderate regime (Fig. 10 and Fig. 6), but for the largest droplets the axisymmetric MGV model becomes unphysical (prolate) and the instability argument takes over. The paper should be careful to distinguish where the humidity effect is a robust model prediction and where it is intertwined with the unresolved axisymmetry issue. This is partly a presentation point, but it affects the interpretation of the abstract's first claim.","section":"§6, Fig. 10"}],"minor_comments":[{"comment":"Units for viscosity are written as 'Pa s−1' in Eqs. (A1) and (A2) and 'mPa s−1' in Table 2; these should be 'Pa s' and 'mPa s', respectively.","section":"Appendix A.2, Table 2 and Eqs. (A1)-(A5)"},{"comment":"The caption says 'where 0 is at the drop of the bottom'; this should be 'at the bottom of the drop'.","section":"Fig. 9 caption"},{"comment":"The mixing rules are said to follow Poling et al. (2000), while the data source in Fig. 15 is given as Poling et al. (2008). Both references are listed, but the text should clarify which edition is used for the mixing rules and which for the pure-species data.","section":"§2.1, Eqs. (2.10)-(2.11)"},{"comment":"The hydrostatic boundary condition is defined by an integral with a minus sign. If p_hydro(z) is intended to be the local hydrostatic pressure, the sign convention should be stated more explicitly, since it is easy to misread.","section":"§3.1, Eq. (3.1)"},{"comment":"The evaporation time is obtained by integrating the quasi-stationary evaporation rate. This is valid only if the quasi-stationary assumption holds over the entire integration range; the paper notes this for take-off, but it would be helpful to state it directly at Eq. (7.10) as well.","section":"§7, Eq. (7.10)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a serious and technically impressive contribution that is likely to be accepted after revision. The main decision is whether the central axisymmetry claim needs a full 3D MGV computation. I think a full 3D MGV run is desirable but perhaps not strictly necessary if the abstract and conclusions are reworded to describe the 3D result as strong qualitative support rather than proof. The overprediction of internal velocities is a separate concern; a sensitivity study with reduced circulation would materially strengthen the paper's quantitative claims. I would not reject: the model synthesis, the humidity/circulation coupling, and the stability analysis are substantial contributions even in their current form."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a carefully built computational study that genuinely shows ambient humidity plus internal circulation changes Leidenfrost evaporation and shape, and it earns serious referee time. The second headline claim—that the prolate-shape discrepancy is 'due to the unrealistic constraint of axisymmetry'—is plausible but only backed by a simplified 3D model, so I would not let the abstract outrun the evidence.\n\nWhat is actually new: the authors run one model across four decades of drop radius, connecting the puddle regime to the take-off regime. They recover the known scalings (h ~ R^{-1/2}, J ~ R^{9/5}, J ~ R) as validation, which is honest. The mixed gas-vapour phase coupled to internal droplet circulation is a real step beyond Chakraborty et al. and Yim et al. The result that the film layer does not dominate evaporation in the large-drop MGV model, and that this improves lifetime predictions against Bleiker & Specht, is a substantive finding. The azimuthal stability cascade m=3→2→1 is also new and qualitatively matches experiment.\n\nThe paper is methodologically careful: material properties from standard sources, no fitted parameters, domain-size and mesh checks, and a quasi-stationary assumption that is explicitly tested. The authors also flag their own limitations, which is credit.\n\nSoft spots, in proportion. The central axisymmetry claim rests on Section 8.2's coupled NS+lubrication model. That model omits the mixed gas-vapour phase and humidity that the rest of the paper shows matter; it uses a patching angle of π/4 and surface-tangent gradients; and in axisymmetric form it produces even more prolate shapes than the full MGV model, so it does not quantitatively reproduce the model whose failure it is meant to explain. The 3D result is suggestive, but it could be an artifact of the omitted physics. The stability analysis itself establishes instability but not the post-critical shape, and the simplified material-property assumptions there are reasonable but not negligible. The model also overpredicts internal velocities by about an order of magnitude, so the mechanism driving the prolate shapes may be quantitatively off. No code or data is provided, which limits direct reproduction, though the method description is unusually detailed.\n\nWho is this for: anyone working on Leidenfrost modelling, droplet evaporation, or heat-transfer applications. It deserves a serious referee. My recommendation: send it to peer review, and ask for either a full 3D MGV simulation of at least one large droplet or a clear quantitative comparison between the simplified model's axisymmetric result and the full model. If that comparison cannot be made convincing, the axisymmetry sentence in the abstract should be softened to a hypothesis.","headline":"Strong computational paper with a real humidity/circulation result; the bolder axisymmetry claim is suggestive but not yet proven.","tokens_in":26092,"tokens_out":1069,"would_cite":true,"duration_ms":14161,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Ambient humidity and internal droplet circulation, not just the vapour film, control Leidenfrost droplet shape and evaporation; large droplets need 3D simulations.","keywords":["Leidenfrost effect","droplet evaporation","ambient humidity","internal circulation","azimuthal instability","axisymmetric simulation","lubrication model","direct numerical simulation"],"falsifier":"Run a full 3D DNS of a roughly 1.6 mm radius Leidenfrost water droplet with mixed gas-vapour phase, internal circulation, and no axisymmetry constraint: if it still yields prolate shapes, the axisymmetry attribution is wrong; if it yields puddle shapes with an azimuthal mode cascade comparable to PIV experiments, the claim is supported. A simpler check is to measure the onset of azimuthal perturbations near R≈0.2 mm and compare the critical radius to the predicted bound.","tokens_in":25252,"feed_emoji":"💧","tokens_out":4599,"duration_ms":52373,"temperature":0.7,"pith_summary":"The paper argues that the standard picture of a Leidenfrost droplet—a pure-vapour film, an isothermal drop, and an axisymmetric shape—misses the two ingredients that set the drying rate and geometry: ambient humidity and the circulation inside the drop. Using one numerical model across four decades of droplet radius, it connects the puddle and take-off regimes and shows that humidity plus internal flow produces evaporative cooling of the drop's top, thinner vapour layers, faster evaporation, and lifetimes closer to experiment. It then shows that axisymmetric simulations of large droplets predict prolate shapes because they forbid azimuthal symmetry breaking; a stability analysis finds a critical radius near 0.2 mm, and a simplified 3D model recovers puddle shapes. If correct, the work implies pure-vapour and isothermal assumptions are inadequate, and large-Leidenfrost-droplet modelling must be three-dimensional.","feed_headline":"Leidenfrost droplets dry faster with humidity and internal flow","feed_subtitle":"Pure-vapour, isothermal models over-predict lifetimes; a single simulation across four decades of radius shows why.","key_machinery":"The central object is a finite-element direct numerical simulation of an axisymmetric Leidenfrost droplet with full Navier-Stokes flow in liquid and gas, species transport of vapour in a mixed gas-vapour phase, a saturation condition at the interface, and evaporative cooling, run quasi-stationarily over droplet radii from roughly 0.05 mm to 5 mm. Circulation is controlled by artificially raising liquid viscosity to isolate its effect. To test whether the prolate shapes are an artefact of geometry, the authors perform azimuthal eigenmode stability analysis, perturbing the axisymmetric base state with wavenumber m, and then a simplified 3D coupled Navier-Stokes/lubrication model with a plane o","core_discovery":"Within one quasi-stationary model, increasing the realism from pure vapour to a mixed gas-vapour phase and allowing internal circulation changes the evaporation mechanism: strong internal flow creates thin thermal boundary layers, cooling the top of the droplet by about 10 K, reducing evaporation at the base, and increasing the contribution from the outer surface, so that the film no longer dominates evaporation. This raises the global evaporation rate and lowers the droplet toward the plate, matching observed lifetimes better. The same model, when constrained to axisymmetry, produces prolate shapes for large droplets; linear stability analysis shows the axisymmetric base state is unstable t","pith_inferences":["This suggests a practical rule of thumb: below roughly 0.2 mm radius, axisymmetric simulations of small and take-off Leidenfrost droplets remain trustworthy; above it, quantitative comparison to experiment needs 3D dynamics.","If the simplified 3D model's conclusion holds, the remaining lifetime gap for large droplets may be closed by full 3D simulations that also resolve humidity inside the gas film, rather than by adding only Marangoni effects.","A testable extension: droplets whose azimuthal instability is suppressed—for example by contamination or confinement—should show prolate shapes closer to the axisymmetric prediction."],"forward_implications":["Pure-vapour and isothermal-droplet assumptions under-predict evaporation rates and over-predict lifetimes; both humidity and circulation must be included.","In large droplets, evaporation is not dominated by the thin vapour film; the outer surface contributes substantially when circulation is present.","Axisymmetric models with internal circulation are unreliable for droplets larger than roughly 0.2 mm in radius, because azimuthal symmetry breaking removes the prolate shapes.","The single model stitches together the puddle regime and the take-off regime, recovering known scalings J~R^{9/5}, J~R, and h~R^{-1/2}."],"fun_headline_variants":["Humidity and flow accelerate Leidenfrost droplet drying","Leidenfrost evaporation: humidity and circulation are key","Axisymmetric models fail for large Leidenfrost droplets","Mixed gas-vapour and internal flow improve Leidenfrost model","Leidenfrost drying rate rises with humidity and internal flow"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The conclusion that axisymmetry, rather than missing physics, causes the prolate shapes rests on the simplified 3D coupled Navier-Stokes/lubrication model being a faithful representation of the full system, including its patched shear coupling and neglect of humidity in the film.","fun_headline_variants_meta":{"raw":{"variants":["Humidity and flow accelerate Leidenfrost droplet drying","Leidenfrost evaporation: humidity and circulation are key","Axisymmetric models fail for large Leidenfrost droplets","Mixed gas-vapour and internal flow improve Leidenfrost model","Leidenfrost drying rate rises with humidity and internal flow"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000377,"raw_usage":{"total_tokens":1825,"prompt_tokens":703,"completion_tokens":1122,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":447,"completion_tokens_details":{"reasoning_tokens":1045}},"tokens_in":447,"tokens_out":1122,"duration_ms":11005,"temperature":1.0,"reasoning_tokens":1045,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T20:06:50.390441+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a full 3D DNS of a roughly 1.6 mm radius Leidenfrost water droplet with mixed gas-vapour phase, internal circulation, and no axisymmetry constraint: if it still yields prolate shapes, the axisymmetry attribution is wrong; if it yields puddle shapes with an azimuthal mode cascade comparable to PIV experiments, the claim is supported. A simpler check is to measure the onset of azimuthal perturbations near R≈0.2 mm and compare the critical radius to the predicted bound.","supporting_citations":[],"review_version":1}