Channelization in porous media driven by erosion and deposition

Modeling the Transition between Localized and Extended Deposition in Flow Networks through Packings of Glass Beads

We use a theoretical model to explore how fluid dynamics, in particular, the pressure gradient and wall shear stress in a channel, affect the deposition of particles flowing in a microfluidic network. Experiments on transport of colloidal particles in pressure-driven systems of packed beads have shown that at lower pressure drop, particles deposit locally at the inlet, while at higher pressure drop, they deposit uniformly along the direction of flow. We develop a mathematical model and use agent-based simulations to capture these essential qualitative features observed in experiments. We explore the deposition profile over a two-dimensional phase diagram defined in terms of the pressure and shear stress threshold, and show that two distinct phases exist. We explain this apparent phase transition by drawing an analogy to simple one-dimensional mass-aggregation models in which the phase transition is calculated analytically.

Temporal Evolution of Erosion in Pore Networks: From Homogenization to Instability

We study the dynamics of flow networks in porous media using two and three dimensional pore-network models. We consider a class of erosion dynamics for a single phase flow with no deposition, chemical reactions, or topology changes assuming a constitutive law depending on flow rate, local velocities, or shear stress at the walls. We show that depending on the erosion law, the flow may become uniform and homogenized or become unstable and develop channels. By defining an order parameter capturing these different behaviors we show that a phase transition occurs depending on the erosion dynamics. Using a simple model, we identify quantitative criteria to distinguish these regimes and correctly predict the fate of the network, and discuss the experimental relevance of our result.

Interpolation for the lattice-Boltzmann method to simulate colloid transport in porous media

The lattice-Boltzmann method is convenient for simulating flow fields in porous media. However, due to its lattice characteristics, the velocity near a solid surface is not accurate, which results in significant errors when simulating colloid transport in porous media. Based on the general properties of a flow field close to a solid surface, we propose an alternative velocity interpolation method in which the velocity at a solid surface is strictly zero. Numerical simulation results show that the proposed method can give more accurate results than the usual bilinear interpolation. In addition, we use this method to simulate the contact efficiency of colloids in porous media and obtain a new power-law form of the contact efficiency.

Multiscale dynamics of colloidal deposition and erosion in porous media

Diverse processes-e.g., environmental pollution, groundwater remediation, oil recovery, filtration, and drug delivery-involve the transport of colloidal particles in porous media. Using confocal microscopy, we directly visualize this process in situ and thereby identify the fundamental mechanisms by which particles are distributed throughout a medium. At high injection pressures, hydrodynamic stresses cause particles to be continually deposited on and eroded from the solid matrix-notably, forcing them to be distributed throughout the entire medium. By contrast, at low injection pressures, the relative influence of erosion is suppressed, causing particles to localize near the inlet of the medium. Unexpectedly, these macroscopic distribution behaviors depend on imposed pressure in similar ways for particles of different charges, although the pore-scale distribution of deposition is sensitive to particle charge. These results reveal how the multiscale interactions between fluid, particles, and the solid matrix control how colloids are distributed in a porous medium.

Fluid flow through packings of elastic shells

Fluid transport in porous materials is commonly studied in geological samples (soil, sediments, etc.) or idealized systems, but the fluid flow through compacted granular materials, consisting of substantially strained granules, remains relatively unexplored. As a step toward filling this gap, we study a model of liquid transport in packings of deformable elastic shells using finite-element and lattice-Boltzmann methods. We find that the fluid flow abruptly vanishes as the porosity of the material falls below a critical value, and the flow obstruction exhibits features of a percolation transition. We further show that the fluid flow can be captured by a simplified permeability model in which the complex porous material is replaced by a collection of disordered capillaries, which are distributed and shaped by the percolation transition. To that end, we numerically explore the divergence of hydraulic tortuosity τ_{H} and the decrease of a hydraulic radius R_{h} as the percolation threshold is approached. We interpret our results in terms of scaling predictions derived from the percolation theory applied to random packings of spheres.

Headward growth and branching in subterranean channels

We investigate the erosive growth of channels in a thin subsurface sedimentary layer driven by hydrodynamic drag toward understanding subterranean networks and their relation to river networks charged by ground water. Building on a model based on experimental observations of fluid-driven evolution of bed porosity, we focus on the characteristics of the channel growth and their bifurcations in a horizontal rectangular domain subject to various fluid source and sink distributions. We find that the erosion front between low- and high-porosity regions becomes unstable, giving rise to branched channel networks, depending on the spatial fluctuations of the fluid flow near the front and the degree to which the flow is above the erodibility threshold of the medium. Focusing on the growth of a network starting from a single channel, and by identifying the channel heads and their branch points, we find that the number of branches increases sublinearly and is affected by the source distribution. The mean angles between branches are found to be systematically lower than river networks in humid climates and depend on the domain geometry.

Mechanism behind Erosive Bursts In Porous Media

Erosion and deposition during flow through porous media can lead to large erosive bursts that manifest as jumps in permeability and pressure loss. Here we reveal that the cause of these bursts is the reopening of clogged pores when the pressure difference between two opposite sites of the pore surpasses a certain threshold. We perform numerical simulations of flow through porous media and compare our predictions to experimental results, recovering with excellent agreement shape and power-law distribution of pressure loss jumps, and the behavior of the permeability jumps as a function of particle concentration. Furthermore, we find that erosive bursts only occur for pressure gradient thresholds within the range of two critical values, independent of how the flow is driven. Our findings provide a better understanding of sudden sand production in oil wells and breakthrough in filtration.