Pith. sign in

REVIEW 2 major objections 4 minor 90 references

Interface-dominated sliding compound drops

T0 review · 2 major / 4 minor · reviewed 2026-08-02 · deepseek-v4-flash

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

Pith's one-line read A two-liquid compound drop slides up to about twice as fast when the slower, flatter liquid is at the front, and the difference traces to where sliding energy is dissipated.

desk verdict Competent mesoscopic study of sliding compound drops with a robust velocity mechanism; the angle-extraction caveat is real but surmountable. read the letter →

arxiv 2603.26601 v1 pith:25B3GRHL submitted 2026-03-27 physics.flu-dyn

Dominik Thy , Jan Diekmann , Uwe Thiele This is my paper · ORCID
classification physics.flu-dyn PACS 47.55.D68.08.Bc
keywords slidingcompounddropsdropconfigurationdynamiccontactanglesNeumannlawYoungdissipationprofilesaddle-nodebifurcationtwo-layerthin-filmmodel
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 aims to establish that, for a compound drop made of two immiscible partially wetting liquids sliding down an incline, the lateral order of the two liquid lobes controls the sliding speed: in every parameter range studied, the configuration with the flatter, slower liquid at the front (2-1) moves up to about twice as fast as the mirror configuration (1-2). Using a mesoscopic two-layer hydrodynamic model in full-curvature form, the authors trace this ordering to the lateral dissipation profile: energy is dissipated most strongly in the Young region where the liquid of smaller equilibrium contact angle meets the solid, and the advancing/receding behavior of that contact angle shortens or lengthens the controlling drop. They further show that stationary sliding compound drops exist only within finite parameter ranges bounded by saddle-node bifurcations, and that beyond these boundaries the drop either overtakes and changes configuration or enters a time-periodic fusion-overtaking-splitting cycle.

What carries the argument

The authors employ a mesoscopic two-layer hydrodynamic thin-film model in full-curvature form: a gradient dynamics for the two interface height profiles (liquid-liquid and liquid-gas), with a mobility matrix derived from the Stokes equations and wetting energies that cover the full macroscopic interface-energy parameter space. This model lets them extract dynamic Young and Neumann angles from the slopes of the height profiles at inflection points (contact lines are smooth transitions rather than sharp boundary conditions), compute stationary sliding states and their stability by numerical continuation, and decompose the total dissipation into per-liquid lateral profiles. The dissipation anal

What would settle it

Form a compound drop of two immiscible partially wetting liquids with different equilibrium contact angles on a smooth inclined solid, in both lateral orderings (1-2 and 2-1), and measure the sliding velocities over a range of inclinations and viscosity ratios: if the configuration with the flatter liquid at the front is not consistently the faster one, or if no saddle-node limit or splitting is observed at high inclination, the central claim fails. A complementary numeric check is to repeat the continuation with different domain sizes to confirm the fusion-overtaking-splitting cycle is not an

Watch

Extended reading notes

Core claim

The central claim is that the sliding velocity of a stationary compound drop is determined by the liquid with the smaller equilibrium contact angle, and that the configuration-dependent speed difference results from the opposite behavior of advancing and receding dynamic contact angles on that liquid. In the 2-1 configuration—where the flatter, slower liquid (liquid 1) is at the front—its leading Young angle is advancing, so it increases with velocity and shortens the controlling drop; in the 1-2 configuration the same liquid trails, its receding angle decreases and lengthens the drop. At equal inclination, viscosity ratio, and volume ratio, the 2-1 configuration is therefore always the fast

Load-bearing premise

The load-bearing assumption is that the mesoscopic two-layer thin-film model faithfully converts substrate inclination and wetting energies into macroscopically correct dynamic Young and Neumann angles during sustained sliding; if that mapping fails under strong driving (e.g., at high inclination), the predicted velocity ordering, dissipation peaks, and saddle-node boundaries would be model artifacts rather than physical predictions.

Editorial extensions

If this is right

  • If the central claim is correct, the sliding speed of a two-liquid compound drop is not a fixed material property but depends on which liquid is placed in front; preparing the desired ordering becomes a way to control transport speed.
  • The finite existence ranges bounded by saddle-node bifurcations imply that there is a maximum inclination (and velocity) beyond which a stationary sliding compound drop cannot persist; applications relying on steady droplet transport must stay below these thresholds or accept periodic splitting/overtaking.
  • The dissipation argument predicts that the speed difference between configurations grows with the contrast in equilibrium contact angles of the two liquids and vanishes when the liquids are identical—a testable scaling.
  • The predicted fusion-overtaking-splitting cycle provides a self-contained periodic droplet-release mechanism on an incline with periodic boundaries; on an open domain it would produce a train of alternating droplets.
  • The observed nonmonotonic velocity of a drop sliding on a liquid substrate of varying thickness offers a concrete route to compare against experiments on adaptive or brush-like substrates.

Reading between the lines

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

  • The mechanism—velocity set by the most dissipative (smallest-contact-angle) component, with front/back placement controlling that component's length—likely generalizes to other sliding composite objects, such as drops containing particles, Janus drops, or soft composites; a natural rule to test is that the slowest part leads to faster overall motion.
  • A direct experimental test could use two immiscible oils on a smooth inclined plate in a microfluidic channel: measuring the velocity ratio between the two orderings as a function of viscosity and volume ratios would verify both the ordering and the predicted factor-of-two bound.
  • The paper reports that the Neumann region rotates almost as a solid body while the Neumann residues stay small; this implies that macroscopic Neumann constructions might remain usable in simpler models of driven compound drops even when angles change strongly.
  • Because the model assumes homogeneous smooth substrates, a next step would be to test whether chemical or topographical patterning—which modifies local Young angles—could extend the existence range or even invert the velocity ordering by strengthening dissipation at the other contact line.
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

2 major / 4 minor

Summary. The paper extends the mesoscopic two-layer thin-film model of Ref. [15] to inclined substrates and studies, in one horizontal dimension, (i) a drop of one liquid sliding on a layer of the other liquid and (ii) stationary sliding compound drops in two configurations (1-2 and 2-1). For the compound drops, the authors map velocity, dynamic Young and Neumann angles, and dissipation profiles as functions of inclination β, viscosity ratio η, and volume ratio ν. They report that the 2-1 configuration slides up to about twice as fast as 1-2 at identical parameters, explain this by the lateral distribution of dissipation and by opposite trends of advancing/receding dynamic contact angles, and show that stable sliding compound drops exist only in finite ranges bounded by saddle-node bifurcations. Beyond those ranges, configuration changes or time-periodic fusion-overtaking-splitting cycles are observed. The results are obtained with path continuation (pde2path) and time simulations (oomph-lib), with data/code made publicly available.

Significance. If the results hold, the paper provides the first systematic mesoscopic characterization of driven compound drops with dynamic Neumann and Young angles, and it offers a falsifiable prediction (2-1 faster than 1-2) together with a mechanistic explanation in terms of dissipation localization. The work builds on a previously published, parameter-free gradient-dynamics model; velocities, angles, and dissipation are computed from the governing equations rather than fitted, which is a strength. The availability of reproducible data and code is a further positive. The main significance is as a benchmark for future sharp-interface, DNS, or experimental studies of sliding compound drops.

major comments (2)
  1. [§4.1, Fig. 4; §2.2] The central mechanistic claim—that the 2-1 configuration is faster because ϑ13 is an increasing advancing angle in 2-1 and a decreasing receding angle in 1-2, and that this 'ultimately results' in the velocity difference—rests entirely on the mesoscopic angle extraction from slopes at inflection points. No independent validation of this mapping is provided for strongly driven compound drops. Table 1 shows ~8% deviations from macroscopic Young angles even at equilibrium, and footnote 2 states that Neumann residues, while small, increase near the bifurcations—precisely the regime where the angle-based mechanism is invoked (e.g., near β_c21 in Fig. 3). The authors should provide a convergence study with respect to domain size, spatial/temporal resolution, and adsorption-layer thickness, and compare at least one representative dynamic angle against an independent method (e.g., sharp-interfac
  2. [§2.2; Figs. 3, 5, 6] No numerical tolerances or convergence studies are reported for either the pde2path continuation or the oomph-lib time simulations. Quantitative claims such as β_c21 ≈ 2.15×10⁻⁵, β_c12 ≈ 2.63×10⁻⁵, η_c21 ≈ 3.68, η_c12 ≈ 0.791, η_c12b ≈ 5.2, and ν_c12 ≈ 121.6 are quoted without any estimate of numerical error. Given that these saddle-node positions are used to support the existence-range claims, the authors should report continuation step sizes, Newton tolerances, finite-element adaptivity thresholds, and a representative mesh/timestep convergence test for at least one bifurcation point.
minor comments (4)
  1. [§2.2] The text states 'we effectively nondimensionalize the dynamic equations by setting η1 = γ12 = h_a1 = 1', but h_a1 is not introduced in the main text. Please define it (it appears to be an equilibrium adsorption-layer thickness) and explain its role in the wetting energy.
  2. [Fig. 2 caption] The text refers to an 'inset of Fig. 2(a)' concerning the nonmonotonic U(¯h12) dependence, but the caption does not describe any inset. Please either label the inset or remove the reference.
  3. [§5] In Fig. 7(b), the space-time diagram is called a 'top view', which is misleading for a one-dimensional substrate plot; a 'side view' or 'space-time diagram' would be clearer. Also, the threshold '1.2' used in the binary thickness visualization is not given units.
  4. [Appendix A] The sentence 'the drop of liquid 1 has a smaller equilibrium contact angle than the drop of liquid 2 ϑ13 < ϑ23' should be punctuated for readability, e.g., '... liquid 2: ϑ13 < ϑ23'.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the reported velocities, angles, and bifurcations are computed outputs of the stated PDE model with fixed parameters, not fitted inputs or self-cited uniqueness constraints.

full rationale

The paper's central claims—the 2-1 configuration sliding faster than 1-2, the dynamic Young/Neumann angle trends, the dissipation profiles, and the saddle-node-limited existence ranges—are all generated by numerically solving the coupled evolution equations (Eq. (1)) with the mobility (Eq. (2)), wetting energy (Eq. (4)), and gravitational potential (Eq. (9)). No parameter is tuned to reproduce the target velocity ordering or angle trends: the interface energies are fixed in Section 2.2 as (γ13, γ23, γs2, γs3) = (1.6, 0.8, 0.85, 1.4), and then β, η, and ν are scanned as independent control parameters. The dynamic angles are extracted from the simulated profiles via 'the slopes at the inflection points of the corresponding numerically obtained profiles' (Section 4), i.e., they are diagnostics, not imposed boundary conditions. The dissipation explanation is a post-hoc analysis of the same computed fields; while D = −βUV follows from energy balance, the spatial distribution of dissipation and its connection to the Young region of liquid 1 is an emergent result, not an input. Citations to the authors' prior work, especially Ref. [15], supply the mesoscopic model and the consistency relations linking meso- and macroscopic energies, but the governing equations are restated in the present paper, and no uniqueness theorem from prior work is invoked to forbid alternative explanations. The lack of an experimental or sharp-interface benchmark for sliding compound drops is a validation/correctness limitation (and is acknowledged indirectly by the ~8% finite-size deviation in Table 1), but it does not amount to a circular reduction of the predictions to the inputs. Hence no circular step is identifiable, and the paper's derivation chain is self-contained; the mild self-citation reliance on the model reference is not load-bearing in a circular sense.

Assumptions & free parameters 2 free parameters · 6 assumptions · 0 invented entities

The paper introduces no new physical entities. Its central claims rest on a previously published mesoscopic model (Ref. [15]), a fixed set of interface-energy parameters, and several domain assumptions: long-wave mobility in a full-curvature setting, neglect of hydrostatic gravity, equal densities, 1D periodic geometry, and the inflection-point definition of contact angles.

free parameters (2)
  • Fixed macroscopic interface energies (gamma13, gamma23, gammas2, gammas3) = (1.6, 0.8, 0.85, 1.4) in nondimensional units
    Chosen by hand in Section 2.2; all results are obtained for this single set of energies, and the stability boundaries and velocity ordering depend on it.
  • Wetting-energy/Hamaker parameters (A1, A2, A3) = approximately (0.5, 0.67, 0.67)
    Resulting from the chosen macroscopic energies and the meso-macro consistency relations of Ref. [15]; they control adsorption-layer thicknesses and contact angles.
assumptions (6)
  • domain assumption The mobility matrix Eq. (2), derived via long-wave approximation from Navier-Stokes with no-slip at the substrate, remains valid in the full-curvature variant used here.
    Section 2.1; the model uses exact metric factors for curvature but keeps the long-wave mobility. If this mixed approximation fails under strong sliding, computed velocities and dissipation are inaccurate.
  • domain assumption The wetting-energy construction Eq. (4) with the consistency relations of Ref. [15] spans the full six-interface-energy parameter space and recovers the macroscopic Neumann and Young laws at equilibrium.
    Section 2.1; the paper's dynamic Young/Neumann angle interpretation depends on this mapping between mesoscopic and macroscopic laws.
  • domain assumption The hydrostatic (cos-alpha) contribution to gravity is negligible; only the effective inclination parameter beta = rho2 g sin(alpha) enters the dynamics.
    Section 2.1: the paper explicitly neglects the cos-alpha flattening terms, justified only for small slowly sliding drops or vertical substrates.
  • domain assumption The two liquids have identical mass densities, so rho_r = rho1/rho2 - 1 = 0.
    Section 2.2; this removes buoyancy-like terms that would otherwise couple the two height fields differently.
  • domain assumption A one-dimensional periodic domain of length L=1000 or 2000 adequately represents sliding compound drops; 2D instabilities such as pearling and satellite-drop formation are not considered.
    Section 2.2 and Section 6; the authors explicitly defer 2D extensions to future work, so the central claims are restricted to the 1D setting.
  • domain assumption Dynamic Neumann and Young angles obtained from slopes at inflection points of the mesoscopic profiles correspond to the macroscopic dynamic contact angles.
    Section 2.2 and Section 4; the interpretation of the dynamic angle trends depends on this correspondence, which is inherited from Ref. [15].

how reviews work

0 comments
Cite this review

Pith. "Pith review of Interface-dominated sliding compound drops." pith.science (2026). https://pith.science/paper/25B3GRHL

@misc{pith2026260326601,
  author       = {Pith},
  title        = {Pith review of: Interface-dominated sliding compound drops},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/25B3GRHL}},
  note         = {Machine review of arXiv:2603.26601}
}
read the original abstract

We investigate compound drops composed of two immiscible nonvolatile partially wetting liquids that slide down an inclined homogeneous smooth solid substrate based on a mesoscopic hydrodynamic two-layer model in full-curvature formulation. First, drops of one liquid stationarily sliding on a layer of the other liquid are briefly investigated with a focus on the dependence of drop velocity and interface profiles on inclination and mean thickness of the adaptive substrate. Then, stationary sliding compound drops are studied with a focus on the dependence of their configuration, velocity, dynamic Young and Neumann angles on three control parameters, namely, the inclination, the volume ratio and the viscosity ratio. The reasons for the encountered dependence of the velocity on configuration are clarified based on a discussion of the lateral dissipation profile. Finally, we briefly consider the time-periodic fusion-overtaking-splitting behavior found outside the existence range of the stationary sliding compound drops as determined by saddle-node bifurcations.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

90 extracted references · 29 canonical work pages

  1. [15]

    Phys Rev Fluids 10:024002

    Diekmann J, Thiele U (2025) Mesoscopic hydrodynamic model for spreading, sliding and coarsening compound drops. Phys Rev Fluids 10:024002. https://doi.org/10.1103/ PhysRevFluids.10.024002

    2025

  2. [1]

    Soft Matter 8:6333–6349

    Amarandei G, Beltrame P, Clancy I, et al (2012) Pattern formation induced by an electric field in a polymer-air-polymer thin film sys- tem. Soft Matter 8:6333–6349. https://doi.org/ 10.1039/c2sm25273b

    work page doi:10.1039/c2sm25273b 2012

  3. [2]

    Annu Rev Fluid Mech 52:285–308

    Andreotti B, Snoeijer JH (2020) Statics and dynamics of soft wetting. Annu Rev Fluid Mech 52:285–308. https://doi.org/10.1146/ annurev-fluid-010719-060147 17 Preprint– contact: u.thiele@uni-muenster.de – www.uwethiele.de – March 30, 2026

    2020

  4. [3]

    Phys Rev E 109:024801

    Areshi M, Tseluiko D, Thiele U, et al (2024) Binding potential and wetting behavior of binary liquid mixtures on surfaces. Phys Rev E 109:024801. https://doi.org/10.1103/PhysRevE. 109.024801

    doi:10.1103/physreve 2024

  5. [4]

    Ind Eng Chem Res 44:1259–1272

    Bandyopadhyay D, Gulabani R, Sharma A (2005) Stability and dynamics of bilayers. Ind Eng Chem Res 44:1259–1272

    2005

  6. [5]

    Langmuir 30:4033–4039

    Ben Said M, Selzer M, Nestler B, et al (2014) A phase-field approach for wetting phe- nomena of multiphase droplets on solid sur- faces. Langmuir 30:4033–4039. https://doi.org/ 10.1021/la500312q

    work page doi:10.1021/la500312q 2014

  7. [6]

    J Fluid Mech 911:A13

    Bertin V , Lee CL, Salez T, et al (2021) Cap- illary levelling of immiscible bilayer films. J Fluid Mech 911:A13. https://doi.org/10.1017/ jfm.2020.1045

    arXiv 2021

  8. [7]

    Eur Phys J E 36:87

    Bommer S, Cartellier F, Jachalski S, et al (2013) Droplets on liquids and their journey into equi- librium. Eur Phys J E 36:87. https://doi.org/10. 1140/epje/i2013-13087-x

    2013

Show all 90 references
  1. [8]

    Colloid Surf A-Physicochem Eng Asp 345:163–165

    Bormashenko E (2009) Young, Boruvka- Neumann, Wenzel and Cassie-Baxter equations as the transversality conditions for the vari- ational problem of wetting. Colloid Surf A-Physicochem Eng Asp 345:163–165. https://doi.org/10.1016/j.colsurfa.2009.04.054

    2009 doi

  2. [9]

    Lang- muir 9:3682–3690

    Brochard-Wyart F, Martin P, Redon C (1993) Liquid/liquid dewetting. Lang- muir 9:3682–3690. https://doi.org/10.1021/ la00036a053

    1993

  3. [10]

    Cox RG (1986) The dynamics of the spread- ing of liquids on a solid surface. Part 1. viscous flow. J Fluid Mech 168:169–194. https://doi.org/ 10.1017/S0022112086000332

    1986 doi

  4. [11]

    J Colloid Interface Sci 303:503–

    Craster RV , Matar OK (2006) On the dynamics of liquid lenses. J Colloid Interface Sci 303:503–

    2006

  5. [12]

    Rev Mod Phys 81:1131–1198

    Craster RV , Matar OK (2009) Dynam- ics and stability of thin liquid films. Rev Mod Phys 81:1131–1198. https: //doi.org/10.1103/RevModPhys.81.1131

    2009 doi

  6. [13]

    The equations of lubrication approximation

    Danov KD, Paunov VN, Alleborn N, et al (1998) Stability of evaporating two-layered liq- uid film in the presence of surfactant - I. The equations of lubrication approximation. Chem Eng Sci 53:2809–2822. https://doi.org/10.1016/ s0009-2509(98)00098-0

    1998

  7. [14]

    Linear analysis

    Danov KD, Paunov VN, Stoyanov SD, et al (1998) Stability of evaporating two-layered liq- uid film in the presence of surfactant - II. Linear analysis. Chem Eng Sci 53:2823–2837. https: //doi.org/10.1016/s0009-2509(98)00099-2

    1998 doi

  8. [16]

    Langmuir 37(8):2575–2585

    Diez JA, Gonz ´alez AG, Garfinkel DA, et al (2021) Simultaneous decomposition and dewet- ting of nanoscale alloys: A comparison of exper- iment and theory. Langmuir 37(8):2575–2585. https://doi.org/10.1021/acs.langmuir.0c02964

    2021 doi

  9. [17]

    Phys Rev Fluids 1:073901

    Engelnkemper S, Wilczek M, Gurevich SV , et al (2016) Morphological transitions of sliding drops - dynamics and bifurca- tions. Phys Rev Fluids 1:073901. https: //doi.org/10.1103/PhysRevFluids.1.073901, http://arxiv.org/abs/1607.05482

    2016 arXiv

  10. [18]

    Macromolecules 54:584–

    Etha SA, Desai PR, Sachar HS, et al (2021) Wet- ting dynamics on solvophilic, soft, porous, and responsive surfaces. Macromolecules 54:584–

    2021

  11. [20]

    J Colloid Interface Sci 291:515–528

    Fisher LS, Golovin AA (2005) Nonlinear sta- bility analysis of a two-layer thin liquid film: Dewetting and autophobic behavior. J Colloid Interface Sci 291:515–528. https://doi.org/10. 1016/j.jcis.2005.05.024 18 Preprint– contact: u.thiele@uni-muenster.de – www.uwethiele.de – ...

    2005

  12. [21]

    Rev Mod Phys 57:827–863

    de Gennes PG (1985) Wetting: Statics and dynamics. Rev Mod Phys 57:827–863. https:// doi.org/10.1103/RevModPhys.57.827

    1985 doi

  13. [22]

    Springer, New York, https://doi.org/10.1007/ 978-0-387-21656-0

    de Gennes PG, Brochard-Wyart F, Qu ´er´e D (2004) Capillarity and Wetting Phe- nomena: Drops, Bubbles, Pearls, Waves. Springer, New York, https://doi.org/10.1007/ 978-0-387-21656-0

    2004

  14. [23]

    Prog Polym Sci 28:261–302

    Geoghegan M, Krausch G (2003) Wet- ting at polymer surfaces and interfaces. Prog Polym Sci 28:261–302. https: //doi.org/10.1016/S0079-6700(02)00080-1

    2003 doi

  15. [24]

    Phys Lett A 353:198–204

    Govor LV , Reiter G, Bauer GH, et al (2006) Bilayer formation in thin films of a binary solu- tion. Phys Lett A 353:198–204

    2006

  16. [25]

    Lang- muir 31:11781–11789

    Guan JH, Wells GG, Xu B, et al (2015) Evaporation of sessile droplets on slippery liquid-infused porous surfaces (slips). Lang- muir 31:11781–11789. https://doi.org/10.1021/ acs.langmuir.5b03240

    2015

  17. [26]

    Lang- muir 40:4001–4021

    Hartmann S, Diekmann J, Greve D, et al (2024) Drops on polymer brushes – Advances in thin- film modelling of adaptive substrates. Lang- muir 40:4001–4021. https://doi.org/10.1021/acs. langmuir.3c03313

    2024 doi

  18. [27]

    In: Bungartz HJ, Sch ¨afer M (eds) Fluid-Structure Interaction: Modelling, Simulation, Optimisa- tion

    Heil M, Hazel AL (2006) Oomph-lib - an object- oriented multi-physics finite-element library. In: Bungartz HJ, Sch ¨afer M (eds) Fluid-Structure Interaction: Modelling, Simulation, Optimisa- tion. Springer, Berlin, Heidelberg, p 19–49, https://doi.org/10.1007/3-540-34596-5 2

    2006 doi

  19. [28]

    Soft Matter 17:10359–10375

    Henkel C, Snoeijer JH, Thiele U (2021) Gradient-dynamics model for liquid drops on elastic substrates. Soft Matter 17:10359–10375. https://doi.org/10.1039/D1SM01032H

    2021 doi

  20. [29]

    Q J Mech Appl Math 36:55–69

    Hocking LM (1983) The spreading of a thin drop by gravity and capillarity. Q J Mech Appl Math 36:55–69. https://doi.org/10.1093/qjmam/ 36.1.55

    1983 doi

  21. [30]

    J Eng Math 94:43–61

    Huth R, Jachalski S, Kitavtsev G, et al (2015) Gradient flow perspective on thin-film bilayer flows. J Eng Math 94:43–61. https://doi.org/10. 1007/s10665-014-9698-1

    2015

  22. [31]

    Langmuir 33:5713–5723

    Iqbal R, Dhiman S, Sen AK, et al (2017) Dynam- ics of a water droplet over a sessile oil droplet: Compound droplets satisfying a Neumann con- dition. Langmuir 33:5713–5723. https://doi.org/ 10.1021/acs.langmuir.6b04621

    2017 doi

  23. [32]

    SIAM J Appl Math 73:1183–1202

    Jachalski S, Huth R, Kitavtsev G, et al (2013) Stationary solutions of liquid two-layer thin- film models. SIAM J Appl Math 73:1183–1202. https://doi.org/10.1137/120886613

    2013 doi

  24. [33]

    J Eng Math 86:9–29

    Jachalski S, Peschka D, M ¨unch A, et al (2014) Impact of interfacial slip on the stability of liquid two-layer polymer films. J Eng Math 86:9–29. https://doi.org/10.1007/s10665-013-9651-8

    2014 doi

  25. [34]

    Proc Natl Acad Sci U S A 113:7403–7407

    Karpitschka S, Pandey A, Lubbers LA, et al (2016) Liquid drops attract or repel by the inverted Cheerios effect. Proc Natl Acad Sci U S A 113:7403–7407. https://doi.org/10.1073/pnas. 1601411113

    2016 doi

  26. [35]

    Soft Matter 13:6981–6987

    Keiser A, Keiser L, Clanet C, et al (2017) Drop friction on liquid-infused materials. Soft Matter 13:6981–6987. https://doi.org/10.1039/ c7sm01226h

    2017

  27. [36]

    J Colloid Interface Sci 247:372–380

    Kim HY , Lee HJ, Kang BH (2002) Sliding of liquid drops down an inclined solid surface. J Colloid Interface Sci 247:372–380

    2002

  28. [37]

    J Nonlinear Sci 35

    Kitavtsev G (2024) Composite solutions to a liq- uid bilayer model. J Nonlinear Sci 35. https: //doi.org/10.1007/s00332-024-10108-5

    2024 doi

  29. [38]

    J Fluid Mech 391:45–65

    Kliakhandler IL (1999) Long interfacial waves in multilayer thin films and coupled Kuramoto- Sivashinsky equations. J Fluid Mech 391:45–65. https://doi.org/10.1017/S0022112099005297

    1999 doi

  30. [39]

    SIAM J Appl Math 64:18–40

    Kriegsmann JJ, Miksis MJ (2003) Steady motion of a drop along a liquid interface. SIAM J Appl Math 64:18–40

    2003

  31. [40]

    J Fluid Mech 541:293–315

    Le Grand N, Daerr A, Limat L (2005) Shape and motion of drops sliding down an inclined plane. J Fluid Mech 541:293–315. https://doi.org/10. 1017/s0022112005006105

    2005

  32. [41]

    Nat Phys 18:713–719

    Li X, Bista P, Stetten AZ, et al (2022) Spon- taneous charging affects the motion of sliding drops. Nat Phys 18:713–719. https://doi.org/10. 19 Preprint– contact: u.thiele@uni-muenster.de – www.uwethiele.de – March 30, 2026 1038/s41567-022-01563-6

    2022

  33. [42]

    J Chem Phys 114:2377–2381

    Lin Z, Kerle T, Baker SM, et al (2001) Electric field induced instabilities at liquid/liquid inter- faces. J Chem Phys 114:2377–2381. https://doi. org/10.1063/1.1338125

    2001 doi

  34. [43]

    Ann Phys (Poggendorf) 137:362–377

    L ¨udtge R (1869) Ueber die Ausbre- itung der Fl ¨ussigkeiten auf einander. Ann Phys (Poggendorf) 137:362–377. https://doi.org/10.1002/andp.18692130703

    doi

  35. [44]

    Phys Rev Appl 7:014017

    Luo JT, Geraldi NR, Guan JH, et al (2017) Slippery liquid-infused porous surfaces and droplet transportation by surface acoustic waves. Phys Rev Appl 7:014017. https://doi.org/10. 1103/PhysRevApplied.7.014017

    2017

  36. [45]

    J Fluid Mech 451:411–420

    Mahadevan L, Adda-Bedia M, Pomeau Y (2002) Four-phase merging in sessile compound drops. J Fluid Mech 451:411–420. https://doi.org/10. 1017/S0022112001007108

    2002

  37. [46]

    Soft Mat- ter 5:3801–3809

    Matar OK, Craster RV (2009) Dynamics of surfactant-assisted spreading. Soft Mat- ter 5:3801–3809. https://doi.org/10.1039/ b908719m

    2009

  38. [47]

    Phys Rev Lett 135

    McHale G, Janahi S, Barrio-Zhang H, et al (2025) Adhesive forces in droplet kinetic fric- tion on liquidlike surfaces. Phys Rev Lett 135. https://doi.org/10.1103/ljbc-k193

    2025 doi

  39. [48]

    J Appl Phys

    Mohammad Karim A (2022) A review of physics of moving contact line dynamics models and its applications in interfacial science. J Appl Phys

    2022

  40. [49]

    Euro- phys Lett 95:44003

    Moradi N, Varnik F, Steinbach I (2011) Con- tact angle dependence of the velocity of slid- ing cylindrical drop on flat substrates. Euro- phys Lett 95:44003. https://doi.org/10.1209/ 0295-5075/95/44003

    2011

  41. [50]

    Soft Matter 11:8717–8740

    Mukherjee R, Sharma A (2015) Instability, self- organization and pattern formation in thin soft films. Soft Matter 11:8717–8740. https://doi.org/ 10.1039/c5sm01724f

    2015 doi

  42. [51]

    Non- linearity 23:1559–1583

    N ´araigh L ´O, Thiffeault JL (2010) Nonlinear dynamics of phase separation in thin films. Non- linearity 23:1559–1583. https://doi.org/10.1088/ 0951-7715/23/7/003

    2010

  43. [52]

    Soft Matter 8:11042– 11050

    Neeson MJ, Tabor RF, Grieser F, et al (2012) Compound sessile drops. Soft Matter 8:11042– 11050. https://doi.org/10.1039/c2sm26637g

    2012 doi

  44. [53]

    Curr Opin Colloid Interface Sci 51:101398

    Nepomnyashchy A (2021) Droplet on a liq- uid substrate: Wetting, dewetting, dynamics, instabilities. Curr Opin Colloid Interface Sci 51:101398. https://doi.org/10.1016/j.cocis.2020. 101398

    2021 doi

  45. [54]

    Phys Rev Lett 102:164501

    Nepomnyashchy A, Simanovskii I (2009) Instabilities and ordered patterns in non- isothermal ultrathin bilayer fluid films. Phys Rev Lett 102:164501. https://doi.org/10.1103/ PhysRevLett.102.164501

    2009

  46. [55]

    Phys Rev Fluids 6:034001

    Nepomnyashchy A, Simanovskii I (2021) Droplets on the liquid substrate: Thermo- capillary oscillatory instability. Phys Rev Fluids 6:034001. https://doi.org/10.1103/ PhysRevFluids.6.034001

    2021

  47. [56]

    Phys Fluids 29:092104

    Nepomnyashchy AA, Simanovskii IB (2017) Novel criteria for the development of monotonic and oscillatory instabilities in a two-layer film. Phys Fluids 29:092104. https://doi.org/10.1063/ 1.5001729

    2017

  48. [57]

    V orlesungen¨uber mathematische Physik

    Neumann FE, Wangerin A (1894) V orlesun- gen ¨uber die Theorie der Capillarit ¨at, gehalten an der Universit ¨at K ¨onigsberg von Franz Neu- mann. V orlesungen¨uber mathematische Physik. (7. Hft.), B. G. Teubner, Leipzig

  49. [58]

    Phys Rev E 55:R1267–R1270

    Oron A, Rosenau P (1997) Evolution and forma- tion of dispersive-dissipative patterns. Phys Rev E 55:R1267–R1270

    1997

  50. [59]

    Non-linear stability analysis

    Paunov VN, Danov KD, Alleborn N, et al (1998) Stability of evaporating two-layered liquid film in the presence of surfactant - III. Non-linear stability analysis. Chem Eng Sci 53:2839–2857. https://doi.org/10.1016/ S0009-2509(98)00100-6

    1998

  51. [60]

    Phys Fluids 30:082115

    Peschka D (2018) Variational approach to dynamic contact angles for thin films. Phys Fluids 30:082115. https://doi.org/10.1063/1. 5040985 20 Preprint– contact: u.thiele@uni-muenster.de – www.uwethiele.de – March 30, 2026

    2018 doi

  52. [61]

    Sci Rep 8:13295

    Peschka D, Bommer S, Jachalski S, et al (2018) Impact of energy dissipation on interface shapes and on rates for dewetting from liquid substrates. Sci Rep 8:13295. https://doi.org/10. 1038/s41598-018-31418-1

    2018

  53. [62]

    Phys Rev Lett 87:036102

    Podgorski T, Flesselles JM, Limat L (2001) Corners, cusps, and pearls in running drops. Phys Rev Lett 87:036102. https://doi.org/10. 1103/PhysRevLett.87.036102

    2001

  54. [63]

    Phys Rev E 70:025201(R)

    Pototsky A, Bestehorn M, Merkt D, et al (2004) Alternative pathways of dewetting for a thin liq- uid two-layer film. Phys Rev E 70:025201(R). https://doi.org/10.1103/PhysRevE.70.025201

    2004 doi

  55. [64]

    J Chem Phys 122:224711

    Pototsky A, Bestehorn M, Merkt D, et al (2005) Morphology changes in the evolution of liquid two-layer films. J Chem Phys 122:224711. https: //doi.org/10.1063/1.1927512

    2005 doi

  56. [65]

    Euro- phys Lett 74:665–671

    Pototsky A, Bestehorn M, Merkt D, et al (2006) Evolution of three-dimensional interface pat- terns in dewetting two-layer liquid films. Euro- phys Lett 74:665–671. https://doi.org/10.1209/ epl/i2006-10026-8

    2006

  57. [66]

    Soft Matter 13:101–110

    Semprebon C, McHale G, Kusumaatmaja H (2017) Apparent contact angle and contact angle hysteresis on liquid infused surfaces. Soft Matter 13:101–110. https://doi.org/10.1039/ c6sm00920d

    2017

  58. [67]

    Annu Rev Fluid Mech 45:269–292

    Snoeijer JH, Andreotti B (2013) Moving contact lines: Scales, regimes, and dynamical transitions. Annu Rev Fluid Mech 45:269–292. https://doi. org/10.1146/annurev-fluid-011212-140734

    2013 doi

  59. [68]

    Colloid Surf A 553:487–495

    Thiele U (2018) Recent advances in and future challenges for mesoscopic hydrodynamic mod- elling of complex wetting. Colloid Surf A 553:487–495. https://doi.org/10.1016/j.colsurfa. 2018.05.049

    2018 doi

  60. [69]

    Thiele U (2026) What is active wetting? http://arxiv.org/abs/2602.10287

    2026 arXiv

  61. [70]

    Phys Rev Lett 111:117801

    Thiele U, Todorova DV , Lopez H (2013) Gradi- ent dynamics description for films of mixtures and suspensions: dewetting triggered by cou- pled film height and concentration fluctuations. Phys Rev Lett 111:117801. https://doi.org/10. 1103/PhysRevLett.111.117801

    2013

  62. [71]

    Phys Rev Fluids 1:083903

    Thiele U, Archer AJ, Pismen LM (2016) Gradi- ent dynamics models for liquid films with solu- ble surfactant. Phys Rev Fluids 1:083903. https: //doi.org/10.1103/PhysRevFluids.1.083903

    2016 doi

  63. [72]

    Bachelor’s thesis, Uni- versity of M¨unster

    Thy D (2024) Modeling of interface-dominated sliding compound drops. Bachelor’s thesis, Uni- versity of M¨unster

    2024

  64. [73]

    Soft Matter 13:3760–

    Tress M, Karpitschka S, Papadopoulos P, et al (2017) Shape of a sessile drop on a flat surface covered with a liquid film. Soft Matter 13:3760–

    2017

  65. [74]

    Phys Fluids 35:122121

    Tsao YH, Liao YC, Tsao HK (2023) Slid- ing motion of highly deformed droplets on smooth and rough surfaces: shape oscillation, chaotic breakage, corner shape, and pearling. Phys Fluids 35:122121. https://doi.org/10.1063/ 5.0181630

    2023

  66. [75]

    Numer Math-Theory Methods Appl 7:58–106

    Uecker H, Wetzel D, Rademacher JDM (2014) pde2path - a Matlab package for continuation and bifurcation in 2D elliptic systems. Numer Math-Theory Methods Appl 7:58–106. https:// doi.org/10.4208/nmtma.2014.1231nm

    2014

  67. [76]

    Langmuir 21:3710–

    Verma R, Sharma A, Kargupta K, et al (2005) Electric field induced instability and pattern for- mation in thin liquid films. Langmuir 21:3710–

    2005

  68. [77]

    Fluid Dyn 11:714–721

    V oinov OV (1976) Hydrodynamics of wetting. Fluid Dyn 11:714–721. https://doi.org/10.1007/ BF01012963

    1976

  69. [78]

    Phys Fluids 23:062105

    Ward MH (2011) Interfacial thin films rup- ture and self-similarity. Phys Fluids 23:062105. https://doi.org/10.1063/1.3604003

    2011 doi

  70. [79]

    Sci Rep 7:14776

    Xu L, Bandyopadhyay D, Reddy PDS, et al (2017) Giant slip induced anomalous dewet- ting of an ultrathin film on a viscous sub- layer. Sci Rep 7:14776. https://doi.org/10.1038/ s41598-017-14861-4

    2017

  71. [80]

    Phys Rev B 85:235446

    Yadavali S, Krishna H, Kalyanaraman R (2012) Morphology transitions in bilayer spinodal dewetting systems. Phys Rev B 85:235446. https://doi.org/10.1103/PhysRevB.85.235446 21 Preprint– contact: u.thiele@uni-muenster.de – www.uwethiele.de – March 30, 2026

    2012 doi

  72. [81]

    Phil Trans R Soc 95:65–87

    Young T (1805) An essay on the cohesion of flu- ids. Phil Trans R Soc 95:65–87. https://doi.org/ 10.1098/rstl.1805.0005

  73. [82]

    Soft Matter 15:6055–

    Yu HT, Kant P, Dyett B, et al (2019) Splitting droplets through coalescence of two different three-phase contact lines. Soft Matter 15:6055–

    2019

  74. [83]

    J Fluid Mech 917:A37

    Zhang CY , Gao P, Li EQ, et al (2021) On the compound sessile drops: Configuration bound- aries and transitions. J Fluid Mech 917:A37. https://doi.org/10.1017/jfm.2021.314

    2021 doi

  75. [84]

    J Colloid Inter- face Sci 462:88–99

    Zhang Y , Chatain D, Anna SL, et al (2016) Stability of a compound sessile drop at the axisymmetric configuration. J Colloid Inter- face Sci 462:88–99. https://doi.org/10.1016/j. jcis.2015.09.043

    2016 doi

  76. [85]

    Adv Mater 18:2311470

    Zhou X, Wang Y , Li X, et al (2024) Thick- ness of nanoscale poly(dimethylsiloxane) layers determines the motion of sliding water drops. Adv Mater 18:2311470. https://doi.org/10.1002/ adma.202311470 22 Preprint– contact: u.thiele@uni-muenster.de – www.uwethiele.de – March 30, 2026

    2024

  77. [132]

    https://doi.org/10.1063/5.0102028

    doi

  78. [516]

    https://doi.org/10.1016/j.jcis.2006.08.009

    2006 doi

  79. [596]

    https://doi.org/10.1021/acs.macromol. 0c02234

    doi

  80. [3721]

    https://doi.org/10.1021/la0472100

    doi

  81. [3767]

    https://doi.org/10.1039/c7sm00437k

    doi

  82. [6061]

    https://doi.org/10.1039/c9sm00638a

    doi

Pith tools

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