We perform simulations of an impulsively started, axisymmetric viscoelastic jet exiting a nozzle and entering a stagnant gas phase using the open-source code Basilisk. This code allows for efficient computations through an adaptively refined volume-of-fluid technique that can accurately capture the deformation of the liquid–gas interface. We use the FENE-P constitutive equation to describe the viscoelasticity of the liquid, and employ the log-conformation transformation, which provides stable solutions for the evolution of the conformation tensor as the jet thins down under the action of interfacial tension. For the first time, the entire jetting and breakup process of a viscoelastic fluid is simulated, including the pre-shearing flow through the nozzle, which results in an inhomogeneous initial radial stress distribution in the fluid thread that affects the subsequent breakup dynamics. The evolution of the velocity field and the elastic stresses in the nozzle are validated against analytical solutions where possible, and the early-stage dynamics of the jet evolution are compared favourably to the predictions of linear stability theory. We study the effect of the flow inside the nozzle on the thinning dynamics of the viscoelastic jet (which develops distinctive ‘beads-on-a-string’ structures) and on the spatio-temporal evolution of the polymeric stresses in order to systematically explore the dependence of the filament thinning and breakup characteristics on the initial axial momentum of the jet and the extensibility of the dissolved polymer chains.

We numerically study the influence of a soluble surfactant on the microjetting mode of the liquid–liquid flow focusing configuration. The surfactant adsorbs on the interface next to the feeding capillary and accumulates in front of the emitted jet, significantly lowering the surface tension there. The resulting Marangoni stress substantially alters the balance of the tangential stresses at the interface but does not modify the interface velocity. The global stability analysis at the minimum flow rate stability limit shows that the Marangoni stress collaborates with soluto-capillarity to stabilize the microjetting mode. Our analysis unveils the noticeable effect of the Marangoni stress associated with the surface tension perturbation. Surfactant diffusion and desorption hardly affect the stability limit. Transient numerical simulations show how subcritical and supercritical base flows respond to a spatially localized initial perturbation. Our parametric study indicates that the minimum flow rate ratio depends on the adsorption constant and the surfactant concentration through the product of these two variables. The surfactant stabilizing effect increases with the outer stream flow rate. We show that surfactants not only stabilize the microemulsion resulting from the jet breakup in hydrodynamic focusing, but also allow for the reduction of droplet size. Our findings advance the fundamental understanding of the complex role of surfactants in tip streaming via hydrodynamic focusing. In particular, our results contradict the common assumption that adding surfactant favours tip streaming simply because it reduces the meniscus tip surface tension.

When subjected to external stress exceeding their yield stress ($\sigma _{Y}$), yield stress fluids (YSFs) undergo transition from solid to liquid. While the flow dynamics of capillary flows – including wetting, coating, spreading, and wicking – can be predicted accurately for fluids without yield stress, only a few asymptotic solutions have been validated against experiments for YSFs. In this study, we explore experimentally and theoretically the entire dynamics and solidifying process during capillary imbibition of YSFs, representing the fundamental aspect of capillarity-driven flow. We observe that during the gradual deceleration of capillary rise, YSFs exhibit apparent flow below $\sigma _{Y}$. This unexpected behaviour is linked to wall slip, manifested as plug flow in particle image velocimetry analyses. We introduce a model with numerical solutions, grounded in the rheological properties of YSFs and slip mechanisms, to accurately describe this liquid–solid transitional behaviour.

We explore the drawing of a shear-thinning or shear-thickening thread with an axisymmetric hole that evolves due to axial drawing, inertia and surface tension effects. The stress is assumed to be proportional to the shear rate raised to the $n$th power. The presence of non-Newtonian rheology and surface tension forces acting on the hole introduces radial pressure gradients that make the derivation of long-wavelength equations significantly more challenging than either a Newtonian thread with a hole or shear-thinning and shear-thickening threads without a hole. In the case of weak surface tension, we determine the steady-state profiles. Our results show that for negligible inertia the hole size at the exit becomes smaller as $n$ is decreased (i.e. strong shear-thinning effects) above a critical draw ratio, but surprisingly the opposite is true below this critical draw ratio. We determine an accurate estimate of the critical draw ratio and also discuss how inertia affects this process. We further show that the dynamics of hole closure is dominated by a different limit, and we determine the asymptotic forms of the hole closure process for shear-thinning and shear-thickening fluids with inertia. A linear instability analysis is conducted to predict the onset of draw resonance. We show that increased shear thinning, surface tension and inlet hole size all act to destabilise the flow. We also show that increasing shear-thinning effects reduce the critical Reynolds number required for unconditional stability. Our study provides valuable insights into the drawing process and its dependence on the physical effects.

Bistable states for a sufficiently large amount of liquid can appear in an eccentric capillary due to the eccentricity effect under zero gravity (J. Fluid Mech, vol. 863, 2019, pp. 364–385). A transverse body force, which can lead to rich physical phenomena of a droplet, may lead to multistable states (bistability, tristability and the likes) of a sufficiently large amount of liquid in a capillary. We theoretically investigate this situation in a circular or annular capillary tube under a transverse body force. The results show that there can be tristable (bistable) states in an annular (circular) capillary tube: an occluding configuration and two (one) non-occluding configurations. In the annular tube, for one of the non-occluding configurations, the gas–liquid interface in the middle cross-section of the droplet meets both the inner and outer walls of the tube (bridging configuration); for the other non-occluding configuration, the gas–liquid interface in the middle cross-section of the droplet does not meet the inner wall (non-bridging configuration). The multistability is dependent on the Bond number, the contact angle and the cross-sectional shape. The multistability cannot occur for a zero or very large Bond number. A hydrophilic condition (the contact angle smaller than 90°) contributes to the non-occluding non-bridging configuration, while the hydrophobic condition (the contact angle larger than 90°) contributes to the non-occluding bridging configuration (only for the annular capillary). For the annular capillary with a not-so-large contact angle, increasing the inner-to-outer radius ratio can lead to a larger range of Bond numbers, in which the multistability occurs.

An interesting resurrection phenomenon (including the initial complete submersion, subsequent resurfacing and final rebounding) of a superhydrophobic sphere impacting onto a liquid bath was observed in experiments and direct numerical simulations by Galeano-Rios et al. (J. Fluid Mech., vol. 912, 2021, A17). We investigate the mechanisms of the liquid entry for a superhydrophobic cylinder in this paper. The superhydrophobic cylinder, commonly employed as supporting legs for insects and robots at the liquid surface, can exhibit liquid-entry mechanisms different from those observed with the sphere. The direct numerical simulation method is applied to the impact of a two-dimensional (2-D) superhydrophobic cylinder (modelled as a pseudo-solid) onto a liquid bath. We find that for the impacting cylinder the resurrection phenomenon can also exist, and the cylinder can either rebound (get detached from the liquid surface) or stay afloat after resurfacing. The cylinder impact behaviour is classified into four regimes, i.e. floating, bouncing, resurrecting (resurrecting-floating and resurrecting-bouncing) and sinking, dependent on the Weber number and the density ratio of the cylinder to the liquid. For the regimes of floating and bouncing, the force analysis indicates that the form drag dominates the motion of the cylinder in the very beginning of the impact, while subsequently the surface tension force also plays a role with the contact line pinning on the horizontal midline of the cylinder. For the critical states of the highlighted resurrecting regime, our numerical results show that the rising height for the completely submerged cylinder of different density ratios remains nearly unchanged. Accordingly, a relation between the maximum ascending velocity and the density ratio is derived to predict whether the completely submerged cylinder can resurface.

Diffuse interface models are an important class of models used to describe multi-phase flows. In the case of incompressible viscous fluids there are a number of different diffuse interface models which have been known for many years. Nevertheless, a model based on mixture theory with a full set of conservation laws for the conservation of linear momentum of each constituent was not yet available. This gap was filled by ten Eikelder et al. (J. Fluid Mech., in press) recently and a first comparison with known models of Navier–Stokes/Cahn–Hilliard type is given. A detailed understanding of the relations between these models remains an important question.

Thin-film equations are utilised in many different areas of fluid dynamics when there exists a direction in which the aspect ratio can be considered small. We consider thin free films with Marangoni effects in the extensional flow regime, where velocity gradients occur predominantly along the film. In practice, because of the local deposition of surfactants or input of energy, asymmetric distributions of surfactants or surface tension more generally, are possible. Such examples include the surface of bubbles and the rupture of thin films. In this study, we consider the asymmetric thin-film equations for extensional flow with Marangoni effects. Concentrating on the case of small Reynolds number $ Re $, we study the deposition of insoluble surfactants on one side of a liquid sheet otherwise at rest and the resulting thinning and rupture of the sheet. The analogous problem with a uniformly thinning liquid sheet is also considered. In addition, the centreline deformation is discussed. In particular, we show analytically that if the surface tension isotherm $\sigma = \sigma (\varGamma )$ is nonlinear (surface tension $\sigma$ varies with surfactant concentration $\varGamma$), then accounting for top–bottom asymmetry leads to slower (faster) thinning and pinching if $\sigma = \sigma (\varGamma )$ is convex (concave). The analytical progress reported in this paper allows us to discuss the production of satellite drops from rupture via Marangoni effects, which, if relevant to surface bubbles, would be an aerosol production mechanism that is distinct from jet drops and film drops.

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.

We explore predictions of two models of one-dimensional capillary rise in rigid and partially saturated porous media. One is an existing one from the literature and the second is a free-boundary model based on Richards’ equation with two moving boundaries of the evolving partially saturated region. Both models involve the specification of saturation-dependent functions for local capillary pressure and permeability and connect to classical models for saturated porous media. Existing capillary-rise experiments show two notable regimes: (i) an early-time regime typically well-described by classical capillary-rise theory in a fully saturated porous media, and (ii) a long-time regime that has anomalous dynamics in which the capillary-rise height may scale with a non-classical power law in time or have more complicated dynamics. We demonstrate that the predictions of both models compare well with experimental capillary-rise data over early- and long-time regimes gathered from three independent studies in the literature. The model predictions also shed light on recent scaling laws that relate the capillary pressure and permeability of the partially saturated media to the capillary-rise height. We use these models to probe computationally observed permeability relationships to capillary-rise height. We demonstrate that a recently proposed permeability scaling for the anomalous capillary-rise regime is indeed realized and is particularly apparent in the lower portion of the partially saturated media. For our free-boundary model we also compute capillary pressure measures and show that these reveal the linear relation between the capillary pressure and capillary-rise height expected for a capillarity–gravity balance in the upper portion of the partially saturated porous media.

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.

Studying liquid jet impacts on a liquid pool is crucial for various engineering and environmental applications. During jet impact, the free surface of the pool deforms and a cavity is generated. Simultaneously, the free surface of the cavity extends radially outward and forms a rim. Eventually the cavity collapses by means of gas inertia and surface tension. Our numerical investigation using an axisymmetric model in Basilisk C explores cavity collapse dynamics under different impact velocities and gas densities. We validate our model against theory and experiments across a previously unexplored parameter range. Our results show two distinct regimes in the cavity collapse mechanism. By considering forces pulling along the interface, we derive scaling arguments for the time of closure and maximum radius of the cavity, based on the Weber number. For jets with uniform constant velocity from tip to tail and $We \leqslant 150$, the cavity closure is capillary-dominated and happens below the surface (deep seal). In contrast, for $We \geqslant 180$ the cavity closure happens above the surface (surface seal) and is dominated by the gas entrainment and the pressure gradient that it causes. Additionally, we monitor gas velocity and pressure throughout the impact process. This analysis reveals three critical moments of maximum gas velocity: before impact, at the instant of cavity collapse and during droplet ejection following cavity collapse. Our results provide information for understanding pollutant transport during droplet impacts on large bodies of water, and other engineering applications, like additive manufacturing, lithography and needle-free injections.

When two fluid drops touch, they coalesce due to surface tension. At early times, there is only a relatively small fluid bridge joining the drops. An asymptotic solution is presented for an inertial regime of early-time coalescence, in which inertial forces balance surface tension at leading order. It is demonstrated that viscosity nevertheless has a leading-order effect. Radial momentum is created at the tightly curved edge of the fluid bridge by the net force $2\gamma$ (per unit length) due to surface tension. This momentum is left behind the radially expanding bridge edge in a thin viscous wake. The divergent volume flux in the wake entrains fluid from above and below the bridge, and drives an inviscid irrotational flow in the drops on the scale of the bridge radius. This flow widens the gap between the drops ahead of the bridge, and the larger gap width results in a lower rate of coalescence. Including viscosity in this way improves the agreement between theory and the available experimental and numerical data.

We study numerically the microjetting mode obtained when a fluid is injected through a tube submerged in a uniaxial extensional flow. The steady solution to the full nonlinear Navier–Stokes equations is calculated. We obtain the linear global modes determining the linear stability of the steady solution. For sufficiently large outer viscosity, the flow remains stable for infinitely small values of the injected flow rate. This implies that jets with vanishing diameters can be produced regardless of the jet viscosity and outer flow strength. For a sufficiently small inner-to-outer viscosity ratio, the microjetting instability is associated only with the flow near the entrance of the jet. The tapering meniscus stretches and adopts a slender quasiconical shape. Consequently, the cone tip is exposed to an intense outer flow, which stabilizes the flow in the cone–jet transition region. This work presents the first evidence that fluid jets with arbitrarily small diameters can be stably produced via tip streaming. The results are related to those of a droplet in a uniaxial extensional flow with its equator pinned to an infinitely thin ring. The pinning of the equator drastically affects the droplet stability and breakup.

To investigate the influence of inertia and slip on the instability of a liquid film on a fibre, a theoretical framework based on the axisymmetric Navier–Stokes equations is proposed via linear instability analysis. The model reveals that slip significantly enhances perturbation growth in viscous film flows, whereas it exerts minimal influence on flows dominated by inertia. Moreover, under no-slip boundary conditions, the dominant instability mode of thin films remains unaltered by inertia, closely aligning with predictions from a no-slip lubrication model. Conversely, when slip is introduced, the dominant wavenumber experiences a noticeable reduction as inertia decreases. This trend is captured by an introduced lubrication model with giant slip. Direct numerical simulations of the Navier–Stokes equations are then performed to further confirm the theoretical findings at the linear stage. For the nonlinear dynamics, no-slip simulations show complex vortical structures within films, driven by fluid inertia near surfaces. Additionally, in scenarios with weak inertia, a reduction in the volume of satellite droplets is observed due to slip, following a power-law relationship.

The rupture of the thin film at the top of a bubble at a liquid–gas interface leads to an axisymmetric collapse of the bubble cavity. We present scaling laws for such a cavity collapse, established from experiments conducted with bubbles spanning a wide range of Bond (${10^{-3}< Bo\leq 1}$) and Ohnesorge numbers (${10^{-3}< Oh<10^{-1}}$), defined with the bubble radius $R$. The cavity collapse is a capillary-driven process, with a dependency on viscosity and gravity, affecting respectively, precursory capillary waves on the cavity boundary and the static bubble shape. The collapse is characterised by the normal interface velocity ($U_n$) and by the tangential wave propagation velocity of the kink ($U_t$), defined by the intersection of the concave cavity boundary formed after the rupture of the thin film with the convex boundary of the bubble cavity. During the collapse, $U_t$ remains constant and is shown to be $U_t=4.5U_c{\mathcal {W}}_R$, where $U_c$ is the capillary velocity and ${\mathcal {W}}_R(Oh,Bo)={(1-\sqrt {Oh {\mathscr {L}}} )^{-1/2}}$ is the wave resistance factor due to the precursory capillary waves, with $\mathscr {L}(Bo)$ being the path correction of the kink motion. The movement of the kink in the normal direction is part of the inward shrinkage of the whole cavity due to the sudden reduction of gas pressure inside the bubble cavity after the thin film rupture. This normal velocity is shown to scale as $U_c$ in the equatorial plane, while at the bottom of the cavity $\bar {U}_{nb}=U_c(Z_c/R)({\mathcal {W}_R}/ {\mathscr {L}})$, where $Z_c(Bo)$ is the static cavity depth. The filling rate of the cavity, which remains a constant throughout the collapse, is shown to be entirely determined by the shrinking velocity and scales as ${Q_T\simeq 2{\rm \pi} R Z_c U_c}$. From $Q_T$ we recover the jet velocity scaling, thereby relating the cavity collapse with the jet velocity scaling.

Keeping a tube from being plugged by a fluid is an important process in applications. An interesting re-entrant phenomenon for the capillary state with the occluding state sandwiching the non-occluding state from both the high- and low-Bond-number regions can appear by inserting a rod into a horizontal tube at an eccentric position (Tan et al., J. Fluid Mech, vol. 946, 2022, A7). Containers with rounded corners are very common. We theoretically investigate a situation for a horizontal open tube with rounded corner(s). The results show that a re-entrant non-occlusion at a contact angle can also appear without the insertion of any object. The competition between the rounded corner wetting/non-wetting effect and gravity effect can lead to a re-entrant non-occlusion. The re-entrant non-occlusion is affected by the shape and orientation of the rounded corner(s). For a tube with only one rounded corner, the re-entrant non-occlusion exists when the rounded corner has a not-so-large corner radius and is not in a landscape orientation. For a tube with two (or more) rounded corners, the corner(s) with the strongest corner effect will determine the existence or non-existence of the re-entrant non-occlusion. This paper provides an effective scheme for designing a high-performance capillary with corners that are not easily occluded by a fluid and removing fluid blockage from a capillary in optofluidic/microfluidic applications.

We explore ${\rm CO}_2$ injection into a layered permeable rock consisting of high permeability reservoir layers, separated by low permeability mudstone, and taking the shape of an anticline within a laterally extensive aquifer. We first show how the storage capacity of the formation depends on the capillary entry pressure of the inter-layer mudstone, so that ${\rm CO}_2$ cannot flow from one layer into the next. We then consider a formation composed of two layers, overlain by a cap rock. For injection into the lowest layer, we show that the injection rate, capillary entry pressure and buoyancy driven flux through the mudstone determine whether the lower or upper layer fills to the spill point first. We also show that at the end of the injection phase, ${\rm CO}_2$ may continue to flow from the lower to the upper layer. This implies that injection should be stopped once the injected volume matches the static capacity of the formation in order to prevent spilling after injection. We present a series of analogue experiments of a two layered system that illustrate some of the principles described by the model, and assess the implications of the results for field scale systems.

We study the capillary rise of viscous liquids into sharp corners formed by two surfaces whose geometry is described by power laws $h_i(x) = c_i x^n$, $i = 1,2$, where $c_2 > c_1$ for $n \geq 1$. Prior investigations of capillary rise in sharp corners have shown that the meniscus altitude increases with time as $t^{1/3}$, a result that is universal, i.e. applies to all corner geometries. The universality of the phenomenon of capillary rise in sharp corners is revisited in this work through the analysis of a partial differential equation for the evolution of a liquid column rising into power-law-shaped corners, which is derived using lubrication theory. Despite the lack of geometric similarity of the liquid column cross-section for $n>1$, there exist a scaling and a similarity transformation that are independent of $c_i$ and $n$, which gives rise to the universal $t^{1/3}$ power law for capillary rise. However, the prefactor, which corresponds to the tip altitude of the self-similar solution, is a function of $n$, and it is shown to be bounded and monotonically decreasing as $n\to \infty$. Accordingly, the profile of the interface radius as a function of altitude is also independent of $c_i$ and exhibits slight variations with $n$. Theoretical results are compared against experimental measurements of the time evolution of the tip altitude and of profiles of the interface radius as a function of altitude.

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.