Liquid droplet dynamics are widely used in biological and engineering applications, which contain complex interfacial instabilities and pattern formation such as droplet merging, splitting and transport. This paper studies a class of mean field control formulations for these droplet dynamics, which can be used to control and manipulate droplets in applications. We first formulate the droplet dynamics as gradient flows of free energies in modified optimal transport metrics with nonlinear mobilities. We then design an optimal control problem for these gradient flows. As an example, a lubrication equation for a thin volatile liquid film laden with an active suspension is developed, with control achieved through its activity field. Lastly, we apply the primal–dual hybrid gradient algorithm with high-order finite-element methods to simulate the proposed mean field control problems. Numerical examples, including droplet formation, bead-up/spreading, transport, and merging/splitting on a two-dimensional spatial domain, demonstrate the effectiveness of the proposed mean field control mechanism.

The confined Stokesian dynamics (CSD) algorithm recently reported equilibrium properties but was missing hydrodynamic functions required for suspension stress and non-equilibrium properties. In this first of a two-part series, we expand the CSD algorithm to model the traceless part of the stress tensor. To obtain quantities needed to solve the integral expressions for the stress, we developed a general method to solve Stokes’ equations in bispherical coordinates. We calculate the traceless stress tensor for arbitrary particle-to-enclosure size ratio. We next compute rheology of a confined suspension by implementing the stresslet hydrodynamic coefficients into CSD, yielding the deviatoric part of the many-body hydrodynamic stresslet. We employed energy methods to relate this stresslet to the high-frequency viscosity of the confined suspension, finding an increase in viscous dissipation with crowding and confinement well beyond the unconfined value. We show that confinement effects on viscosity are dominated by near-field interactions between the particles that reside very near the cavity wall (rather than particle–wall interactions). Surprisingly, this near-field effect is stronger than the viscosity of an unconfined suspension, showing that entropic exclusion driven by the wall sets up many lubrication interactions that then generate strong viscous dissipation. The limiting case of a particle near a flat wall reveals a correction to prior literature. The theory presented in this work can be expanded to study the Brownian contribution to the viscosity of confined suspensions in and away from equilibrium. In part 2, we report the osmotic pressure, via the trace of the stress tensor.

The nonlinear waves in a sheared liquid film on a horizontal plate at small Reynolds numbers are examined by theoretical and numerical approaches. The analysis employs the long-wave approximation along with finite difference schemes. The results show that the surface tension can suppress disturbances and prevent the occurrence of singularities. While the film flow is driven by the shear stress on the interface, its instability highly depends on the magnitude and direction of gravity. Specifically, when the direction of gravity is opposite to the wall-normal direction, perturbations are stabilized by gravity. In contrast, when these two directions are the same, the gravitational force is destabilizing, and stationary travelling waves can exist if a balance is reached between the effects of gravity and surface tension. For the steady solitary waves, there are quasi-periodic oscillations occurring between two stationary points, indicating the presence of heteroclinic trajectories. For periodic waves, the evolutions are sensitive to several parameters and initial disturbances, while one steady-state wave exhibits a sine function-like behaviour.

Diffusion-driven flow is a boundary layer flow arising from the interplay of gravity and diffusion in density-stratified fluids when a gravitational field is non-parallel to an impermeable solid boundary. This study investigates diffusion-driven flow within a nonlinearly density-stratified fluid confined between two tilted parallel walls. We introduce an asymptotic expansion inspired by the centre manifold theory, where quantities are expanded in terms of derivatives of the cross-sectional averaged stratified scalar (such as salinity or temperature). This technique provides accurate approximations for velocity, density and pressure fields. Furthermore, we derive an evolution equation describing the cross-sectional averaged stratified scalar. This equation takes the form of the traditional diffusion equation but replaces the constant diffusion coefficient with a positive-definite function dependent on the solution's derivative. Numerical simulations validate the accuracy of our approximations. Our investigation of the effective equation reveals that the density profile depends on a non-dimensional parameter denoted as $\gamma$ representing the flow strength. In the large $\gamma$ limit, the system is approximated by a diffusion process with an augmented diffusion coefficient of $1+\cot ^{2}\theta$, where $\theta$ signifies the inclination angle of the channel domain. This parameter regime is where diffusion-driven flow exhibits its strongest mixing ability. Conversely, in the small $\gamma$ regime, the density field behaves like pure diffusion with distorted isopycnals. Lastly, we show that the classical thin film equation aligns with the results obtained using the proposed expansion in the small $\gamma$ regime but fails to accurately describe the dynamics of the density field for large $\gamma$.

We investigate theoretically the steady incompressible viscoelastic flow in a rigid axisymmetric cylindrical pipe with varying cross-section. We use the Oldroyd-B viscoelastic constitutive equation to model the fluid viscoelasticity. First, we derive exact general formulae: for the total average pressure-drop as a function of the wall shear rate and the viscoelastic axial normal extra-stress; for the viscoelastic extra-stress tensor and the Trouton ratio as functions of the fluid velocity on the axis of symmetry; and for the viscoelastic extra-stress tensor along the wall in terms of the shear rate at the wall. Then we exploit the classic lubrication approximation, valid for small values of the square of the aspect ratio of the pipe, to simplify the original governing equations. The final equations are solved analytically using a regular perturbation scheme in terms of the Deborah number, De, up to eighth order in De. For a hyperbolically shaped pipe, we reveal that the reduced pressure-drop and the Trouton ratio can be recast in terms of a modified Deborah number, Dem, and the polymer viscosity ratio, η, only. Furthermore, we enhance the convergence and accuracy of the eighth-order solutions by deriving transformed analytical formulae using Padé diagonal approximants. The results show the decrease of the pressure drop and the enhancement of the Trouton ratio with increasing Dem and/or increasing η. Comparison of the transformed solutions with numerical simulations of the lubrication equations using pseudospectral methods shows excellent agreement between the results, even for high values of Dem and all values of η, revealing the robustness, validity and efficiency of the theoretical methods and techniques developed in this work. Last, it is shown that the exact solution for the Trouton ratio gives a well-defined and finite solution for any value of Dem and reveals the reason for the failure of the corresponding high-order perturbation series for Dem > 1/2.

We present an experimental investigation aimed at understanding the effects of surface roughness on the time-mean drag coefficient ($\bar {C}_{D}$) of finite-span cylinders ($\text {span/diameter} = \text {aspect ratio}$, $0.51 \le AR \le 6.02$) freely rolling without slip on an inclined plane. While lubrication theory predicts an infinite drag force for ideally smooth cylinders in contact with a smooth surface, experiments yield finite drag coefficients. We propose that surface roughness introduces an effective gap ($G_{eff}$) resulting in a finite drag force while allowing physical contact between the cylinder and the plane. This study combines measurements of surface roughness for both the cylinder and the plane panel to determine a total relative roughness ($\xi$) that can effectively describe $G_{eff}$ at the point of contact. It is observed that the measured $\bar {C}_{D}$ increases as $\xi$ decreases, aligning with predictions of lubrication theory. Furthermore, the measured $\bar {C}_{D}$ approximately matches combined analytical and numerical predictions for a smooth cylinder and plane when the imposed gap is approximately equal to the mean peak roughness ($R_p$) for rough cylinders, and one standard deviation peak roughness ($R_{p, 1\sigma }$) for relatively smooth cylinders. As the time-mean Reynolds number ($\overline {Re}$) increases, the influence of surface roughness on $\bar {C}_{D}$ decreases, indicating that wake drag becomes dominant at higher $\overline {Re}$. The cylinder aspect ratio ($AR$) is found to have only a minor effect on $\bar {C}_{D}$. Flow visualisations are also conducted to identify critical flow transitions and these are compared with visualisations previously obtained numerically. Variations in $\xi$ have little effect on the cylinder wake. Instead, $AR$ was observed to have a more pronounced effect on the flow structures observed. The Strouhal number ($St$) associated with the cylinder wake shedding was observed to increase with $\overline {Re}$, while demonstrating little dependence on $AR$.

We derive a generalised asymptotic model for the flow of a thin fluid film over an arbitrarily parameterised non-axisymmetric curved substrate surface based on the lubrication approximation. In addition to surface tension, gravity and centrifugal force, our model incorporates the effects of the Coriolis force and disjoining pressure, together with a non-uniform initial condition, which have not been widely considered in existing literature. We use this model to investigate the impact of the Coriolis force and fingering instability on the spreading of a non-axisymmetric spin-coated film at a range of substrate angular velocities, first on a flat substrate, and then on parabolic cylinder- and saddle-shaped curved substrates. We show that, on flat substrates, the Coriolis force has a negligible impact at low angular velocities, and at high angular velocities results in a small deflection of fingers formed at the contact line against the direction of substrate rotation. On curved substrates, we demonstrate that, as the angular velocity is increased, spin-coated films transition from being dominated by gravitational drainage with no fingering to spreading and fingering in the direction with the greatest component of centrifugal force tangent to the substrate surface. For both curved substrates and all angular velocities considered, we show that the film thickness and total wetted substrate area remain similar over time to those on a flat substrate, with the key difference being the shape of the spreading droplet.

A phenomenological model is proposed to estimate the initial thickness of the liquid microlayer forming beneath a vapour bubble growing on a solid surface upon nucleate boiling. The model employs an analogy between the microlayer formation and the classic plate withdrawal problem. It calculates the microlayer thickness by considering it as a Landau–Levich film, where the thickness is a function of the meniscus speed and radius of curvature. Given the nearly hemispherical shape of the bubble during the early growth stage when the microlayer is first deposited, we assume that the meniscus speed can be approximated by the bubble expansion rate, and estimate the meniscus curvature using the Rayleigh equations. Unlike previous theories that assume that the bubble radius growth is proportional to the square root of time, the proposed model does not rely on any specific law of growth for vapour bubbles. The model is validated for predicting the microlayer thickness in water and ethanol, showing good agreement with experimental measurements and empirical correlations. Subsequent analyses of the microlayer interface profile address inconsistent reports – some described a wedge-like shape, whereas others reported a slight outward curvature with decreasing thickness in the outer region. This discrepancy is attributed to a reduction in the expansion rate of the microlayer's outer edge, particularly when the bubble reaches its maximum width. Our model provides insights into microlayer dynamics, essential to boiling heat transfer, as the evaporative heat flux through the microlayer is very sensitive to its initial thickness.

A viscous, lubrication-like response can be triggered in a thin film of fluid squeezed between a rigid flat surface and the tip of an incoming projectile. We develop a scaling for this viscous approach stage of fluid-mediated normal impact, applicable to soft impactors. Under the assumption of mediating fluid being incompressible, the impacting solid displays two limit regimes: one dominated by elasticity, and the other by inertia. The transition between the two is predicted by a dimensionless parameter, which can be interpreted as the ratio between two time scales that are the time that it takes for the surface waves to warn the leading edge of the impactor of the forthcoming impact, and the characteristic duration of the final viscous phase of the approach. Additionally, we elucidate why nearly incompressible solids feature (a) substantial ‘gliding’ prior to contact at the transition between regimes, (b) the largest size of entrapped bubble between the deformed tip of the impactor and the flat surface, and (c) a sudden drop in entrapped bubble radius past the transition between regimes. Finally, we argue that the above time scale ratio (a dimensionless number) can govern the different dynamics reported experimentally for a fluid droplet as a function of its viscosity and surface tension.

The interaction between deformable surfaces and oscillatory driving is known to produce complex secondary time-averaged flows due to inertial and elastic nonlinearities. Here, we revisit the problem of oscillatory flow in a cylindrical tube with a deformable wall, and analyse it under a long-wave theory for small deformations, but for arbitrary Womersley numbers. We find that the oscillatory pressure does not vary linearly along the length of a deformable channel, but instead decays exponentially with spatial oscillations. We show that this decay occurs over an elasto-visco-inertial length scale that depends on the material properties of the fluid and the elastic walls, the geometry of the system, and the frequency of the oscillatory flow, but is independent of the amplitude of deformation. Inertial and geometric nonlinearities associated with the elastic deformation of the channel drive a time-averaged secondary flow. We quantify the flow using numerical solutions of the perturbation theory, and gain insight by developing analytic approximations. The theory identifies a complex non-monotonic dependence of the time-averaged flux on the elastic compliance and inertia, including a reversal of the flow. Finally, we show that our analytic theory is in excellent quantitative agreement with the three-dimensional direct numerical simulations of Pande et al. (Phys. Rev. Fluids, vol. 8, no. 12, 2023, 124102).

The flow of a nematic liquid crystal in a Hele-Shaw cell with an electrically controlled viscous obstruction is investigated using both a theoretical model and physical experiments. The viscous obstruction is created by temporarily electrically altering the viscosity of the nematic in a region of the cell across which an electric field is applied. The theoretical model is validated experimentally for a circular cylindrical obstruction, demonstrating user-controlled flow manipulation of an anisotropic liquid within a heterogeneous single-phase microfluidic device.

Instability and rupture dynamics of a liquid nano-thread, subjected to external hydrodynamic perturbations, are captured by a stochastic lubrication equation (SLE) incorporating thermal fluctuations via Gaussian white noise. Linear instability analysis of the SLE is conducted to derive the spectra and distribution functions of thermal capillary waves influenced by external perturbations and thermal fluctuations. The SLE is also solved numerically using a second-order finite difference method with a correlated noise model. Both theoretical and numerical solutions, validated through molecular dynamics, indicate that surface tension forces due to specific external perturbations overcome the random effects of thermal fluctuations, determining both the thermal capillary waves and the evolution of perturbation growth. The results also show two distinct regimes: (i) the hydrodynamic regime, where external perturbations dominate, leading to uniform ruptures, and (ii) the thermal-fluctuation regime, where external perturbations are surpassed by thermal fluctuations, resulting in non-uniform ruptures. The transition between these regimes, modelled by a criterion developed from linear instability theory, exhibits a strong dependence on the amplitudes and wavenumbers of the external perturbations.

We study the dynamics of hydraulic fracturing of an elastic solid in a Hele-Shaw cell. Compared with hydraulic fractures in an infinite elastic bulk, the viscous resistance comes mainly from the drag by the two parallel plates that forms the Hele-Shaw cell rather than by the fluid–solid interface. Such a feature leads to a different nonlinear differential–integral system that describes the coupled evolution of the fracture shape and pressure field. Our theory leads to hydraulic fractures of cusp shapes in the neighbourhood of the fracture tip, which is consistent with recent experimental observations. Accordingly, there exists no pressure singularity at the location of the fracture tip, which is also fundamentally different from our previous understandings of hydraulic fracturing of elastic solids.

Immersed nonlinear elements are prevalent in biological systems that require a preferential flow direction, such as the venous and the lymphatic system. We investigate here a certain class of models where the fluid is driven by peristaltic pumping and the nonlinear elements are ideal valves that completely suppress backflow. This highly nonlinear system produces discontinuous solutions that are difficult to study. We show that, as the density of valves increases, the pressure and flow are well approximated by a continuum of valves which can be analytically treated, and we demonstrate through numeric simulation that the approximation works well even for intermediate valve densities. We find that the induced flow is linear in the peristaltic amplitude for small peristaltic forces and, in the case of sinusoidal peristalsis, is independent of pumping direction. Despite the continuum approximation used, the physical valve density is accounted for by modifying the resistance of the fluid appropriately. The suppression of backflow causes a net benefit in adding valves when the valve density is low, but once the density is high enough, valves predominately suppress forward flow, suggesting there is an optimum number of valves per wavelength. The continuum model for peristaltic pumping through an array of valves presented in this work can eventually provide insights about the design and operating principles of complex flow networks with a broad class of nonlinear elements.

When a viscous fluid partially fills a Hele-Shaw channel, and is pushed by a pressure difference, the fluid interface is unstable due to the Saffman–Taylor instability. We consider the evolution of a fluid region of finite extent, bounded between two interfaces, in the limit that the interfaces are close, that is, when the fluid region is a thin liquid filament separating two gases of different pressure. In this limit, we derive a second-order ‘thin-filament’ model that describes the normal velocity of the filament centreline, and evolution of the filament thickness, as functions of the thickness, centreline curvature and their derivatives. We show that the second-order terms in this model, that include the effect of transverse flow along the filament, are necessary to regularise the instability. Numerical simulation of the thin-filament model is shown to be in accordance with level-set computations of the complete two-interface model. Solutions ultimately evolve to form a bubble of rapidly increasing radius and decreasing thickness.

Lubrication theory is adapted to incorporate the large normal stresses that occur for order-one Deborah numbers, $De$, the ratio of the relaxation time to the residence time. Comparing with the pressure drop for a Newtonian viscous fluid with a viscosity equal to that of an Oldroyd-B fluid in steady simple shear, we find numerically a reduced pressure drop through a contraction and an increased pressure drop through an expansion, both changing linearly with $De$ at high $De$. For a constriction, there is a smaller pressure drop that plateaus at high $De$. For a contraction, much of the change in pressure drop occurs in the stress relaxation in a long exit channel. An asymptotic analysis for high $De$, based on the idea that normal stresses are stretched by an accelerating flow in proportion to the square of the velocity, reveals that the large linear changes in pressure drop are due to higher normal stresses pulling the fluid through the narrowest gap. A secondary cause of the reduction is that the elastic shear stresses do not have time to build up to their steady-state equilibrium value while they accelerate through a contraction. We find for a contraction or expansion that the high $De$ analysis works well for $De>0.4$.

When a sphere rotates near a rigid boundary coated with a thin layer of viscous liquid, ‘tracks’ are generated both behind and over the sphere. This paper describes a theory for the simplest one-track configuration which can occur under particular experimental conditions. The theoretical predictions are in good agreement with experimental observations.

Rain drops form in clouds by collision of submillimetric droplets falling under gravity: larger drops fall faster than smaller ones and collect them on their path. The puzzling stability of fogs and non-precipitating warm clouds with respect to this avalanche mechanism has been a longstanding problem. How can droplets of diameter around 10 $\mathrm {\mu }$m have a low collision probability, inhibiting the cascade towards larger and larger drops? Here we review the dynamical mechanisms that have been proposed in the literature and quantitatively investigate the frequency of drop collisions induced by Brownian diffusion, electrostatics and gravity, using an open-source Monte Carlo code taking all of them into account. Inertia dominates over aerodynamic forces for large drops, when the Stokes number is larger than $1$. Thermal diffusion dominates over aerodynamic forces for small drops, when the Péclet number is smaller than $1$. We show that there exists a range of size (typically 3–30 $\mathrm {\mu }$m for water drops in air) where neither inertia nor Brownian diffusion are significant, leading to a gap in the collision rate. The effect is particularly important, due to the lubrication film forming between the drops immediately before collision, and secondarily to the long-range aerodynamic interaction. Two different mechanisms regularise the divergence of the lubrication force at vanishing separation: the transition to a non-continuum regime in the lubrication film, when the separation is comparable to the mean free path of air, and the induction of a flow inside the drops due to shear at their surfaces. In the gap between inertia-dominated and diffusion-dominated regimes, dipole–dipole electrostatic interactions becomes the major effect controlling the efficiency of drop collisions.

We investigate theoretically and numerically the impact of an elastic sphere on a rigid wall in a viscous fluid. Our focus is on the dynamics of the contact, employing the soft lubrication model in which the sphere is separated from the wall by a thin liquid film. In the limit of large sphere inertia, the sphere bounces and the dynamics is close to the Hertz theory. Remarkably, the film thickness separating the sphere from the wall exhibits non-trivial self-similar properties that vary during the spreading and retraction phases. Leveraging these self-similar properties, we establish the energy budget and predict the coefficient of restitution for the sphere. The general framework derived here opens many perspectives to study the lubrication film in impact problems.

Magnetic soap films under the forcing of an inhomogeneous magnetic field are governed by a wide range of interconnected physics. The study of magnetic soap films requires the development of comprehensive models to support experimental observations. In this study, the thin film approximation is applied to the Navier–Stokes equations to derive a model for the film thickness of magnetic soap films that incorporates the effects of interfacial mobility, surfactant transport and magnetite nanoparticle (NP) transport. This derived model consists of a coupled system of equations for the film thickness, interfacial velocity, interfacial surfactant concentration and magnetite NP concentration. Simulations are performed for both soap films and magnetic soap films by solving the system of equations using the Galerkin finite element method, and results are compared with experiments. Simulation results highlight that interfacial flows can dominate the rate of film thinning and that accounting for the dependence of the magnetisation on the local magnetite NP concentration can influence the predicted speed of magnetically driven flows. Furthermore, simulation results demonstrate that the model is able to predict marginal regeneration in qualitative agreement with the experiments for soap films; the model also predicts the same flow pattern as seen in the experiments for magnetic soap films. Overall, this study advances the state of soap film and magnetic soap film modelling and will contribute to acquiring control over the drainage and stability of magnetic soap films in the long term.