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.

Both experimental and theoretical studies of fast and microscale physical phenomena occurring during the growth of vapour bubbles in nucleate pool boiling are reported. The focus is on the liquid film of micrometric thickness (a ‘microlayer’) that can form between the heater and the liquid–vapour interface of a bubble. The microlayer strongly affects the macroscale heat transfer and is thus important to be understood. The microlayer appears as a result of the inertial forces that cause the hemispherical bubble shape. It is shown that the microlayer can be seen as the Landau–Levich film deposited by the bubble foot edge during its receding. Paradoxically, the deposition is controlled by viscosity and surface tension. The microlayer profile measured with white-light interferometry, the temperature distribution over the heater, and the bubble shape are observed with synchronised high-speed cameras. According to the numerical simulations, the microlayer consists of two regions: a dewetting ridge near the contact line, and a longer and flatter bumped part. It is shown that the ridge cannot be measured by interferometry because of its intrinsic limitation on the interface slope. The ridge growth is linked to the contact line receding. The simulated dynamics of both the bumped part and the contact line agrees with the experiment. The physical origin of the bump in the flatter part of microlayer is explained.

Surface tension and wetting are dominating physical effects in microscale and nanoscale flows. We present an efficient and reliable model of surface tension and equilibrium contact angles in smoothed particle hydrodynamics for free-surface problems. We demonstrate its robustness and accuracy by simulating several three-dimensional free-surface flow problems driven by interfacial tension.

Droplet spreading is ubiquitous and plays a significant role in liquid-based energy systems, thermal management devices and microfluidics. While the spreading of non-volatile droplets is quantitatively understood, the spreading and flow transition in volatile droplets remains elusive due to the complexity added by interfacial phase change and non-equilibrium thermal transport. Here we show, using both mathematical modelling and experiments, that the wetting dynamics of volatile droplets can be scaled by the spatial–temporal interplay between capillary, evaporation and thermal Marangoni effects. We elucidate and quantify these complex interactions using phase diagrams based on systematic theoretical and experimental investigations. A spreading law of evaporative droplets is derived by extending Tanner's law (valid for non-volatile liquids) to a full range of liquids with saturation vapour pressure spanning from $10^1$ to $10^4$ Pa and on substrates with thermal conductivity from $10^{-1}$ to $10^3\ {\rm W}\ {\rm m}^{-1}\ {\rm K}^{-1}$. In addition to its importance in fluid-based industries, the conclusions also enable a unifying explanation to a series of individual works including the criterion of flow reversal and the state of dynamic wetting, making it possible to control liquid transport in diverse application scenarios.

Using as a starting point conservation of momentum, a multiphase mechanical energy balance equation is derived that accounts for multiple material phases and interfaces present within a moving control volume. This balance is applied to a control volume that is anchored to a three-phase contact line as it advances continuously over the surface of a rough and chemically homogeneous and inert solid. Using semi-quantitative models for the material behaviour occurring within the control volume, an order of magnitude analysis is performed to neglect insignificant terms, producing an equation for predicting contact-angle hysteresis from a knowledge of the interface dynamics occurring around the three-phase contact line. It is shown that the viscous energy dissipation that occurs during the ‘stick–slip’ motion of the three-phase contact line, being the cause of contact-angle hysteresis on rough surfaces, can be calculated from changes in intermediate equilibrium interface states. The balance is applied to the Wenzel, Cassie–Baxter and Fakir (super-hydrophobic) wetting states, showing for the Fakir case that significant dissipation occurs during both interface advance and recede, and relating these dissipations to interfacial area changes that occur around the ‘stick–slip’ events.

The classical Cox–Voinov theory of contact line motion provides a relation between the macroscopically observable contact angle, and the microscopic wetting angle as a function of contact-line velocity. Here, we investigate how viscoelasticity, specifically the normal stress effect, modifies the wetting dynamics. Using the thin film equation for the second-order fluid, it is found that the normal stress effect is dominant at small scales yet can significantly affect macroscopic motion. We show that the effect can be incorporated in the Cox–Voinov theory through an apparent microscopic angle, which differs from the true microscopic angle. The theory is applied to the classical problems of drop spreading and dip coating, which shows how normal stress facilitates (inhibits) the motion of advancing (receding) contact lines. For rapid advancing motion, the apparent microscopic angle can tend to zero, in which case the dynamics is described by a regime that was already anticipated in Boudaoud (Eur. Phys. J. E, vol. 22, 2007, pp. 107–109).

Capillary folding is the process of folding planar objects into three-dimensional (3-D) structures using capillary force. It has many important industrial applications, e.g. the fabrication of microelectromechanical systems. In this work, we propose a 3-D model for the capillary folding of thin elastic sheets with pinned contact lines. The energy of the system consists of interfacial energies between the different phases and the elastic energy of the sheet. The latter is described by the nonlinear Koiter's model, which can accommodate large deformations of the sheet. From the energy, we derive the governing equations for the static system using a variational approach. We then develop a numerical method to find equilibrium solutions via a relaxation dynamics. At each time step, we evolve the sheet by using the subdivision element method, and update the droplet surface by minimizing a squared area functional using linear finite elements. Qualitatively, numerical solutions for the equilibrium configurations of the sheet–droplet system agree well with those observed in experiments. Furthermore, we identify the critical bendabilities and present bifurcation diagrams for the folding of a triangular sheet. The results exhibit rich and fully 3-D behaviours that were not captured by previous two-dimensional models. Our results provide new insights into the nonlinear process of capillary folding, and may contribute to the advancement of microfabrication techniques.

From industrial-scale production to small-scale fabrication of functional films, the blade coating method is used widely to apply a uniform thin liquid film on a moving substrate. However, conventional hydrodynamic models are inadequate for laboratory-scale low-speed blade coating, where capillary forces dominate. In this study, the low-speed blade coating of non-evaporative Newtonian fluids was investigated in experimental, computational and analytical approaches. The transient free boundary problem was solved utilizing a two-dimensional finite element method, and a simple viscocapillary model was developed to describe the viscous stress and capillary forces within the puddle, and predict the thickness of the wet film as a function of the speed of the substrate. Comparing the predicted film thickness with the computational results, and a flow visualization experiment of blade coating with silicone oil, respectively, confirmed the model's validity. The study indicates that the proposed model may be a useful tool for optimizing laboratory coating processes, as it provides a greater understanding of the low-speed blade coating system on a laboratory scale.

We investigate jumping of sessile droplets from a solid surface in ambient oil using modulated electrowetting actuation. We focus on the case in which the electrowetting effect is activated to cause droplet spreading and then deactivated exactly at the moment the droplet reaches its maximum deformation. By systematically varying the control parameters such as the droplet radius, liquid viscosity and applied voltage, we provide detailed characterisation of the resulting behaviours including a comprehensive phase diagram separating detachment from non-detachment behaviours, as well as how the detach velocity and detach time, i.e. duration leading to detachment, depend on the control parameters. We then construct a theoretical model predicting the detachment condition using energy conservation principles. We finally validate our theoretical analysis by experimental data obtained in the explored ranges of the control parameters.

Evaporating sessile droplets are critical to many industrial applications and are also ubiquitous in nature. Two predominant evaporation models have emerged in the literature, one-sided and diffusion-limited, with differing assumptions on the evaporation process. Both models are used widely, and their predictions can differ greatly from each other, but the physical mechanisms responsible for these differences are not yet well understood. Here, we develop a lubrication-theory-based model of a thin evaporating sessile droplet, and compare predictions from both evaporation models to elucidate the origins of the differences in their predictions. For the one-sided model, we derive expressions for the droplet lifetime, show that in certain parameter regimes the total evaporation rate is proportional to the droplet surface area, and demonstrate that the contact line is always warmer than the bulk of the droplet. Furthermore, we show that differences in the structures of the evaporation models near the contact line lead to qualitatively different behaviour of the apparent contact angles and interface temperature profiles. The fundamental understanding gained from this work is expected to be helpful in determining which evaporation model is most appropriate for describing experimental observations.

We study the stability and collapse of holes at the wall in liquid layers on circular bounded containers with various wettabilities. Three distinct wetting modes of the hole are observed, which are related to the wettability of the container: when the substrate and the inner wall of the container are superhydrophobic, a stable hole remains as the liquid volume is continuously increased until the liquid layer covers the entire substrate; when the substrate and the inner wall are hydrophobic, an eye-shaped hole remains stable as the projected area of the hole exceeds a critical value $A_c$, however, the hole collapses if the liquid volume is further increased; when the substrate is superhydrophobic but the wall is hydrophilic, on increasing the liquid volume, the hole suddenly transfers into a circular hole and is pushed against the wall, leaving the hole dwelling around the centre of the container. Theoretical analyses and numerical simulations are conducted to establish the phase diagram for different wetting modes. It is found that, in the second mode, $A_c$ increases with the size of the container but decreases with the contact angle of the substrate and the wall. Furthermore, we experimentally investigate the dynamics of the hole. The time evolution of the area of the hole obeys a scaling relationship $A \sim (t_0 - t)^{1.1}$ after the hole collapses at time $t_0$.

A thin liquid droplet spreads on a soft viscoelastic substrate with arbitrary rheology. Lubrication theory is applied to the governing field equations in the liquid and solid domains, which are coupled through the free boundary at the solid–liquid interface, to derive a set of reduced equations that describe the spreading dynamics. Fourier transform techniques and the finite difference method are used to construct a solution for the dynamic liquid–gas and solid–liquid interface shapes, as well as the macroscopic contact angle. Substrate properties affect the spreading dynamics through the contact angle and internal droplet flow fields, and these mechanisms are revealed. Increased substrate softness increases the spreading rate, whereas increased viscoelasticity decreases the spreading rate. For the case of a purely elastic substrate, the spreading power-law exponent recovers Tanner's law in the rigid limit and increases with substrate softness.

We study experimentally the dynamics of a water droplet on a tilted and vertically oscillating rigid fibre. As we vary the frequency and amplitude of the oscillations the droplet transitions between different modes: harmonic pumping, subharmonic pumping, a combination of rocking and pumping modes, and a combination of pumping and swinging modes. We characterize these responses and report how they affect the sliding speed of the droplet along the fibre. The swinging mode is explained by a minimal model making an analogy between the droplet and a forced elastic pendulum.

We investigate numerically the spreading dynamics of a bubble coming into contact with a smooth solid substrate in a viscous liquid. The substrate is partially wettable, and the singularity of the moving contact line is relieved by adopting the Navier-slip model. The Stokes equations are solved by employing a boundary element method coupled with adaptive mesh refinement. This allows us to realize sufficiently small slip lengths down to $O(10^{-5})$ in dimensionless form, which is crucial to resolve the local interface structures at the early stage of spreading. The results show that the early-stage spreading of the bubble is always characterized by the growth and propulsion of a dewetting liquid rim close to the contact line, while the macroscopic interface remains unchanged. The evolution of the contact line and the morphology of the rim depends on the wettability and the slip length, and a parametric investigation is performed. Based on mass conservation, a relation between the rim size and the spreading radius is established. We also propose an analytical prediction of the temporal variation in the contact line radius at the early stage of spreading, which is found to follow a logarithmically corrected linear relation rather than a pure power law. Moreover, the early stages of bubble spreading are qualitatively similar for two-dimensional and axisymmetric configurations.

The dynamics of compound drops impacting on a flat substrate is numerically investigated using a ternary-fluid diffuse-interface method, with the aim of assessing the effect of a density difference between the inner and outer droplets (denoted by $\lambda$) on the evolution of the interfaces. With the help of numerical simulations, we find that, at the intermediate stage of drop impact, the inner droplet exhibits a self-similar deformation at $\lambda =1$ and relatively high Weber number, and experiences more or less a uniform acceleration for various $\lambda$. In particular, the acceleration magnitude at $\lambda \ne 1$ can be correlated with the acceleration at $\lambda =1$ and the Atwood number. When the inner droplet is denser than the outer one, a lamella occurs at the spreading front of the inner droplet. We present a scaling analysis of the thickness of the lamella, and the resultant theoretical prediction is in good agreement with numerical results. At the maximal spreading of the compound drop, a bulging structure is formed around the symmetry axis due to the presence of the inner droplet, thereby effectively reducing the liquid supply to the spreading front and leading to a decrease of maximal spreading ratio $\beta _{max}$ as compared with a pure drop. We proposed a corrected Weber number $We^*_\lambda$ by taking account of the combined effects of $\lambda$, volume fraction of the inner droplet, Weber number and morphology of the compound drop. Integrating $We^*_\lambda$ with the universal model of $\beta _{max}$ for impacting pure drops, we successfully build up a new model for predicting the maximal spreading ratio of impacting compound drops with various $\lambda$.

The deformation, movement and breakup of a wall-attached droplet subject to Couette flow are systematically investigated using an enhanced lattice Boltzmann colour-gradient model, which accounts for not only the viscoelasticity (described by the Oldroyd-B constitutive equation) of either droplet (V/N) or matrix fluid (N/V) but also the surface wettability. We first focus on the steady-state deformation of a sliding droplet for varying values of capillary number ($Ca$), Weissenberg number ($Wi$) and solvent viscosity ratio ($\beta$). Results show that the relative wetting area $A_r$ in the N/V system is increased by either increasing $Ca$, or by increasing $Wi$ or decreasing $\beta$, where the former is attributed to the increased viscous force and the latter to the enhanced elastic effects. In the V/N system, however, $A_r$ is restrained by the droplet elasticity, especially at higher $Wi$ or lower $\beta$, and the inhibiting effect strengthens with an increase of $Ca$. Decreasing $\beta$ always reduces droplet deformation when either fluid is viscoelastic. The steady-state droplet motion is quantified by the contact-line capillary number $Ca_{cl}$, and a force balance is established to successfully predict the variations of $Ca_{cl}/Ca$ with $\beta$ for each two-phase viscosity ratio in both N/V and V/N systems. The droplet breakup is then studied for varying $Wi$. The critical capillary number of droplet breakup monotonically increases with $Wi$ in the N/V system, while it first increases, then decreases and finally reaches a plateau in the V/N system.

Hydraulic fracturing for oil and gas production from shale formations, as well as natural geological phenomena, involve the propagation of thin viscous films within elastic media. For viscous fluids, stress diverges as the thickness of the film tends to zero, arresting the propagation of the film, and thus implying the contact line paradox. For free-surface films, this paradox is resolved by considering a precursor film, leading to Tanner's law. This approach was extended recently for viscous films between a thin elastic plate and a rigid solid, allowing calculation of the film propagation rate. In this work, we examine the effect of a pre-wetting layer on the rate of propagation of a viscous flow within an infinitely deep and long domain. We analyse the linear and nonlinear dynamic problems, and perform a self-similarity analysis. We find that peeling front propagation scales as time to the power of $1/9$ and $1/3$ for thin and thick pre-wetting layer limits, respectively. Our results contribute to the understanding of the contact line paradox in elastic media and the crucial role of the pre-wetting layer in resolving it.