We apply control-based continuation, an experimental method used to explore the bifurcation structure of nonlinear systems, to a spatially extended system by examining the equilibrium configurations of a bubble placed at the center of a straining flow. Feedback control is delivered by injecting/withdrawing fluid at actuation points and is chosen according to real-time observations of the bubble as viewed from above. The amplitude of actuation is used to locate its unstable equilibria, and thus enable their observation. We map the unstable steady shapes of the bubble as a function of straining flow rate and demonstrate excellent agreement of the resulting experimental bifurcation diagram with numerical predictions.

Liquid films coating vertical cylinders can form annular liquid collars which translate downwards under gravity. We investigate the dynamics of a thin viscoplastic liquid film coating the interior or exterior of a vertical cylindrical tube, quantifying how the yield stress modifies both the Rayleigh-Plateau instability leading to collar formation and the translation of collars down the tube. We use thin-film theory to derive an evolution equation for the layer thickness, which we solve numerically to examine the nonlinear dynamics. Instability and collar formation occur when gravity is sufficiently strong to make the fluid yield initially. We use matched asymptotics to derive a model describing the quasi-steady translation of a slender liquid collar when the Bond number is small. The structure of the asymptotic solution for a viscoplastic collar shares some features with the Newtonian version, but there are several novel asymptotic regions that emerge at the two ends of the collar. The global force balance, which determines the collar's speed, is modified by a leading-order contribution from viscous drag in the collar when the liquid is viscoplastic. We use the asymptotic model to describe slow changes in collar volume when the film thicknesses ahead of, and behind, the collar are unequal. When the film thickness ahead of the collar is less than a critical value that we determine, viscoplastic collars adjust their volume and reach a steadily-translating state. This contrasts with the Newtonian problem, where the only state in which steady translation occurs is unstable to small changes in the film thickness.

Airway constriction and blockage in obstructive lung diseases cause ventilation heterogeneity and create barriers to effective drug deposition. Established computational particle-deposition models have not accounted for these impacts of disease. We present a new particle-deposition model that calculates ventilation based on the resistance of each airway, such that ventilation responds to airway constriction. The model incorporates distal airway constrictions representative of cystic fibrosis, allowing us to investigate the resulting impact on patterns of deposition. Unlike previous models, our model predicts how constrictions affect deposition in airways throughout the lungs, not just in the constricted airways. Deposition is reduced in airways directly distal and proximal to constrictions. When constrictions are clustered together, central-airways deposition can increase significantly in regions away from constrictions, but distal-airways deposition in those regions remains largely unchanged. We use our model to calculate lung clearance index (LCI), a clinical measure of ventilation heterogeneity, after applying constrictions of varying severities in one lobe. We find an increase in LCI coinciding with significantly reduced deposition in the affected lobe. Our results show how the model provides a framework for development of computational tools that capture the impacts of airway disease, which could significantly affect predictions of regional dosing.

To assess how the presence of surfactant in lung airways alters the flow of mucus that leads to plug formation and airway closure, we investigate the effect of insoluble surfactant on the instability of a viscoplastic liquid coating the interior of a cylindrical tube. Evolution equations for the layer thickness using thin-film and long-wave approximations are derived that incorporate yield-stress effects and capillary and Marangoni forces. Using numerical simulations and asymptotic analysis of the thin-film system, we quantify how the presence of surfactant slows growth of the Rayleigh-Plateau instability, increases the size of initial perturbation required to trigger instability and decreases the final peak height of the layer. When the surfactant strength is large, the thin-film dynamics coincide with the dynamics of a surfactant-free layer but with time slowed by a factor of four and the capillary Bingham number, a parameter proportional to the yield stress, exactly doubled. By solving the long-wave equations numerically, we quantify how increasing surfactant strength can increase the critical layer thickness for plug formation to occur and delay plugging. The previously established effect of the yield stress in suppressing plug formation [Shemilt et al., J. Fluid Mech., 2022, vol. 944, A22] is shown to be amplified by introducing surfactant. We discuss the implications of these results for understanding the impact of surfactant deficiency and increased mucus yield stress in obstructive lung diseases.

One mechanism for airway closure in the lung is the surface-tension-driven instability of the mucus layer which lines the airway wall. We study the instability of an axisymmetric layer of viscoplastic Bingham liquid coating the interior of a rigid tube, which is a simple model for an airway that takes into account the yield stress of mucus. An evolution equation for the thickness of the liquid layer is derived using long-wave theory, from which we also derive a simpler thin-film evolution equation. In the thin-film case we show that two branches of marginally yielded static solutions of the evolution equation can be used to both predict the size of the initial perturbation required to trigger instability and quantify how increasing the capillary Bingham number (a parameter measuring yield stress relative to surface tension) reduces the final deformation of the layer. Using numerical solutions of the long-wave evolution equation, we quantify how the critical layer thickness required to form a liquid plug in the tube increases as the capillary Bingham number is increased. We discuss the significance of these findings for modelling airway closure in obstructive conditions such as cystic fibrosis, where the mucus layer is often thicker and has a higher yield stress.

We study the dynamics of two air bubbles driven by the motion of a suspending viscous fluid in a Hele-Shaw channel with a small elevation along its centreline via physical experiment and numerical simulation of a depth-averaged model. For a single-bubble system we establish that, in general, the bubble propagation speed monotonically increases with bubble volume so that two bubbles of different sizes, in the absence of any hydrodynamic interactions, will either coalesce or separate in a finite time. However, our experiments indicate that the bubbles interact and that an unstable two-bubble state is responsible for the eventual dynamical outcome: coalescence or separation. These results motivate us to develop an edge-tracking routine and to calculate these weakly unstable two-bubble steady states from the governing equations. The steady states consist of pairs of ‘aligned’ bubbles that appear on the same side of the centreline with the larger bubble leading. We also discover, through time-dependent simulations and physical experiment, another class of two-bubble states that, surprisingly, are stable. In contrast to the ‘aligned’ steady states, these bubbles appear on either side of the centreline and are ‘offset’ from each other. We calculate the bifurcation structures of both classes of steady states as the flow rate and bubble volume ratio are varied. We find that they exhibit intriguing similarities to the single-bubble bifurcation structure, which has implications for the existence of $n$ -bubble steady states.

The steady propagation of air bubbles through a Hele-Shaw channel with either a rectangular or partially occluded cross section is known to exhibit solution multiplicity for steadily propagating bubbles, along with complicated transient behaviour where the bubble may visit several edge states or even change topology several times, before typically reaching its final propagation mode. Many of these phenomena can be observed both in experimental realisations and in numerical simulations based on simple Darcy models of flow and bubble propagation in a Hele-Shaw cell. In this paper, we investigate the corresponding problem for the propagation of a viscous drop (with viscosity $$\nu $$ relative to the surrounding fluid) using a Darcy model. We explore the effect of drop viscosity on the steady solution structure for drops in rectangular channels or with imposed height variations. Under the Darcy model in a uniform channel, steady solutions for bubbles map directly on to those for drops with any internal viscosity $$\nu \ne 1$$ . Hence, the solution multiplicity predicted for bubbles also occurs for drops, although for $$\nu >1$$ , the interface shape is reversed with inflection points appearing at the rear rather than the front of the drop. The equivalence between bubbles and drops breaks down for transient behaviour, at the introduction of any height variation, for multiple bodies of different viscosity ratios and for more detailed models which produce a more complicated flow in the interior of the drop. We show that the introduction of topography variations affects bubbles and drops differently, with very viscous drops preferentially moving towards more constricted regions of the channel. Both bubbles and drops can undergo transient behaviour which involves breakup into two almost equal bodies, which then symmetry break before either recombining or separating indefinitely.

Motivated by the desire to understand complex transient behaviour in fluid flows, we study the dynamics of an air bubble driven by the steady motion of a suspending viscous fluid within a Hele-Shaw channel with a centred depth perturbation. Using both experiments and numerical simulations of a depth-averaged model, we investigate the evolution of an initially centred bubble of prescribed volume as a function of flow rate and initial shape. The experiments exhibit a rich variety of organised transient dynamics, involving bubble break up as well as aggregation and coalescence of interacting neighbouring bubbles. The long-term outcome is either a single bubble or multiple separating bubbles, positioned along the channel in order of increasing velocity. Up to moderate flow rates, the life and fate of the bubble are reproducible and can be categorised by a small number of characteristic behaviours that occur in simply-connected regions of the parameter plane. Increasing the flow rate leads to less reproducible time evolutions with increasing sensitivity to initial conditions and perturbations in the channel. Time-dependent numerical simulations that allow for break up and coalescence are found to reproduce most of the dynamical behaviour observed experimentally including enhanced sensitivity at high flow rate. An unusual feature of this system is that the set of steady and periodic solutions can change during temporal evolution because both the number of bubbles and their size distribution evolve due to break up and coalescence events. Calculation of stable and unstable solutions in the single- and two-bubble cases reveals that the transient dynamics are orchestrated by weakly-unstable solutions of the system termed edge states that can appear and disappear as the number of bubbles changes.

Place a droplet of mineral oil on water and the oil will spread to cover the water surface in a thin film. In this paper we study the everted problem: an aqueous droplet deposited onto a deep layer of silicone oil. As it is energetically favourable for the oil phase to spread to cover the droplet surface completely, the droplet is ultimately engulfed in the oil layer. We present a detailed study of engulfment dynamics, from the instant the droplet impacts the oil surface until it finally sediments into the less dense oil. We study a broad range of droplet sizes (micrometric to millimetric) and oil kinematic viscosities (10(2)-10(5) cSt), corresponding to a viscosity-dominated parameter regime. We find that droplet engulfment involves the rapid submersion of the droplet driven by capillary forces in the oil surface, followed by the much slower peeling of the droplet from the interface, to which it is adhered via a thin cloaking layer of oil formed during the earlier stage. The later peeling stage is driven by a combination of geometric constraints at the apparent contact line and gravity pulling on the droplet. Gravitational effects are therefore essential to complete engulfment, even for micrometric droplets. Furthermore, the opposing effects of geometry and gravity result in the longest engulfment times for droplets of intermediate size. Experiments at fixed droplet size reveal a power law dependence of engulfment time on oil kinematic viscosity, which we argue reflects the dynamical formation of the oil cloaking layer.

Accurate low-dimensional models for the dynamics of falling liquid films subject to localized or time-varying heating are essential for applications that involve patterning or control. However, existing modelling methodologies either fail to respect fundamental thermodynamic properties or else do not accurately capture the effects of advection and diffusion on the temperature profile. We argue that the best-performing long-wave models are those that give the surface temperature implicitly as the solution of an evolution equation in which the wall temperature alone (and none of its derivatives) appears as a source term. We show that, for both flat and non-uniform films, such a model can be rationally derived by expanding the temperature field about its free-surface values. We test this model in linear and nonlinear regimes, and show that its predictions are in remarkable quantitative agreement with full Navier–Stokes calculations regarding the surface temperature, the internal temperature field and the surface displacement that would result from temperature-induced Marangoni stresses.

The fabrication of a high-quality POLED (Polymeric Organic Light Emitting Diode) display requires the deposition of identical, uniform fluid films into a large number of shallow recessed regions that form a regular array of pixels on a display backplane. We determine the protocols required to achieve continuous liquid coverage of the entire pixel area for the case where equally-spaced fluid droplets are sequentially printed along a straight line within a stadium-shaped pixel, and explore how these protocols depend on the wetting properties of the pixel surface. Our investigation uses a combination of experiments and numerical modelling, based on the assumption of fluid redistribution via capillary spreading according to a Cox-Voinov law. We show that the model is able to predict quantitatively the evolution of the liquid deposited in a pixel and provides a computationally inexpensive design tool to determine efficient printing strategies that account for uncertainties arising from imperfect substrate preparation or printhead dysfunction.

We hypothesize that dynamical systems concepts used to study the transition to turbulence in shear flows are applicable to other transition phenomena in fluid mechanics. In this paper, we consider a finite air bubble that propagates within a Hele-Shaw channel containing a depth-perturbation. Recent experiments revealed that the bubble shape becomes more complex, quantified by an increasing number of transient bubble tips, with increasing flow rate. Eventually, the bubble changes topology, breaking into multiple distinct entities with non-trivial dynamics. We demonstrate that qualitatively similar behaviour to the experiments is exhibited by a previously established, depth-averaged mathematical model and arises from the model's intricate solution structure. For the bubble volumes studied, a stable asymmetric bubble exists for all flow rates of interest, while a second stable solution branch develops above a critical flow rate and transitions between symmetric and asymmetric shapes. The region of bistability is bounded by two Hopf bifurcations on the second branch. By developing a method for a numerical weakly nonlinear stability analysis we show that unstable periodic orbits (UPOs) emanate from the first Hopf bifurcation. Moreover, as has been found in shear flows, the UPOs are edge states that influence the transient behaviour of the system.

We study the propagation of finite bubbles in a Hele-Shaw channel, where a centred occlusion (termed a rail) is introduced to provide a small axially uniform depth constriction. For bubbles wide enough to span the channel, the system’s behaviour is similar to that of semi-infinite fingers and a symmetric static solution is stable. Here, we focus on smaller bubbles, in which case the symmetric static solution is unstable and the static bubble is displaced towards one of the deeper regions of the channel on either side of the rail. Using a combination of experiments and numerical simulations of a depth-averaged model, we show that a bubble propagating axially due to a small imposed flow rate can be stabilised in a steady symmetric mode centred on the rail through a subtle interaction between stabilising viscous forces and destabilising surface tension forces. However, for sufficiently large capillary numbers Ca, the ratio of viscous to surface tension forces, viscous forces in turn become destabilising thus returning the bubble to an off-centred propagation regime. With decreasing bubble size, the range of Ca for which steady centred propagation is stable decreases, and eventually vanishes through the coalescence of two supercritical pitchfork bifurcations. The depth-averaged model is found to accurately predict all the steady modes of propagation observed experimentally, and provides a comprehensive picture of the underlying steady bifurcation structure. However, for sufficiently large imposed flow rates, we find that initially centred bubbles do not converge onto a steady mode of propagation. Instead they transiently explore weakly unstable steady modes, an evolution which results in their break-up and eventual settling into a steady propagating state of changed topology.

We use a combination of experiments and numerical modelling to investigate the influence of physico-chemical-patterned substrates on the spreading of fluid deposited as a partially overlapping sequence of droplets via inkjet printing. Our investigation is motivated by the manufacture of polymeric organic light-emitting-diode displays, where the substrate is textured with a regular array of shallow recessed regions (pixels) that are highly wetting compared to the remainder of the substrate. We examine the roles of topography and wettability patterning separately and in combination. On a substrate with uniform wettability, we find that the presence of bounding side walls enhances the local spreading and facilitates fluid coverage of the entire recessed region, but containment within the pixel is not guaranteed. In contrast, wettability patterning alone leads to robust containment of the fluid within the wetting region, but fluid coverage is reduced in the absence of side walls. Our theoretical calculations use a simplified numerical model of fluid redistribution via purely capillary effects, augmented by a Cox-Voinov spreading law. The neglect of fluid viscosity in this model means that, after an initial period of agreement, the predicted evolution is faster than in the experiments. Nonetheless, the simplified model achieves excellent predictions both for the liquid morphologies and for the conditions required for successful pixel filling on substrates with topographical and wettability variations.

We present a novel experimental system that can be used to study the dynamics of picoliter droplet spreading over substrates with topographic variations. We concentrate on spreading of a droplet within a recessed stadium-shaped pixel, with applications to the manufacture of POLED displays, and find that the sloping side wall of the pixel can either locally enhance or hinder spreading depending on whether the topography gradient ahead of the contact line is positive or negative. Locally enhanced spreading occurs via the formation of thin pointed rivulets along the side walls of the pixel through a mechanism similar to capillary rise in sharp corners. We demonstrate that a thin-film model combined with an experimentally measured spreading law, which relates the speed of the contact line to the contact angle, provides excellent predictions of the evolving liquid morphologies. We also show that the spreading can be adequately described by a Cox-Voinov law for the majority of the evolution. The model does not include viscous effects and hence, the timescales for the propagation of the thin pointed rivulets are not captured. Nonetheless, this simple model can be used very effectively to predict the areas covered by the liquid and may serve as a useful design tool for systems that require precise control of liquid on substrates.