Bubbles breaking the wall: Two-dimensional stress and stability analysis

Frictional fluid instabilities shaped by viscous forces

Multiphase flows involving granular materials are complex and prone to pattern formation caused by competing mechanical and hydrodynamic interactions. Here we study the interplay between granular bulldozing and the stabilising effect of viscous pressure gradients in the invading fluid. Injection of aqueous solutions into layers of dry, hydrophobic grains represent a viscously stable scenario where we observe a transition from growth of a single frictional finger to simultaneous growth of multiple fingers as viscous forces are increased. The pattern is made more compact by the internal viscous pressure gradient, ultimately resulting in a fully stabilised front of frictional fingers advancing as a radial spoke pattern.

Pressure evolution and deformation of confined granular media during pneumatic fracturing

By means of digital image correlation, we experimentally characterize the deformation of a dry granular medium confined inside a Hele-Shaw cell due to air injection at a constant overpressure high enough to deform it (from 50 to 250 kPa). Air injection at these overpressures leads to the formation of so-called pneumatic fractures, i.e., channels empty of beads, and we discuss the typical deformations of the medium surrounding these structures. In addition we simulate the diffusion of the fluid overpressure into the medium, comparing it with the Laplacian solution over time and relating pressure gradients with corresponding granular displacements. In the compacting medium we show that the diffusing pressure field becomes similar to the Laplace solution on the order of a characteristic time given by the properties of the pore fluid, the granular medium, and the system size. However, before the diffusing pressure approaches the Laplace solution on the system scale, we find that it resembles the Laplacian field near the channels, with the highest pressure gradients on the most advanced channel tips and a screened pressure gradient behind them. We show that the granular displacements more or less always move in the direction against the local pressure gradients, and when comparing granular velocities with pressure gradients in the zone ahead of channels, we observe a Bingham type of rheology for the granular paste (the mix of air and beads), with an effective viscosity μ_{B} and displacement thresholds ∇[over ⃗]P_{c} evolving during mobilization and compaction of the medium. Such a rheology, with disorder in the displacement thresholds, could be responsible for placing the pattern growth at moderate injection pressures in a universality class like the dielectric breakdown model with η=2, where fractal dimensions are found between 1.5 and 1.6 for the patterns.

Pneumatic fractures in confined granular media

We perform experiments where air is injected at a constant overpressure P_{in}, ranging from 5 to 250 kPa, into a dry granular medium confined within a horizontal linear Hele-Shaw cell. The setup allows us to explore compacted configurations by preventing decompaction at the outer boundary, i.e., the cell outlet has a semipermeable filter such that beads are stopped while air can pass. We study the emerging patterns and dynamic growth of channels in the granular media due to fluid flow, by analyzing images captured with a high speed camera (1000 images/s). We identify four qualitatively different flow regimes, depending on the imposed overpressure, ranging from no channel formation for P_{in} below 10 kPa, to large thick channels formed by erosion and fingers merging for high P_{in} around 200 kPa. The flow regimes where channels form are characterized by typical finger thickness, final depth into the medium, and growth dynamics. The shape of the finger tips during growth is studied by looking at the finger width w as function of distance d from the tip. The tip profile is found to follow w(d)∝d^{β}, where β=0.68 is a typical value for all experiments, also over time. This indicates a singularity in the curvature d^{2}d/dw^{2}∼κ∼d^{1-2β}, but not of the slope dw/dd∼d^{β-1}, i.e., more rounded tips rather than pointy cusps, as they would be for the case β>1. For increasing P_{in}, the channels generally grow faster and deeper into the medium. We show that the channel length along the flow direction has a linear growth with time initially, followed by a power-law decay of growth velocity with time as the channel approaches its final length. A closer look reveals that the initial growth velocity v_{0} is found to scale with injection pressure as v_{0}∝P_{in}^{3/2}, while at a critical time t_{c} there is a cross-over to the behavior v(t)∝t^{-α}, where α is close to 2.5 for all experiments. Finally, we explore the fractal dimension of the fully developed patterns. For example, for patterns resulting from intermediate P_{in} around 100-150 kPa, we find that the box-counting dimensions lie within the range D_{B}∈[1.53,1.62], similar to viscous fingering fractals in porous media.

Frictional Fluid Dynamics and Plug Formation in Multiphase Millifluidic Flow

We study experimentally the flow and patterning of a granular suspension displaced by air inside a narrow tube. The invading air-liquid interface accumulates a plug of granular material that clogs the tube due to friction with the confining walls. The gas percolates through the static plug once the gas pressure exceeds the pore capillary entry pressure of the packed grains, and a moving accumulation front is reestablished at the far side of the plug. The process repeats, such that the advancing interface leaves a trail of plugs in its wake. Further, we show that the system undergoes a fluidization transition-and complete evacuation of the granular suspension-when the liquid withdrawal rate increases beyond a critical value. An analytical model of the stability condition for the granular accumulation predicts the flow regime.

Evolution of Porosity and Channelization of an Erosive Medium Driven by Fluid Flow

We demonstrate that a homogeneous porous medium composed of sedimentary particles develops channels due to curvature driven growth of fluid flow coupled with an increase in porosity. While the flux is increased linearly, the evolution of porosity is observed to be intermittent with erosion occurring at the boundaries between low and high porosity regions. Calculating the spatial distribution of the flow within the medium and the fluid stress given by the product of the fluid flux and the volume fraction of the particles, we find that the system organizes itself to be locally near the threshold needed to erode the weakest particles. A statistical model simulating the coupling of the erosion, transport, and deposition of the particles to the local fluid flow and porosity is found to capture the overall development of the observed channels.

Numerical approach to frictional fingers

Experiments on confined two-phase flow systems, involving air and a dense suspension, have revealed a diverse set of flow morphologies. As the air displaces the suspension, the beads that make up the suspension can accumulate along the interface. The dynamics can generate "frictional fingers" of air coated by densely packed grains. We present here a simplified model for the dynamics together with a new numerical strategy for simulating the frictional finger behavior. The model is based on the yield stress criterion of the interface. The discretization scheme allows for simulating a larger range of structures than previous approaches. We further make theoretical predictions for the characteristic width associated with the frictional fingers, based on the yield stress criterion, and compare these to experimental results. The agreement between theory and experiments validates our model and allows us to estimate the unknown parameter in the yield stress criterion, which we use in the simulations.