Dispersive entrainment into gravity currents in porous mediaS

</jats:p

Self-similar dynamics of two-phase flows injected into a confined porous layer

We study the dynamics of two-phase flows injected into a confined porous layer. A model is derived to describe the evolution of the fluid–fluid interface, where the effective saturation of the injected fluid is zero. The flow is driven by pressure gradients due to injection, the buoyancy due to density contrasts and the interfacial tension between the injected and ambient fluids. The saturation field is then computed after the interface evolution is obtained. The results demonstrate that the flow behaviour evolves from early-time unconfined to late-time confined behaviours. In particular, at early times, the influence of capillary forces drives fluid flow and produces a new self-similar spreading behaviour in the unconfined limit that is distinct from the gravity current solution. At late times, we obtain two new similarity solutions, a modified shock solution and a compound wave solution, in addition to the rarefaction and shock solutions in the sharp-interface limit. A schematic regime diagram is also provided, which summarises all possible similarity solutions and the time transitions between them for the partially saturating flows resulting from fluid injection into a confined porous layer. Three dimensionless control parameters are identified and their influence on the fluid flow is also discussed, including the viscosity ratio, the pore-size distribution and the relative contributions of capillary and buoyancy forces. To underline the relevance of our results, we also briefly describe the implications of the two-phase flow model to the geological storage of $\text{CO}_{2}$, using representative geological parameters from the Sleipner and In Salah sites.NER

Static and dynamic fluid-driven fracturing of adhered elastica

The transient spreading of a viscous fluid beneath an elastic sheet adhered to the substrate is controlled by the dynamics at the tip where the divergence of viscous stresses necessitates the formation of a vapour tip separating the fluid front and fracture front. The model for elastic-plated currents is extended for an axisymmetric geometry with analysis showing that adhesion gives rise to the possibility of static, elastic droplets and to two dynamical regimes of spreading; viscosity dominant spreading controlled by flow of viscous fluid into the vapour tip, and adhesion dominant spreading. Constant flux experiments using clear, PDMS elastic sheets enable new, direct measurements of the vapour tip, and confirm the existence of spreading regimes controlled by viscosity and adhesion. The theory and experiments thereby provide an important test coupling the dynamics of flow with elastic deformation and have implications in fluid-driven fracturing of elastic media more generally.J. A. N. is supported by a Royal Society University Research Fellowship

Deformation of an elastic beam on a winkler foundation

AbstractWe present a simple model for geophysical systems involving sources of deformation, such as magmatic intrusions, subglacial lakes, and the subsurface storage of CO2. We consider the idealized system of a uniform elastic layer overlying a localized region of constant pressure that is surrounded by a Winkler foundation composed of springs. We investigate the effect of source depth and foundation stiffness on the resulting displacement profiles at both the surface and the level of the source. The system is characterized by three key features: the maximum uplift, the maximum subsidence, and the distance to the point of zero displacement. For each of these, we determine asymptotic scaling behavior in the limits of a thin/thick layer and a soft/stiff foundation and form composite curves that allow specific parameter values to be determined from field data. Both two-dimensional and axisymmetric pressure patches are considered, and in the thin-layer limit we derive analytical solutions.</jats:p

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

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

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

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