Indentation of a floating elastic sheet: Geometry versus applied tension

The localized loading of an elastic sheet floating on a liquid bath occurs at scales from a frog sitting on a lily pad to a volcano supported by the Earth's tectonic plates. The load is supported by a combination of the stresses within the sheet (which may include applied tensions from, for example, surface tension) and the hydrostatic pressure in the liquid. At the same time, the sheet deforms, and may wrinkle, because of the load. We study this problem in terms of the (relatively weak) applied tension and the indentation depth. For small indentation depths, we find that the force--indentation curve is linear with a stiffness that we characterize in terms of the applied tension and bending stiffness of the sheet. At larger indentations the force--indentation curve becomes nonlinear and the sheet is subject to a wrinkling instability. We study this wrinkling instability close to the buckling threshold and calculate both the number of wrinkles at onset and the indentation depth at onset, comparing our theoretical results with experiments. Finally, we contrast our results with those previously reported for very thin, highly bendable membranes.Comment: 24 pages, revised version submitted to Proc. R. Soc.

Convective shutdown in a porous medium at high Rayleigh number

Convection in a closed domain driven by a dense buoyancy source along the upper boundary soon starts to wane owing to the increase of the average interior density. In this paper, theoretical and numerical models are developed of the subsequent long period of shutdown of convection in a two-dimensional porous medium at high Rayleigh number Ra\mathit{Ra}. The aims of this paper are twofold. Firstly, the relationship between this slowly evolving ‘one-sided’ shutdown system and the statistically steady ‘two-sided’ Rayleigh–Bénard (RB) cell is investigated. Numerical measurements of the Nusselt number Nu\mathit{Nu} from an RB cell (Hewitt et al., Phys. Rev. Lett., vol. 108, 2012, 224503) are very well described by the simple parametrization Nu=2.75+0.0069Ra\mathit{Nu}= 2. 75+ 0. 0069\mathit{Ra}. This parametrization is used in theoretical box models of the one-sided shutdown system and found to give excellent agreement with high-resolution numerical simulations of this system. The dynamical structure of shutdown can also be accurately predicted by measurements from an RB cell. Results are presented for a general power-law equation of state. Secondly, these ideas are extended to model more complex physical systems, which comprise two fluid layers with an equation of state such that the solution that forms at the (moving) interface is more dense than either layer. The two fluids are either immiscible or miscible. Theoretical box models compare well with numerical simulations in the case of a flat interface between the fluids. Experimental results from a Hele-Shaw cell and numerical simulations both show that interfacial deformation can dramatically enhance the convective flux. The applicability of these results to the convective dissolution of geologically sequestered CO2{\mathrm{CO} }_{2} in a saline aquifer is discussed

Heat production and tidally driven fluid flow in the permeable core of Enceladus

Author Posting. © American Geophysical Union, 2020. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Planets 125(9), (2020): e2019JE006209, doi:10.1029/2019JE006209Saturn's moon Enceladus has a global subsurface ocean and a porous rocky core in which water‐rock reactions likely occur; it is thus regarded as a potentially habitable environment. For icy moons like Enceladus, tidal heating is considered to be the main heating mechanism, which has generally been modeled using viscoelastic solid rheologies in existing studies. Here we provide a new framework for calculating tidal heating based on a poroviscoelastic model in which the porous solid and interstitial fluid deformation are coupled. We show that the total heating rate predicted for a poroviscoelastic core is significantly larger than that predicted using a classical viscoelastic model for intermediate to large (>1014 Pa·s) rock viscosities. The periodic deformation of the porous rock matrix is accompanied by interstitial pore fluid flow, and the combined effects through viscous dissipation result in high heat fluxes particularly at the poles. The heat generated in the rock matrix is also enhanced due to the high compressibility of the porous matrix structure. For a sufficiently compressible core and high permeability, the total heat production can exceed 10 GW—a large fraction of the moon's total heat budget—without requiring unrealistically low solid viscosities. The partitioning of heating between rock and fluid constituents depends most sensitively on the viscosity of the rock matrix. As the core of Enceladus warms and weakens over time, pore fluid motion likely shifts from pressure‐driven local oscillations to buoyancy‐driven global hydrothermal convection, and the core transitions from fluid‐dominated to rock‐dominated heating.2021-01-2

Dispersive entrainment into gravity currents in porous media

he effects of dispersion acting on gravity currents propagating through porous media are considered theoretically and experimentally. We exploit the large aspect ratio of these currents to formulate a depth-averaged model of the evolution of the mass and buoyancy. Dispersion, acting predominantly at the interface between the current and the ambient, is velocity dependent and acts to entrain fluid into the gravity current, in direct analogy to turbulent mixing. Here, we show that when the gravity current is fed by a constant buoyancy and mass flux the buoyancy of the current is self-similar and recovers the classical solution for gravity currents in porous media. In contrast, the profile and the depth-averaged concentration of the current evolve in a non-self-similar manner. The total volume of the current increases with time as due to this dispersive entrainment. We test our theoretical predictions using a suite of laboratory experiments in which the evolution of the concentration within the current was mapped using a dye-attenuation technique. These experimental results show good agreement with the early-time limits of our theoretical model, and in particular accurately predict the evolution of the depth-averaged concentration profile. These results suggest that mixing within porous media may be modelled using an effective dispersive entrainment, the magnitude of which may be set by the underlying structure of the porous medium

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 vapor 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 vapor tip, and adhesion dominant spreading. Constant flux experiments using clear, PDMS elastic sheets enable new, direct measurements of the vapor 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

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

Crystal settling and convection in the Shiant Isles Main Sill

The 168 m-thick Shiant Isles Main Sill is a composite body, dominated by an early, 24 m-thick, picrite sill formed by the intrusion of a highly olivine-phyric magma, and a later 135 m-thick intrusion of olivine-phyric magma that split the earlier picrite into a 22 m-thick lower part and a 2 m-thick upper part, forming the picrodolerite/crinanite unit (PCU). The high crystal load in the early picrite prevented effective settling of the olivine crystals, which retain their initial stratigraphic distribution. In contrast, the position of the most evolved rocks of the PCU at a level ~80% of its total height point to significant accumulation of crystals on the floor, as evident by the high olivine mode at the base of the PCU. Crystal accumulation on the PCU floor occurred in two stages. During the first, most of the crystal load settled to the floor to form a modally and size-sorted accumulation dominated by olivine, leaving only the very smallest olivine grains still in suspension. The second stage is recorded by the coarsening-upwards of individual olivine grains in the picrodolerite, and their amalgamation into clusters which become both larger and better sintered with increasing stratigraphic height. Large clusters of olivine are present at the roof, forming a foreshortened mirror image of the coarsening-upwards component of the floor accumulation. The coarsening-upwards sequence records the growth of olivine crystals while in suspension in a convecting magma, and their aggregation into clusters, followed by settling over a prolonged period (with limited trapping at the roof). As olivine was progressively lost from the convecting magma, crystal accumulation on the (contemporaneous) floor of the PCU was increasingly dominated by plagioclase, most likely forming clusters and aggregates with augite and olivine, both of which form large poikilitic grains in the crinanite. While the PCU is unusual in being underlain by an earlier, still hot, intrusion that would have enhanced any driving force for convection, we conclude from comparison with microstructures in other sills that convection is likely in tabular bodies >100 m thickness

Convective shutdown in a porous medium at high Rayleigh number

Convection in a closed domain driven by a dense buoyancy source along the upper boundary soon starts to wane owing to the increase of the average interior density. In this paper, theoretical and numerical models are developed of the subsequent long period of shutdown of convection in a two-dimensional porous medium at high Rayleigh number Ra\mathit{Ra}. The aims of this paper are twofold. Firstly, the relationship between this slowly evolving ‘one-sided’ shutdown system and the statistically steady ‘two-sided’ Rayleigh–Bénard (RB) cell is investigated. Numerical measurements of the Nusselt number Nu\mathit{Nu} from an RB cell (Hewitt et al., Phys. Rev. Lett., vol. 108, 2012, 224503) are very well described by the simple parametrization Nu=2.75+0.0069Ra\mathit{Nu}= 2. 75+ 0. 0069\mathit{Ra}. This parametrization is used in theoretical box models of the one-sided shutdown system and found to give excellent agreement with high-resolution numerical simulations of this system. The dynamical structure of shutdown can also be accurately predicted by measurements from an RB cell. Results are presented for a general power-law equation of state. Secondly, these ideas are extended to model more complex physical systems, which comprise two fluid layers with an equation of state such that the solution that forms at the (moving) interface is more dense than either layer. The two fluids are either immiscible or miscible. Theoretical box models compare well with numerical simulations in the case of a flat interface between the fluids. Experimental results from a Hele-Shaw cell and numerical simulations both show that interfacial deformation can dramatically enhance the convective flux. The applicability of these results to the convective dissolution of geologically sequestered CO2{\mathrm{CO} }_{2} in a saline aquifer is discussed