We developed a two-phase lattice Boltzmann model by coupling the entropic multiple-relaxation-time (EMRT or KBC) collision operator enabling low fluid viscosity, with a source term (Wang et al. 2022, Phys. Rev. E vol. 105, no 4) to independently adjust surface tension. The coupling is implemented via the exact difference method (EDM), which allows full consideration of external-force effects on the entropic stabiliser in KBC, in contrast to the recent work of Wang et al. (2022 Phys. Rev. E vol. 105) and Xu et al. (2024 Comput. Math. Appl. vol. 159, 92–101). More importantly, we address a major drawback of the EDM by explicitly demonstrating how its high-order error terms influence the pressure tensor and surface tension. Using the developed model, we investigated droplet impact and splashing on a thin liquid film at a remarkably high Weber number of ${\textit{We}} = 5000$ and Reynolds number of ${\textit{Re}} = 5000$. Droplet impact and splashing on flat surfaces and mesh structures at very high ${\textit{Re}}$ (15 200) and ${\textit{We}}$ (1020) are also studied after validating four representative cases against experiments. For droplet impact on flat surfaces, hydrophobicity promotes the growth of peripheral instabilities, leading to fingering splashing. Corona splashing transitions to fingering splashing as the liquid–gas viscosity ratio increases. For droplet impact on mesh structures, large openings promote liquid penetration, whereas small openings enhance spreading. As the solid ratio increases, the maximum spreading ratio increases monotonically but nonlinearly, whereas the maximum penetrated liquid pillar length first rises and then drops. These simulations demonstrate the proposed model offers significant advantages for accurately capturing and elucidating complex droplet impact and splashing dynamics at high ${\textit{Re}}$ and ${\textit{We}}$.

A lattice Boltzmann method is adopted to investigate the breakup of surfactant-free and surfactant-laden droplets in both regular and irregular T-junction microchannels. During droplet neck contraction, the neck thinning shifts from inertia dominated to interfacial tension dominated, causing spontaneous rapid neck collapse due to Rayleigh–Plateau instability. For the regular rectangular microchannels, we find that the prerequisite for the spontaneous breakup of a surfactant-free droplet is that the local capillary pressure in the triggering area exceeds the Laplace pressure difference between the inside and outside of the droplet neck. Results show that the critical neck thickness $\delta _\textit{cr}^{*}$ for the droplet spontaneous breakup increases with increasing height-to-width ratio $\chi$ of the microchannel in both surfactant-free and surfactant-laden systems. The presence of surfactants decreases $\delta _\textit{cr}^{*}$ at the identified $\chi$, while the surfactant effects are gradually enhanced as $\chi$ increases. Subsequently, a constriction section is incorporated into the upper microchannel wall to establish an irregular microchannel. As constriction depth (length) increases, $\delta _\textit{cr}^{*}$ linearly decreases (increases) in the surfactant-free system, while $\delta _\textit{cr}^{*}$ exponentially decreases (linearly increases) in the surfactant-laden system. Four empirical formulas are proposed to predict the values of $\delta _\textit{cr}^{*}$ under varying constriction depths and lengths in the two systems.

In this work, we propose a lattice Boltzmann model (LBM) to simulate diverse particle deposition patterns induced by isothermal droplet evaporation. The model is composed of two distributions, for the multiphase flow with phase change, and the particle transport with deposition, coupled with a contact angle hysteresis model for the contact line stick-slip dynamics. The model is validated by two benchmarks, and our simulations agree well with the theoretical solutions or experimental results. With the validated LBM, we first reproduced diverse deposition patterns, ranging from the coffee ring, uniform, to mountain-type patterns in single and multiple symmetrical/unsymmetrical forms. Then a parametric study is conducted to investigate how the solvent/particle/substrate properties affect the evaporation dynamics and resultant deposition patterns. Afterwards, we apply the average ratio ($r_{\phi ,a}$) of particles deposited at the droplet periphery and the centre to quantitatively classify the diverse emerging patterns. We show that $r_{\phi ,a}$ is controlled by the competition between the capillary transport and particle diffusion, leading to a linear dependence on the average Péclet number $\textit{Pe}_{a}$. Finally, we validated the scaling by lattice Boltzmann simulations with the proposed $\textit{Pe}_{a}$ spanning over three orders of magnitude, supplemented by discussions from the aspect of the particle transport equation.

We consider the problem of a cylindrical (quasi-two-dimensional) droplet impacting on a hard surface. Cylindrical droplet impact can be engineered in the laboratory, and a theoretical model of the system can also be used to shed light on various complex experiments involving the impact of liquid sheets. We formulate a rim-lamella model for the droplet-impact problem. Using Gronwall’s inequality applied to the model, we establish theoretical bounds for the maximum spreading radius $\mathcal{R}_{\textit{max}}$ in droplet impact, specifically $k_1 {\textit{Re}}^{1/3}-k_2(1-\cos \vartheta _a)^{1/2}({\textit{Re}}/{\textit{We}})^{1/2}\leq \mathcal{R}_{\textit{max}}/R_0\leq k_1{\textit{Re}}^{1/3}$, valid for ${\textit{Re}}$ and ${\textit{We}}$ sufficiently large. Here, ${\textit{Re}}$ and ${\textit{We}}$ are the Reynolds and Weber number based on the droplet’s pre-impact velocity and radius $R_0$, $\vartheta _a$ is the advancing contact angle (assumed constant in our simplified analysis) and $k_1$ and $k_2$ are constants. We perform several campaigns of simulations using the volume of fluid method to model the droplet impact, and we find that the simulation results fall within the theoretical bounds.

The hydrodynamics of wetting involves a singularity of viscous stress, and its microscopic regularisation ultimately determines the speed at which contact lines move over a surface. In a recent paper, Luo & Gao (J. Fluid Mech., vol. 1019, 2025, A52) explore a new analytical solution, based on which they construct a model for ‘slippery wedge flow’. This lucid approach provides an accurate description of viscous wetting flows in the presence of slip, without the usual restriction to small contact angles, and offers a quantitative multiscale formalism for slippery contact lines.

Solidification of droplets is of great importance to various technological applications, drawing considerable attention from scientists aiming to unravel the fundamental physical mechanisms. In the case of multicomponent droplets undergoing solidification, the emergence of concentration gradients may trigger significant interfacial flows that dominate the freezing dynamics. Here, we experimentally investigate the fascinating interfacial freezing dynamics of supercooled ethanol–water droplets, accompanied with the migration and growth of massive ice particles. We reveal that this unique freezing dynamics is driven by solidification-induced solutal Marangoni flow within the droplets. Our model, which incorporates the temperature- and concentration-dependent properties of the ethanol–water mixture, quantitatively predicts both the migration velocity and the growth rate of the ice particles. The former is determined by the solutal Marangoni flow velocity, while the latter is governed by a balance between the latent heat release and the enhanced thermal dissipation by the Marangoni flow. Moreover, we show that the final wrapping state of droplets can be modulated by the concentration of ethanol. Our findings may pave the way for novel insights into the physicochemical hydrodynamics of multicomponent liquids undergoing phase transitions.

In this paper, we numerically investigate the orbit dynamics of three-dimensional symmetric Janus drops in shear flow using an improved ternary-fluids phase field method, focusing on how drop deformation and initial orientation affect the orbit drift of two configurations of Janus drops: dumbbell-shaped and near-spherical. We find that the motion of dumbbell-shaped drops eventually evolves into tumbling, while near-spherical drops attain stable spinning. We attribute this bifurcation in orbit drift to contrasting deformation dynamics and shape-dependent hydrodynamics of the two configurations. Specifically, the drift bifurcation is closely related to the aspect ratio of Janus drops at equilibrium, giving rise to two distinct mechanisms: (1) coupling between outer interface deformation and the surrounding flow field; and (2) interplay between inner interface deformation and vortices enclosed within the drop. In addition, we observe that for the dumbbell-shaped Janus drops with different aspect ratios, their tumbling dynamics resembles ellipsoids in shear flow. Moreover, the trajectories of the dumbbell-shaped Janus drops during orbit drift collapse onto a universal curve, independent of their initial orientations, and significant deformation and inertia accelerate the orbit transition. To quantitatively evaluate the effect of drop deformation on the orbit drift of the dumbbell-shaped Janus drops, we propose an effective aspect ratio model based on the drop shapes at equilibrium and at the maximum elongation. By incorporating the effective aspect ratio into Jeffery’s theory for solid particles, we accurately predict the rotation period and angular velocity of Janus drops in the tumbling regime and during the orbit drift, especially for drops with linear deformation. Moreover, the orbit parameter $C$ is found to vary exponentially with time for drops with linear deformation, while the time variation of $C$ transits from one exponential function to another for drops with nonlinear deformation.

The thermal interactions of liquid droplets impacting a moving substrate are investigated, combining theoretical modelling with experimental validation. An analytical model is developed to predict the time-evolving contact temperature and heat flux at the droplet–substrate interface. Accounting for the convective heat transport induced by the impacting drop, the model incorporates a finite thermal contact resistance, which is a critical parameter that was often neglected in earlier studies for drop impact. High-speed, spatially resolved infrared thermography is used to record the two-dimensional, transient temperature evolution at the droplet–substrate interface during drop impact on a rotating disc. Measured temperature maps are used for numerical simulations to reconstruct local interfacial heat fluxes. The model is validated for different droplet diameters, substrate velocities and thermal conditions. The findings demonstrate that the substrate velocity and droplet diameter have negligible influence on the thermal behaviour within the tested parameter space.

The deposition of droplets onto a swollen polymer network induces the formation of a wetting ridge at the contact line. Current models typically consider either viscoelastic effects or poroelastic effects, while polymeric gels often exhibit both properties. In this study, we investigate the growth of the wetting ridge using a comprehensive large-deformation theory that integrates both dissipative mechanisms – viscoelasticity and poroelasticity. In the purely poroelastic case, following an initial instantaneous incompressible deformation, the growth dynamics exhibits scale-free behaviour, independent of the elastocapillary length or system size. A boundary layer of solvent imbibition between the solid surface (in contact with the reservoir) and the region of minimal chemical potential is created. At later times, the ridge equilibrates on the diffusion time scale given by the elastocapillary length. When viscoelastic properties are incorporated, our findings show that, during the early stages (prior to the viscoelastic relaxation time scale), viscoelastic effects dominate the growth dynamics of the ridge and solvent transport is significantly suppressed. Beyond the relaxation time, the late-time dynamics closely resembles that of the purely poroelastic case. These findings are discussed in light of recent experiments, showing how our approach offers a new interpretation framework for wetting of polymer networks of increasing complexity.

A thin, evaporating sessile droplet with a pinned contact line containing inert particles is considered. In the limit in which the liquid flow decouples from the particle transport, we discuss the interplay between particle advection, diffusion and adsorption onto the solid substrate on which the droplet sits. We perform an asymptotic analysis in the physically relevant regime in which the Péclet number is large, i.e. ${\textit{Pe}}\gg 1$, so that advection dominates diffusion in the droplet except in a boundary layer near the contact line, and in which the ratio of the particle velocities due to substrate adsorption and diffusion is at most of order unity as ${\textit{Pe}}\rightarrow \infty$. We use the asymptotic model alongside numerical simulations to demonstrate that substrate adsorption leads to a different leading-order distribution of particle mass compared with cases with negligible substrate adsorption, with a significant reduction of the mass in the suspension – the nascent coffee ring reported in Moore et al. (J. Fluid Mech., vol. 920, 2021, A54). The redistribution leads to an extension of the validity of the dilute suspension assumption, albeit at the cost of breakdown due to the growth of the deposited layer, which are important considerations for future models that seek to accurately model the porous deposit regions.

We highlight the complete transition from liquid-wall-film instability of an annular gas–liquid flow inside a nozzle to spray formation at the trailing edge, aiming to identify two distinct flow regimes of ripple waves and disturbance waves and to clarify their distinct fragmentation mechanisms. Experiments conducted under strictly controlled boundary conditions support our theoretical analysis, revealing that the onset of disturbance waves coincides with the liquid-film Weber number (${\textit{We}}$) of unity, marking a significant change in following fragmentation dynamics. For ${\textit{We}}\lt 0.5$, the liquid wall film forms three-dimensional ripple waves driven by the superposition of Kelvin–Helmholtz and Rayleigh–Taylor (RT) instabilities, with no disturbance waves present. At the trailing edge, the liquid film temporarily accumulates, extends into isolated ligaments along the axial direction via RT instability, and subsequently fragments into droplets through Plateau–Rayleigh instability, displaying a weak coupling between ripple wave dynamics and fragmentation. In contrast, for ${\textit{We}}\gt 0.5$, disturbance waves with long wavelengths and large amplitudes become prominent, superimposed on the base ripple waves. As these disturbance waves reach the trailing edge, they are spontaneously ejected as liquid sheets at the same frequency, forming transverse rims through RT instability and rapidly disintegrating into fine droplets. This regime demonstrates a direct coupling between disturbance-wave dynamics and fragmentation.

The forced breakup of liquid jets in ambient gas surroundings is studied systematically through numerical simulations and theoretical analyses, with particular emphasis on characterising the response modes of jet breakup across wide ranges of perturbation frequency and amplitude. Simulations reveal that the breakup of liquid jet can be effectively synchronised with external actuation within a frequency range encompassing the natural breakup frequency, thereby enabling the generation of highly uniform droplets. As the perturbation frequency exceeds an upper critical value, the external perturbation cannot dominate the jet breakup, while below a lower critical frequency, the jet breaks up with multiple droplets generated within one period. A high perturbation amplitude can result in liquid accumulation, leading to the formation of a pancake-shaped jet configuration. Through spectrum analyses, the development of jet interface perturbations under different response modes is elucidated, revealing the competition between the natural frequency and the external frequency. A linear instability analysis of a liquid jet is performed, which successfully predicts the synchronised frequency range by comparing the breakup time between the free liquid jet and the actuated jet, along with the variation tendencies of jet breakup length with varying perturbation frequency, amplitude and jet velocity. Quantitative numerical results demonstrate that in the case of multiple droplet generation under low perturbation frequency, the rear droplet maintains a higher velocity than its leading counterpart and the emergence of a high-pressure zone at the leading edge of a droplet train facilitates the droplet coalescence. Furthermore, the study introduces an innovative approach by superimposing periodic pulses onto the sinusoidal perturbation waveform, enabling active modulation of multiple droplet merging dynamics. This fundamental study is intended to offer valuable guidance for the on-demand generation of droplets in various industrial applications.

The present study deals with the electrophoresis of a non-polarizable droplet with irreversibly adsorbed ionic surfactants suspended in monovalent or multivalent electrolyte solutions. The impact of the non-uniform surface charge density, governed by the interfacial surfactant concentration, along with Marangoni, hydrodynamic and Maxwell stresses on droplet electrophoresis is analysed. At a large ionic concentration, the hydrodynamic steric interactions and correlations among finite-sized ions manifest. In this case the viscosity of the medium rises as the local volume fraction of the finite-sized ions is increased. The governing equations, incorporating these short-range effects, are solved numerically based on the regular linear perturbation analysis under a weak applied electric field consideration. We find that the electrophoretic velocity consistently decreases with an increase in the droplet-to-electrolyte viscosity ratio due to the Marangoni stress caused by inhomogeneous surfactant distribution. This monotonic relationship with droplet viscosity is absent for the case of constant surface charge density, where a low-viscosity droplet may exhibit a lower mobility than a high-viscosity droplet. In the presence of ionic surfactant, a continuous variation of mobility with surfactant concentration is found. For a monovalent electrolyte, mobility decreases significantly at an elevated ionic concentration due to the short-range effects described above. Reversal in mobility is observed in multivalent electrolytes due to the correlations among finite-sized ions, attributed to overscreening of surface charge and formation of a coion-rich layer within the electric double layer. In this case a toroidal vortex develops adjacent to the droplet and the reversed mobility enhances as the Marangoni number is increased. This mobility reversal is delayed for low-viscosity droplets.

Robust surfaces capable of reducing flow drag, controlling heat and mass transfer, and resisting fouling in fluid flows are important for various applications. In this context, textured surfaces impregnated with a liquid lubricant show promise due to their ability to sustain a liquid–liquid interface that induces slippage. However, theoretical and numerical studies suggest that the slippage can be compromised by surfactants in the overlying fluid, which contaminate the liquid–liquid interface and generate Marangoni stresses. In this study, we use Doppler-optical coherence tomography, an interferometric imaging technique, combined with numerical simulations to investigate how surfactants influence the slip length of lubricant-infused surfaces with longitudinal grooves in a laminar flow. Surfactants are endogenously present in the contrast agent (milk) which is added to the working fluid (water). Local measurements of slip length at the liquid–liquid interface are significantly smaller than theoretical predictions for clean interfaces (Schönecker & Hardt 2013). In contrast, measurements are in good agreement with numerical simulations of fully immobilized interfaces, indicating that milk surfactants adsorbed at the interface are responsible for the reduction in slippage. This work provides the first experimental evidence that liquid–liquid interfaces within textured surfaces can become immobilised in the presence of surfactants and flow.

This study investigates droplet impact on elastic plates using a two-phase lattice Boltzmann method in both two-dimensional (2-D) and three-dimensional (3-D) configurations, with a focus on rebound dynamics and contact time. The 2-D simulations reveal three distinct rebound modes – conventional bounce, early bounce and rim rising – driven by fluid–structure interaction. Among them, the early bounce mode uniquely achieves a significant reduction in contact time, occurring only at moderate plate oscillation frequency. Momentum analysis shows a non-monotonic relationship between vertical momentum transfer and rebound efficiency: increased momentum does not necessarily promote rebound if it concentrates in a central jet, which contributes minimally to lift-off. This introduces a novel rebound mechanism governed by momentum distribution morphology rather than total magnitude. A theoretical model treating the droplet–plate system as coupled oscillators is developed to predict contact time in the early bounce regime, showing good agreement with numerical results. The mechanism and model are further validated through fully 3-D simulations, confirming the robustness of the findings.

Pendant drops appear in many engineering applications, such as inkjet printing and optical tensiometry, and they have also been the subject of studies of droplet–particle interaction. While the hydrostatics of pendant drops has been studied extensively, the influence of external flow disturbances has received limited attention. This research aims to incorporate aerodynamic factors into the understanding of pendant drop behaviour. Employing a simplified model, an irrotational flow aligned with the drop’s axis is derived from a distribution of singularity elements within the drop. The drop’s equilibrium shape is then determined using a numerical model that couples the flow field with the Young–Laplace equation. The model’s predictions are compared to droplet images captured via high-speed shadowgraph in a vertical wind tunnel, showing good agreement with the experimentally observed shapes. Additionally, under certain flow conditions, the drop exhibits instability in the form of periodic pendulum-like motion. This instability was linked to two distinct critical drop heights, and the corresponding stability criterion was mathematically derived from the numerical model. Our theoretical and experimental findings provide the first quantitative description of the equilibrium shape and stability criterion of pendant drops under the influence of external flow.

We investigate the shape of a tin sheet formed from a droplet struck by a nanosecond laser pulse. Specifically, we examine the dynamics of the process as a function of laser beam properties, focusing on the outstanding puzzle of curvature inversion: tin sheets produced in experiments and state-of-the-art extreme ultraviolet (EUV) nanolithography light sources curve in a direction opposite to previous theoretical predictions. We resolve this discrepancy by combining direct numerical simulations with experimental data, demonstrating that curvature inversion can be explained by an instantaneous pressure impulse with low kurtosis. Specifically, we parametrise a dimensionless pressure width, $ W$, using a raised cosine function and successfully reproduce the experimentally observed curvature over a wide range of laser-to-droplet diameter ratios, $ 0.3 \lt d/D_0 \lt 0.8$. The simulation process described in this work has applications in the EUV nanolithography industry, where a laser pulse deforms a droplet into a sheet, which is subsequently ionised by a second pulse to produce EUV-emitting plasma.

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.

We study the rebound of drops impacting non-wetting substrates at low Weber number ($\textit{We}$) through experiment, direct numerical simulation and reduced-order modelling. Submillimetre-sized drops are normally impacted onto glass slides coated with a thin viscous film that allows them to rebound without contact line formation. Experiments are performed with various drop viscosities, sizes and impact velocities, and we directly measure metrics pertinent to spreading, retraction and rebound using high-speed imaging. We complement experiments with direct numerical simulation and a fully predictive reduced-order model that applies natural geometric and kinematic constraints to simulate the drop shape and dynamics using a spectral method. At low $\textit{We}$, drop rebound is characterised by a weaker dependence of the coefficient of restitution on $\textit{We}$ than in the more commonly studied high-$\textit{We}$ regime, with nearly $\textit{We}$-independent rebound in the inertio-capillary limit, and an increasing contact time as $\textit{We}$ decreases. Drops with higher viscosity or size interact with the substrate longer, have a lower coefficient of restitution and stop bouncing sooner, in good quantitative agreement with our reduced-order model. In the inertio-capillary limit, low-$\textit{We}$ rebound has nearly symmetric spreading and retraction phases and a coefficient of restitution near unity. Increasing $\textit{We}$ or viscosity breaks this symmetry, coinciding with a drop in the coefficient of restitution and an increased dependence on $\textit{We}$. Lastly, the maximum drop deformation and spreading are related through energy arguments, providing a comprehensive framework for drop impact and rebound at low $\textit{We}$.

We must address active matter in the context of soft boundaries to bridge the gap between our understanding of active matter and the dynamics of biological systems (represented as active matter) under natural conditions. However, the physics of such active drops (matter) in contact with a soft and deformable surface has remained elusive. In this paper, we attempt to fill this gap and develop a theory for soft, active wetting. Our theory, which accounts for the various free energies for passive substrate and active drops as well as the active stresses, provides an equilibrium description of (active) particle orientation inside the drop and an equilibrium shape of the drop–soft-solid system. We obtain an analytical equation relating the activity to the internal pressure of an active drop. The equilibrium calculation further yields an ordered state of the polarisation field inside the drop. As compared to the non-active drops, the active drops with extensile activity press more into the soft surface, while the active drops with contractile activity either rise out of the soft surface (for smaller magnitude of negative activity) or make the soft surface bulge (for larger magnitude of negative activity). Finally, the three-phase contact line undergoes a rotation that depends on the strength of activity. These findings shed light on the manner in which the active stresses interact with surface tension and elasticity at the fundamental level.