{"id":"590d9400-2721-4880-b0da-011fa0a2f5c5","arxiv_id":"2607.02730","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":5,"one_line_summary":"Extending single-drop lamella theory to simultaneous pair impact yields closed-form sheet fields in the lamella-fed stage, an empirical post-lamella continuation, and a Rayleigh–Plateau height bound that together explain Weber-number crossover and pinch-off.","lead":"A semi-analytical model predicts the shape, thickness, and height of the free-standing liquid sheet that rises when two drops hit a dry surface at once. The model explains why reported Weber-number scalings differ and bounds the height when the sheet pinches off, which matters for spray coating and cooling.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"No significant objection identified beyond the already-flagged empirical post-lamella inlet.","rationale":"The paper’s strongest claim—that the framework supplies the full 3-D sheet fields, an explicit centreline height that explains the literature’s differing We exponents as a crossover, and a Rayleigh–Plateau bound at high We—rests on two stages. The lamella-fed stage is closed-form once the single-drop solution is granted; the post-lamella stage and the subsequent RP cut-off inherit the empirical inlet of Appendix B. That is precisely the weakest assumption identified by the Reader. All other modelling choices (ballistic characteristics, free-rim balances with α=β=1, neglect of gravity during the rise, quasi-steady RP growth) are either standard or shown to be consistent with the data inside the rising phase. Because the paper already flags the limited range of the fit and supplies independent validation inside that range, the concern does not warrant a change of verdict. The recommended concrete test simply quantifies how far the present inlet can be extrapolated, which is the natural next check already implied by the Reader’s CONDITIONAL assessment.","tokens_in":21101,"tokens_out":589,"duration_ms":5973,"concrete_test":"Extract the post-lamella inlet (T and ū_PL along the collision line) from an independent DNS at a point outside the present calibration window (e.g. We=100, a=1.5 or Oh=0.005) and recompute H_m from Eqs. (28)–(30) and H_RP from Eqs. (31)–(33). If the new inlet changes either height by more than ~20 % relative to the original fit, the quantitative reach of the post-lamella and pinch-off claims is limited exactly as the Reader states.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The Reader correctly isolates the load-bearing soft spot: after t_ℓ the inlet is prescribed by empirical forms (ū_PL = C/t^ζ with ζ=1.2 fixed and a cubic T_P(t) whose four coefficients are taken from the same DNS suite; Appendix B and Eq. (1)). Those forms are stated to be valid only inside the calibrated (We,a) window. Outside that window the post-lamella height prediction (Eqs. 28–30) and the Rayleigh–Plateau cut-off (Eqs. 31–33) lose their quantitative basis. Inside the window the paper already supplies independent experimental and DNS support for both the full sheet shape and the centreline H(t), and the lamella-fed stage itself remains essentially parameter-free once the Gordillo et al. single-drop solution is accepted. No deeper internal inconsistency or hidden assumption that would overturn the central claims was found.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The paper develops a semi-analytical model for the free-standing vertical sheet formed when two low-Oh drops impact a dry substrate simultaneously. Extending the single-drop lamella theory of Gordillo et al. (2019), the authors couple thin-film sheet flow to colliding-lamella inlet conditions and a capillary-retarded rim, then solve the sheet interior along ballistic characteristics in a lamella-fed (L) stage (closed-form velocity and thickness) and a post-lamella (PL) stage (empirical inlet from DNS). The framework yields three-dimensional sheet fields and the full side-view shape; on the centreline it supplies explicit apex height and thickness, attributing literature scatter in We exponents to an L/PL crossover rather than a single power law. At high We a linear Rayleigh–Plateau analysis, driven by the model’s time-dependent jet diameter and deceleration, bounds the maximum height and closes the pinch-off regime. Predictions are compared with Goswami–Hardalupas experiments and the authors’ Basilisk DNS over We, a and modest Oh.","tokens_in":21357,"tokens_out":1779,"duration_ms":20745,"significance":"If the results hold, the work supplies the first predictive, largely first-principles description of the unsteady central sheet that has been the defining feature of simultaneous two-drop impact since Barnes et al. The L-stage construction is essentially parameter-free once the Gordillo lamella solution is accepted, yields closed-form centreline fields (Eqs. 20–25), and cleanly explains why reported We exponents differ. The full side-view profiles (Figs. 3–4) and H(t) trajectories (Fig. 5) match independent experiments as well as DNS, and the RP cut-off (Eqs. 31–33) is a falsifiable, literature-threshold closure with no additional tuning. These strengths—analytic L-stage fields, explicit crossover interpretation, and a parameter-light capillary bound—make the paper a useful advance for multi-drop impact, spray cooling and related applications, even though the PL inlet remains empirical.","major_comments":[{"comment":"Appendix B and Eq. (1): the post-lamella inlet (ū_PL = C/t^ζ with ζ = 1.2 fixed, and the four-coefficient cubic T_P(t)) is extracted from the same Basilisk DNS suite later used for validation of H_m and sheet shape in the L+PL regime (Figs. 5–7). Consequently part of the ‘prediction’ of H_m via Eqs. (28)–(30) and of H_RP via Eqs. (31)–(33) is conditioned on those data. The manuscript should (i) state explicitly which comparisons are independent of the PL fit (experiments of Goswami & Hardalupas; pure L-stage cases) versus which reuse the calibration DNS, (ii) report a sensitivity study of H_m and H_RP to modest variations of ζ and the T_P coefficients, and (iii) if possible hold out at least one (We, a) case from the fit and show it a posteriori. Without this, the quantitative claim for the PL and pinch-off regimes is weaker than the abstract suggests.","section":"Appendix B, Eq. (1), §§IV B–V"},{"comment":"§IV A, Eqs. (22)–(25) and §IV B, Eqs. (28)–(30): the central claim that ‘different Weber-number exponents reported in the literature arise from a crossover rather than from a single universal scaling law’ is attractive but only weakly demonstrated. The L-stage asymptote H_m ∼ We^{1/(1+ϕ_L)} is given, yet the paper never overlays the literature exponents (or the empirical scalings of Goswami & Hardalupas and Zhang et al.) on Fig. 6 or 7, nor does it show where each dataset sits relative to t_m ≶ t_ℓ,0. A short quantitative comparison—listing the published exponents against the model’s local effective slope d ln H_m / d ln We in the L-only, crossover and PL-dominated windows—would make the claim load-bearing rather than interpretive.","section":"§IV, Figs. 6–7"},{"comment":"§V A, Eq. (32): the Rayleigh–Plateau cut-off adopts ln(r_jet/ε_0) ≈ 12 as a ‘standard value’. While the range 8–15 is cited, H_RP is exponentially sensitive to this threshold. The manuscript should either (a) show that H_RP remains within the experimental scatter for the conventional interval 8–15, or (b) calibrate the constant once against a single high-We pinch-off event and then predict the remaining cases. As written, the cut-off is not fully parameter-free and the high-We bound in Fig. 6a is less robust than claimed.","section":"§V A, Eq. (32), Fig. 6"}],"minor_comments":[{"comment":"§II and Appendix A: ϕ_L is introduced as a slowly varying function of (r,t) yet is thereafter treated as a single constant read from numerics for each (We, Oh). A one-sentence statement of the relative variation of ϕ across the lamella at fixed Oh (the text already says <10% in We) would clarify how much error the constant-ϕ_L approximation introduces into Eqs. (20) and (25).","section":"§II, Appendix A"},{"comment":"Figure 2c: the regime map is drawn only for water. Adding a second Oh contour (or a brief remark that the L/PL boundary t_ℓ,0 = a²/3 is Oh-independent while s_max is not) would help readers place the glycerol–water cases of Appendix D.","section":"Fig. 2c"},{"comment":"§III A: gravity is neglected on the basis of We^{1/2} Fr^{-2} ≲ 0.1, yet late-time descent is visible in both experiment and DNS (Figs. 3–5). A short estimate of the cumulative gravitational deceleration over t ∼ 1–5 would quantify when the ballistic assumption begins to fail and would justify the rising-phase-only comparisons.","section":"§III A"},{"comment":"Notation: lower-case (ū, h, s, b) for lamella quantities and upper-case (W̄, T, H, B) for sheet quantities is helpful, but T_0, W̄_0, V̄_0 in Eq. (1) and T_b, T_r later are easy to confuse with the sheet thickness field T(y,z,t). A compact symbol table in §III A would reduce cognitive load.","section":"§III A, Eq. (1)"},{"comment":"Appendix D, Fig. 9: the Oh = 0.0141 comparison is valuable; stating explicitly that the same PL inlet coefficients (no re-fit) were used would strengthen the claim that the model is not retuned outside the water window.","section":"Appendix D"},{"comment":"Typos / style: ‘§,III’ (p. 3) should be ‘§III’; ‘theycomponents’ / ‘they–zplane’ missing spaces appear in several places; ‘DeepSeek for assistance with spell-checking’ in the acknowledgements is unusual for a journal and may be better omitted or rephrased.","section":"Throughout / Acknowledgements"}],"recommendation":"minor_revision","confidential_remarks":"The paper is a solid JFM-level contribution. The empirical PL inlet is the only real soft spot and the authors already flag it in §VI; requiring a sensitivity/hold-out analysis and a clearer literature-exponent comparison should be enough for acceptance. I do not see a load-bearing error that would justify major revision or reject. Fit to the journal is good."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"This is a solid, usable piece of multi-drop impact modelling. The real advance is the extension of the Gordillo–Riboux single-drop lamella solution to the colliding-lamella geometry: ballistic characteristics give closed-form velocity and thickness fields in the lamella-fed stage, an algebraic centreline stall condition that recovers H_m ~ We^{1/(1+ϕ_L)}, and a clean demonstration that the scatter in reported We exponents is just the L-to-PL crossover rather than competing universal laws. The Rayleigh–Plateau cut-off that bounds the high-We pinch-off regime is parameter-free once the model supplies the jet diameter and deceleration. Side-view profiles and H(t) match both the authors’ Basilisk DNS and the Goswami–Hardalupas experiments across We and a (Figs. 3–5), and the rim balances are standard free-sheet balances. That is real progress for spray-cooling and coating people who need a reduced description of the central sheet.\n\nThe soft spot is exactly the one the reader flags and the authors themselves state: after t_ℓ the inlet is an empirical fit (ū_PL = C/t^{1.2}, cubic T_P) taken from the same DNS suite later used for validation. Inside the calibrated window the comparisons still hold and the L-stage plus RP cut-off remain independent of that fit; outside it the PL height prediction loses its quantitative basis. Gravity is neglected, so late-time descent is missed, but that is minor for the rising phase they care about. Free parameters are few and openly listed; citation pattern is appropriate.\n\nThis is for anyone working on multi-drop impact or needing a reduced model of the rising sheet. It deserves a serious referee. I would engage with it and cite the L-stage results and the crossover explanation.","headline":"Clean, usable extension of Gordillo–Riboux to the two-drop rising sheet; L-stage is essentially closed-form and the We-exponent scatter is explained as a stage crossover, with the only real soft spot being the openly empirical post-lamella inlet.","tokens_in":21998,"tokens_out":496,"would_cite":true,"duration_ms":5118,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"A semi-analytical model predicts the full shape and height of the free-standing sheet that rises when two drops collide on a dry surface, and shows that reported Weber-number scalings are a crossover, not a single law.","keywords":["droplet-pair impact","rising sheet","lamella collision","ballistic characteristics","Weber number scaling","Rayleigh–Plateau pinch-off","semi-analytical model","thin-film flow"],"falsifier":"Measure maximum sheet height versus Weber number at fixed half-spacing for a liquid whose Ohnesorge number lies inside the low-viscosity range but whose post-lamella inlet velocity and thickness deviate measurably from the fitted power-law and cubic forms; if the measured heights still follow the model’s inertial and Rayleigh–Plateau curves, the empirical inlet is robust, otherwise the post-lamella prediction fails.","tokens_in":21962,"feed_emoji":"💧","tokens_out":1077,"duration_ms":9622,"temperature":0.7,"pith_summary":"When two identical low-viscosity drops hit a dry surface at the same time and close enough together, their spreading films meet and throw up a free-standing vertical sheet bounded by a retracting rim. Existing descriptions of how high that sheet climbs rest on empirical fits or energy balances calibrated to simulations. This paper extends an established single-drop lamella theory to the two-drop geometry, couples the thin-film flow inside the sheet to the colliding base and the capillary rim, and solves the interior along ballistic particle paths in two stages. In the early lamella-fed stage the velocity and thickness fields are closed-form; later they are continued with inlet data taken once from simulations. The resulting centreline formulae for apex height and thickness show that the different Weber-number exponents in the literature are successive regimes of the same continuous description rather than competing universal laws. At high Weber number the same fields feed a Rayleigh–Plateau calculation that caps the height once the apex pinches off. The model therefore supplies a predictive, largely analytical account of the sheet’s shape, rise, and breakup limit over the low-Ohnesorge regime of practical spray and coating flows.","feed_headline":"Two drops throw up a sheet whose height follows from ballistic paths","feed_subtitle":"Closed-form centreline formulae show reported Weber exponents are a crossover, not a universal law","key_machinery":"Ballistic characteristic construction of the free-sheet interior: each fluid parcel injected at the collision line keeps constant velocity while thickness evolves by mass conservation; the construction is algebraic in the lamella-fed stage and continued numerically with an empirical post-lamella inlet, then closed by rim balances and a Rayleigh–Plateau cut-off.","core_discovery":"By matching the colliding single-drop lamellae to a free sheet whose fluid parcels follow ballistic characteristics and whose rim obeys mass–momentum balances, one obtains the three-dimensional velocity and thickness fields of the rising sheet and, on the centreline, explicit algebraic relations for apex height and thickness. Those relations demonstrate that the scattered Weber-number exponents reported experimentally arise from a smooth crossover between a lamella-fed regime and a post-lamella regime rather than from any single power law. A linear Rayleigh–Plateau analysis driven by the model’s own time-dependent jet diameter and deceleration then supplies an upper bound on attainable heigh","pith_inferences":["Because the post-lamella inlet is the sole empirical input, any future closed-form description of the residual bulk convergence after the lamella has passed would render the entire sheet dynamics analytical.","The framework already separates geometry (ballistic paths, rim balances) from feeding; the same separation should apply directly to multi-drop arrays once pairwise collision lines are identified.","The observed crossover in Weber scaling suggests that earlier energy-balance models may have been sampling different stages of the same continuous process rather than capturing distinct physics."],"forward_implications":["Different published Weber exponents for maximum sheet height can be reconciled as successive segments of one continuous centreline solution rather than as competing scalings.","The full three-dimensional sheet shape, not only apex height, becomes available from a single rim integration once the inlet is known.","At high Weber number the Rayleigh–Plateau cut-off supplies a parameter-free upper bound on height once the apex jet diameter and deceleration are taken from the sheet model.","The same characteristic-plus-rim construction can be reused for non-simultaneous or unequal-drop impacts once the appropriate collision-line inlet is supplied."],"fun_headline_variants":["Ballistic model yields sheet shape from colliding drop lamellae","Drop-pair sheet height solved by ballistic characteristics","Weber exponents for apex height arise from regime crossover","Closed-form centreline apex of free sheet from two-drop impact","Model bounds pinch-off height of rising sheet from drop-pair collision"],"cache_read_input_tokens":16512,"weakest_assumption_plain":"After the spreading film no longer reaches the collision line, the velocity and thickness that feed the sheet are taken from simple empirical fits calibrated to the same simulations the model is later compared against; outside that fitted window the post-lamella height and the pinch-off bound lose their quantitative footing.","fun_headline_variants_meta":{"raw":{"variants":["Ballistic model yields sheet shape from colliding drop lamellae","Drop-pair sheet height solved by ballistic characteristics","Weber exponents for apex height arise from regime crossover","Closed-form centreline apex of free sheet from two-drop impact","Model bounds pinch-off height of rising sheet from drop-pair collision"]},"model":"grok-4.5","effort":"low","cost_usd":0.005566,"raw_usage":{"total_tokens":1531,"prompt_tokens":857,"num_sources_used":0,"completion_tokens":86,"cost_in_usd_ticks":55660000,"prompt_tokens_details":{"text_tokens":857,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":588,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":857,"tokens_out":86,"duration_ms":5548,"temperature":1.0,"reasoning_tokens":588,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-12T07:26:08.246067+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Measure maximum sheet height versus Weber number at fixed half-spacing for a liquid whose Ohnesorge number lies inside the low-viscosity range but whose post-lamella inlet velocity and thickness deviate measurably from the fitted power-law and cubic forms; if the measured heights still follow the model’s inertial and Rayleigh–Plateau curves, the empirical inlet is robust, otherwise the post-lamella prediction fails.","supporting_citations":[],"review_version":1}