REVIEW 4 major objections 5 minor 40 references
Inelastic spreading of viscoelastic drops
T0 review · 4 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read Polymer elasticity does not change how far an impacting drop spreads, because the elastic energy stored during spreading is negligible compared with viscous dissipation.
desk verdict A careful experimental study with a plausible energetic criterion, but the missing supplementary derivation and dataset placeholder make the central collapse impossible to audit. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the dimensionless elastic-to-viscous damping ratio $\Gamma$ (Eq. 2), derived from an energy balance in which the impact kinetic energy is split into surface energy, viscous dissipation, and elastic storage. The elastic storage is deliberately overestimated by treating the drop as a Hookean solid in Hertzian contact, giving an elastic term scaling as $G \beta_{\max}^5$; any weaker polymer-constitutive model (Neo-Hookean, Oldroyd-B, FENE-P) would give even less resistance. The kinematic closure $D_{\max} = t_s V_s$ and the leading-order average-velocity closure $V \approx \sqrt{E}$ turn the balance into a closed expression $\beta_{\max}(\Gamma)$, with $\Gamma = \frac{8}{5\pi} \frac{G D_0}{\mu_i V_0 \sqrt{E}}$. The paper also introduces the impact Carreau number $Cu_i = \dot{\gamma}_s / \dot{\gamma}^*_{\mathrm{rheo}}$ to separate Newtonian-adherent Boger fluids from shear-thinning concentrated solutions.
What would settle it
The predicted elastic regime is currently untested: prepare a constant-viscosity elastic fluid with $\Gamma \gtrsim 1$ by the paper's own formula (for example a low-viscosity Boger fluid with $G \approx 2 \mu_i V_0 / D_0$), impact it, and measure $\beta_{\max}$. The paper predicts a measurable reduction below the Newtonian baseline in that regime; if none appears, the criterion fails.
Extended reading notes
Core claim
The central discovery is that the maximum spreading ratio $\beta_{\max} = D_{\max}/D_0$ of Boger-fluid drops (PEO and PAAM, up to 1000 ppm) is indistinguishable from the Newtonian baseline for the same impact parameters, with deviations below the typical 10% scatter. The authors formulate an energy balance with inertial, capillary, viscous, and elastic terms, and define $\Gamma$ as the ratio of elastic to viscous energy, $\Gamma = \frac{8}{5\pi} \frac{G D_0}{\mu_i V_0 \sqrt{E}}$. Evaluating $\Gamma$ with two independent relaxation-time estimates (Zimm lower bound and CaBER upper bound) gives $\Gamma \ll 1$ for every impact tested, so the elastic term is too small to shift the spreading. They conclude that Weissenberg and Deborah numbers, while large, do not predict spreading because they compare microscopic polymer relaxation to the deformation rate rather than the macroscopic energy budget. The observed enhanced spreading of concentrated shear-thinning fluids is attributed to viscosity reduction, not elastic resistance, and concentration actually lowers $G/\mu_i$ and $\Gamma$.
Load-bearing premise
Everything rests on the Newtonian baseline from the authors' companion model being accurate for the tested conditions; if that baseline is off by more than the experimental scatter, an elasticity-induced change in spreading could be hidden inside the roughly 10% noise.
Editorial extensions
If this is right
- For Boger-fluid drops up to 1000 ppm, the maximum spreading ratio falls on the pure Newtonian master curve within experimental scatter, so adding elasticity without changing viscosity does not alter how far a drop spreads.
- The dimensionless criterion $\Gamma$, not $Wi$ or $De$, sets the detectability of elastic resistance: only when $\Gamma \approx 1$ should a measurable (greater than 10%) reduction below the Newtonian baseline appear.
- The material threshold for elastic resistance is approximately $G \approx 2 \mu_i V_0 / D_0$, which impact conditions do not reach because the viscous stress at impact shear rates overwhelms the polymer network's elastic modulus.
- In concentrated polymer drops, the extra spreading comes from shear thinning, not elasticity, and polymer concentration lowers $G/\mu_i$ and therefore $\Gamma$, making elastic resistance even less relevant.
Reading between the lines
- The authors do not test the predicted elastic regime, since all data sit at $\Gamma \ll 1$; a deliberate attempt to reach $\Gamma \approx 1$ (for instance a low-viscosity solvent with a high-modulus polymer) would provide the sharpest test of the criterion.
- The same energy-budget logic could be applied to other rapid deformations where $Wi > 1$ is used to infer elasticity, such as splashing thresholds or atomization, replacing relaxation-time criteria with an elastic-to-viscous energy ratio.
- If the Newtonian baseline's kinematic closures fail at the extremes of the tested Weber and Ohnesorge ranges, part of the claimed collapse could be masking an elasticity effect inside the noise; direct measurements of spreading time and average velocity would reveal it.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports an experimental study of the maximum spreading ratio of viscoelastic drops, using Boger fluids (PEO and PAAM in Newtonian solvents) to decouple elasticity from shear-thinning over a wide parameter range (We = 155–763, Re = 18–11,500, Oh = 0.002–1). The authors derive a dimensionless elasto-viscocapillary energy balance, introduce the ratio Gamma of elastic to viscous energy (Eq. 2), and show that Gamma < 1 for all tested conditions, concluding that elasticity is energetically negligible during spreading and that the maximum spreading ratio is identical to the Newtonian prediction, contrary to the expectations from Wi and De. The central evidence is the collapse of experimental data onto a Newtonian master curve in Fig. 2b, together with the theoretical sensitivity analysis in Fig. 3a.
Significance. If the central claim holds, the paper resolves a long-standing paradox by showing that polymer elasticity does not resist the expansion phase of drop impact, so maximum spreading remains governed by inertia, capillarity, and viscous dissipation alone. The study has clear strengths: the use of constant-viscosity Boger fluids is the right experimental strategy to separate elasticity from shear thinning; the parameter range is broad; the Hertzian-model choice is explicitly framed as a conservative upper bound; and the Gamma criterion is a parameter-free, falsifiable prediction rather than a fit to the observed spreading data. However, the verification is currently incomplete: the derivation of Eq. (3), the rheological characterization, and the dataset are all deferred to a Supplementary Material that is not included, and the Newtonian baseline is taken from the authors' own prior model. These gaps are significant because the null result is defined as a collapse onto that baseline.
major comments (4)
- [Eq. (3) and Fig. 2b] The central comparison in Fig. 2b rests on Eq. (3), but the derivation of Eq. (3)—including the kinematic closure Dmax = ts Vs, the average-velocity closure V ≈ sqrt(E), and the harmonic-mean interpolation—is stated to be in the Supplementary Material, which is not included in this preprint. The data-availability link is also the placeholder "XXX". Without the derivation and the dataset, the claimed collapse onto the Newtonian master curve cannot be independently checked. Please provide the full derivation and make the data available at a permanent repository before the central claim can be assessed.
- [Newtonian baseline, ref. [4]] Eq. (3) with Γ = 0 recovers the Newtonian baseline βmax(Γ = 0) of ref. [4], which is the authors' own prior model; neither the derivation nor the raw dataset of that baseline is shown here. The statement that "this collapse further validates the Newtonian baseline in [4]" uses the same comparison that is supposed to establish the null result. If the baseline is inaccurate for the tested range, elasticity-induced deviations of order the quoted 10% scatter could be hidden. Please provide an independent validation of the baseline, for example Newtonian solvent data measured with the same setup, or the complete baseline dataset of ref. [4].
- [Boger-fluid characterization, near Eq. (2)] The decoupling from shear thinning assumes µi = µ0 for the Boger and quasi-Boger fluids at impact shear rates much greater than 1000 s^-1. The rheological evidence for constant viscosity in this range is said to be in the Supplementary Material, which is absent. If these fluids shear-thin at impact rates, enhanced spreading from thinning could offset elastic resistance and mimic the null result. Please provide high-shear rheometry for all fluids, or an explicit estimate of µi at the characteristic impact shear rate γdot_s defined in the text.
- [Fig. 3a and experimental scatter] The detectability argument in Fig. 3a defines the observable elastic regime as deviations exceeding approximately 10% scatter, and the text states that the ratio of experimental to theoretical spreading yields an average deviation below 10%. However, no error analysis, per-condition statistics, or table of βmax values is presented in the manuscript. Since the 10% noise floor is a load-bearing threshold for the conclusion that Γ ≪ 1 produces no observable effect, please report the actual deviations and their uncertainties.
minor comments (5)
- [Fig. 1 caption] The word "captublack" in the Fig. 1 caption should be "captured".
- [Data Availability] The data availability statement contains the placeholder "link XXX"; a working permanent repository link is needed.
- [Main text, first paragraph of Section 2] There are missing spaces in phrases such as "diffusionandadsorptiontimescalesaresignificantly"; please correct the line-breaking artifacts throughout.
- [Fig. 2b axis label] The axis label "Eq. [4]" is ambiguous; it should refer to Eq. (3) of this paper or explicitly to the reference, not to a citation number.
- [Fluid classification] The term "quasi-Boger fluids" is used in Fig. 2b and Fig. 3 but is not defined in the main text; please provide a definition or a criterion for this category.
Circularity Check
Minor self-citation burden: the Newtonian reference curve is the authors' own prior model (Eq. 3 with Γ=0, ref. [4]), but the experimental collapse provides independent support; no fitted-to-spreading parameter is involved.
-
self citation load bearing
[Section 'To evaluate deviations caused by fluid elasticity...' and Eq. (3), Fig. 2b]
"To evaluate deviations caused by fluid elasticity, we must establish a purely Newtonian baseline. We follow a similar energy balance presented in [4]. ... As Γ→0, the equation recovers the purely Newtonian spreading ratio βmax(Γ = 0) limit [4]. ... This collapse further validates the Newtonian baseline in [4]."
The reference curve used to establish the central claim that elastic properties do not alter maximum spreading is not an independent empirical Newtonian master curve presented in this paper; it is the Γ→0 limit of the authors' own prior model, ref. [4] (Abbot and Bonn, arXiv:2605.11797). The paper explicitly says it follows 'a similar energy balance presented in [4]' and that Eq. (3) recovers the Newtonian baseline from [4]. The key closures of that baseline (Dmax≡tsVs and V≈√E) are stated but their derivation is not shown here, being relegated to a missing Supplementary Material. The paper then uses the collapse of its own new data to 'further validate' this same self-cited baseline.
full rationale
The central null result is primarily an experimental observation: Boger-fluid and quasi-Boger-fluid drops up to 1000 ppm collapse onto the Newtonian baseline in Fig. 2b, with deviations below the quoted ~10% scatter. Because no parameter of Eq. (3) is fitted to the observed spreading values (We, Re, Oh, G, and τ are all independently measured), the collapse is genuine evidence and the paper does not commit a fitted-input-called-prediction circularity. The Γ-criterion is model-based, but it is a conservative upper bound (Hertzian β^5 scaling) and is not defined in terms of the observed βmax, so it is not self-definitional in the strict sense. The only notable circularity-adjacent burden is the reference curve itself: the Newtonian baseline βmax(Γ=0) is taken from the authors' own prior model (ref. [4]), with the kinematic closure Dmax≡tsVs and average-velocity closure V≈√E relegated to a missing Supplementary Material. The paper then uses the collapse of its own data to 'validate' that self-cited baseline. This makes the comparison dependent on an unshown, self-cited derivation, but the experimental data provide independent content, so the severity is low (score 2, not a forced circle). Missing data link and Supplementary Material are reproducibility concerns rather than circularity.
Assumptions & free parameters
assumptions (7)
- domain assumption The Newtonian maximum-spreading baseline beta_max(Gamma=0) from ref. 4 is accurate for the tested fluids and substrates.
- domain assumption Elastic energy storage during spreading is bounded by a linear elastic Hertzian contact model with scaling of order G beta^5.
- domain assumption Surface tension during the spreading phase equals the pure-solvent plateau value (about 72 mN/m) because polymer adsorption and diffusion timescales exceed impact timescales.
- domain assumption The Boger fluids have essentially constant viscosity at impact shear rates, so mu_i = mu_0 and shear thinning is absent.
- domain assumption The kinematic closure Dmax = ts Vs and the leading-order average velocity closure V approximately sqrt(E) are valid.
- domain assumption The characteristic impact shear rate is gamma_dot_s approximately (3 V0 / 2 D0) V beta_max^2, from volume conservation with h_min approximately 2 D0^3 / (3 Dmax^2).
- domain assumption Zimm relaxation times give a lower bound on tau (and therefore an upper bound on G), while CaBER relaxation times give an upper bound on tau (and therefore a lower bound on G).
Cite this review
Pith. "Pith review of Inelastic spreading of viscoelastic drops." pith.science (2026). https://pith.science/paper/XMBVTPHF
@misc{pith2026260810851,
author = {Pith},
title = {Pith review of: Inelastic spreading of viscoelastic drops},
year = {2026},
howpublished = {\url{https://pith.science/paper/XMBVTPHF}},
note = {Machine review of arXiv:2608.10851}
}
abstract
When a liquid drop impacts a solid surface, rebound and splashing are known to be suppressed by additives that induce elastic effects; however, the influence of elasticity on the maximum spreading radius remains debated. The difficulty lies in isolating elastic resistance from the enhanced spreading caused by viscous shear thinning. Here, we experimentally decouple these effects by investigating Boger fluid drops (constant-viscosity solutions with high elasticity) of polyethylene oxide (PEO) and polyacrylamide (PAAM). We demonstrate that elastic properties do not alter the maximum spreading ratio, even at high concentrations up to 1000 ppm. We formulate an energy balance that accounts for the inertial, capillary, viscous, and elastic contributions, and yields a dimensionless criterion ($\Gamma$) that estimates the degree of elastic effects. We show that $\Gamma \ll 1$ for all impact conditions tested, demonstrating that elastic effects are energetically negligible during droplet spreading, opposite to what Weissenberg and Deborah numbers would predict.
Figures
Figures from the paper (1 more)
Reference graph
Works this paper leans on
-
[4]
Kinematic Closure of Drop Impact
M. Abbot, D. Bonn, Kinematic closure of drop impact, arXiv preprint arXiv:2605.11797 (2026)
work page Pith review arXiv 2026
-
[1]
C. Josserand, S. T. Thoroddsen, Drop impact on a solid surface, Annual review of fluid me- chanics 48 (1) (2016) 365–391
work page 2016
-
[2]
A. L. Yarin, Drop impact dynamics: splashing, spreading, receding, bouncing..., Annu. Rev. Fluid Mech. 38 (1) (2006) 159–192
work page 2006
-
[3]
X. Cheng, T.-P. Sun, L. Gordillo, Drop impact dynamics: Impact force and stress distributions, Annual Review of Fluid Mechanics 54 (1) (2022) 57–81
work page 2022
-
[5]
N. Laan, K. G. De Bruin, D. Bartolo, C. Josserand, D. Bonn, Maximum diameter of impacting liquid droplets, Physical Review Applied 2 (4) (2014) 044018
work page 2014
-
[6]
J. B. Lee, N. Laan, K. G. de Bruin, G. Skantzaris, N. Shahidzadeh, D. Derome, J. Carmeliet, D. Bonn, Universal rescaling of drop impact on smooth and rough surfaces, Journal of Fluid Mechanics 786 (2016) R4
work page 2016
-
[7]
T. C. de Goede, K. G. de Bruin, N. Shahidzadeh, D. Bonn, Predicting the maximum spreading of a liquid drop impacting on a solid surface: effect of surface tension and entrapped air layer, Physical review fluids 4 (5) (2019) 053602
work page 2019
-
[8]
T. Varkevisser, D. Bonn, Capillary-viscous retraction dynamics of droplets: The role of the dynamic contact angle, Physical Review Fluids 11 (1) (2026) 013601
work page 2026
Show all 40 references
-
[9]
Abbot, M
M. Abbot, M. H. Iqbal, L. Liu, E. Koos, I. V. Roisman, J. Hussong, A. A. Castrejón-Pita, J. R. Castrejón-Pita, Nanoparticles do not influence droplet break-up, spreading, or splashing, Journal of Colloid and Interface Science 693 (2025) 137570
2025
-
[10]
R. G. Larson, The rheology of dilute solutions of flexible polymers: Progress and problems, Journal of Rheology 49 (1) (2005) 1–70
2005
-
[11]
Gaillard, M
A. Gaillard, M. A. Herrada, A. Deblais, J. Eggers, D. Bonn, Beware of caber: Filament thinning rheometry does not always give ‘the’relaxation time of polymer solutions, Physical Review Fluids 9 (7) (2024) 073302. 7
2024
-
[12]
Vadillo, T
D. Vadillo, T. Tuladhar, A. Mulji, M. Mackley, The rheological characterization of linear vis- coelasticityforinkjetfluidsusingpiezoaxialvibratorandtorsionresonatorrheometers, Journal of Rheology 54 (4) (2010) 781–795
2010
-
[13]
Clasen, J
C. Clasen, J. Plog, W.-M. Kulicke, M. Owens, C. Macosko, L. Scriven, M. Verani, G. H. McKinley, How dilute are dilute solutions in extensional flows?, Journal of Rheology 50 (6) (2006) 849–881
2006
-
[14]
Bergeron, D
V. Bergeron, D. Bonn, J. Y. Martin, L. Vovelle, Controlling droplet deposition with polymer additives, Nature 405 (6788) (2000) 772–775
2000
-
[15]
Rozhkov, B
A. Rozhkov, B. Prunet-Foch, M. Vignes-Adler, Impact of drops of polymer solutions on small targets, Physics of Fluids 15 (7) (2003) 2006–2019
2003
-
[16]
Bartolo, A
D. Bartolo, A. Boudaoud, G. Narcy, D. Bonn, Dynamics of non-newtonian droplets, Physical review letters 99 (17) (2007) 174502
2007
-
[17]
Kansal, V
M. Kansal, V. Bertin, C. Datt, J. Eggers, J. H. Snoeijer, Viscoelastic wetting: Cox–voinov theory with normal stress effects, Journal of fluid mechanics 985 (2024) A17
2024
-
[18]
P. Dhar, S. R. Mishra, D. Samanta, Onset of rebound suppression in non-newtonian droplets post-impact on superhydrophobic surfaces, Physical Review Fluids 4 (10) (2019) 103303
2019
-
[19]
Varkevisser, On droplets and sprays, Ph.d
T. Varkevisser, On droplets and sprays, Ph.d. thesis, Universiteit van Amsterdam (2026)
2026
-
[20]
Singh, S
M. Singh, S. Basu, D. Samanta, Effect of substrate flexibility on non-newtonian droplet impact dynamics on flexible superhydrophobic surfaces, Journal of Colloid and Interface Science 679 (2025) 100–112. doi:10.1016/j.jcis.2025.02.012
2025 doi
-
[21]
Gaillard, R
A. Gaillard, R. Sijs, D. Bonn, What determines the drop size in sprays of polymer solutions?, Journal of Non-Newtonian Fluid Mechanics 305 (2022) 104813
2022
-
[22]
Varkevisser, D
T. Varkevisser, D. Bonn, Effect of high wind speeds on droplet formation in sprays of dilute polymer solutions, Journal of Non-Newtonian Fluid Mechanics (2026) 105618
2026
-
[23]
E. Vega, A. Castrejón-Pita, Suppressing prompt splash with polymer additives, Experiments in Fluids 58 (5) (2017) 57
2017
-
[24]
Mobaseri, S
A. Mobaseri, S. Kumar, X. Cheng, Maximum spreading of impacting shear-thinning and shear-thickening drops, Proceedings of the National Academy of Sciences 122 (22) (2025) e2500163122
2025
-
[25]
Abbot, M
M. Abbot, M. Lannert, A. Kiran, S. Bakshi, J. Hussong, I. V. Roisman, Spreading of a viscous drop after impact onto a spherical target, Journal of Fluid Mechanics 996 (2024) A10
2024
-
[26]
Abbot, F
M. Abbot, F. Yin, M. Piskunov, I. V. Roisman, J. Hussong, Elucidating on the drop impact dynamics of water-in-oil microemulsions via advanced rheometry, Colloids and Surfaces A: Physicochemical and Engineering Aspects 720 (2025) 136952
2025
-
[27]
Abbot, F
M. Abbot, F. Yin, A. L. Azevedo Costa, I. V. Roisman, J. Hussong, Minimal addition of fumed nanoparticles suppresses droplet splashing, Applied Physics Letters 127 (8) (2025)
2025
-
[28]
S. Tang, X. Li, W. Wang, W. Zhao, Y. Zhou, S. Wang, Y. Li, P. Yu, X. Wen, G. Hu, et al., Droplet impact-based microliter viscometry, Analytical Chemistry 97 (25) (2025) 13076–13085. 8
2025
-
[29]
Piskunov, X
M. Piskunov, X. Zhang, W. Yuan, F. Chen, L. Chen, M. Abbot, Spreading of water-in-oil emulsion drop on oleophilic surface, Physics of Fluids 37 (5) (2025)
2025
-
[30]
O. Avni, D. Wang, M. Ravisankar, R. Zenit, Maximal spreading of impacting viscoelastic droplets, arXiv preprint arXiv:2601.15246 (2026)
2026 arXiv
-
[31]
Shende, E
T. Shende, E. Moeendarbary, I. Eames, Forces and energy of viscoelastic droplets impacting a wall, Journal of Rheology 70 (5) (2026) 995–1021. doi:10.1122/8.0001202
2026 doi
-
[32]
K. L. Johnson, Contact mechanics, Cambridge university press, 1987
1987
-
[33]
S. Jana, J. Kolinski, D. Lohse, V. Sanjay, Impacting spheres: from liquid drops to elastic beads, Soft matter 22 (11) (2026) 2226–2236
2026
-
[34]
Richard, C
D. Richard, C. Clanet, D. Quéré, Contact time of a bouncing drop, Nature 417 (6891) (2002) 811–811
2002
-
[35]
Mousavi, S
M. Mousavi, S. A. Faroughi, Spreading and rebound of viscoelastic droplets on surfaces with hybrid wettability, Physics of Fluids 37 (1) (2025)
2025
-
[36]
R. A. Figueiredo, C. M. Oishi, J. A. Cuminato, J. C. Azevedo, A. M. Afonso, M. A. Alves, Numerical investigation of three dimensional viscoelastic free surface flows: impacting drop problem, in: 11th World Congress on Computational Mechanics, WCCM, 2014
2014
-
[37]
Zhang, Q
H. Zhang, Q. Luan, W. Yuan, F. Chen, B. Meng, Impact and spread dynamics of a viscoelastic droplet on an inclined hydrophilic surface, Physics of Fluids 36 (8) (2024)
2024
-
[38]
Chowdhury, S
A. Chowdhury, S. Mitra, S. K. Mitra, Bridging liquid and elastic solid impact regimes using flexible hydrogels, Langmuir (2026)
2026
-
[39]
I. V. Roisman, Inertia dominated drop collisions. ii. an analytical solution of the navier–stokes equations for a spreading viscous film, Physics of Fluids 21 (5) (2009)
2009
-
[40]
Campo-Deano, C
L. Campo-Deano, C. Clasen, The slow retraction method (srm) for the determination of ultra- short relaxation times in capillary breakup extensional rheometry experiments, Journal of Non- Newtonian Fluid Mechanics 165 (23-24) (2010) 1688–1699. 9
2010
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.