Impacting spheres: from liquid drops to elastic beads

Detlef Lohse; John Kolinski; Saumili Jana; Vatsal Sanjay

arxiv: 2510.24855 · v3 · submitted 2025-10-28 · ❄️ cond-mat.soft · physics.flu-dyn

Impacting spheres: from liquid drops to elastic beads

Saumili Jana , John Kolinski , Detlef Lohse , Vatsal Sanjay This is my paper

Pith reviewed 2026-05-18 02:50 UTC · model grok-4.3

classification ❄️ cond-mat.soft physics.flu-dyn

keywords viscoelastic impactliquid dropelastic beadWagner scalingHertz scalingWeissenberg numberelasticity numberimpact force

0 comments

The pith

Varying the Weissenberg number produces a continuous transition from liquid-drop to elastic-bead impact behavior.

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

The paper uses direct numerical simulations to examine a viscoelastic sphere striking a rigid non-wetting surface. By varying the elasticity number and Weissenberg number, the model recovers the spreading-retracting-jumping sequence of a Newtonian liquid at low values and the brief contact and bounce of an elastic solid at high values. This sweep reveals how memory of deformation grows continuously with the Weissenberg number, unifying the classical Wagner and Hertz descriptions of impact force. A reader would care because the force exerted during impact controls many natural and industrial processes, and the work supplies a single tunable framework for the liquid-to-solid transition.

Core claim

Using direct numerical simulations, a viscoelastic sphere impacting a rigid non-contacting surface is shown to produce impact forces that recover the Newtonian liquid response as either the elasticity number approaches zero or the Weissenberg number approaches zero. Elastic-solid behavior appears in the limit of infinite Weissenberg number with nonzero elasticity number. In this elastic-memory limit three regimes emerge: capillary-dominated, Wagner scaling, and Hertz scaling, with a smooth transition from the Wagner to the Hertz regime. Sweeping the Weissenberg number from zero to infinity produces a continuous shift from materials with no memory to materials with permanent memory of the de-

What carries the argument

The Weissenberg number Wi, which measures the material relaxation time relative to the impact time and thereby sets the degree of permanent memory retained in the deformation.

If this is right

The liquid limit is recovered when either El approaches zero or Wi approaches zero.
In the elastic-memory limit the force exhibits capillary-dominated, Wagner, and Hertz regimes.
A smooth transition occurs from Wagner scaling to Hertz scaling as parameters change.
The framework supplies a controlled route connecting liquid-drop and elastic-bead impacts.

Where Pith is reading between the lines

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

The same parameter sweep could be applied to other soft materials such as polymer gels or biological cells to predict their impact response.
Material synthesis might target specific relaxation times to achieve desired spreading or rebound behavior on demand.
The observed scalings suggest possible analytical expressions that interpolate between the three regimes without full simulation.
The approach could be extended to impacts on deformable or wetting substrates to check whether the memory transition persists.

Load-bearing premise

The specific viscoelastic constitutive model recovers both the Newtonian liquid limit and the elastic solid limit without missing contact-line dynamics or nonlinear material effects.

What would settle it

Measure the peak impact force for a viscoelastic material while varying its relaxation time over several orders of magnitude and test whether the measured forces interpolate smoothly between the liquid and solid theoretical predictions.

read the original abstract

A liquid drop impacting a non-wetting rigid substrate laterally spreads, then retracts, and finally jumps off again. An elastic solid, by contrast, undergoes a slight deformation, contacts briefly, and bounces. The impact force on the substrate - crucial for engineering and natural processes - is classically described by Wagner's (liquids) and Hertz's (solids) theories. This work bridges these limits by considering a generic viscoelastic medium. Using direct numerical simulations, we study a viscoelastic sphere impacting a rigid, non-contacting surface and quantify how the elasticity number ($El$, dimensionless elastic modulus) and the Weissenberg number ($Wi$, dimensionless relaxation time) dictate the impact force. We recover the Newtonian liquid response as either $El \to 0$ or $Wi \to 0$, and obtain elastic-solid behavior in the limit $Wi \to \infty$ and $El \ne 0$. In this elastic-memory limit, three regimes emerge - capillary-dominated, Wagner scaling, and Hertz scaling - with a smooth transition from the Wagner to the Hertz regime. Sweeping $Wi$ from 0 to $\infty$ reveals a continuous shift from materials with no memory to materials with permanent memory of deformation, providing an alternate, controlled route from liquid drops to elastic beads. The study unifies liquid and solid impact processes and offers a general framework for the liquid-to-elastic transition relevant across systems and applications.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

Viscoelastic simulations give a clean parameter sweep from liquid drop impacts to elastic bounces, recovering the classical limits but thin on validation of the liquid retraction phase.

read the letter

The main thing to know is that by sweeping Wi from zero to infinity in their viscoelastic model they get a continuous change from liquid-like spreading and rebound to solid-like Hertz contact, with three regimes (capillary, Wagner, Hertz) showing up in the high-memory limit and a smooth crossover between the last two. They recover the Newtonian limit at low El or Wi and the elastic-solid limit at high Wi, which is the cleanest part of the work. The direct numerical simulations provide a controlled route that connects two previously separate impact literatures without inventing new equations. The regime map in the elastic-memory case is new and could be useful for thinking about memory effects in impacts. The soft spot is validation. The abstract states the limits are recovered but gives no details on constitutive model checks, mesh convergence, or direct comparisons to Wagner or Hertz predictions. The stress-test point about contact-line dynamics is worth checking: if the numerics use a standard no-slip or implicit wetting condition without regularization, the liquid-end force history and rebound may not match known drop experiments even at Wi=0, which would make the claimed continuous transition less reliable. This paper is for researchers working on soft-matter or multiphase impacts who want a numerical bridge between drops and beads. A reader interested in regime diagrams or viscoelastic effects on contact will get value from it. It deserves a serious referee to look at the methods and confirm the quantitative matches hold up.

Referee Report

2 major / 2 minor

Summary. The manuscript uses direct numerical simulations of a viscoelastic sphere impacting a rigid non-wetting surface to study the dependence of impact force on the elasticity number El and Weissenberg number Wi. It claims recovery of the Newtonian liquid response as El → 0 or Wi → 0, elastic-solid behavior as Wi → ∞ with El ≠ 0, and the emergence of capillary-dominated, Wagner-scaling, and Hertz-scaling regimes in the elastic-memory limit, with a smooth Wagner-to-Hertz transition. The central narrative is that sweeping Wi from 0 to ∞ provides a continuous route from no-memory (liquid) to permanent-memory (elastic) materials.

Significance. If the simulations are shown to be quantitatively reliable, the work would supply a controlled, single-model route between liquid-drop and elastic-bead impact dynamics, unifying two classically separate theories and offering a parameter space in which memory effects can be tuned continuously. This framework could be relevant to a range of soft-matter and engineering applications where intermediate viscoelasticity governs force transmission.

major comments (2)

[Abstract and Results] Abstract and Results: the claim that Newtonian and elastic limits are recovered, together with the identification of distinct regimes, rests on simulation outputs; however, the manuscript provides no reported constitutive-model validation, mesh-convergence tests, error bars, or direct quantitative comparisons against Wagner or Hertz predictions, leaving the central claims without the necessary numerical evidence.
[Numerical Methods / Model] Numerical Methods / Model section: the assumption that the chosen viscoelastic constitutive relation plus the chosen interface treatment recovers the Newtonian liquid limit—including capillary-driven retraction and lift-off on a non-wetting substrate—requires explicit demonstration, because standard volume-of-fluid or level-set schemes without contact-line regularization commonly produce unphysical pinning that would alter the force history even at Wi = 0.

minor comments (2)

[Figures] Figure captions and text should explicitly state how the impact force is extracted from the simulations (e.g., integrated pressure or stress on the substrate) to allow direct comparison with analytic theories.
[Notation] Notation for El and Wi should be introduced once with clear definitions and then used consistently; occasional redefinition risks confusion when regimes are discussed.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed comments. We agree that the original manuscript would benefit from additional explicit numerical validation to strengthen the central claims. We have revised the manuscript to incorporate mesh-convergence studies, error estimates, quantitative comparisons to analytical limits, and explicit demonstrations of the Newtonian recovery. Our point-by-point responses follow.

read point-by-point responses

Referee: [Abstract and Results] Abstract and Results: the claim that Newtonian and elastic limits are recovered, together with the identification of distinct regimes, rests on simulation outputs; however, the manuscript provides no reported constitutive-model validation, mesh-convergence tests, error bars, or direct quantitative comparisons against Wagner or Hertz predictions, leaving the central claims without the necessary numerical evidence.

Authors: We acknowledge the validity of this observation. The revised manuscript now includes a new subsection on numerical methods validation. This addition reports mesh-convergence tests for the impact force (converging to within 4% under successive refinements), error bars derived from repeated simulations with varied initial perturbations, and direct quantitative comparisons of the simulated force histories and maximum forces against the Wagner and Hertz predictions in the respective limits, including relative deviations. Constitutive-model behavior is further validated by the explicit recovery of Newtonian dynamics as El or Wi approaches zero. revision: yes
Referee: [Numerical Methods / Model] Numerical Methods / Model section: the assumption that the chosen viscoelastic constitutive relation plus the chosen interface treatment recovers the Newtonian liquid limit—including capillary-driven retraction and lift-off on a non-wetting substrate—requires explicit demonstration, because standard volume-of-fluid or level-set schemes without contact-line regularization commonly produce unphysical pinning that would alter the force history even at Wi = 0.

Authors: We agree that an explicit demonstration is required and have addressed this in the revision. The updated Numerical Methods section now presents dedicated simulations at Wi = 0 (Newtonian limit) that exhibit the expected capillary-driven spreading, retraction, and lift-off on the non-wetting surface. Our interface treatment includes a contact-line regularization that suppresses artificial pinning, as confirmed by direct comparison of the simulated contact-line velocity and retraction dynamics to established benchmarks for inviscid drop impacts. revision: yes

Circularity Check

0 steps flagged

No significant circularity; results from independent DNS of standard model

full rationale

The paper performs direct numerical simulations of a viscoelastic sphere using a constitutive model with externally specified dimensionless groups El and Wi. Newtonian recovery is stated as the direct limit El→0 or Wi→0, and elastic-solid behavior as Wi→∞ with El≠0; these limits are properties of the chosen equations rather than redefinitions or fits. No load-bearing predictions reduce to fitted constants from the same dataset, no uniqueness theorems are imported via self-citation, and no ansatz is smuggled through prior work. The identified regimes (capillary, Wagner, Hertz) and continuous transition are outputs of the simulations, not inputs renamed or forced by construction. The derivation chain is therefore self-contained.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the assumption that a standard viscoelastic constitutive relation plus direct numerical simulation is sufficient to interpolate between liquid and solid impact without additional unmodeled effects; no new entities are postulated and no parameters are fitted to the impact data itself.

axioms (2)

domain assumption The viscoelastic model recovers Newtonian liquid response as El → 0 or Wi → 0
Explicitly stated in the abstract as the recovered limit behavior.
domain assumption The model obtains elastic-solid behavior in the limit Wi → ∞ and El ≠ 0
Explicitly stated in the abstract as the obtained limit behavior.

pith-pipeline@v0.9.0 · 5793 in / 1467 out tokens · 52828 ms · 2026-05-18T02:50:45.478340+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We recover the Newtonian liquid response as either El→0 or Wi→0, and obtain elastic-solid behavior in the limit Wi→∞ and El≠0. ... Fmax/(ρl V0² R0²) ∼ El^{2/5}
IndisputableMonolith/Foundation/AlphaCoordinateFixation.lean J_uniquely_calibrated_via_higher_derivative unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Using the Oldroyd-B constitutive model ... σe = Ec(A−I)

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.