Locomotion with a wavy cylindrical filament in a yield-stress fluid

A yield stress is added to Taylor’s (1952, Proc. Royal Soc. A, 211, 225-239) model of a microscopic organism with a wavy cylindrical tail swimming through a viscous fluid. Viscoplastic slender-body theory is employed for the task, generalising existing results for Bingham fluid to the Herschel–Bulkley constitutive model. Numerical solutions are provided over a range of the two key parameters of the problem: the wave amplitude relative to the wavelength, and a Bingham number which describes the strength of the yield stress. Numerical solutions are supplemented with discussions of various limits of the problem in which analytical progress is possible. If the wave amplitude is sufficiently small, the yield stress of the material inevitably dominates the flow; the resulting ‘plastic locomotion’ results in swimming speeds that depend strongly on the swimming gait, and can, in some cases, even be negative. Conversely, when the yield stress is large, swimming becomes possible at the wave speed, with the swimmer sliding or burrowing along its centreline with a relatively high efficiency

Taylor's swimming sheet in a yield-stress fluid

A yield stress is added to Taylor’s (Proc. R. Soc. Lond. A, vol. 209, 1951, pp. 447–461) model of a two-dimensional flexible sheet swimming through a viscous fluid. Both transverse waves along the sheet, as in Taylor’s original model, and longitudinal waves are considered as means of locomotion. In each case, numerical solutions are provided over a range of the two key parameters of the problem: the wave amplitude relative to the wavelength and a Bingham number which describes the strength of the yield stress. The numerical solutions are supplemented with discussions of various limits of the problem in which analytical progress is possible. When the yield stress is large, the swimming speed for low wave amplitude is exactly double that for a Newtonian fluid, for either type of wave

Viscoplastic slender-body theory

The theory of slow viscous flow around a slender body is generalized to the situation where the ambient fluid has a yield stress. The local flow around a cylinder that is moving along or perpendicular to its axis, and rotating, provides a first step in this theory. Unlike for a Newtonian fluid, the nonlinearity associated with the viscoplastic constitutive law precludes one from linearly superposing solutions corresponding to each independent component of motion, and instead demands a full numerical approach to the problem. This is accomplished for the case of a Bingham fluid, along with a consideration of some asymptotic limits in which analytical progress is possible. Since the yield stress of the fluid strongly localizes the flow around the body, the leading-order slender-body approximation is rendered significantly more accurate than the equivalent Newtonian problem. The theory is applied to the sedimentation of inclined cylinders, bent rods and helices, and compared with some experimental data. Finally, the theory is applied to the locomotion of a cylindrical filament driven by helical waves through a viscoplastic fluid

The dynamics of miscible viscous fingering from onset to shutdown

We examine the full ‘life cycle’ of miscible viscous fingering from onset to shutdown with the aid of high-resolution numerical simulations. We study the injection of one fluid into a planar two-dimensional porous medium containing another, more viscous fluid. We find that the dynamics are distinguished by three regimes: an early-time linearly unstable regime, an intermediate-time nonlinear regime and a late-time single-finger exchange-flow regime. In the first regime, the flow can be linearly unstable to perturbations that grow exponentially. We identify, using linear stability theory and numerical simulations, a critical Péclet number below which the flow remains stable for all times. In the second regime, the flow is dominated by the nonlinear coalescence of fingers which form a mixing zone in which we observe that the convective mixing rate, characterized by a convective Nusselt number, exhibits power-law growth. In this second regime we derive a model for the transversely averaged concentration which shows good agreement with our numerical experiments and extends previous empirical models. Finally, we identify a new final exchange-flow regime in which a pair of counter-propagating diffusive fingers slow exponentially. We derive an analytic solution for this single-finger state which agrees well with numerical simulations. We demonstrate that the flow always evolves to this regime, irrespective of the viscosity ratio and Péclet number, in contrast to previous suggestions

Shear dilation of subglacial till results in time-dependent sliding laws

The dynamics of glacial sliding over water-saturated tills are poorly constrained and difficult to capture realistically in large-scale models. Experiments characterize till as a plastic material with a pressure-dependent yield stress, but the subglacial water pressure may fluctuate on annual to daily timescales, leading to transient adjustment of the till. We construct a continuum two-phase model of coupled fluid and solid deformation, describing the movement of water through the pore space of a till that is itself dilating and deforming. By forcing the model with time-dependent effective pressure at the ice-till interface, we infer the resulting relationships between basal traction, solid fraction and rate of deformation. We find that shear dilation introduces internal pressure variations and transient dilatant strengthening emerges, leading to hysteretic behaviour in low-permeability materials. The result is a time-dependent effective sliding law, with permeability-dependent lag between changes in effective pressure and the slidingspeed. This deviation from traditional steady-state sliding laws may play an important role in a wide range of transient ice-sheet phenomena, from glacier surges to the tidal response of ice streams

Stable and unstable miscible displacements in layered porous media

The effect of permeability heterogeneities and viscosity variations on miscible displacement processes in porous media is examined using high-resolution numerical simulations and reduced theoretical modelling. The planar injection of one fluid into a fluid-saturated, two-dimensional porous medium with a permeability that varies perpendicular to the flow direction is studied. Three cases are considered, in which the injected fluid is equally viscous, more viscous or less viscous than the ambient fluid. In general it is found that the flow in each case evolves through three regimes. At early times, the flow exhibits the concentration evolves diffusively, independent of both the permeability structure and the viscosity ratio. At intermediate times, the flow exhibits different dynamics including channelling and fingering, depending on whether the injected fluid is more or less viscous than the ambient fluid, and depending on the relative magnitude of the viscosity and permeability variations. Finally, at late times, the flow becomes independent of the viscosity ratio and dominated by shear-enhanced (Taylor) dispersion. For each of the regimes identified above, we develop reduced-order models for the evolution of the transversely averaged concentration and compare them to the full numerical simulations

Building on Oldroyd’s viscoplastic legacy: Perspectives and new developments

The decade following the second world war heralded the publication of a collection of important papers on non-Newtonian fluid mechanics; Oldroyd’s work featured heavily in this collection. Not only did these articles establish important results, but Oldroyd’s style and methods set the scene for subsequent work in the area, exploiting mathematical analysis to formulate problems, establish results and guide further research. While Oldroyd’s name will forever be linked with the study of elastic fluids, the purpose of the present paper is to offer a modern perspective on a number of Oldroyd’s papers on viscoplastic fluids from 1947–1951 [1], [2], [3], [4], [5], [6], [7], [8]. Along the way, we sprinkle in a brief review of some of the subsequent developments stemming from Oldroyd’s advances, together with a few new results guided by his work. Following the approach of most of Oldroyd’s original papers, we focus on unidirectional flow down conduits. In an Appendix, we complement this discussion with a lubrication analysis, extending, clarifying and correcting the important original analysis of Walton and Bittleston (1991) [9]; although lubrication theory was not directly utilized by Oldroyd, the methodology aligns with his philosophy of using asymptotic and analytical approaches

The influence of a poroelastic till on rapid subglacial flooding and cavity formation

We develop a model of the rapid propagation of water at the contact between elastic glacial ice and a poroelastic subglacial till, motivated by observations of the rapid drainage of supraglacial lakes in Greenland. By treating the ice as an elastic bending beam, the fluid dynamics of contact with the subglacial hydrological network, which is modelled as a saturated poroelastic till, can be examined in detail. The model describes the formation and dynamics of an axisymmetric subglacial cavity, and the spread of pore pressure, in response to injection of fluid. A combination of numerical simulation and asymptotic analysis is used to describe these dynamics for both a rigid and a deformable porous till, and for both laminar and turbulent fluid flow. For constant injection rates and laminar flow, the cavity is isostatic and its spread is controlled by bending of the ice and suction of pore water in the vicinity of the ice–till contact. For a deformable till, this control can be modified: generically, a flexural wave that is initially trapped in advance of the contact point relaxes over time by diffusion of pore pressure ahead of the cavity. While the dynamics are found to be relatively insensitive to the properties of the subglacial till during injection with a constant flux, significant dependence on the till properties is manifest during the subsequent spread of a constant volume. A simple hybrid turbulent–laminar model is presented to account for fast injection rates of water: in this case, self-similar turbulent propagation can initially control the spread of the cavity, but there is a transition to laminar control in the vicinity of the ice–till contact point as the flow slows. Finally, the model results are compared with recent geophysical observations of the rapid drainage of supraglacial lakes in Greenland; the comparison provides qualitative agreement and raises suggestions for future quantitative comparison.</jats:p

Tidal Grounding-Line Migration Modulated by Subglacial Hydrology

We present a mathematical model of the hydrology of grounding-line migration on tidal timescales, in which the ice acts elastically, overlying a connected hydrological network, with the ocean tides modelled by an oscillating far-field fluid height. The upstream grounding-line migration is driven by a fluid pressure gradient through the grounding zone, while the downstream migration is limited by fluid drainage through the till. The two processes are described using separate travelling-wave solutions, based on a model of fluid flow under an elastic sheet. The asymmetry between the up- and downstream motion allows the grounding line to act as a non-linear filter on the tidal forcing as the pressure signal propagates upstream, and this frequency modulation is discussed in the context of velocity data from ice streams across Antarctica to provide a novel constraint on till permeability

Translating and squirming cylinders in a viscoplastic fluid

Three related problems of viscoplastic flow around cylinders are considered. First, translating cylinders with no-slip surfaces appear to generate adjacent rotating plugs in the limit where the translation speed becomes vanishingly small. In this plastic limit, analytical results are available from plasticity theory (slipline theory) which indicate that no such plugs should exist. Using a combination of numerical computations and asymptotic analysis, we show that the plugs of the viscoplastic theory actually disappear in the plastic limit, albeit very slowly. Second, when the boundary condition on the cylinder is replaced by one that permits sliding, the plastic limit corresponds to a partially rough cylinder. In this case, no plasticity solution has been previously established; we provide evidence from numerical computations and slipline theory that a previously proposed upper bound (Martin & Randolph, Geotechnique, vol. 56, 2006, pp. 141–145) is actually the true plastic solution. Third, we consider how a prescribed surface velocity field can propel cylindrical squirmers through a viscoplastic fluid. We determine swimming speeds and contrast the results with those from the corresponding Newtonian problem