How fluid-mechanical erosion creates anisotropic porous media

Abstract

Using a Cauchy integral formulation of the boundary integral equations, we simulate the erosion a porous medium comprised of up to 100 solid bodies embedded in a Stokes flow. The grains of the medium are resolved individually and erode under the action of surface shear stress. Through nonlinear feedback with the surrounding flow fields, microscopic changes in grain morphology give way to larger-scale features in the medium such as channelization. The Cauchy-integral formulation and associated quadrature formulas enable us to resolve dense configurations of nearly contacting bodies. We observe substantial anisotropy to develop over the course of erosion; that is, the configurations that result from erosion generally permit flow in the longitudinal direction more easily than in the transverse direction by up to a factor of six. These results suggest that the erosion of solid material from groundwater flows may contribute to previously observed anisotropy of natural porous media.

References in corpus (7)

• Multiscale dynamics of colloidal deposition and erosion in porous media
• Anomalous convective flows carve pinnacles and scallops in melting ice
• Shape and dynamics of seepage erosion in a horizontal granular bed
• Channelization in Porous Media driven by Erosion and Deposition
• Morphological attractors in natural convective dissolution
• Evolution of an eroding cylinder in single and lattice arrangements
• Alcove formation in dissolving cliffs driven by density inversion instability