Pith. sign in

REVIEW 4 major objections 5 minor 32 references

An SPH model with physically prescribed parameters for droplet dynamics on complex surfaces

T0 review · 4 major / 5 minor · reviewed 2026-08-01 · deepseek-v4-flash

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

Pith's one-line read SPH droplet forces now follow from surface tension and adhesion energy alone.

desk verdict A clean first-moment theorem gives SPH cohesion forces from macroscopic surface tension; the contact-angle validation is self-consistency, and dynamic claims need convergence support before they are quantitative. read the letter →

arxiv 2607.16670 v1 pith:EH2WK7O5 submitted 2026-07-18 physics.flu-dyn math-phmath.MP

Zhonghua Qiao , Yifan Wei , Xianmin Xu This is my paper · ORCID
classification physics.flu-dynmath-phmath.MP MSC 76M2876D4565M9976T10
keywords smoothedparticlehydrodynamicssurfacetensioncontactanglewettingdropletimpactYoung–Dupréequationpairwiseinteractionforceinterfacialenergy
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper proposes a single-phase SPH model for droplet dynamics in which the strength of every pairwise interfacial force is computed directly from two macroscopic quantities: the liquid surface tension γl and the liquid–solid work of adhesion Wsl. The key step is a theorem stating that the interfacial energy density across a planar interface equals one quarter of the first absolute moment of the pair potential. Inverting this relation gives the kernel amplitude Φll = -8γl / M1[W], so no empirical fitting of interaction strength is needed; the equilibrium contact angle follows from the Young–Dupré equation. Static wetting tests confirm cos θC = 2αadh - 1, and dynamic tests reproduce Tanner spreading, rolling rebound on patterned surfaces, and coalescence-induced jumping.

What carries the argument

Theorem 3.1, the moment representation: for two half-spaces separated by a planar interface, the interfacial potential energy density JQ[φ] equals π∫r³φ(r)dr = ¼ M1[φ], the first absolute moment of the potential over R³. This identity carries the argument: it converts the continuum surface energy into a kernel amplitude, giving Eq. (4.10) and (4.14). The rest of the model is standard SPH with a single-phase treatment: long-range attraction via the kernel, short-range repulsion absorbed into pressure.

What would settle it

Run the static-droplet test of Section 5.1 at several kernel support radii Rcho (e.g., 0.1, 0.2, 0.4 mm at fixed h/Δx) and measure cos θC at fixed αadh. If the measured angle shifts away from 2αadh - 1 as curvature-to-Rcho ratio increases, the planar moment calibration does not generalize to curved interfaces.

Watch

Extended reading notes

Core claim

The central claim is that macroscopic surface tension and work of adhesion can be prescribed exactly at the level of SPH pair potentials through the first absolute moment of the smoothing kernel. For a liquid–liquid interaction, the pair potential is set as -8γl W(r)/M1[W], and for liquid–solid adhesion as -4Wsl W(r)/M1[W]. These choices make the interfacial energy per unit area match γl and Wsl and make the equilibrium contact angle satisfy the Young–Dupré relation cos θC = Wsl/γl - 1. The paper validates this parameter-free calibration on static droplets and on three dynamic problems.

Load-bearing premise

The contact-angle calibration assumes the interface is locally planar with a constant number density; if the kernel support radius is not much smaller than the local curvature radius, the moment theorem no longer predicts the surface energy correctly, and the paper provides no convergence study in Rcho or resolution for the curved, dynamic interfaces it simulates.

Editorial extensions

If this is right

  • Static contact angles can be set by choosing one scalar, the adhesion ratio αadh = Wsl/2γl, without tuning force parameters.
  • The same calibrated kernel preserves the liquid surface tension γl and the work of adhesion Wsl, so dynamic wetting problems inherit the macroscopic energetics.
  • Single-phase treatment lowers cost because only the liquid is discretized; solid boundary effects enter through the calibrated adhesive force and a kernel-weighted hydrostatic pressure correction.
  • Reproduces the t^1/10 Tanner spreading law, rolling rebound on a patterned-wettability surface, and coalescence-induced jumping, suggesting the calibration transfers to transient, curved interfaces.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The moment theorem is derived for a planar interface; a convergence study in kernel radius Rcho near highly curved interfaces would determine how small Rcho must be relative to local curvature for the calibration to remain exact.
  • The same first-moment argument could be extended to liquid–liquid adhesion between two droplet phases, offering a parameter-free route for multi-component SPH simulations.
  • Because the calibration is energetic rather than force-curvature based, it may be more robust than CSF-type surface tension models for topological changes such as coalescence and pinch-off, but this advantage remains to be demonstrated at higher resolution.
  • The model's reliance on the first absolute moment suggests a testable extension where the kernel is chosen to match higher moments as well, which would control curvature-dependent corrections beyond the planar limit.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 5 minor

Summary. The manuscript introduces a single-phase SPH formulation for droplet dynamics on solid surfaces. It models liquid-liquid and liquid-solid interactions through effective kernel potentials whose prefactors are fixed by a first-moment relation derived for planar, constant-density interfaces (Theorem 3.1). In particular, Eq. (4.10) prescribes the liquid-liquid cohesion potential from the macroscopic surface tension coefficient γ_l and the kernel's first absolute moment, and Eq. (4.14) prescribes the liquid-solid adhesion from the work of adhesion W_sl. The Young-Dupré relation is rewritten in terms of an adhesion coefficient α_adh = W_sl/(2γ_l), giving cos θ_C = 2α_adh - 1. The paper validates the method with static contact-angle simulations, a Tanner-spreading test, droplet impact and rolling rebound on a patterned wettability surface, and coalescence-induced jumping.

Significance. If the proposed mapping remains valid at finite SPH resolution and on curved interfaces, the model would be a useful step toward reducing empirical parameter calibration in pairwise-force SPH. The proof of Theorem 3.1 is correct and transparent, and the derivation from the pair-potential energy through Eq. (3.24) is clean. The public availability of the source code and validation videos is also a strength. However, the current evidence does not yet establish the central claim that γ_l is physically prescribed rather than effectively tuned: the static contact-angle test is a self-consistency check, and the dynamic tests are qualitative or use a prescribed fitting exponent without convergence studies.

major comments (4)
  1. [§5.1, Eq. (4.22)] The static contact-angle validation is a self-consistency test, not an independent check of the parameter mapping. α_adh is defined as W_sl/(2γ_l), and W_sl is prescribed from the target angle through Young-Dupré; Eq. (4.22) is therefore an identity. The measured slope 2.081 vs 2 and intercept -1.104 vs -1 in Table 1 show a systematic deviation, yet no error bars, repeated-realization statistics, or resolution study are reported. Because cos θ_C depends only on the ratio W_sl/γ_l, uniform errors in the effective surface tension are invisible to this test. An independent measure of γ_l (e.g., Laplace pressure, capillary wave frequency, droplet oscillation) and a resolution study are needed.
  2. [Theorem 3.1 / Remark 3.2, Eq. (4.10)] The calibration fixes the first absolute moment for a planar, constant-density half-space interface. In the impact simulations the droplet radius is about 2.5 mm while R_cho = 0.20 mm, and during spreading/recoil the local interface curvature is much sharper; the film can become comparable to or thinner than R_cho. Remark 3.2 merely asserts the asymptotic regime. No R_cho- or h-convergence study is supplied. Since the first-moment matching does not control curvature corrections of order R_cho/R, the effective surface tension produced by Eq. (4.16) may differ from the prescribed γ_l at the resolutions used. Please quantify this with a convergence study or provide a curved-interface correction/error estimate.
  3. [§5.2, Eq. (5.2)] The Tanner-law test prescribes the exponent 0.1 in the fitting function R(t) = 2(t - t_min)^0.1 rather than fitting it. The text states that 'the fitted exponent' is 0.1, but the only free parameter is t_min. With t_min = 15 ms and no sensitivity or uncertainty analysis, the data do not provide evidence for the 1/10 exponent. Report a free-exponent fit with confidence intervals, or compare models with different fixed exponents using an appropriate information criterion.
  4. [§5.3–5.4] The comparisons for rolling rebound and coalescence jumping are qualitative: visual shape agreement is claimed, but no quantitative metrics are reported (e.g., spreading factor, contact-line position time series, jumping velocity, rebound angle), and there are no error estimates or grid-convergence checks. Given the transient, strongly curved interfaces in these tests, they cannot substitute for the planar-interface validation and do not yet support the claim of 'good agreement' with the experimental references [2,31].
minor comments (5)
  1. [Eq. (4.10)] Specify how M1[W_cs_{Rcho}] is evaluated—analytically or by numerical quadrature—and report its value or the quadrature error.
  2. [§5.1] Please state how contact angles are extracted from simulated droplet profiles (e.g., circular fit, tangent at the contact line) and report the measurement uncertainty.
  3. [§4.2, Eq. (4.16)] The momentum equation uses the cubic-spline kernel for interaction forces and the Wendland kernel for SPH interpolation. Clarify why these two kernels are used and whether the choice affects the moment-calibration consistency.
  4. [Abstract / Conclusion] The phrase 'physically prescribed parameters' is stronger than what is demonstrated, since R_cho, h, Δx, and the artificial viscosity coefficient remain numerical parameters. Consider softening to 'reduced empirical calibration of interaction prefactors.'
  5. [Figure 6] Show error bars and define how the spreading radius is measured. Log-log plots with a fitted t_min can be misleading without these details.

Circularity Check

3 steps flagged · score 6.0 of 10

Static contact-angle "validation" is Young–Dupré with the input W_sl renamed; Tanner-law fit hard-codes the exponent; the first-moment potential mapping itself is non-circular.

  1. self definitional [§3.2 Eq. (3.22); §4.2 Eqs. (4.21)–(4.22); §5.1]
    "α adh := Wsl 2γl . Substituting this definition into Eq. (3.22) gives (4.22) cosθ C = 2αadh −1. ... For each prescribed value of αadh, the equilibrium droplet profile is measured and the corresponding contact angle is extracted."

    The "theoretical relation" is obtained by substituting the definition of α_adh into the Young–Dupré equation. The simulations then prescribe W_sl (equivalently α_adh) to target specific angles, e.g. W_sl = (1+√3/2)γ_l for 30°, and measure the resulting angle. Agreement with Eq. (4.22) is therefore a self-consistency check of the prescribed adhesive force, not an independent prediction. It also cannot test the absolute liquid surface tension, since cosθ_C depends only on the ratio W_sl/γ_l.

  2. self definitional [§5.1, Eq. (5.1)]
    "For this symmetric configuration, Eqs. (3.16)–(3.17) simplify to cosθ eq = Wl1l2/(2γl1) −1. ... When W l1l2 = γ l1, the equilibrium contact angle is θ eq = 2π/3."

    This two-droplet test is the same construction: W_l1l2 is an input chosen to produce the reported contact angles, and Eq. (5.1) is just the Young–Dupré-type relation written in terms of that input. The simulated equilibrium shape consequently reproduces the prescribed angle by design, so it provides no independent validation of the model's contact-angle physics.

1 more flagged steps
  1. fitted input called prediction [§5.2, Eq. (5.2) and following paragraph]
    "the numerical data were fitted using the function (5.2) R(t) = 2×(t−tmin)^0.1 ... The fitted exponent of 0.1 agrees well with the theoretical exponent predicted by Tanner’s law."

    The fitting function hard-codes the Tanner exponent 0.1; only the time offset t_min is adjusted. The subsequent statement that 'the fitted exponent of 0.1 agrees' with Tanner's law is therefore true by construction. No independent exponent is estimated from the data, so the exercise cannot confirm or refute the t^{1/10} scaling.

full rationale

The core derivation—Theorem 3.1 and Eq. (4.10)—is self-contained. For a planar half-space interface, J_Q[φ] = (1/4)M_1[φ] is proven. Combining this with W_ll = 2γ_l and W_sl = -e_φsl yields Eq. (4.10), which sets the kernel amplitude so that the first absolute moment of the liquid–liquid kernel equals -8γ_l. This is a legitimate constitutive construction: it prescribes the SPH force from the macroscopic surface tension and does not depend on contact-angle data. It is not itself a circular prediction, even though the relation is by construction. The circularity lies in the validation claims. α_adh is defined as W_sl/(2γ_l), so Eq. (4.22) is an algebraic identity with the Young–Dupré equation; the static-wetting simulations prescribe W_sl to target exactly the reported angles (30°, 90°, 150°). The agreement in §5.1 is thus a self-consistency test, not an independent prediction. The two-droplet test in Eq. (5.1) is the same tautology. Section 5.2's Tanner-law 'verification' is vacuous because the exponent 0.1 is imposed in the fitting function and then reported as confirmed. These issues do not infect the first-moment mapping, which is mathematically independent of the contact-angle validation. No load-bearing self-citation or imported uniqueness theorem is present; the cited prior work by the authors is not used as the basis for the main result. Overall, the central parameter-mapping derivation is non-circular, but several announced validations reduce to their prescribed inputs, so the paper is partially circular.

Assumptions & free parameters 6 free parameters · 4 assumptions · 1 invented entities

The model's parameter economy comes from replacing the microscopic potential's first moment with γl and Wsl. What remains free or assumed is the kernel support Rcho, the SPH resolution, and especially the prescribed Wsl (or α_adh) that actually sets the contact angle. The contact-angle validation is therefore a self-consistency check rather than an independent prediction.

free parameters (6)
  • R_cho (kernel support radius for interaction potential) = 0.20 mm
    Chosen cutoff for the cohesive/adhesive kernel in §5.1; not determined by physics and no sensitivity study is provided.
  • smoothing length h = 0.11 mm
    SPH smoothing length chosen in §5.1; numerical resolution parameter.
  • particle spacing Δx = 0.05 mm
    Resolution parameter in §5.1.
  • liquid–solid adhesion work W_sl (or α_adh) = varied per case (e.g., W_sl = γl, (1±√3/2)γl, 2γl)
    Prescribed to target desired contact angles in §5.1 and used in the dynamic cases; this input determines the contact angle, so it is the key empirical knob.
  • t_min in Tanner fit = 15 ms
    Time offset fitted to the spreading data in §5.2.
  • artificial viscosity coefficient α = not stated
    Monaghan artificial viscosity in Eq. (4.18); value not reported, so an unspecified numerical parameter.
assumptions (4)
  • domain assumption The two phases are homogeneous with constant number densities n1, n2, and the interface is planar on the scale of the kernel support.
    Used to derive Eq. (3.14) and the first-moment calibration; Remark 3.2 extends it to curved interfaces only when the local curvature radius is much larger than Rcho.
  • ad hoc to paper The SPH kernel represents the long-range attractive tail of the intermolecular potential; the short-range repulsive part contributes negligibly to the first absolute moment and is represented by pressure.
    Stated in Remark 4.1 and §4.2; if pressure does not reproduce the short-range structure, the surface tension calibration may break.
  • domain assumption Young–Dupré and Dupré relations hold, and adsorption at solid–gas and liquid–gas interfaces is neglected so that γ_s ≈ γ_sg and γ_l ≈ γ_lg.
    Used in §3.2 to derive Eq. (3.24). Standard surface-chemistry approximation.
  • domain assumption Single-phase representation: the gas phase has no dynamical role; liquid–gas surface tension is produced entirely by liquid–liquid cohesion.
    Adopted for efficiency; ignores air entrainment/cushioning that can matter in droplet impact and coalescence jumping.
invented entities (1)
  • Effective kernel-based cohesion/adhesion potentials Φll and Φls
    purpose: Replace microscopic molecular pair potentials in SPH; amplitude fixed by surface tension γl and adhesion work Wsl via first-moment matching.
    These are numerical modeling constructs, not independently measurable physical entities; they have no falsifiable handle outside the SPH simulation.

how reviews work

0 comments
Cite this review

Pith. "Pith review of An SPH model with physically prescribed parameters for droplet dynamics on complex surfaces." pith.science (2026). https://pith.science/paper/EH2WK7O5

@misc{pith2026260716670,
  author       = {Pith},
  title        = {Pith review of: An SPH model with physically prescribed parameters for droplet dynamics on complex surfaces},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/EH2WK7O5}},
  note         = {Machine review of arXiv:2607.16670}
}
read the original abstract

Numerical simulation of droplet dynamics on complex surfaces with varying wettability is of great significance to both engineering applications and fundamental research. However, existing numerical methods still face challenges in accurately capturing interfacial interactions while preserving physical consistency and computational efficiency. In this work, a physically grounded and efficient smoothed particle hydrodynamics (SPH) model is developed for droplet dynamics simulation. To reduce computational cost, a single-phase droplet modeling strategy is employed. At the interface, long-range interactions are approximated using the SPH kernel function, whereas short-range interactions are represented through pressure. Based on this treatment, an explicit relationship between the intermolecular potential energy and the macroscopic surface tension coefficient is further established, thereby reducing reliance on empirical parameter calibration. The proposed method is first validated through static wetting simulations, where the equilibrium contact angles agree well with the Young--Dupr\'e equation. Further simulations of wetting and droplet impact demonstrate that the method is capable of capturing complex dynamic wetting behaviors.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

32 extracted references

  1. [1]

    Contact angle, wetting, and adhesion: a critical review , author=. J. Adhes. Sci. Technol. , volume=. 1992 , publisher=

    1992

  2. [2]

    Monaghan, J. J. , journal=. 1992 , publisher=

    1992

  3. [3]

    and Rogers, B

    Gomez-Gesteira, M. and Rogers, B. D. and Dalrymple, R. A. and Crespo, A. J. C. , journal=. State-of-the-art of classical. 2010 , publisher=

    2010

  4. [4]

    The limit of droplet rebound angle , author=. Nat. Commun. , volume=. 2025 , publisher=

    2025

  5. [5]

    The spreading of silicone oil drops on horizontal surfaces , author=. J. Phys. D: Appl. Phys. , volume=

  6. [6]

    Role of volatility and thermal properties in droplet spreading: a generalisation to

    Wang, Zhenying and Karapetsas, George and Valluri, Prashant and Inoue, Chihiro , journal=. Role of volatility and thermal properties in droplet spreading: a generalisation to. 2024 , publisher=

    2024

  7. [7]

    Interface Sci

    The apolar and polar properties of liquid water and other condensed-phase materials , author=. Interface Sci. Technol. , volume=. 2008 , publisher=

    2008

  8. [8]

    Self-propelled dropwise condensate on superhydrophobic surfaces , author=. Phys. Rev. Lett. , volume=. 2009 , publisher=

    2009

Show all 32 references
  1. [9]

    Pairwise force smoothed particle hydrodynamics model for multiphase flow: surface tension and contact line dynamics , author=. J. Comput. Phys. , volume=. 2016 , publisher=

    2016

  2. [10]

    A continuum method for modeling surface tension , author=. J. Comput. Phys. , volume=. 1992 , publisher=

    1992

  3. [11]

    Modeling of surface tension and contact angles with smoothed particle hydrodynamics , author=. Phys. Rev. E , volume=. 2005 , publisher=

    2005

  4. [12]

    Wetting: intrinsically robust hydrophobicity , author=. Nat. Mater. , volume=. 2013 , publisher=

    2013

  5. [13]

    and Oger, G

    Vergnaud, A. and Oger, G. and Le Touz. Comput. Methods Appl. Mech. Engrg. , volume=. 2022 , publisher=

    2022

  6. [14]

    Drop dynamics in the inkjet printing process , author=. Curr. Opin. Colloid Interface Sci. , volume=. 2018 , publisher=

    2018

  7. [15]

    Fundamental fluid dynamics challenges in inkjet printing , author=. Annu. Rev. Fluid Mech. , volume=. 2022 , publisher=

    2022

  8. [16]

    Nanotechnol

    Micro/nano manufacturing aircraft surface with anti-icing and deicing performances: an overview , author=. Nanotechnol. Rev. , volume=. 2023 , publisher=

    2023

  9. [17]

    An icephobic coating for aircraft anti-icing system , author=. J. Mater. Sci. , volume=. 2023 , publisher=

    2023

  10. [18]

    Agronomy , volume=

    Drop size measurement techniques for agricultural sprays: a state-of-the-art review , author=. Agronomy , volume=. 2023 , publisher=

    2023

  11. [19]

    Global Challenges , volume=

    Sustained agricultural spraying: from leaf wettability to dynamic droplet impact behavior , author=. Global Challenges , volume=. 2023 , publisher=

    2023

  12. [20]

    Review of the role of surfactant dynamics in drop microfluidics , author=. Adv. Colloid Interface Sci. , volume=. 2023 , publisher=

    2023

  13. [21]

    Droplet-based microfluidics , author=. Nat. Rev. Methods Primers , volume=. 2023 , publisher=

    2023

  14. [22]

    A smoothed particle hydrodynamics model for droplet and film flow on smooth and rough fracture surfaces , author=. Adv. Water Resour. , volume=. 2013 , publisher=

    2013

  15. [23]

    and Menon, Dileep and Govind, Rahul and Sanka, Suraj and Sebastian, Amal S

    Ramachandran, Prabhu and Bhosale, Aditya and Puri, Kunal and Negi, Pawan and Muta, Abhinav and Dinesh, A. and Menon, Dileep and Govind, Rahul and Sanka, Suraj and Sebastian, Amal S. and Sen, Ananyo and Kaushik, Rohan and Kumar, Anshuman and Kurapati, Vikas and Patil, Mrinalgou...

    2021

  16. [24]

    Physics of droplet impact on various substrates and its current advancements in interfacial science: a review , author=. J. Appl. Phys. , volume=. 2023 , publisher=

    2023

  17. [25]

    Multicomponent droplet evaporation in a geometric volume-of-fluid framework , author=. J. Comput. Phys. , volume=. 2024 , publisher=

    2024

  18. [26]

    A unified variational framework and lattice

    Qiao, Zhonghua and Xu, Xianmin and Yang, Xuguang and Zhang, Yuze , journal=. A unified variational framework and lattice. 2026 , publisher=

    2026

  19. [27]

    Effect of the surface peak-valley features on droplet splash dynamics , author=. Phys. Fluids , volume=. 2023 , publisher=

    2023

  20. [28]

    and Fox, Patrick J

    Morris, Joseph P. and Fox, Patrick J. and Zhu, Yi , journal=. Modeling low. 1997 , publisher=

    1997

  21. [29]

    Smoothed particle hydrodynamics , author=. Rep. Progr. Phys. , volume=. 2005 , publisher=

    2005

  22. [30]

    Monaghan, J. J. , journal=. Simulating free surface flows with. 1994 , publisher=

    1994

  23. [31]

    A generalized wall boundary condition for smoothed particle hydrodynamics , author=. J. Comput. Phys. , volume=. 2012 , publisher=

    2012

  24. [32]

    A high-order

    Qin, Jian and Wu, Jie and Qin, Qiushuo , journal=. A high-order. 2026 , volume =

    2026

Pith tools

Reviewed August 1, 2026 · model on record in the stance chip above.