A boundary-integral framework to simulate viscous erosion of a porous medium

Abstract

We develop numerical methods to simulate the fluid-mechanical erosion of many bodies in two-dimensional Stokes flow. The broad aim is to simulate the erosion of a porous medium (e.g. groundwater flow) with grain-scale resolution. Our fluid solver is based on a second-kind boundary integral formulation of the Stokes equations that is discretized with a spectrally-accurate Nystrom method and solved with fast-multipole-accelerated GMRES. The fluid solver provides the surface shear stress which is used to advance solid boundaries. We regularize interface evolution via curvature penalization using the $\theta$-$L$ formulation, which affords numerically stable treatment of stiff terms and therefore permits large time steps. The overall accuracy of our method is spectral in space and second-order in time. The method is computationally efficient, with the fluid solver requiring $\mathcal{O}(N)$ operations per GMRES iteration, a mesh-independent number of GMRES iterations, and a one-time $\mathcal{O}(N^2)$ computation to compute the shear stress. We benchmark single-body results against analytical predictions for the limiting morphology and vanishing rate. Multibody simulations reveal the spontaneous formation of channels between bodies of close initial proximity. The channelization is associated with a dramatic reduction in the resistance of the porous medium, much more than would be expected from the reduction in grain size alone.

arXiv comment10 figures, 3 tables, 28 pages

References in corpus (17)

• Quadrature by Expansion: A New Method for the Evaluation of Layer Potentials
• How river rocks round: resolving the shape-size paradox
• Dissolution instability and roughening transition
• Shape and dynamics of seepage erosion in a horizontal granular bed
• A fast integral equation method for solid particles in viscous flow using quadrature by expansion
• High-volume fraction simulations of two-dimensional vesicle suspensions
• An integral equation formulation for rigid bodies in Stokes flow in three dimensions
• Geometry of Valley Growth
• Erosion patterns in a sediment layer
• Boundary integral equation analysis for suspension of spheres in Stokes flow
• Integral equation methods for elastance and mobility problems in two dimensions
• A generalized traction integral equation for Stokes flow, with applications to near-wall particle mobility and viscous erosion
• Asymmetric collapse by dissolution or melting in a uniform flow
• Contact-aware simulations of particulate Stokesian suspensions
• Low-resolution simulations of vesicle suspensions in 2D
• Evolution of an eroding cylinder in single and lattice arrangements
• The pear-shaped fate of an ice melting front

Cited by in corpus (9)

• Morphological attractors in natural convective dissolution
• Viscous Transport in Eroding Porous Media
• Predicting convection configurations in coupled fluid-porous systems
• Parallel Algorithms for Successive Convolution
• Multiphase modeling of precipitation-induced membrane formation
• How fluid-mechanical erosion creates anisotropic porous media
• Rigorous criteria for anomalous waves induced by abrupt depth change using truncated KdV statistical mechanics
• Boundary Integral Analysis for the Non-homogeneous 3D Stokes Equation
• Convection in a coupled fluid-porous media system