{"id":"028f2dbf-9d6e-41a0-a4cf-428be2af005f","arxiv_id":"2505.22379","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A lubrication model with vapor absorption and Marangoni effects predicts that absorption can trigger droplet coalescence and regime transitions in films flowing down vertical fibers.","lead":"This paper builds a mathematical model of a water-absorbing silicone oil film flowing down a vertical fiber. The model predicts that vapor absorption, together with surface tension changes, can make droplets merge, and it maps the conditions where merging happens.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The coalescence regime boundary hinges on a non-quantified J0; the paper itself admits the absorption scale is uncertain, yet the central phase diagram is computed only for an unvalidated Λ range.","rationale":"The reader's weakest_assumption is exactly the constitutive law J = J0(c − cs) and the scale of J0, and the paper itself admits Λ is hard to estimate. The stress-test identifies the same gap as the most load-bearing concern: the entire claim of regime transitions and coalescence (Figures 9–12) is made in a parameter range Λ ∈ [0,5], but there is no physical measurement or independent bound on Λ, so the predicted coalescence window may not exist at experimentally accessible conditions. I do not find a separate internal inconsistency strong enough to shift the verdict to REJECT. The linear stability analysis is consistently derived from the model once the frozen-time approximation is accepted; the unphysical cs = 0.2 is a deliberate and disclosed modeling choice, so it weakens the quantitative prediction but does not invalidate the theoretical construction. The missing code/data is a genuine reproducibility concern, but not the central claim's logical weak point. The paper credits prior literature, the derivation is internally coherent, and the qualitative mechanism—a spatially decaying concentration profile combined with a Marangoni-driven speed increase can bring droplets together—is plausible. The crucial missing piece is the physical scale of the absorption rate; the regime map is not anchored to the physical system until J0 is measured or inferred. A quasi-static absorption measurement of Dow XX-8810 is the natural decisive check. In the meantime, an update to the manuscript that justifies the range of Λ from data, or reframes the results as a parametric study with a clear experimental target, would be the minimal condition for stronger acceptance.","tokens_in":26100,"tokens_out":2063,"duration_ms":21553,"concrete_test":"Set up a quasi-static saturation experiment for Dow XX-8810 at a controlled humidity (e.g., 80% RH) and measure the absorbed mass as a function of time for a thin film of known thickness; fit J(c) to infer J0. If the inferred J0 yields Λ = LJ0/(HV) in the range 0.001–0.5, the regime map is physically relevant. Independently, rerun the linear stability calculation for cs = 0.86 and a realistic initial concentration, using the full 2×2 dispersion relation from (4.5b)–(4.5c) without dropping the bc/bh coupling term; if the effective growth rate deviates from (4.8) by more than 20% for the parameter range explored in Figure 12, the phase-boundary interpretation is not supported by the linear analysis.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's central claim is that vapor absorption, characterized by Λ, triggers a transition from Regime I to Regime II, with the phase boundary in Figure 12 as a key prediction. But Λ is not measured: Λ = LJ0/(HV), and in §3.2 the authors state that “Λ is more difficult to estimate due to the uncertainty in the J0 scale.” They then arbitrarily choose Λ ∈ [0,5]. The entire regime map is computed for this guessed range. If the true J0 for Dow XX-8810 is an order of magnitude smaller or larger, the values Λ = 0.0083 and Λ = 0.32 may be off by orders of magnitude, so the existence of a coalescence window at physically realizable conditions is unsupported. In addition, the stability analysis in §4.2 intentionally uses cs = 0.2 rather than the estimated cs ≈ 0.86 from §3.2. The paper justifies this as needed to make absorption persist, but it means the linear-stability argument is made at an unphysical saturation level: for the actual cs ≈ 0.86, the base-state concentration decays much more slowly, altering the neutral curves and the quasi-static profile C(z) in §5.1. The regime boundaries in Figure 12, and their dependence on Ma, are therefore only established for the model when artificially loaded with a nonphysical low cs. The O(χ) analysis also drops the term proportional to bc/bh exp(λc − λh) after asserting that λc,r < 0, but the approximation in (4.8) has not been tested against the full coupled dispersion relation, so the effective growth rate underlying the dynamic regime interpretation remains an uncontrolled simplification. The most load-bearing issue, however, is the unquantified Λ: the parameter map in Figure 12 cannot be converted into a prediction for the physical silicone‑oil system until J0 is measured or bounded.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a one-sided lubrication model for a thin film of water-absorbing silicone oil flowing down a vertical cylindrical fibre. The model couples an evolution equation for the film thickness h with an advection-diffusion equation for the oil concentration c, and includes gravity, surface tension, Marangoni transport, and a non-mass-conserving absorption flux J(c)=c-cs. The authors perform a linear stability analysis of the coupled system around time-dependent base states using a frozen-time approximation, and then carry out extensive numerical simulations with Dirichlet inlet and Neumann outlet boundary conditions. The central results are the identification of two droplet regimes (Regime I: no coalescence; Regime II: coalescence) and a phase diagram in the (Ma, Λ) plane, with a critical Marangoni number Mac ≈ 8.28, together with an approximate mass-balance model that rationalizes the dependence of the averaged liquid mass on Λ.","tokens_in":26467,"tokens_out":6440,"duration_ms":63616,"significance":"If the model is faithful, it extends classical fibre-coating theory to a non-mass-conserving setting and proposes a concrete mechanism, vapor absorption, for triggering droplet coalescence and regime transitions. The derivation is coherent and reduces to known mass-conserving models in the limits Ma=Λ=0; the numerical exploration is systematic and the quasi-static concentration comparison in Figure 8b is a useful check. The paper provides testable predictions (the shape of the regime boundaries and the coalescence thresholds), although the primary axis of the phase diagram, Λ, is not calibrated to a measured mass-transfer coefficient. The presentation is generally clear, and the authors are candid about the uncertainty in the absorption scale.","major_comments":[{"comment":"The regime-transition thresholds and the phase diagram in Figure 12 depend on Λ, yet §3.2 states that 'Λ is more difficult to estimate due to the uncertainty in the J0 scale' and the simulations simply set Λ ∈ [0,5]. Consequently, the reported values ΛI→II ≈ 0.0083 and ΛII→I ≈ 0.32 for Ma=20, and the shaded Regime II region in Figure 12, are not tied to the Dow XX-8810 system; an order-of-magnitude change in J0 would shift these thresholds correspondingly. The central quantitative claim should be reframed as a model prediction over a hypothesized Λ range, or supplemented by a sensitivity discussion of how the regime boundaries depend on J0 / Λ.","section":"§3.2 and Figure 12"},{"comment":"The linear stability analysis is deliberately performed with cs=0.2 while the PDE simulations in §5 use the estimated cs≈0.86, and the authors acknowledge this choice is made to make absorption persist. This changes the base-state dynamics in (4.5a): for cs≈0.86 the concentration decays much more slowly and the effective growth rate (4.8) and neutral curve (4.9) would differ. Since the stability results are invoked to interpret the regime transitions, the paper should either repeat the stability calculation at cs=0.86 or explicitly demonstrate (e.g., with a supplementary plot) that the conclusions about absorption-driven instability are insensitive to the choice of cs.","section":"§4.2, Eqs. (4.8)-(4.9)"},{"comment":"The effective film-thickness growth rate (4.8) is obtained by dropping the Marangoni coupling term proportional to (bc/bh)exp(λc−λh) in (4.5b) and the corresponding exponential term in Γ, based on the assertion that λc,r<0. This is an uncontrolled truncation: it has not been tested against the full coupled linear system (4.5b)-(4.5c), and during the initial transient both exponents are small. The comparison in Figure 6b tests the final prediction but does not isolate the error introduced by this neglect. A direct comparison of (4.8) with numerical solutions of the coupled linearized equations for representative parameters is needed to validate the approximation.","section":"§4.2, Eqs. (4.5b)-(4.8)"},{"comment":"The mass approximation (5.18) fixes hmin=0.455 'based on numerical observations' and then uses that same approximation to explain the trends in ⟨Ml⟩T. This is a post-hoc calibration, not a predictive derivation, so the sentence 'this figure concludes that the change in total mass can be estimated by equation (5.18)' overstates the closure. This issue does not affect the phase diagram, but the claim should be softened or an independent estimate of hmin should be provided.","section":"§5.3, Eqs. (5.17)-(5.18)"}],"minor_comments":[{"comment":"The text uses 'Raleigh-Plateau regime' but the correct spelling is 'Rayleigh-Plateau'; please correct this typo.","section":"§3 (before boundary conditions)"},{"comment":"In the sentence defining km, 'absoption parameter' should be 'absorption parameter'.","section":"§4.2"},{"comment":"The reference to 'Burelbachet al. 1988' in the frozen-time discussion is a typo and should read 'Burelbach et al. 1988'.","section":"§4.2"},{"comment":"The classification of Regime I versus Regime II in Figure 12 appears to be based on visual inspection of coalescence events; the paper would benefit from stating an explicit quantitative criterion (e.g., a threshold in spacing variance or number of peak crossings) so that the phase boundaries are reproducible.","section":"§5.2, Eq. (5.12) and Fig. 12"},{"comment":"The notation for time averages in (5.6) (⟨X⟩T) and spatiotemporal averages in (5.12) (⟨X⟩) is similar and could be confused; consider using a different symbol, such as an overbar or double bracket, for the double average.","section":"§5.3"}],"recommendation":"major_revision","confidential_remarks":"The paper fits the scope of the journal and the central idea is attractive, but the quantitative predictive claims (especially Figure 12) rest on an uncalibrated Λ and a stability analysis performed at a nonphysical cs. These issues are addressable by reframing the claims and adding the requested sensitivity checks. The mass-approximation concern in §5.3 is secondary but should be tidied."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know: this paper gives a careful derivation of a coupled thickness-concentration thin-film model for a silicone-oil film on a vertical fibre with vapor absorption, including Marangoni effects and a non-mass-conserving absorption flux. The derivation is clean and reduces correctly to the usual Craster-Matar and Kalliadasis-Chang limits. The numerical finding that absorption, for intermediate Λ and Marangoni number above about 8.28, widens a coalescing regime is a genuinely new result in fiber-coating theory.\n\nWhat the paper does well: the model is well-posed in the sense that total silicone-oil mass is conserved while total liquid mass grows, and the base-state ODEs track the PDE averages nicely in the linear stage. The quasi-static concentration profile in (5.11) is elegant and compares favorably to full simulations. The authors also cite the fiber-coating and volatile-film literature thoroughly, and they are honest about several limitations in the conclusions.\n\nThe soft spots, in proportion: the biggest one is the absorption parameter Λ. It is defined as LJ0/(HV), and the paper admits in §3.2 that J0 is uncertain. The whole regime map in Figure 12 is computed on the guessed range Λ ∈ [0,5]. If the true mass-transfer coefficient is an order of magnitude off, the thresholds ΛI→II and ΛII→I move out of the plotted window, so the predictive claim for the physical Dow XX-8810 system is not yet supported. This is an acknowledged but unresolved gap. Second, the linear stability analysis uses cs=0.2 instead of the estimated cs≈0.86, because absorption dies too fast at the physical value. That is an honest statement, but it means the stability results in §4 are for an artificially long-lived absorption regime; the paper goes back to cs=0.86 only in the numerics. Third, the dropping of the Marangoni coupling term in (4.8) rests on λc,r<0, which is true, but the magnitude of the dropped term is never tested against the full coupled dispersion relation. The authors could check this in a few lines; they don't. Fourth, the mass approximation in §5.3 fixes hmin=0.455 from the numerics and then uses that approximation to explain the mass trend—a data-informed fit, labeled as such, but not a derivation. None of these are fatal, but together they mean the paper is a strong model-formulation paper with a plausible phase diagram, not a validated prediction.\n\nWho gets value from this: people working on fiber coating, thin films with phase change, and droplet dynamics in dehumidification/water-harvesting contexts. It deserves a serious referee, and I would send it out. I would ask the authors for code or detailed convergence data, a sensitivity study on Λ with realistic J0 bounds, and a direct check of the (4.8) approximation against the full linear system. With those, it could be a solid JFM-type contribution.","headline":"A coherent new lubrication model for absorbing films on fibres; the droplet-coalescence phase diagram is plausible but rests on an unquantified absorption parameter and a nonphysical stability analysis.","tokens_in":27028,"tokens_out":2043,"would_cite":true,"duration_ms":27247,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["76A20","76D45","76E17"],"pacs":[],"model":"deepseek-v4-flash","headline":"Vapor absorption can force film droplets on a vertical fibre to merge.","keywords":["thin liquid films","vertical fibre","vapor absorption","droplet coalescence","Marangoni effect","Rayleigh-Plateau instability","non-conservative film flow","lubrication model"],"falsifier":"Run a controlled fibre experiment with a water-absorbing silicone oil at fixed fibre radius, flow rate, and inlet concentration, varying only the surrounding humidity (which sets the absorption parameter), and track droplet peak positions over time. The model predicts that at a Marangoni number of 20 the droplet train remains non-coalescing for absorption parameters below about 0.0083 and above about 0.32, but coalesces in between; observing no coalescence across a continuous humidity sweep, or coalescence at all humidities, would falsify the regime diagram. A direct check of the constitutive assumption is also possible: measure the absorption flux as a function of concentration difference for films of different thickness, since the predicted regime structure depends on the flux being linear.","tokens_in":2064,"feed_emoji":"💧","tokens_out":5689,"duration_ms":119061,"temperature":0.7,"pith_summary":"This paper proposes a lubrication model for a water-absorbing silicone oil film flowing down a vertical fibre, in which mass is not conserved because the film takes up water vapor. The central claim is that vapor absorption, acting through a non-mass-conserving source and a concentration-gradient Marangoni force, can destabilize the usual Rayleigh-Plateau droplet train and trigger droplet coalescence. The model identifies two regimes, a non-coalescing Regime I and a coalescing Regime II, with coalescence occurring only for intermediate absorption rates when the Marangoni number is above a critical value near 8.28. If correct, this gives a mechanism by which ambient humidity alone can control droplet size, spacing, and collision behaviour in fibre-based dehumidification and water-harvesting devices.","feed_headline":"Vapor absorption can force film droplets on a vertical fibre to merge","feed_subtitle":"The model predicts a coalescence window bounded by absorption rates 0.008 and 0.32, once the Marangoni number exceeds 8.28.","key_machinery":"The carrying object is the coupled PDE system (3.9): a thickness equation with non-conservative source $\\Lambda(1+\\alpha h)(c-c_s)$, and a concentration equation of advection-diffusion form, linked by the flow rate $q$ that contains gravity, surface tension with destabilizing azimuthal curvature $\\alpha/[\\varsigma(1+\\alpha h)]$ and stabilizing streamwise curvature $-h_{zz}$, and a Marangoni term proportional to $Ma\\, h^2\\psi(\\alpha h)c_z$. The analytical workhorse is the effective linear growth rate (4.8) obtained by freezing a slowly absorbing base state: absorption adds a positive $\\Lambda\\Gamma$ term to the classical Rayleigh-Plateau growth rate and shifts the critical wavenumber. The regime boundaries are read off from this growth-rate structure and from the numerical phase diagram (Figure 12), where the two threshold curves $\\Lambda_{I\\to II}(Ma)$ and $\\Lambda_{II\\to I}(Ma)$ meet at $(Ma,\\Lambda)\\approx(8.28,0.063)$.","core_discovery":"The paper's central discovery is a coupled system of two nonlinear fourth-order PDEs for film thickness $h(z,t)$ and silicone-oil concentration $c(z,t)$, equation (3.9), that extends the classical mass-conserving fibre-coating equation to include water vapor absorption through the flux $J(c)=c-c_s$ and Marangoni effects through the concentration gradient. For weak absorption, the paper derives a frozen-time linear stability result showing that absorption enlarges the unstable wavenumber band and raises the effective growth rate of interfacial perturbations. Numerical simulations with realistic inlet conditions then show that, for sufficiently strong Marangoni effects, increasing the absorption parameter moves the film from a non-coalescing Regime I into a coalescing Regime II, and back into Regime I at still larger values. The paper also derives a quasi-static logistic concentration profile and an approximate droplet-mass formula that captures the non-monotone dependence of total liquid mass on absorption rate.","pith_inferences":["Beyond the paper: if the linear absorption law is replaced by a thickness-dependent or Langmuir-type kinetics, the predicted coalescence window would likely shift and could widen or close; measuring the absorption flux of the actual sorbent as a function of film thickness would tell which.","Beyond the paper: because the regime transition is controlled by humidity through the absorption parameter, a fibre dehumidifier could in principle switch between coalescing and non-coalescing modes by adjusting ambient vapour pressure; the phase diagram gives the control margins.","Beyond the paper: the prediction that no coalescence occurs for Marangoni numbers below about 8.28 suggests a simple experimental check: vary the initial oil concentration (which changes the surface-tension difference) at fixed humidity, and see whether coalescence turns on only above a threshold concentration difference.","Beyond the paper: the same coupled thickness-concentration structure may apply to other non-conservative fibre coatings, such as reactive or evaporating films, where the source term in the thickness equation has a different functional form; the regime-boundary machinery is a testable generic prediction."],"forward_implications":["At fixed Marangoni number above the critical value, there is an intermediate absorption window in which droplet trains cannot remain in the non-coalescing regime; device operation would need to avoid or deliberately enter that window depending on the purpose.","Stronger absorption within the coalescing regime moves the collision point upstream toward the nozzle, so the location of coalescence can be used as a readout of absorption rate.","Below the critical Marangoni number (about 8.28), absorption alone is not enough to trigger coalescence; both concentration-gradient forcing and non-conservative growth are needed.","The total liquid mass held in a fixed downstream window is non-monotone in the absorption rate, and the paper's parabolic droplet approximation gives a formula that reproduces this trend, connecting droplet spacing and height to mass capture.","The quasi-static logistic concentration profile predicts where along the fibre saturation sets in and where droplet compression begins, so it furnishes a design relation between humidity, sorbent concentration, and droplet spacing."],"supporting_citations":[{"why":"Supplies the vertical-fibre lubrication model, the scaling $\\epsilon = We^{-1/3}$, and the mass-conserving base equation that the new model extends.","marker":"Ji et al. (2019)"},{"why":"Provides the classical viscous beads-on-fibre model and the treatment of destabilizing azimuthal curvature with retained streamwise curvature.","marker":"Craster & Matar (2006)"},{"why":"Gives the drop-formation equation recovered in the $\\alpha h\\to 0$ limit and frames the Rayleigh-Plateau instability context.","marker":"Kalliadasis & Chang (1994)"},{"why":"One-sided evaporating/condensing film model whose frozen-time stability approach is adopted for weak absorption.","marker":"Burelbach et al. (1988)"},{"why":"Frozen-time quasi-static base-state stability method used in the weak-absorption linear analysis.","marker":"Shklyaev & Fried (2007)"},{"why":"Thermally-driven coalescence on a fibre; provides the boundary-condition setup and the Regime I versus Regime II comparison for coalescence.","marker":"Ji et al. (2021)"},{"why":"Fibre-film modelling and regime classification; used for the Nusselt-thickness parameter estimation.","marker":"Ruyer-Quil et al. (2008)"},{"why":"Supplies physical properties of the water-absorbing silicone oil, including the saturated water uptake that sets $c_s\\approx 0.86$.","marker":"Ahn et al. (2017)"}],"fun_headline_variants":["Absorption rate tunes droplet merging on vertical fibres","Vapor uptake drives coalescence in fibre film flows","Marangoni and absorption switch film droplet regimes","Absorption windows control droplet coalescence on fibres","Thin film on fibre: absorption flips coalescence on/off"],"cache_read_input_tokens":29056,"weakest_assumption_plain":"The model assumes the absorption rate into the film is simply proportional to how far the oil concentration is above its saturation value, with a fixed proportionality constant; if real sorption is slower, nonlinear, or depends on film thickness, the predicted coalescence window and regime boundaries change.","fun_headline_variants_meta":{"raw":{"variants":["Absorption rate tunes droplet merging on vertical fibres","Vapor uptake drives coalescence in fibre film flows","Marangoni and absorption switch film droplet regimes","Absorption windows control droplet coalescence on fibres","Thin film on fibre: absorption flips coalescence on/off"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000179,"raw_usage":{"total_tokens":1293,"prompt_tokens":933,"completion_tokens":360,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":549,"completion_tokens_details":{"reasoning_tokens":283}},"tokens_in":549,"tokens_out":360,"duration_ms":4347,"temperature":1.0,"reasoning_tokens":283,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T13:09:33.456684+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a controlled fibre experiment with a water-absorbing silicone oil at fixed fibre radius, flow rate, and inlet concentration, varying only the surrounding humidity (which sets the absorption parameter), and track droplet peak positions over time. The model predicts that at a Marangoni number of 20 the droplet train remains non-coalescing for absorption parameters below about 0.0083 and above about 0.32, but coalesces in between; observing no coalescence across a continuous humidity sweep, or coalescence at all humidities, would falsify the regime diagram. A direct check of the constitutive assumption is also possible: measure the absorption flux as a function of concentration difference for films of different thickness, since the predicted regime structure depends on the flux being linear.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the classical viscous beads-on-fibre model and the treatment of destabilizing azimuthal curvature with retained streamwise curvature."},{"cited_title":"Journal of Fluid Mechanics 261 , 135--168","cited_arxiv_id":null,"evidence_quote":"Gives the drop-formation equation recovered in the $\\alpha h\\to 0$ limit and frames the Rayleigh-Plateau instability context."},{"cited_title":"Journal of Fluid Mechanics 195 , 463--494","cited_arxiv_id":null,"evidence_quote":"One-sided evaporating/condensing film model whose frozen-time stability approach is adopted for weak absorption."},{"cited_title":"Journal of Fluid Mechanics 584 , 157--183","cited_arxiv_id":null,"evidence_quote":"Frozen-time quasi-static base-state stability method used in the weak-absorption linear analysis."},{"cited_title":"Journal of Fluid Mechanics 603 , 431--462","cited_arxiv_id":null,"evidence_quote":"Fibre-film modelling and regime classification; used for the Nusselt-thickness parameter estimation."},{"cited_title":"US Patent 9,731,245","cited_arxiv_id":null,"evidence_quote":"Supplies physical properties of the water-absorbing silicone oil, including the saturated water uptake that sets $c_s\\approx 0.86$."}],"review_version":1}