Pith. sign in

REVIEW 3 major objections 4 minor 5 references

Fluid elasticity alone—not viscosity or surface tension—is claimed to control droplet-pool impacts, with up to 30–40% of impact energy stored as polymer stretch and Worthington-jet filaments thinning exponentially.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review

2026-08-04 00:49 UTC pith:KEOK57HN

load-bearing objection A visually rich, useful experimental study of elastic droplet-pool impacts whose headline energy-storage and jet-thinning claims are fitted residuals, not predictions. the 3 major comments →

arxiv 2608.00274 v1 pith:KEOK57HN submitted 2026-07-31 physics.flu-dyn cond-mat.softphysics.app-phphysics.pop-ph

Elasto-hydrodynamics of droplet-pool-interactions

Md Sultan , Purbarun Dhar This is my paper
classification physics.flu-dyn cond-mat.softphysics.app-phphysics.pop-ph
keywords droplet-pool impactBoger fluidselastic energy storageWorthington jetelasto-capillary thinningFENE-P modelDeborah numbercavity dynamics
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper argues that when a droplet of an elastic (Boger) fluid falls into a pool of the same fluid, the impact is governed by fluid elasticity, not by shear viscosity or surface tension. It reports new morphological regimes—trapezoidal cavity reversal, suppressed crown breakup, beads-on-a-string filaments, and no-pinch-off jets—and quantifies that 30–40% of the droplet's kinetic impact energy can be stored as elastic energy in stretched polymer chains during cavity expansion. It also derives, from an energy balance and the FENE-P constitutive model, that the Worthington jet filament radius decays exponentially as exp(-t/2λ), in contrast to the classical exp(-t/3λ) capillary-thinning law. A sympathetic reader would care because the paper presents elastic droplet-pool impact as a distinct elastohydrodynamic phenomenon, with direct consequences for predicting and controlling splashing, jetting, and aerosol formation in polymer-laden fluids.

Core claim

With experiments on dilute polyethylene-oxide (PEO) Boger fluids, the paper claims that droplet-pool impact enters a regime where elasticity alone sets the outcome. The evidence is a matched-viscosity comparison: a water-glycerol Newtonian pool and a 2000 ppm PEO pool have essentially the same shear viscosity and surface tension, yet only the PEO pool produces a trapezoidal cavity, a flattened floor, delayed capillary waves, and suppressed crown breakup. Quantitatively, the paper claims that during cavity expansion roughly 30–40% of the droplet's kinetic impact energy is stored as elastic energy in stretched polymer chains, estimated by fitting an energy-conservation balance to measured cavi

What carries the argument

The argument is carried by two theoretical objects. First, an energy-conservation ODE for the hemispherical cavity radius (Eq. 4.13) with two calibration factors: β, which absorbs all neglected losses (crown formation, shape deviations, minor viscous dissipation) and is fixed by the measured maximum cavity radius, and α, which calibrates the potential-flow kinetic-energy estimate. The elastic-energy claim is the residual difference in β between water and polymer systems. Second, a one-dimensional elasto-capillary thinning model of the Worthington jet built on the FENE-P constitutive equation, assuming a constant axial tensile force and polymer recoil to equilibrium at jet birth; this yields

Load-bearing premise

The load-bearing premise is that the two calibration factors in the cavity energy balance—which are supposed to absorb crown energy, shape deviations, droplet deformation, and minor viscous losses—take the same values for water and polymer fluids; the 30–40% elastic-storage figure is the leftover difference in one of those factors, so if the hidden losses differ between fluids, the number is not actually elastic energy.

What would settle it

Recompute the same energy balance using the measured trapezoidal cavity geometry and directly measured crown energy, or measure polymer stretch along the cavity wall by birefringence; if the water-versus-polymer gap in the fitted factor β vanishes, the claimed 30–40% elastic energy is an artifact of unmodeled losses rather than stored polymer stress.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • If 30–40% of impact energy is stored elastically, then maximum cavity size and Worthington-jet height should decrease monotonically with increasing Deborah number, as observed.
  • Elasticity should delay or suppress capillary pinch-off, shifting post-impact outcomes from pinch-off to satellite-filament to no-pinch-off as pool or droplet elasticity rises.
  • The exponential r(t) ~ exp(-t/2λ) law gives a one-line diagnostic: the thinning slope of the Worthington-jet filament can be used to extract an apparent relaxation time from simple droplet-pool experiments.
  • Beads-on-a-string formation should be governed by whether the Rayleigh-instability growth rate exceeds the inverse polymer relaxation time, with early polymer stretching suppressing bead formation.
  • Cavity reversal morphology changing from hemispherical to trapezoidal is a purely elastic signature, since matched-viscosity Newtonian fluids retain hemispherical cavities.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the 30–40% elastic-energy share is real, droplet-pool impact becomes a passive energy-storage event: polymer chains act as transient springs, so elasticity could be tuned to suppress splash and aerosol generation in industrial sprays without changing viscosity or surface tension.
  • The 1/2λ exponent (vs the standard 1/3λ) emerges from a 0+1-dimensional constant-tension reduction with a recoiled initial polymer state; a full spatially resolved simulation with finite extensibility and axial stress variation could test whether the exponent is robust or an artifact of that reduction.
  • The paper itself notes that the fitted relaxation time is about two orders of magnitude larger than rheometric values; this suggests a single-mode FENE-P description is incomplete, and a multi-mode or polydisperse relaxation spectrum should be tested against the same experiments.
  • The β-difference method could be independently verified by measuring stored elastic energy directly, for example through flow birefringence near the cavity wall; without such a check, the 30–40% figure remains a fitted residual.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. This manuscript reports a combined experimental, theoretical, and computational study of droplet impact onto deep liquid pools for elastic (Boger) fluids made from aqueous PEO solutions. The authors document several new morphological regimes — trapezoidal cavity reversal, suppression of crown breakup, beads-on-a-string filaments, and no-pinch-off states — and organize them in Weber-number/Deborah-number regime maps. They propose an energy-conservation model for the cavity expansion in which two correction factors α and β are introduced, and they interpret differences in β between water and PEO cases as elastic energy stored in stretched polymer chains, claiming that 30–40% of the droplet's initial kinetic energy is stored elastically. They also develop a FENE-P-based one-dimensional thinning model for the Worthington jet filament and claim an exponential decay r(t) ~ exp(-t/(2λ)) in the elasto-capillary regime. The simulations couple the Cahn-Hilliard phase-field method with FENE-P viscoelasticity and are validated qualitatively against the experiments.

Significance. If the central quantitative claims were supported, this would be a substantial contribution: it would identify a distinct class of elastohydrodynamic droplet-pool interactions and provide both a regime classification and a predictive framework for elastic-energy partitioning and jet thinning. The experimental database is extensive and carefully controlled — the matched-viscosity water-glycerol comparison is a good control, the parameter ranges in We and De are broad, and the regime maps are a useful organizing device. The FENE-P/phase-field simulations are novel for this configuration and reproduce many observed morphologies. However, the two headline quantitative results — the 30–40% elastic-energy fraction and the exp(-t/2λ) thinning law — are currently resting on fitted parameters and are not independently verified. The significance of the paper therefore depends on whether these identifications can be strengthened.

major comments (3)
  1. [§4.2, Eqs. (4.13)–(4.14)] The elastic-energy storage fraction is not measured; it is a difference of fitted β values. β is determined from the measured maximum cavity radius via Eq. (4.14) and α is tuned to reproduce R*(t). The identification of Δβ with stored elastic energy assumes that all unmodeled losses — crown formation, non-hemispherical cavity geometry, droplet deformation, altered momentum transfer — are identical in the water and PEO cases. The authors' own numbers contradict this: the pool-only reduction is Δβ=0.27, the droplet-only reduction is Δβ=0.43, but the combined case gives Δβ=0.35, not ~0.70. This non-additivity shows that β absorbs configuration-specific non-elastic losses, so the headline claim of 30–40% elastic energy storage (Abstract, §4.2, §6) is unsupported as stated.
  2. [§4.4, Eq. (4.35), Figs. 15–16] The claimed exponential thinning law is not a parameter-free prediction. The authors state that Eq. (4.35) predicts a relaxation time of about 29 ms, whereas the independently measured extensional relaxation time of the 500 ppm PEO solution is about 0.14 ms — a discrepancy of roughly a factor of 200. Since λ appears explicitly in the exponent of Eq. (4.35), the agreement in Fig. 15 is achieved only by using an effective relaxation time inferred from the same thinning data. This means the exponential decay is a fitted functional form, not a validation of the FENE-P model, and it undermines the conclusion in §6 that the measurement can be used as a simple method for estimating the extensional relaxation time.
  3. [§5, Figs. 17–19] The simulations are the natural place to test the elastic-energy partition independently, but they do not. Only qualitative morphological comparisons, elastic-stress profiles, and velocity fields are reported. No integrated polymer elastic energy is computed or compared with the β-based estimate. Without such a diagnostic, the simulations cannot corroborate the '30–40% stored energy' claim; they only show that elastic stresses are present and localized. This should be remedied by post-processing the existing simulations to compute the volume-integrated polymeric energy during cavity expansion.
minor comments (4)
  1. [§4.3] Cross-referencing errors: the text says 'as shown in figure 12' for the normalized jet height evolution, but the data appear in figure 13. Similarly, §3.4 says the regime maps are 'presented in figure 9' when they are in figure 8. Please correct the figure numbering throughout.
  2. [Table I] The μ_p/μ_s column appears inconsistent with the total viscosity column. For 1000 ppm, μ=1.15 mPa·s and μ_s=1.0 mPa·s implies μ_p/μ_s=0.15, not 0.281; for 2000 ppm the ratio should be 0.38, not 0.533. Clarify whether μ_p denotes the polymeric contribution or the zero-shear polymer viscosity, and ensure the table entries are consistent.
  3. [§4.1, Eqs. (4.6)–(4.8)] Please verify the numerical prefactors in the velocity potential and kinetic-energy expressions. Direct differentiation of Eq. (4.6) gives u_r = R^2 \dot R / r^2; the form shown in Eq. (4.7) appears to have an extra factor of 2, which would propagate into Eq. (4.8). If the factor is intentional (e.g., due to volume averaging), state the definition explicitly.
  4. [References] Sen et al. (2022) is cited as an arXiv preprint. If a peer-reviewed version has appeared, it should be cited to allow readers to verify the rheological data. Also, there are duplicate entries in the reference list (e.g., Bazilevsky et al. 1990 appears twice).

Circularity Check

2 steps flagged

Headline elastic-energy fraction is a fitted β-residual, and the jet exponential decay uses a relaxation time inferred from the same thinning data.

specific steps
  1. fitted input called prediction [§4.2, Eqs. (4.13)-(4.14), discussion of Fig. 10]
    "To account for the kinetic energy of the cavity and the initial droplet impact, two correction factors, α and β, are introduced. ... The β compensates for energy losses that are neglected in the model, such as minor viscous dissipation, and the energy associated with crown formation. ... By comparing the value of β corresponding to the 0D-0P case, with those of the other droplet-pool combinations, we quantify the amount of elastic energy stored by the stretching of the long-chain polymer molecules during cavity formation."

    Eq. (4.14) fixes β from the measured R*_max, and α is then tuned to the measured R*(t); so the 'model prediction' in Eq. (4.13) is a fit to the cavity-radius data, not an independent prediction. The 30-40% elastic-storage number is simply 0.80 minus the fitted β of the polymer case, with no independent measurement of stored elastic energy. The claimed segregation is also internally inconsistent: the separately attributed reductions (0.27 pool + 0.43 droplet) do not sum to the 0.35 reduction in the combined case, showing that β absorbs configuration-specific crown/shape/viscous losses rather than a unique elastic energy fraction.

  2. fitted input called prediction [§4.4, Eq. (4.35), discussion after Fig. 16]
    "At time 40t > ms (figure 16b), the filament radius of the elastic Worthington jet starts decreasing exponentially with time as ( ) exp( / 2 )r t t λ − , in agreement with the elasto-capillary thinning model. ... The relaxation time predicted by Eq. (4.35) is approximately 29 ms, whereas the experimentally measured relaxation time of a 500 ppm aqueous PEO solution, obtained from extensional rheometry, is ~0.14 ms (Sen et al. 2022)."

    The exp(-t/2λ) law is presented as a derived prediction and confirmed by data, but the λ used to match the thinning curves is an effective 29 ms inferred from those same curves; the independently measured λ is 0.14 ms, which the paper admits would miss by two orders of magnitude. Thus the comparison does not test the elasto-capillary balance; it only reports an exponential fit with a fitted parameter, and the 'prediction' collapses when the independently measured input is used.

full rationale

The paper contains genuinely independent observational content: matched-viscosity water-glycerol vs PEO controls, regime maps, morphology tracking, and FENE-P simulations that are validated against experiments. Those components are not circular. The circularity is concentrated in the quantitative headline claims. The energy-balance model is calibrated to the measured cavity radius through α and β and then the fitted β difference is relabeled as stored elastic energy; no independent elastic-energy measurement is supplied. The additivity failure (0.27 + 0.43 ≠ 0.35) shows that β is absorbing more than elastic storage. Similarly, the jet-thinning 'prediction' r(t) ~ exp(-t/2λ) uses an effective relaxation time obtained from the same thinning data rather than the independently measured 0.14 ms, so the agreement is a fit, not a first-principles forecast. Self-citations in the paper (e.g., Dhar et al. 2019) are not load-bearing, and key rheological inputs come from Sen et al. 2022, an external source. Overall, the paper is partially circular in its central quantitative inferences, meriting a 6 rather than a lower score.

Axiom & Free-Parameter Ledger

3 free parameters · 4 axioms · 0 invented entities

No new physical entities are introduced. The central quantitative claims rest on fitted correction factors (α, β) and an effective relaxation time that differs by two orders of magnitude from the independently measured value. The modeling axioms are standard for viscoelastic filament thinning but are applied here to a spatially evolving Worthington jet without direct validation.

free parameters (3)
  • β (cavity energy-loss correction factor) = 0.80 (0D-0P), 0.45 (500D-500P), 0.53 (0D-500P), 0.37 (500D-0P)
    Computed from the experimental maximum cavity radius via Eq. 4.14; it 'compensates for energy losses that are neglected in the model'. The residual differences in β are interpreted as elastic energy storage, so the headline 30–40% number is a difference of fitted values.
  • α (cavity kinetic-energy correction factor) = not stated
    The paper says 'α is adjusted to achieve the best agreement between the model predictions and experimentally measured time evolution of the mean cavity radius' (§4.2). It is a free fit parameter used to make the energy-balance model match R*(t).
  • effective jet relaxation time λ_eff = ~29 ms (inferred from thinning data)
    The paper notes Eq. 4.35 'overestimates the relaxation time by approximately two orders of magnitude' relative to the measured 0.14 ms. The exponential decay exp(-t/2λ) therefore only matches data if λ is effectively extracted from the same filament-thinning curves rather than taken from independent rheometry.
axioms (4)
  • domain assumption Potential-flow, irrotational, incompressible flow around a hemispherical cavity, with viscous dissipation neglected
    Invoked before Eq. 4.1 in §4.2 to build the cavity energy balance; supported by the authors' claim that Oh<0.0046 makes viscosity negligible, but the hemispherical shape is explicitly violated by their own observed trapezoidal cavities.
  • ad hoc to paper Single-relaxation-mode FENE-P with f≈1, A_zz >> A_rr, and A_zz << b (finite extensibility not reached)
    Assumed in §4.4 to reduce the FENE-P equations to Eqs. 4.30-4.31 and obtain the exponential decay law. If finite extensibility is approached, the exponential form changes.
  • domain assumption Polymer chains fully recoil to equilibrium when the Worthington jet emerges, so A_zz(0)=1 and the initial tensile force is T=σ r0
    Used in §4.4 (Eqs. 4.15-4.16 and the initial condition for Eq. 4.31) to fix the integration constant. The paper offers indirect support via De-independent initial jet velocity, but no direct measurement of polymer conformation at jet birth.
  • domain assumption The tensile force is constant along the jet length
    Standard slender-filament assumption adopted from Clasen et al. 2009, used to equate the elastic stress contribution in Eq. 4.17.

pith-pipeline@v1.3.0-alltime-deepseek · 27315 in / 13115 out tokens · 124581 ms · 2026-08-04T00:49:40.240307+00:00 · methodology

0 comments
read the original abstract

In Newtonian fluids, impact of a droplet on a liquid pool births a cavity, crown, capillary waves, and Worthington jet. The corresponding hydrodynamic events for elastic or Boger fluids, however, remain an uncharted domain of comprehension and exploration. We thoroughly investigate, via experiments, theory, and simulations, how elastic energy storage, fluid relaxation, and competitive inertio elasto capillarity govern the spatio temporal evolution of the cavity, the crown, and the ensuing Worthington jet in polymeric elastic fluids. The events are systematically explored over a wide range of impact Weber and Deborah numbers, considering varied Newtonian and elastic fluid droplet pool combinations, and revealing new, and distinct morphological regimes compared to Newtonian counterparts. We illustrate that these new findings are purely driven by fluid elasticity, and not by viscosity or interfacial tension. We derive a theory for cavity radius evolution, using energy conservation within potential-flow framework. We show that 30-40 % of the droplets kinetic impact energy may be stored as elastic energy by the stretching polymer chains during cavity expansion. Appealing to the FENE P model, we derive a theory for the temporal evolution of the radius of the elongated Worthington jet. We show that in elasto capillary regime, competitive elastic and capillary stresses lead to exponential decay of the jet radius. The role of elastic stresses and the local velocity field in governing cavity evolution, morphology, and jet formation are further elucidated through computer simulations. Our findings significantly advance the uncharted paradigm of interplay between inertia, capillarity, and elasticity in droplet-pool interaction elastohydrodynamics.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

5 extracted references · 1 linked inside Pith

  1. [1]

    Amarouchene, Y ., Bonn, D., Meunier, J., & Kellay, H. (2001). Inhibition of the finite-time singularity during droplet fission of a polymeric fluid. Physical Review Letters, 86(16),

    2001

  2. [110]

    Eggers, J. (1997). Nonlinear dynamics and breakup of free -surface flows. Reviews of modern physics, 69(3),

    1997

  3. [865]

    Eggers, J., & Villermaux, E. (2008). Physics of liquid jets. Reports on progress in physics , 71(3), 036601. Eggers, J., Herrada, M. A., & Snoeijer, J. H. (2020). Self -similar breakup of polymeric threads as described by the Oldroyd-B model. Journal of fluid mechanics, 887, A19. Engel, O. G. (1966). Crater depth in fluid impacts. Journal of Applied Physi...

    2008

  4. [3558]

    D., & Basaran, O

    Ambravaneswaran, B., Wilkes, E. D., & Basaran, O. A. (2002). Drop formation from a capillary tube: Comparison of one -dimensional and two -dimensional analyses and occurrence of satellite drops. Physics of Fluids, 14(8), 2606-2621. Anna, S. L., McKinley, G. H., Nguyen, D. A., Sridhar, T., Muller, S. J., Huang, J., & James, D. F. (2001). An interlaboratory...

    2002

  5. [6083]

    S., & Samanta, D

    Kamaluddin, K., Dhar, P., Sharma, C. S., & Samanta, D. (2025). Quantification of sub -surface capillary pressure distribution within unstructured superhydrophobic surfaces. Colloids and Surfaces A: Physicochemical and Engineering Aspects, 138498. Kamaluddin, K., Pal, G. C., Dhar, P., Sharma, C. S., & Samanta, D. (2024). Parallel governing criteria for non...

    Pith/arXiv arXiv 2025

This paper was first reviewed by deepseek-v4-flash on August 4, 2026.