Gravity currents in a porous medium at an inclined plane

We consider the release from a point source of relatively heavy fluid into a porous saturated medium above an impermeable slope. We consider the case where the volume of the resulting gravity current increases with time like $t^\alpha$ and show that for $\alpha<3$, at short times the current spreads axisymmetrically, with radius $r\sim t^{(\alpha+1)/4}$, while at long times it spreads predominantly downslope. In particular, for long times the downslope position of the current scales like $t$ while the current extends a distance $t^{\alpha/3}$ across the slope. For $\alpha>3$, this situation is reversed with spreading occurring predominantly downslope for short times. The governing equations admit similarity solutions whose scaling behaviour we determine, with the full similarity form being evaluated by numerical computations of the governing partial differential equation. We find that the results of these analyses are in good quantitative agreement with a series of laboratory experiments. Finally, we briefly discuss the implications of our work for the sequestration of carbon dioxide in aquifers with a sloping, impermeable cap.Comment: 10 pages, 6 figures - revised versio

The competition between gravity and flow focusing in two-layered porous media

The gravitationally driven flow of a dense fluid within a two-layered porous media is examined experimentally and theoretically. We find that in systems with two horizontal layers of differing permeability a competition between gravity driven flow and flow focusing along high-permeability routes can lead to two distinct flow regimes. When the lower layer is more permeable than the upper layer, gravity acts along high-permeability pathways and the flow is enhanced in the lower layer. Alternatively, when the upper layer is more permeable than the lower layer, we find that for a sufficiently small input flux the flow is confined to the lower layer. However, above a critical flux fluid preferentially spreads horizontally within the upper layer before ultimately draining back down into the lower layer. This later regime, in which the fluid overrides the low-permeability lower layer, is important because it enhances the mixing of the two fluids. We show that the critical flux which separates these two regimes can be characterized by a simple power law. Finally, we briefly discuss the relevance of this work to the geological sequestration of carbon dioxide and other industrial and natural flows in porous media

Topographic controls on gravity currents in porous media

We present a theoretical and experimental study of the propagation of gravity currents in porous media with variations in the topography over which they flow, motivated in part by the sequestration of carbon dioxide in saline aquifers. We consider cases where the height of the topography slopes upwards in the direction of the flow and is proportional to the nth power of the horizontal distance from a line or point source of a constant volumetric flux. In two-dimensional cases with n>1/2, the current evolves from a self-similar form at early times, when the effects of variations in topography are negligible, towards a late-time regime that has an approximately horizontal upper surface and whose evolution is dictated entirely by the geometry of the topography. For n<1/2, the transition between these flow regimes is reversed. We compare our theoretical results in the case n=1 with data from a series of laboratory experiments in which viscous glycerine is injected into an inclined Hele-Shaw cell, obtaining good agreement between the theoretical results and the experimental data. In the case of axisymmetric topography, all topographic exponents n>0 result in a transition from an early-time similarity solution towards a topographically controlled regime that has an approximately horizontal free surface. We also analyse the evolution over topography that can vary with different curvatures and topographic exponents between the two horizontal dimensions, finding that the flow transitions towards a horizontally topped regime at a rate which depends strongly on the ratio of the curvatures along the principle axes. Finally, we apply our mathematical solutions to the geophysical setting at the Sleipner field, concluding that topographic influence is unlikely to explain the observed non-axisymmetric flow

Interaction of viscous free-surface flows with topography

The interaction of gravitationally driven, free-surface flows of viscous fluid with topographic features is investigated theoretically. The motion is studied in the regime where the depth of the flow is much smaller than the streamwise extent of the topography. A lubrication model of the motion is developed, integrated numerically and analysed asymptotically. For small mounds, it is shown that the flow surmounts the obstacles, but for larger mounds the flow is deflected around it and can form dry zones in its wake into which fluid does not flow, as well as forming deeper ponded regions upstream. Which of these phenomena prevails is shown to depend upon the amplitude of the mound height and the thickness of the oncoming flow relative to the streamwise length scale over which the topography varies. By using numerical and asymptotic results, we demonstrate that relatively wide mounds lead to the development of deep ponds of material upstream, which may lead to flow overtopping if the mound is not sufficiently high. These insights can be used to inform the design of barriers that defend built infrastructures from lava flows, and it is shown how this model can also provide an upper bound on the force exerted by the flow on them.</jats:p

Shallow free-surface Stokes flow around a corner.

The steady lateral spreading of a free-surface viscous flow down an inclined plane around a vertex from which the channel width increases linearly with downstream distance is investigated analytically, numerically and experimentally. From the vertex the channel wall opens by an angle α to the downslope direction and the viscous fluid spreads laterally along it before detaching. The motion is modelled using lubrication theory and the distance at which the flow detaches is computed as a function of α using analytical and numerical methods. Far downslope after detachment, it is shown that the motion is accurately modelled in terms of a similarity solution. Moreover, the detachment point is well approximated by a simple expression for a broad range of opening angles. The results are corroborated through a series of laboratory experiments and the implication for the design of barriers to divert lava flows are discussed. This article is part of the theme issue 'Stokes at 200 (Part 1)'

Capillary pinning and blunting of immiscible gravity currents in porous media

Gravity-driven flows in the subsurface have attracted recent interest in the context of geological carbon dioxide (CO2) storage, where supercritical CO2 is captured from the flue gas of power plants and injected underground into deep saline aquifers. After injection, the CO2 will spread and migrate as a buoyant gravity current relative to the denser, ambient brine. Although the CO2 and the brine are immiscible, the impact of capillarity on CO2 spreading and migration is poorly understood. We previously studied the early time evolution of an immiscible gravity current, showing that capillary pressure hysteresis pins a portion of the macroscopic fluid-fluid interface and that this can eventually stop the flow. Here we study the full lifetime of such a gravity current. Using tabletop experiments in packings of glass beads, we show that the horizontal extent of the pinned region grows with time and that this is ultimately responsible for limiting the migration of the current to a finite distance. We also find that capillarity blunts the leading edge of the current, which contributes to further limiting the migration distance. Using experiments in etched micromodels, we show that the thickness of the blunted nose is controlled by the distribution of pore-throat sizes and the strength of capillarity relative to buoyancy. We develop a theoretical model that captures the evolution of immiscible gravity currents and predicts the maximum migration distance. By applying this model to representative aquifers, we show that capillary pinning and blunting can exert an important control on gravity currents in the context of geological CO2 storage