We consider shear-driven longitudinal flow of an exterior fluid over a periodic array of rectangular grooves filled with an immiscible interior fluid (the "lubricant"), the grooves being formed by infinitely thin ridges protruding from a flat substrate. The ratio lambda of the effective slip length to the semiperiod is a function of the ratio mu of the interior to exterior viscosities and the ratio h of the groove depth to the semiperiod. We focus on the limit mu << 1, which is singular for that geometry. We find that the viscous resistance to the imposed shear is dominated by a boundary layer of exponentially small extent about the ridge tips, resulting in the effective slip length scaling as mu(-1/2)-not mu(-1) as implied by intuitive arguments overlooking the tip contributions (and by proposed approximations in the literature). Analyzing that exponential region in conjunction with an integral force balance, we find the simple asymptotic approximation lambda approximate to mu(-1/2); using conformal mappings, we also calculate the leading-order correction to that result, which introduces a dependence upon h. The ensuing asymptotic expansion breaks down for h=O(mu(1/2)), upon transitioning to a lubrication geometry. We accordingly conduct a companion asymptotic analysis in the distinguished limit of small mu and fixed H=h/mu(1/2), which gives lambda approximate to mu H-1/2/(1+H) as well as a closed-form leading-order correction to that approximation; the intuitive mu(-1) scaling is accordingly only relevant to the regime H << 1 corresponding to extremely shallow grooves. We demonstrate excellent agreement between our predictions and numerical solutions constructed using a boundary-integral formulation.

The charge on a conducting cylinder becomes uniformly distributed at large aspect ratios, albeit very slowly and with nonuniformity persisting near the ends. We show that the charge density and field magnitude are enhanced near the ends by a factor scaling as the square-root of the logarithm of the aspect ratio. This scaling is obtained by locally resumming a perturbation-series solution to the integral equation of slender-body theory. The same scaling applies to a broad class of "truncated" slender bodies – including cylinders and shapes tapered near the ends on the cross-sectional scale – as well as other physical setups. We validate this scaling through boundary-element simulations for cylinders with flat and hemispherical caps, and demonstrate its applicability to diffusion, plasmonics and Stokes flow.

Particles in a sufficiently strong electric field spontaneously rotate, provided that charge relaxation is slower in the particle than in the suspending fluid. It has long been known that drops also exhibit such "Quincke rotation," with the electrohydrodynamic flow induced by electrical shear stresses at the interface leading to an increased critical field. However, the hysteretic onset of this instability, observed for sufficiently low-viscosity drops, has so far eluded theoretical understanding-including simulations that have struggled in this regime owing to charge-density steepening driven by surface convection. Here, we conduct a numerical study of the leaky-dielectric model in a simplified two-dimensional setting involving a circular drop, considering arbitrary viscosity ratios and field strengths. As the viscosity of the drop is decreased relative to the suspending fluid, the pitchfork bifurcation marking the onset of drop rotation is found to transition from supercritical to subcritical, giving rise to a field-strength interval of bistability. In this subcritical regime, the critical field is always large enough that, at the bifurcation, the symmetric base-state solution exhibits equatorial charge-density blowup singularities of the type recently described by Peng et al. [Phys. Rev. Fluids 9, 083701 (2024)]. As the rotation speed increases along the initially unstable solution branch from the bifurcation, the singularities gradually shift from the equator and ultimately disperse once the rotational component of the flow is strong enough to eliminate the surface stagnation points.

We investigate flow-induced choking in soft Hele-Shaw cells comprising a fluid-filled gap in between a rigid plate and a confined block of elastomer. Fluid injected from the centre of the circular rigid plate flows radially outwards, causing the elastomeric block to deform, before exiting through the cell rim. The pressure in the fluid deforms the elastomer, increasing the size of the gap near the inlet, and decreasing the gap near the cell rim, because of volume conservation of the solid. At a critical injection flow rate, the magnitude of the deformation becomes large enough that the flow is occluded entirely at the rim. Here, we explore the influence of elastomer geometry on flow-induced choking and, in particular, the case of a thick block with radius smaller than its depth. We show that choking can still occur with small-aspect-ratio elastomers, even though the confining influence of the back wall that bounds the elastomer becomes negligible; in this case, the deformation length scale is set by the radial size of the cell rather than the depth of the block. Additionally, we reveal a distinction between flow-induced choking in flow-rate-controlled flows and flow-rate-limiting behaviour in pressure-controlled flows.

We revisit the model problem of Squires & Brady ( Phys. Fluids , vol. 17, 2005, 073101), where a Brownian probe is dragged through a dilute dispersion of Brownian bath particles. In this problem, the microrheology due to excluded-volume interactions is represented by an effective viscosity, with the nonlinearity in the driving force entering via the dependence of the viscosity increment (relative to the viscosity of a pure solvent) upon the deformation of the dispersion microstructure. Our interest is in the limit of large Péclet numbers, $ P{\kern-1pt}e\gg 1$ , where the microstructural deformation adopts the form of a boundary layer about the upstream hemisphere of the probe. We show that the boundary-layer solution breaks down at the equator of the probe and identify a transition region about the equator, connecting the layer to a downstream wake. The microstructural deformation in this region is governed by a universal boundary-value problem in a semi-bounded two-dimensional domain. The equatorial region continues downstream as a transition layer, which separates the wake of the probe from the undisturbed ambient; in that layer, the microstructure is governed by a one-dimensional heat-like equation. Accounting for the combined contributions from the respective asymptotic provinces we find the approximation $ ({1}/{2})[1+ (\ln P{\kern-1pt}e + 1.046)/ P{\kern-1pt}e]$ for the ratio of the large- $ P{\kern-1pt}e$ viscosity increment to the corresponding linear-response increment. Our asymptotic approximation is in excellent agreement with the increment predicted by a finite-difference numerical calculation of the microstructure deformation, tailored to the large- $ P{\kern-1pt}e$ topology.

We illuminate effects of surface-charge convection intrinsic to leaky-dielectric electrohydrodynamics by analyzing the symmetric steady state of a circular drop in an external field at arbitrary electric Reynolds number ReE. In formulating the problem, we identify an exact factorization that reduces the number of dimensionless parameters from four—ReE and the conductivity, permittivity and viscosity ratios—to two: a modified electric Reynolds number Rẽ and a charging parameter ϖ. In the case ϖ<0, where charge relaxation in the drop phase is slower than in the suspending phase, and, as a consequence, the interface polarizes antiparallel to the external field, we find that above a critical Rẽ value the solution exhibits a blowup singularity such that the surface-charge density diverges antisymmetrically with the −1/3 power of distance from the equator. We use local analysis to uncover the structure of that blowup singularity, wherein surface charges are convected by a locally induced flow towards the equator where they annihilate. To study the blowup regime, we devise a numerical scheme encoding that local structure where the blowup prefactor is determined by a global charging-annihilation balance. We also employ asymptotic analysis to construct a universal problem governing the blowup solutions in the regime Rẽ≫1, far beyond the blowup threshold. In the case ϖ>0, where charge relaxation is faster in the drop phase and the interface polarizes parallel to the external field, we numerically observe and asymptotically characterize the formation at large Rẽ of stagnant, perfectly conducting surface-charge caps about the drop poles. The cap size grows and the cap voltage decreases monotonically with increasing conductivity or decreasing permittivity of the drop phase relative to the suspending phase. The flow in this scenario is nonlinearly driven by electrical shear stresses at the complement of the caps. In both polarization scenarios, the flow at large Rẽ scales linearly with the magnitude of the external field, contrasting the familiar quadratic scaling under weak fields. Published by the American Physical Society 2024

While the problem governing Stokes flow about a single particle that is subject to an external force is ill posed in two dimensions (the ‘Stokes paradox’), the related problem of two mutually repellent particles is well posed. Motivated by self-assembly phenomena in thin viscous membranes, we consider this problem in the limit of remote particles. Such limits are typically handled in the literature using reflection techniques, which provide successive approximations to the mutual hydrodynamic interactions. Since their starting point is a single particle in an unbounded fluid domain, these techniques are futile in the present two-dimensional problem. We show how this apparent contradiction is resolved via use of singular perturbations. We obtain a two-term approximation for the velocity acquired by circular disks, considering both rigid and free particle surfaces. We also illustrate our perturbation scheme for elliptic disks, deriving a renormalised single-particle velocity. The utility of our asymptotic scheme is illustrated in the general problem of hydrodynamic interaction between a cluster of remote disks.

We consider phoretic self-propulsion of a chemically active colloid where solute is produced on the colloid surface (with a spatially varying rate) and consumed in the bulk solution (or vice versa). Assuming first-order kinetics, the dimensionless transport problem is governed by the surface Damkohler number S and the bulk Damkohler number B. The dimensionless colloid velocity U, normalized by a self-phoretic scale, is a nonlinear function of these two parameters. In the limit of small S, the solute flux is effectively prescribed by the surface activity distribution, resulting in an explicit expression for U that is proportional to S. In the limit of large B, the deviations of solute concentration from the equilibrium value are restricted to a narrow layer about the active portion of the colloid boundary. The associated boundary-layer analysis yields another explicit expression for U. Both asymptotic predictions are corroborated by an eigenfunction expansion solution of the exact problem for the cases when all physical parameters are held fixed except for a varying colloid size (resulting in S proportional to B1/2) or a varying solute diffusivity (resulting in S proportional to B). The boundary-layer structure breaks down near the transition between the active and inactive portions of the boundary. The local solution in the transition region partially resembles the classical Sommerfeld solution of wave diffraction from an edge.

We develop a reduced model for the slow unsteady dynamics of an isotropic chemically active particle near the threshold for spontaneous motion. Building on the steady theory developed in part I of this series, we match a weakly nonlinear expansion valid on the particle scale with a leading-order approximation in a larger-scale unsteady remote region, where the particle acts as a moving point source of diffusing concentration. The resulting amplitude equation for the velocity of the particle includes a term representing the interaction of the particle with its own concentration wake in the remote region, which can be expressed as a time integral over the history of the particle motion, allowing efficient simulation and theoretical analysis of fully three-dimensional unsteady problems. To illustrate how to use the model, we study the effects of a weak force acting on the particle, including the stability of the steady states and how the velocity vector realigns towards the stable one, neither of which previous axisymmetric and steady models were able to capture. This unsteady formulation could also be applied to most of the other perturbation scenarios studied in part I as well as the dynamics of interacting active particles.

A circular Hele-Shaw cell bounded by a volumetrically confined elastic solid can act as a fluidic fuse: during radially outward fluid flow, the solid deforms in response to the viscous pressure field such that the gap expands near the inlet (at the centre) and contracts near the outlet (around the rim). If the flow rate exceeds a critical value, the gap at the outlet can close completely, interrupting/choking the flow. Here, we consider the injection of gas into such a soft-walled Hele-Shaw cell filled with viscous liquid. Our theoretical model and numerical simulations for axisymmetric flow driven by the injection of an expanding gas bubble show that the bubble increases the critical flow rate of choking via two mechanisms. Firstly, as the interface approaches the rim, it reduces the length over which the viscous pressure gradient deforms the solid, which increases the critical flow rate above which choking occurs. Secondly, compression of the gas reduces the outlet flow rate relative to the inlet flow rate. As a consequence, for large injection rates, a near-choking regime is established in which the outlet flow rate becomes independent of the injection rate and instead depends only on the instantaneous position of the interface. Our travelling-wave model for the advancement of the bubble front will enable future reduced-order modelling of non-axisymmetric problems, such as viscous fingering.

Viscous flow in the narrow gap between a rigid plate and a confined elastic solid has been observed to ???choke??? at high flow rates, due to the deforming solid making contact with the plate and sealing the gap. When the viscous flow is driven by injection of a gas bubble, the advancing meniscus is susceptible to the viscous-fingering instability. By comparing fingering experiments with axisymmetric numerical simulations, we demonstrate that, depending on the width of the fingers, the fingering instability can either promote or suppress choking, i.e., cause the system to choke when an axisymmetric system would not, or vice versa.

Hydraulic fracturing for production of oil and gas from shale formations releases fluid waste, by-products that must be managed carefully to avoid significant harm to human health and the environment. These fluids are presumed to result from a variety of fracture relaxation processes, and are commonly referred to as ‘flowback’ and ‘produced water’, depending primarily on the time scale of their appearance. Here, a model is presented for investigating the dynamics of backflows caused by the elastic relaxation of a pre-strained medium, namely a single fracture and two model fracture network systems: a single bifurcated channel and its generalization for $n$ bifurcated fracture generations. Early- and late-time asymptotic solutions are obtained for the model problems and agree well with numerical solutions. In the late-time period, the fracture apertures and backflow rates exhibit a time dependence of $t^{-1/3}$ and $t^{-4/3}$ , respectively. In addition, the pressure distributions collapse to universal curves when scaled by the maximum pressure in the system, which we calculate as a function of $n$ . The pressure gradient along the network is steepest near the outlet while the bulk of the network serves as a ‘reservoir’. Fracture networks with larger $n$ are less efficient at evicting fluids, manifested through a longer time required for a given fractional reduction of the initial volume. The developed framework may be useful for informing engineering design and environmental regulations.

Drop coalescence occurs through the rapid growth of a liquid bridge that connects the two drops. At early times after contact, the bridge dynamics is typically self-similar, with details depending on the geometry and viscosity of the liquid. In this paper we analyse the coalescence of two-dimensional viscous drops that float on a quiescent deep pool; such drops are called liquid lenses. The analysis is based on the thin-sheet equations, which were recently shown to accurately capture experiments of liquid lens coalescence. It is found that the bridge dynamics follows a self-similar solution at leading order, but, depending on the large-scale boundary conditions on the drop, significant corrections may arise to this solution. This dynamics is studied in detail using numerical simulations and through matched asymptotics. We show that the liquid lens coalescence can involve a global translation of the drops, a feature that is confirmed experimentally.

We study the spreading of viscous fluid injected under an elastic sheet, which is driven by gravity and by elastic bending and tension forces and resisted by viscous forces. The injected fluid forms a large blister and spreads outwards analogously to a viscous gravity current or a capillary droplet. The relative strengths of the three driving forces are determined by how the horizontal length scales of the system compare with three key transition length scales. Bending is dominant on small length scales, tension is dominant on intermediate length scales and gravity is dominant on large length scales. We show how to use the method of matched asymptotic expansions to predict the spreading rate and thickness profile of the blister of fluid in the seven possible asymptotic regimes, for both two-dimensional and axisymmetric geometries. Consideration of different physical effects at the fluid front increases the number of regimes yet further.

Thin, roughly horizontal low-permeability layers are a common form of large-scale heterogeneity in geological porous formations. In this paper, the dynamics of a buoyancy-driven plume in a two-dimensional layered porous medium is studied theoretically, with the aid of high-resolution numerical simulations. The medium is uniform apart from a thin, horizontal layer of a much lower permeability, located a dimensionless distance L >> 1 below the dense plume source. If the dimensionless thickness 2 epsilon L and permeability Pi of the low-permeability layer are small, the effect of the layer is found to be well parameterized by its impedance Omega = 2 epsilon L/Pi. Five different regimes of flow are identified and characterized. For Omega << L-1/3, the layer has no effect on the plume, but as Omega is increased the plume widens and spreads over the layer as a gravity current. For still larger Omega, the flow becomes destabilized by convective instabilities both below and above the layer, until, for Omega >> L, the spread of the plume is dominated by convective mixing and buoyancy is transported across the layer by diffusion alone. Analytical models for the spread of the plume over the layer in the various different regimes are presented.

Microfluidic systems are usually fabricated with soft materials that deform due to the fluid stresses. Recent experimental and theoretical studies on the steady flow in shallow deformable microchannels have shown that the flow rate is a nonlinear function of the pressure drop due to the deformation of the upper soft wall. Here, we extend the steady theory of Christov et al. (J. Fluid Mech., vol. 841, 2018, pp. 267-286) by considering the start-up flow from rest, both in pressure-controlled and in flow-rate-controlled configurations. The characteristic scales and relevant parameters governing the transient flow are first identified, followed by the development of an unsteady lubrication theory assuming that the inertia of the fluid is negligible, and that the upper wall can be modelled as an elastic plate under pure bending satisfying the Kirchhoff-Love equation. The model is governed by two non-geometrical dimensionless numbers: a compliance parameter fi, which compares the characteristic displacement of the upper wall with the undeformed channel height, and a parameter that compares the inertia of the solid with its flexural rigidity. In the limit of negligible solid inertia, ! 0, a quasi-steady model is developed, whereby the fluid pressure satisfies a nonlinear diffusion equation, with fi as the only parameter, which admits a self-similar solution under pressure-controlled conditions. This simplified lubrication description is validated with coupled three-dimensional numerical simulations of the Navier equations for the elastic solid and the Navier-Stokes equations for the fluid. The agreement is very good when the hypotheses behind the model are satisfied. Unexpectedly, we find fair agreement even in cases where the solid and liquid inertia cannot be neglected.

After centuries of striving for structural rigidity, engineers and scientists alike are increasingly looking to harness the deformation, buckling, and failure of soft materials for functionality. In fluidic devices, soft deformable components that respond to the flow have the advantage of being passive; they do not require external actuation. Harnessing flow-induced deformation for passive functionality provides a means of developing flow analogs of electronic circuit components such as fluidic diodes and capacitors. The electronic component that has so far been overlooked in the microfluidics literature-the fuse-is a passive safety device that relies on a controlled failure mechanism (melting) to protect a circuit from overcurrent. Here, we describe how a compliant Hele-Shaw cell behaves in a manner analogous to the electrical fuse; above a critical flux, the flow-induced deformation of the cell blocks the outflow, interrupting (choking) the flow. In particular, the pressure distribution within the fluid applies a spatially variant normal force to the soft boundary, which causes nonuniform deformation. As a consequence of lateral confinement and incompressibility of the soft material, this flow-induced elastic deformation manifests as bulging near the cell outflow; bulges that come into contact with the rigid cell roof interrupt the flow. We identify two nondimensional parameters that govern the central deflection and the choking of the cell, respectively. This study therefore provides the mechanical foundations for engineering passive-flow limiters into fluidic devices.

We study the mechanisms affecting the viscous-fingering instability in an elastic-walled Hele-Shaw cell by considering the stability of steady states of unidirectional peeling-by-pulling and peeling-by-bending. We demonstrate that the elasticity of the wall influences the steady base state but has a negligible direct effect on the behaviour of linear perturbations, which thus behave like in the ‘printer’s instability’ with rigid walls. Moreover, the geometry of the cell can be very well approximated as a triangular wedge in the stability analysis. We identify four distinct mechanisms – surface tension acting on the horizontal and the vertical interfacial curvatures, kinematic compression in the longitudinal base flow, and the films deposited on the cell walls – that each contribute to stabilizing the system. The vertical curvature is the dominant stabilizing mechanism for small capillary numbers, but all four mechanisms have a significant effect in a large region of parameter space.

We develop a model for predicting the flow resulting from the relaxation of pre-strained, fluid-filled, elastic network structures. This model may be useful for understanding relaxation processes in various systems, e.g. deformable microfluidic systems or by-products from hydraulic fracturing operations. The analysis is aimed at elucidating features that may provide insight on the rate of fluid drainage from fracturing operations. The model structure is a bifurcating network made of fractures with uniform length and elastic modulus, which allows for general self-similar branching and variation in fracture length and rigidity between fractures along the flow path. A late-time $t<^>{-1/3}$ power law is attained and the physical behaviour can be classified into four distinct regimes that describe the late-time dynamics based on the location of the bulk of the fluid volume (which shifts away from the outlet as branching is increased) and pressure drop (which shifts away from the outlet as rigidity is increased upstream) along the network. We develop asymptotic solutions for each of the regimes, predicting the late-time flux and evolution of the pressure distribution. The effects of the various parameters on the outlet flux and the network's drainage efficiency are investigated and show that added branching and a decrease in rigidity upstream tend to increase drainage time.