{"id":"fd545352-aa73-4b8b-a273-273798c4842f","arxiv_id":"2507.05947","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A generalized Hadamard-Rybczynski terminal-velocity formula is derived for a spherical drop with partial slip at the liquid-liquid interface, recovering the classic no-slip result and the slip-modified Stokes law as limits.","lead":"A theoretical paper extends a classic 1911 formula for how fast a liquid drop falls through another liquid, adding a slip effect at the interface where the two liquids meet. The new formula includes the old one as a special case and could help interpret emulsion experiments where droplets move faster or slower than expected.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (74) is internally consistent, but the Fig. 3 example uses slip lengths (λ=0.108 mm, λ'=-0.159 mm) that violate the paper's own relation Eq. (16); the flow-field illustration is therefore not a solution of the model.","rationale":"The reader's conditional verdict is appropriate, but the specific load-bearing concern I find is not the lack of molecular evidence for the slip law; it is a concrete numerical inconsistency in the paper's own illustration. Eq. (74) is a correct algebraic consequence of the stated boundary conditions, and its HRE and solid-slip limits are correct. The reported slip lengths in Fig. 3, however, cannot both come from Eqs. (93)-(94) because they violate Eq. (16) by more than an order of magnitude. This directly affects the paper's experimental interpretation and its claim that the continuous-stress model is a special case of the generalized Navier model. Because the flaw is confined to the worked example rather than to the general formula, I would keep the reader's CONDITIONAL verdict rather than moving to ACCEPT or REJECT.","tokens_in":26207,"tokens_out":35927,"duration_ms":391376,"concrete_test":"Recompute λ and λ' from Eqs. (93)-(94) for the castor/silicate parameters and R=1 mm, then check whether the pair satisfies λ'/λ = -η'/η. Regenerate Fig. 3(b) using a self-consistent pair, e.g., λ'=-η'λ/η with the reported λ. If the reported pair fails Eq. (16), the plotted internal velocity field is not part of the generalized Navier solution and the figure must be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The general derivation of Eq. (74) is standard and reduces correctly to HRE and to the slip-modified Stokes limit, so I do not object to the central formula as a mathematical result. The load-bearing problem is in the worked example used to interpret the experiment. In Fig. 3 the paper reports λ=0.108 mm (external) and λ'=-0.159 mm (internal) for castor oil (η=0.693 Pa·s), silicate oil (η'=0.0232 Pa·s), and R=1 mm. These values violate the paper's own consistency condition Eq. (16), η/λ = -η'/λ'. Starting from λ=0.108 mm, Eq. (16) requires λ'=-0.0036 mm, not -0.159 mm. Starting from λ'=-0.159 mm, Eq. (16) requires λ=4.75 mm, not 0.108 mm. Consequently, the velocity profiles in Fig. 3(b), and the conclusion that the internal droplet velocity is an order of magnitude larger in the partial-slip model, are not generated by the generalized Navier model of Eqs. (14)-(16). This is a purely arithmetic defect that can be fixed, but as published it undermines the experimental illustration and the claim that Eqs. (93)-(94) realize the continuous-stress model. The central Eq. (74) itself is not invalidated.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives an analytical solution for the slow motion of a spherical liquid drop in another immiscible liquid under a generalized Navier partial-slip boundary condition at the liquid-liquid interface. The main result, Eq. (74), is a generalized Hadamard-Rybczynski equation that reduces to the classical HRE for zero slip length and to the slip-modified Stokes law for a solid sphere in the infinite-drop-viscosity limit. The paper also proposes a specific 'continuous viscous stress tensor' model in Appendix A, which yields closed-form slip lengths (93)-(94) and a terminal velocity formula (95), and applies this model to published experimental data for silicate-oil drops falling in castor oil, presenting velocity profiles and streamlines. The central formula (74) appears mathematically consistent and correctly reproduces the stated limits; however, the special model and the experimental application contain serious internal inconsistencies that undermine the paper's secondary claims.","tokens_in":26557,"tokens_out":10132,"duration_ms":100757,"significance":"If Eq. (74) were the only contribution, the paper would constitute a useful and conceptually interesting extension of the Hadamard-Rybczynski theory, with the strength of an explicit solution and clear asymptotic reductions. The idea of introducing a liquid-liquid slip length and linking it to the Boussinesq surface-viscosity model is appealing. However, the paper's additional claims—the continuous-stress-tensor closure, the specific slip-length formulas, and the experimental interpretation—are currently not sound because the derived slip lengths violate the consistency relation (16) that the same model establishes. These errors are load-bearing for the experimental sections and for the claim that a particular partial-slip model explains deviations from HRE. The central generalization itself is not invalidated, but the manuscript as a whole requires major revision before the secondary claims can be accepted.","major_comments":[{"comment":"The slip lengths derived in the Appendix do not satisfy the consistency condition (16) that the main model imposes. From λ = R(1-η'/η)/3 and λ' = R(1-η/η')/3, one obtains η/λ = 3η²/[R(η-η')] and -η'/λ' = 3η'²/[R(η-η')], which are equal only when η=η'. The statement in Section 6 that 'obviously, Eqs. (93) and (94) satisfy the condition given by Eq. (16)' is therefore false. More seriously, substituting λ from Eq. (93) into the main formula (74) gives V₀ = (8/15)(ρ-ρ')gR²/η, whereas the Appendix's Eq. (95) gives V₀ = (4/15)(ρ-ρ')gR²/η; these results differ by a factor of two. The Appendix model is thus not a special case of the generalized HRE, contradicting the paper's central narrative.","section":"Appendix A, Eqs. (93)-(94) and Eq. (16)"},{"comment":"The numerical values used in Fig. 3, λ = 0.108 mm and λ' = -0.159 mm for R = 1 mm, η = 0.693 Pa·s, η' = 0.0232 Pa·s, violate the model's own consistency relation. Using Eq. (16), η/λ ≈ 6417 Pa·s/m while -η'/λ' ≈ 146 Pa·s/m; starting from λ = 0.108 mm requires λ' = -0.0036 mm, and starting from λ' = -0.159 mm requires λ = 4.75 mm. Moreover, these values do not even follow from Eqs. (93)-(94), which for the stated parameters give λ ≈ 0.322 mm and λ' ≈ -9.62 mm. Consequently, the velocity profiles in Fig. 3(b) and the streamlines in Fig. 4(b) are not solutions of the generalized Navier model, and the reported order-of-magnitude enhancement of the internal droplet velocity is not a prediction of the model.","section":"Fig. 3 and Eq. (16)"},{"comment":"The additional boundary conditions (91)-(92), which require continuity of the diagonal stress components σ_θθ and σ_φφ, are not physically appropriate for a fluid-fluid interface. For two immiscible fluids, the correct interfacial conditions are continuity of the traction vector (with the normal-stress jump balanced by surface tension) and, in the absence of surfactant gradients, continuity of tangential stress. There is no physical mechanism requiring σ_θθ and σ_φφ to be continuous across the interface. This over-constraint is the source of the unphysical behavior seen in the Appendix, including a terminal velocity independent of the internal viscosity η' (Eq. (95)) and slip lengths that diverge as η'/η→0 or η/η'→0. This is a modeling error rather than a mere presentation issue.","section":"Appendix A, boundary conditions (91)-(92)"},{"comment":"The experimental comparison in Fig. 2 uses the Appendix's Eq. (95) as the 'partial slip model' curve. Since Eq. (95) is not a consequence of the generalized HRE (74) with any slip length that satisfies Eq. (16), the orange curve in Fig. 2 does not represent the partial-slip model that the paper proposes. The Discussion's claim that the partial-slip model gives a deviation from the experiment comparable to HRE is therefore unsupported. The experimental evidence presented does not validate the generalized HRE or the continuous-stress closure; at most, it illustrates a family of curves that could be obtained if λ were fit freely, but no such fit is performed.","section":"Fig. 2 and Discussion"}],"minor_comments":[{"comment":"The algebraic forms of Eqs. (72)-(75) are hard to read because of spacing and missing parentheses; please re-typeset these formulas with clear denominators and brackets, and verify that Eq. (74) is exactly equivalent to Eq. (73) after substituting α = η'/η.","section":"Eqs. (72)-(75) and text around them"},{"comment":"There is a typo: 'on can obtain' should be 'one can obtain'. Also, in Eq. (86), the drag coefficient expression should be checked for dimensional consistency, as C_D = β/(π ρ R² V₀) appears to have units of Pa·s/(kg·m⁻³·m·m/s) = 1, but the stated formula may have an error in the prefactor.","section":"Sec. 5, after Eq. (81)"},{"comment":"The reference to the Navier boundary condition is numbered inconsistently: [5-6] in the introduction but [14-15] in the Introduction and Section 6. Please standardize the citation numbering throughout.","section":"References"},{"comment":"The caption contains a typo: 'liquids-liquid interface' should be 'liquid-liquid interface'.","section":"Fig. 2 caption"},{"comment":"The rearrangement from Eq. (A3) to Eq. (A4) is not immediately transparent; the derivation of Eq. (A5) appears to contain typographical errors in the Legendre-function identities and should be carefully checked.","section":"Appendix A, Eq. (A5)"}],"recommendation":"major_revision","confidential_remarks":"The main formula (74) is a valuable generalization and the derivation up to that point is largely sound. However, the Appendix's special model is internally inconsistent with the main model's own relation (16), and the numerical example in Fig. 3 uses slip lengths that satisfy neither Eq. (16) nor the Appendix's formulas. These are not cosmetic errors; they invalidate the paper's claimed special case and its experimental interpretation. The author should be asked to either remove the continuous-stress-tensor model and the associated experimental claims, or re-derive them as a genuine special case of the generalized HRE with a consistent closure. Given the centrality of these issues to the paper's claimed novelty, major revision is appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The formula in Eq. (74) is worth knowing. The algebra is straightforward but the reduction to the classical HRE at λ=0 and to the slip-modified Stokes drag for a solid sphere in the infinite-viscosity limit is done carefully, and the result is a compact one-parameter extension that could be useful for emulsion sedimentation problems. Compared to the Boussinesq surface-viscosity model, the slip-length picture is at least as natural, and the paper makes that point without overselling it. The derivation itself is standard l=1 Stokes flow, and the checks against known limits are the right checks. That part is solid.\n\nThe soft spots are real and load-bearing for the paper's own claims. First, the numerical values in Fig. 3 do not satisfy the paper's own relation Eq. (16). For castor oil and silicate oil, Eq. (16) demands λ' = -0.0036 mm if λ = 0.108 mm, or λ = 4.75 mm if λ' = -0.159 mm. The values quoted in the text, λ = 0.108 mm and λ' = -0.159 mm, violate that condition by a factor of roughly 40. So the velocity profiles and the claim of an order-of-magnitude internal flow are not solutions of the model they are supposed to illustrate. This is an arithmetic defect, but it sits in the central example of the paper.\n\nSecond, the special stress-continuity model in the Appendix is questionable. Imposing continuity of all viscous stress components leads to a terminal velocity (Eq. A42) that is independent of the internal viscosity and does not reduce to Stokes drag for an infinitely viscous droplet. The paper never flags this missing limit. Since that model is the one used for the experimental comparison in Fig. 2, the comparison is not on firm ground.\n\nThird, the physical justification for a constant Navier slip length at a liquid-liquid interface is thin. The paper is honest that this is a model, but it also implies that emulsions should exhibit slip without offering any direct molecular or experimental evidence. That is a softer concern than the arithmetic one, but it limits the paper's reach. The experimental data are adapted without error bars, so the agreement in Fig. 2 is suggestive at best.\n\nWho gets value from this? A researcher working on low-Reynolds-number multiphase flow who wants a compact slip-modified drop velocity formula. The formula can be cited as a mathematical result, though the experimental and physical conclusions need correction and external validation.\n\nRecommendation: deserve a serious referee, not a desk reject, but it needs major revision before publication. The arithmetic inconsistency and the missing Stokes limit in the Appendix model should be fixed or explicitly discussed, and the experimental claims toned down.","headline":"The generalized Hadamard-Rybczynski formula (Eq. 74) is a legitimate and clean mathematical extension, but the paper's own experimental illustration uses slip lengths that violate its consistency condition, and the special stress-continuity model lacks a Stokes limit.","tokens_in":27035,"tokens_out":2824,"would_cite":true,"duration_ms":34375,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["47.15.G-","47.55.Dr"],"model":"deepseek-v4-flash","headline":"The paper derives a generalized Hadamard–Rybczynski equation in which a single slip length at the liquid–liquid interface controls the terminal velocity of a spherical drop, recovering the classical no-slip result at zero slip and the…","keywords":["Hadamard–Rybczynski equation","partial slip","Navier boundary condition","liquid–liquid interface","Stokes flow","slip length","drop sedimentation","emulsions"],"falsifier":"Measure terminal velocities for drops of liquid A in liquid B and liquid B in liquid A at the same radius and extract the two slip lengths from Eqs. (109)–(110); if $\\eta_A/\\lambda_A + \\eta_B/\\lambda_B$ is not zero within experimental error, the generalized Navier condition with linked slip lengths is falsified for that pair.","tokens_in":25944,"feed_emoji":"💧","tokens_out":10408,"duration_ms":67314,"temperature":0.7,"pith_summary":"The paper tries to establish that the slow steady motion of a spherical liquid drop in another immiscible liquid is governed by a one-parameter generalization of the Hadamard–Rybczynski equation, in which the no-slip condition is replaced by a linear partial-slip law at the liquid–liquid interface. The extra parameter is a slip length $\\lambda$, and the terminal velocity takes the closed form of Eq. (74). This matters because the standard HRE often deviates from emulsion experiments, and the usual explanations invoke surfactants or interfacial viscosity; the paper offers a purely hydrodynamic alternative that can be tested by measuring slip lengths directly. If the equation is right, the terminal velocity depends on $\\lambda$ in addition to the two viscosities, densities, radius, and gravity, with the classical HRE and Stokes law recovered as limiting cases.","feed_headline":"One extra parameter generalizes the classic drop-speed equation","feed_subtitle":"One slip length links the no-slip drop law, partial-slip drops, and solid-sphere drag; a swap experiment can measure it.","key_machinery":"The load-bearing object is the generalized Navier partial-slip boundary condition at the liquid–liquid interface, Eqs. (14)–(15): the tangential viscous stress of each liquid is proportional to the relative tangential velocity between the two fluids across the interface, with slip lengths $\\lambda$ (external) and $\\lambda'$ (internal) linked by $\\eta/\\lambda = -\\eta'/\\lambda'$ through stress continuity. This replaces no-slip with a linear friction law. The derivation then uses the standard $l=1$ axisymmetric Stokes solution in spherical coordinates, enforces vanishing radial velocity at the drop surface, shear-stress continuity, and the force balance with the Archimedes force, which fixes all coefficients and yields Eq. (74). In the appendix, requiring continuity of all components of the viscous stress tensor fixes the slip length in terms of the viscosities and radius, producing the special parameter-free formula.","core_discovery":"The central claim is that partial slip at a liquid–liquid interface, written as a generalized Navier condition, produces a closed-form terminal velocity $$V_0 = \\frac{2(\\rho-\\rho')$gR^{2}$}{3\\eta}\\,\\frac{\\$\\alpha$ + 1 + 3\\$\\lambda$/R}{2\\$\\alpha$ + 3 + 6\\$\\lambda$/R},\\tag{74}$$ with $\\alpha=\\eta'/\\eta$, which reduces to the Hadamard–Rybczynski result when $\\lambda=0$ and to the slip-modified Stokes drag law for a solid sphere when the droplet viscosity is infinite. Choosing the slip length $\\lambda = \\frac{R}{3}(\\eta/\\eta'-1)$ enforces continuity of all components of the viscous stress tensor at the interface and gives the parameter-free speed $V_0 = 4(\\rho-\\rho')gR^2/(15\\eta)$, independent of the internal viscosity. The paper applies the model to silicate-oil drops falling in castor oil: the experimental points sit between the HRE curve and the stress-continuous partial-slip curve, and the two models differ by an order of magnitude in the predicted speed of the internal toroidal flow.","pith_inferences":["Beyond the paper: the cleanest test of the generalized Navier law is a micro-particle-image-velocimetry measurement of the tangential velocity jump at a settling drop's interface; if the jump is absent while the terminal velocity still deviates from HRE, the missing physics is not linear partial slip.","Beyond the paper: the structural similarity between the generalized HRE and the Boussinesq surface-viscosity formula suggests that a measured slip length may absorb surfactant and Marangoni effects; comparing $\\lambda$ extracted from sedimentation with independent surface-viscosity measurements would probe whether slip is a genuine material property.","Beyond the paper: because $\\lambda$ enters through $\\lambda/R$, measuring the terminal velocity of the same liquid pair over a range of drop radii would give a direct radius-dependence signature; a constant $\\lambda$ would support the model, while $\\lambda$ growing with $R$ would point to a surface-viscosity interpretation instead."],"forward_implications":["At $\\lambda=0$ the generalized equation reduces exactly to the Hadamard–Rybczynski equation, so the classical result is contained as a special case.","For a droplet of infinite viscosity the generalized equation becomes the known partial-slip correction to Stokes drag for a solid sphere; at $\\lambda=0$ this is the ordinary Stokes law.","If the slip mechanism is real, swapping the roles of the two liquids gives two measured velocities from which slip lengths can be extracted, and the model requires $\\eta_A/\\lambda_A = -\\eta_B/\\lambda_B$ for the same pair.","The stress-continuous version predicts a drop terminal velocity $4(\\rho-\\rho')gR^2/(15\\eta)$ that is independent of the droplet viscosity and exceeds the solid-sphere Stokes value by a factor of $6/5$, providing a direct quantitative target for experiment.","In the silicate-oil/castor-oil comparison, the partial-slip model and HRE bracket the measured fall speeds; the decisive difference is that the partial-slip model predicts internal droplet circulation roughly an order of magnitude faster than HRE."],"supporting_citations":[{"why":"defines the Hadamard–Rybczynski equation that this paper generalizes.","marker":"[1-2]"},{"why":"supplies the standard treatment of drop motion, stress boundary conditions, and the surfactant-based explanation for deviations from HRE.","marker":"[3]"},{"why":"is the reference for the no-slip boundary condition and its partial-slip generalizations for solid surfaces.","marker":"[4]"},{"why":"documents the experimental basis of the Navier partial-slip condition at solid–liquid interfaces that motivates the liquid–liquid extension.","marker":"[5-6]"},{"why":"gives the previously obtained slip-modified Stokes drag law for a solid sphere, which the generalized HRE must reproduce in the infinite-viscosity limit.","marker":"[9]"},{"why":"provides the axisymmetric Stokes-flow solution in spherical coordinates used to solve the boundary-value problem.","marker":"[15-16]"},{"why":"introduces the Boussinesq surface-viscosity model whose structure parallels the generalized HRE.","marker":"[21]"},{"why":"provides the experimental falling-drop data for silicate oil in castor oil used for comparison.","marker":"[22]"}],"fun_headline_variants":["Drop speed law gets a slip-length generalization","Slip length extends Hadamard–Rybczynski to partial-slip drops","One slip length tunes drop speed from no-slip to stress-free","A slip length turns Hadamard–Rybczynski into Stokes drag for solids","Generalized drop speed matches experiment with measured slip length"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the liquid–liquid interface obeys a linear partial-slip law with constant slip lengths, linked by $\\eta/\\lambda = -\\eta'/\\lambda'$, and with no surfactant effects, surface viscosity, or dependence of slip on the local flow.","fun_headline_variants_meta":{"raw":{"variants":["Drop speed law gets a slip-length generalization","Slip length extends Hadamard–Rybczynski to partial-slip drops","One slip length tunes drop speed from no-slip to stress-free","A slip length turns Hadamard–Rybczynski into Stokes drag for solids","Generalized drop speed matches experiment with measured slip length"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000882,"raw_usage":{"total_tokens":3882,"prompt_tokens":1089,"completion_tokens":2793,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":705,"completion_tokens_details":{"reasoning_tokens":2702}},"tokens_in":705,"tokens_out":2793,"duration_ms":22438,"temperature":1.0,"reasoning_tokens":2702,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T19:15:31.162381+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure terminal velocities for drops of liquid A in liquid B and liquid B in liquid A at the same radius and extract the two slip lengths from Eqs. (109)–(110); if $\\eta_A/\\lambda_A + \\eta_B/\\lambda_B$ is not zero within experimental error, the generalized Navier condition with linked slip lengths is falsified for that pair.","supporting_citations":[{"cited_title":"(66) Stream function is (Table 1):   12 24 2( ') 6 2, sin 36 ' gR r r R rr R R R r R                                         ","cited_arxiv_id":null,"evidence_quote":"supplies the standard treatment of drop motion, stress boundary conditions, and the surfactant-based explanation for deviations from HRE."},{"cited_title":"(70) Stream function is (Table 1): 1 4 24 2( ') 6 2' sin 36 ' ' gR r r R R R                               ","cited_arxiv_id":null,"evidence_quote":"is the reference for the no-slip boundary condition and its partial-slip generalizations for solid surfaces."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"gives the previously obtained slip-modified Stokes drag law for a solid sphere, which the generalized HRE must reproduce in the infinite-viscosity limit."},{"cited_title":"The Fall of Mercury Droplets in a Viscous Medium","cited_arxiv_id":null,"evidence_quote":"introduces the Boussinesq surface-viscosity model whose structure parallels the generalized HRE."},{"cited_title":"Boinovich","cited_arxiv_id":null,"evidence_quote":"provides the experimental falling-drop data for silicate oil in castor oil used for comparison."}],"review_version":1}