We investigate the motion of a thin liquid drop on a pre-stretched, highly bendable elastic sheet. Under the lubrication approximation, we derive a system of fourth-order partial differential equations, along with appropriate boundary and contact line conditions, to describe the evolution of the fluid interface and the elastic sheet. Extending the classical analysis of Cox and Voinov, we perform a four-region matched asymptotic analysis of the model in the limit of small slip length. The central result is an asymptotic relation for the contact line speed in terms of the apparent contact angles. We validate the relation through numerical simulations. A key implication of this result is that a soft substrate retards drop spreading but enhances receding, compared to the dynamics on a rigid substrate. The relation remains valid across a wide range of bending modulus, despite the distinguished limit assumed in the analysis.

Floating particles deform the liquid–gas interface, which may lead to capillary repulsion or attraction and aggregation of nearby particles (e.g. the Cheerios effect). Previous studies employed the superposition of capillary multipoles to model interfacial deformation for circular or ellipsoidal particles. However, the induced interfacial deformation depends on the shape of the particle and becomes more complex as the geometric complexity of the particle increases. This study presents a generalised solution for the liquid–gas interface near complex anisotropic particles using the domain perturbations approach. This method enables a closed-form solution for interfacial deformation near particles with an anisotropic shape, as well as the varying height of the pinned liquid–gas contact line. We verified the model via experiments performed with fixed particles held at the water level with shapes such as a circle, hexagon and square, which have either flat or sinusoidal pinned contact lines. Although in this study we concentrate on the equilibrium configuration of the liquid–gas interface in the vicinity of particles placed at fixed positions, our methodology paves the way to explore the interactions among multiple floating anisotropic particles and, thus, the role of particle geometry in self-assembly processes of floating particles.

The interface shape near a moving contact line is described by the Cox–Voinov theory, which contains a constant term that is not trivially obtained. In this work, an approximate expression of this term in explicit form is derived under the condition of a Navier slip. Introducing the approximation of a local slippery wedge flow, we first propose a novel form of the generalised lubrication equation. A matched asymptotic analysis of this equation yields the Cox–Voinov relation with the constant term expressed in elementary functions. For various viscosity ratios and contact angles, the theoretical predictions are rigorously validated against full numerical solutions of the Stokes equations and available asymptotic results.

The study of the shape of droplets on surfaces is an important problem in the physics of fluids and has applications in multiple industries, from agrichemical spraying to microfluidic devices. Motivated by these real-world applications, computational predictions for droplet shapes on complex substrates – rough and chemically heterogeneous surfaces – are desired. Grid-based discretisations in axisymmetric coordinates form the basis of well-established numerical solution methods in this area, but when the problem is not axisymmetric, the shape of the contact line and the distribution of the contact angle around it are unknown. Recently, particle methods, such as pairwise-force smoothed particle hydrodynamics (PF-SPH), have been used to conveniently forego explicit enforcement of the contact angle. The pairwise-force model, however, is far from mature, and there is no consensus in the literature on the choice of pairwise-force profile. We propose a new pair of polynomial force profiles with a simple motivation and validate the PF-SPH model in both static and dynamic tests. We demonstrate its capabilities by computing droplet shapes on a physically structured surface, a surface with a hydrophilic stripe and a virtual wheat leaf with both micro-scale roughness and variable wettability. We anticipate that this model can be extended to dynamic scenarios, such as droplet spreading or impaction, in the future.

In gas evolving electrolysis, bubbles grow at electrodes due to a diffusive influx from oversaturation generated locally in the electrolyte by the electrode reaction. When considering electrodes of micrometre size resembling catalytic islands, direct numerical simulations show that bubbles may approach dynamic equilibrium states at which they neither grow nor shrink. These are found in undersaturated and saturated bulk electrolytes during both pinning and expanding wetting regimes of the bubbles. The equilibrium is based on the balance of local influx near the bubble foot and global outflux. To identify the parameter regions of bubble growth, dissolution and dynamic equilibrium by analytical means, we extend the solution of Zhang & Lohse (2023 J. Fluid Mech. vol. 975, R3) by taking into account modified gas fluxes across the bubble interface, which result from a non-uniform distribution of dissolved gas. The Damköhler numbers at equilibrium are found to range from small to intermediate values. Unlike pinned nanobubbles studied earlier, for micrometre-sized bubbles the Laplace pressure plays only a minor role. With respect to the stability of the dynamic equilibrium states, we extend the methodology of Lohse & Zhang (2015a Phys. Rev. E vol. 91, 031003(R)) by additionally taking into account the electrode reaction. Under contact line pinning, the equilibrium states are found to be stable for flat nanobubbles and for microbubbles in general. For unpinned bubbles, the equilibrium states are always stable. Finally, we draw conclusions on how to possibly enhance the efficiency of electrolysis.

Underwater capillary tubes fill rapidly with the surrounding liquid. Capillary and hydrostatic pressures push the liquid into the tube, causing the air to exit as bubbles at the other end. We study the natural filling process of a vertical capillary tube immersed in water during several bubble formation events. A theoretical model is proposed that captures the dynamics of the meniscus inside the capillary tube as it fills with water. We find good agreement with the experimental data that describe this special case of spontaneous flow using a dynamic contact angle model based on molecular kinetic theory.

The breakup dynamics of viscous liquid bridges on solid surfaces is studied experimentally. It is found that the dynamics bears similarities to the breakup of free liquid bridges in the viscous regime. Nevertheless, the dynamics is significantly influenced by the wettability of the solid substrate. Therefore, it is essential to take into account the interaction between the solid and the liquid, especially at the three-phase contact line. It is shown that when the breakup velocity is low and the solid surface is hydrophobic, the dominant channel of energy dissipation is likely due to thermally activated jumping of molecules, as described by the molecular kinetic theory. Nevertheless, the viscous dissipation in the bulk due to axial flow along the bridge can be of importance for long bridges. In view of this, a scaling relation for the time dependence of the minimum width of the liquid bridge is derived. For high viscosities, the scaling relation captures the time evolution of the minimum width very well. Furthermore, it is found that external geometrical constraints alter the dynamic behaviour of low and high viscosity liquid bridges in a different fashion. This discrepancy is explained by considering the dominant forces in each regime. Lastly, the morphology of the satellite droplets deposited on the surface is qualitatively compared with that of free liquid bridges.

In gas evolving electrolysis, bubbles grow at electrodes due to a diffusive influx from oversaturation generated locally in the electrolyte by the electrode reaction. When considering electrodes of micrometre size resembling catalytic islands, direct numerical simulations show that bubbles may approach dynamic equilibrium states at which they neither grow nor shrink. These are found in under- and saturated bulk electrolytes during both pinning and expanding wetting regimes of the bubbles. The equilibrium is based on the balance of local influx near the bubble foot and global outflux. To identify the parameter regions of bubble growth, dissolution and dynamic equilibrium by analytical means, we extend the solution of Zhang & Lohse (2023) J. Fluid Mech. 975, R3, by taking into account modified gas fluxes across the bubble interface, that result from a non-uniform distribution of dissolved gas. The Damköhler numbers at equilibrium are found to range from small to intermediate values. Unlike pinned nano-bubbles studied earlier, for micrometre-sized bubbles the Laplace pressure plays only a minor role. With respect to the stability of the dynamic equilibrium states, we extend the methodology of Lohse & Zhang (2015a) Phys. Rev. E 91 (3), 031003(R), by additionally taking into account the electrode reaction. Under contact line pinning, the equilibrium states are found to be stable for flat nano-bubbles and for micro-bubbles in general. For unpinned bubbles, the equilibrium states are always stable. Finally, we draw conclusions on how to possibly enhance the efficiency of electrolysis.

We consider the dynamics of a liquid film with a pinned contact line (for example, a drop), as described by the one-dimensional, surface-tension-driven thin-film equation $h_t + (h^n h_{xxx})_x = 0$, where $h(x,t)$ is the thickness of the film. The case $n=3$ corresponds to a film on a solid substrate. We derive an evolution equation for the contact angle $\theta (t)$, which couples to the shape of the film. Starting from a regular initial condition $h_0(x)$, we investigate the dynamics of the drop both analytically and numerically, focusing on the contact angle. For short times $t\ll 1$, and if $n\ne 3$, the contact angle changes according to a power law $\displaystyle t^{\frac {n-2}{4-n}}$. In the critical case $n=3$, the dynamics become non-local, and $\dot {\theta }$ is now of order $\displaystyle {\rm{e}}^{-3/(2t^{1/3})}$. This implies that, for $n=3$, the standard contact line problem with prescribed contact angle is ill posed. In the long time limit, the solution relaxes exponentially towards equilibrium.

We study the dynamics of a thin liquid sheet that flows upwards along the sides of a vertically aligned, impacting plate. Upon impact of the vertical solid plate onto a liquid pool, the liquid film is ejected and subsequently continues to flow over the solid surface while the plate enters the water. With increasing impact velocity, the liquid film is observed to rise up faster and higher. We focus on the time evolution of the liquid film height and the thickness of its upper rim and discuss their dynamics in detail. Similar to findings in previous literature on sheet fragmentation during drop impact, we find the rim thickness to be governed by the local instantaneous capillary number based on gravity and the deceleration of the liquid sheet, showing that the retraction of the rim is primarily due to capillarity. In contrast, for the liquid film height, we demonstrate that the viscous dissipation in the thin boundary layer is the dominant factor for the vertical deceleration of the liquid sheet, by modelling the time evolution of the film height and showing that the influences of capillarity, gravity and deceleration due to the air phase are all negligible compared with the viscous term. Finally, we introduce characteristic viscous time and length scales based on the initial rim thickness and show that the maximum height of the film and the corresponding time can be determined from these viscous scales.

Spin coating is the process of generating a uniform coating film on a substrate by centrifugal forces during rotation. In the framework of lubrication theory, we investigate the axisymmetric film evolution and contact-line dynamics in spin coating on a partially wetting substrate. The contact-line singularity is regularized by imposing a Navier slip model. The interface morphology and the contact-line movement are obtained by numerical solution and asymptotic analysis of the lubrication equation. The results show that the evolution of the liquid film can be classified into two modes, depending on the rotational speed. At lower speeds, the film eventually reaches an equilibrium state, and we provide a theoretical description of how the equilibrium state can be approached through matched asymptotic expansions. At higher speeds, the film exhibits two or three distinct regions: a uniform thinning film in the central region, an annular ridge near the contact line, and a possible Landau–Levich–Derjaguin-type (LLD-type) film in between that has not been reported previously. In particular, the LLD-type film occurs only at speeds slightly higher than the critical value for the existence of the equilibrium state, and leads to the decoupling of the uniform film and the ridge. It is found that the evolution of the ridge can be well described by a two-dimensional quasi-steady analysis. As a result, the ridge volume approaches a constant and cannot be neglected to predict the variation of the contact-line radius. The long-time behaviours of the film thickness and the contact radius agree with derived asymptotic solutions.

A model of imbibition dynamics in a channel of flattened triangular cross-section is presented, taking into account the liquid film flow in the corners of the channel. The quasi-analytical solutions are derived on the basis of a lubrication approximation. The analysis encompasses two imbibition scenarios corresponding to a constant flow rate or constant pressure imposed in the wetting fluid at the inlet of the channel. In the former case, the process starts with a liquid film flow regime in the corners that is followed by a bulk and corner film flow regime characterised by a triple point advancing (far) ahead of the bulk meniscus after its entrance in the channel. In the latter case, the occurrence of the bulk and corner film flow regime is conditioned by an imposed pressure yielding a capillary pressure at the inlet smaller than a threshold capillary pressure. Above this threshold, the liquid film regime remains. For both imbibition scenarios under concern, important features are highlighted, including (i) the time scalings of the dynamics of both the triple point and apex of the bulk meniscus (when it exists), (ii) the contrast in the positions of these two points showing that the classical Washburn approach, which neglects the effect of the corner films, overpredicts the dynamics of the bulk meniscus. The important consequence is an early wetting fluid breakthrough at the channel outlet much before the bulk meniscus arrival. Comparisons with experimental data available in the literature are provided, validating the approach proposed in this work.

The flow near a moving contact line depends on the dynamic contact angle, viscosity ratio and capillary number. We report experiments involving immersing a plate into a liquid bath, concurrently measuring the interface shape, interfacial velocity and fluid flow using digital image processing and particle image velocimetry. All experiments were performed at low plate speeds to maintain small Reynolds and capillary numbers for comparison with viscous theories. The dynamic contact angle, measured in the viscous phase, was kept below $90^{\circ }$ and the viscosity ratio, $\lambda < 1$. This region of parameter space is largely unexplored for advancing contact lines. An important aim of the present study is to provide new experimental data against which new contact line models can be developed. The flow field is directly compared against the prediction from the viscous theory of Huh & Scriven (J. Colloid Interface Sci., vol. 35, issue 1, 1971, pp. 85–101) but with a slight modification involving the curved interface. Remarkable agreement is found between experiments and theory across a wide parameter range. The prediction for interfacial speed from Huh & Scriven is also in excellent agreement with experiments except in the vicinity of the contact line. Material points along the interface were found to rapidly slow down near the contact line, thus alleviating the singularity at the moving contact line. To the best of our knowledge, such a detailed test of theoretical models has not been performed before and we hope the present study will spur new modelling efforts in the field.

In spray cleaning, a multitude of small drops, violently accelerated by a high-speed gas stream, strike a dirty surface. This process is extremely effective: very little dirt can resist it. This is true even for dirt particles whose characteristic size is less than 100 nm. Spray cleaning is classically modelled by balancing particle adhesion with either inertial stress or viscous shear near the surface, the latter being calculated using droplet size and velocity as the characteristic length and velocity. This results in dimensionless numbers that are often well below one, suggesting that the mechanical stress exerted on the surface by the drop impact that detaches the particle is not well captured. Using quantitative nanoscale measurements, we show that the remarkable efficiency of spray cleaning results from the forced spreading of each droplet on the surface, which generates an unsteady and inhomogeneous shear confined to a boundary layer entrained in the wake of the liquid–solid contact line. In the very first moments of impact, the boundary layer is extremely thin, yielding a gigantic stress: the contact line of the spreading droplets sweeps all the surface particles away. We propose a quantitative model of spray cleaning based on this unsteady surface stress, which agrees well with (i) experimental data obtained with spray droplets of $34\ \mathrm {\mu }$m mean radius impacting the surface to be cleaned at a mean velocity ranging between 30 and 70 m s$^{-1}$ and contamination by nanoparticles of varying nature and shape and (ii) data in the literature on spray cleaning.

During oscillatory wetting, a phase retardation emerges between contact angle variation and contact line velocity, presenting as a hysteresis loop in their correlation – an effect we term dynamic hysteresis. This phenomenon is found to be tunable by modifying the surface with different molecular layers. A comparative analysis of dynamic hysteresis, static hysteresis and contact line friction coefficients across diverse substrates reveals that dynamic hysteresis is not a result of dissipative effects but is instead proportionally linked to the static hysteresis of the surface. In the quest for appropriate conditions to model oscillatory contact line motion, we identify the generalized Hocking's linear law and modified generalized Navier boundary condition as alternative options for predicting realistic dynamic hysteresis.

The nanoscale is the new frontier of fluid dynamics and its phenomenology can echo at the macroscale as in the canonical example of drop impact on a planar substrate. Unprecedented advances in measurement technology have recently equipped fluid dynamicists with the ability to probe nanoscale effects. The paper by Li et al. (J. Fluid Mech., vol. 785, 2015, R2) uses ultrafast imaging at the hundreds of nanoseconds scale to resolve the first contact between the drop and the substrate and thereby reveal the effect of prescribed nano-roughness on contact line motion.

Superhydrophobic (SHPo) surfaces can capture a thin layer of air called a plastron under water to reduce skin friction. Although a ~30 % drag reduction has been recently reported with longitudinal micro-trench SHPo surfaces under a boat and in a towing tank, the results lacked the consistency to establish a clear trend. Designed based on Yu et al. (J. Fluid Mech, vol. 962, 2023, A9), this work develops and tests a series of high-performance SHPo surface coupons that can sustain a pinned plastron underneath a passenger motorboat revamped to reach 14 knots. Importantly, plastrons in a pinned state, not just their existence, are confirmed during flow experiments for the first time. All the drag-reduction data measured on different longitudinal micro-trenches are found to collapse if plotted against slip length in wall units. In comparison, aligned posts and transverse trenches show less and little drag reduction, respectively, confirming the adverse effect of the spanwise slip in turbulent flows. This report not only verifies SHPo surfaces can provide a consistent drag reduction at high speeds in open sea but also shows that one may predict the amount of drag reduction in turbulent flows using the physical slip length obtained for Stokes flows.

When placed at the surface of a volatile liquid, a sphere of hot dense non-volatile material remains suspended until it cools sufficiently. The duration of this ‘inverse Leidenfrost’ phenomenon depends on the Nusselt number $Nu$ of the sphere, itself determined by flow in the film of vapour separating particle and liquid. It is shown that provided the Nusselt number is large, it can be calculated numerically using only the Laplace relation and the equations governing the thin film; patching to a solution for the outer thick film is not necessary. This method is demonstrated by using it to determine $Nu$ for a sphere sufficiently small that in the governing equations, the acceleration due to gravity is negligible except where multiplied by the density of the sphere. Numerical results giving $Nu$ as a function of a dimensionless measure of sphere weight are supplemented with analysis showing that, when the weight is of the order of the maximum supportable by surface tension alone, the film consists of a spherical bubble cap bounded by its contact rim. The solutions for these regions are coupled: although the apparent contact angle $\chi$ for the cap is determined within the rim, its value depends on the flow rate arriving from the cap as well as on the additional evaporation from the rim. The latter acts to reduce $\chi$ from the value it would otherwise have, thereby reducing the thickness of the entire cap. For the example treated here, the value of $Nu$ is doubled by this mechanism.

Waves are formed on the surface of a sessile drop driven through substrate vibrations oriented at a slanting angle from the normal. A mathematical model is derived, which leads to an infinite system of coupled Mathieu equations governing the wave dynamics that are solved using Floquet theory. The spatial structure of the waves is described by the mode number pair $[\ell,m]$, where $\ell$ and $m$ are the polar and azimuthal mode numbers, respectively. Limiting cases corresponding to horizontal and vertical vibrations are discussed with predictions agreeing well with prior literature. We focus our results on three drop motions – (1) harmonic $[1,1]$ rocking mode, (2) harmonic $[2,0]$ pumping mode, and (3) subharmonic rocking $[1,1]$ mode – as they depend upon the slanting angle, static contact angle, and contact-line conditions, which we assume to be either pinned or freely moving with fixed contact angle. New theoretical predictions are tested through experiments over a range of parameters, showing good agreement.

Depinning of liquid droplets on substrates by flow of a surrounding immiscible fluid is central to applications such as cross-flow microemulsification, oil recovery and waste cleanup. Surface roughness, either natural or engineered, can cause droplet pinning, so it is of both fundamental and practical interest to determine the flow strength of the surrounding fluid required for droplet depinning on rough substrates. Here, we develop a lubrication-theory-based model for droplet depinning on a substrate with topographical defects by flow of a surrounding immiscible fluid. The droplet and surrounding fluid are in a rectangular channel, a pressure gradient is imposed to drive flow and the defects are modelled as Gaussian-shaped bumps. Using a precursor-film/disjoining-pressure approach to capture contact-line motion, a nonlinear evolution equation is derived describing the droplet thickness as a function of distance along the channel and time. Numerical solutions of the evolution equation are used to investigate how the critical pressure gradient for droplet depinning depends on the viscosity ratio, surface wettability and droplet volume. Simple analytical models are able to account for many of the features observed in the numerical simulations. The influence of defect height is also investigated, and it is found that, when the maximum defect slope is larger than the receding contact angle of the droplet, smaller residual droplets are left behind at the defect after the original droplet depins and slides away. The model presented here yields considerably more information than commonly used models based on simple force balances, and provides a framework that can readily be extended to study more complicated situations involving chemical heterogeneity and three-dimensional effects.