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

DOI10.1016/J.JCP.2018.07.037MaRDI QIDQ2002225zbMATH OpenOpenAlexWikidataFDO

Authors Bryan D. Quaife, M. Nicholas J. Moore

Publication date 11 July 2019

Published in Journal of Computational Physics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1804.06016

zbMATH Keywords

\(\theta\)-\(L\) formulationboundary integral equationserosionporous mediumshear stressStokes flow

Mathematics Subject Classification ID

Numerical methods for integral equations (65R20) Fluid-solid interactions (including aero- and hydro-elasticity, porosity, etc.) (74F10) Stokes and related (Oseen, etc.) flows (76D07) Flows in porous media; filtration; seepage (76S05)

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 $heta$-$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 $mathcalO(N)$ operations per GMRES iteration, a mesh-independent number of GMRES iterations, and a one-time $mathcalO(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.

Recommendations

• How fluid-mechanical erosion creates anisotropic porous media
• Simulating flow over and through porous media with application to erosion of particulate deposits
• Viscous transport in eroding porous media
• Numerical modelling of interfacial soil erosion with viscous incompressible flows
• Flows in porous media with erosion of the solid matrix

Cites work

• A fast algorithm for particle simulations
• A fast and stable solver for singular integral equations on piecewise smooth curves
• A fast integral equation method for solid particles in viscous flow using quadrature by expansion
• A generalized traction integral equation for Stokes flow, with applications to near-wall particle mobility and viscous erosion
• A kernel-independent adaptive fast multipole algorithm in two and three dimensions
• A method for computing nearly singular integrals
• A simplified technique for the efficient and highly accurate discretization of boundary integral equations in 2D on domains with corners
• A unified integral equation scheme for doubly periodic Laplace and Stokes boundary value problems in two dimensions
• Adaptive time stepping for vesicle suspensions
• An integral equation formulation for rigid bodies in Stokes flow in three dimensions
• Asymmetric collapse by dissolution or melting in a uniform flow
• Boundary Integral and Singularity Methods for Linearized Viscous Flow
• Boundary integral equation analysis for suspension of spheres in Stokes flow
• Channelization of plumes beneath ice shelves
• Contact-aware simulations of particulate Stokesian suspensions
• Convergence estimates for solution of integral equations with GMRES
• Coupled Stokes-Darcy model with Beavers-Joseph interface boundary condition
• Dissolution instability and roughening transition
• Dynamics of fluids in porous media.
• Evaluation of layer potentials close to the boundary for Laplace and Helmholtz problems on analytic planar domains
• Geometry of valley growth
• GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems
• High-volume fraction simulations of two-dimensional vesicle suspensions
• scientific article; zbMATH DE number 976282 (Why is no real title available?)
• scientific article; zbMATH DE number 3196902 (Why is no real title available?)
• Inclusion of lubrication forces in dynamic simulations
• Integral equation methods for elastance and mobility problems in two dimensions
• Integral equation methods for Stokes flow and isotropic elasticity in the plane
• Low-resolution simulations of vesicle suspensions in 2D
• On optimum profiles in Stokes flow
• On preconditioners for the Laplace double-layer in 2D.
• On the evaluation of layer potentials close to their sources
• On the numerical solution of the biharmonic equation in the plane
• On the Nyström discretization of integral equations on planar curves with corners
• On the solution of elliptic partial differential equations on regions with corners
• Quadrature by expansion: a new method for the evaluation of layer potentials
• Quadrature methods for periodic singular and weakly singular Fredholm integral equations
• Removing the stiffness from interfacial flows with surface tension
• Riemann-Hilbert problems for the shapes formed by bodies dissolving, melting, and eroding in fluid flows
• Second Kind Integral Equation Formulation of Stokes’ Flows Past a Particle of Arbitrary Shape
• Self-similar evolution of a body eroding in a fluid flow
• Slender-body theory for particles of arbitrary cross-section in Stokes flow
• Spectrally accurate quadratures for evaluation of layer potentials close to the boundary for the 2D Stokes and Laplace equations
• Steadily translating parabolic dissolution fingers
• Stratified flows with vertical layering of density: experimental and theoretical study of flow configurations and their stability
• The completed double layer boundary integral equation method for two-dimensional Stokes flow
• The exponentially convergent trapezoidal rule
• The inverse fast multipole method: using a fast approximate direct dolver as a preconditioner for dense linear systems

Cited in

(18)

• Simulating flow over and through porous media with application to erosion of particulate deposits
• Boundary integral analysis for non-homogeneous, incompressible Stokes flows
• Parallel algorithms for successive convolution
• Mathematical modeling and experimental study of erosion-deposition process in deformable porous media
• Soil erosion in the boundary layer flow along a slope: a theoretical study
• How fluid-mechanical erosion creates anisotropic porous media
• Self-similar evolution of a body eroding in a fluid flow
• Viscous transport in eroding porous media
• Convection in a Coupled Free Flow-Porous Media System
• Predicting convection configurations in coupled fluid–porous systems
• Dissolution of a cylindrical disk in Hele-Shaw flow: a conformal-mapping approach
• Multiphase modelling of precipitation-induced membrane formation
• An SDE model for deterioration of rock surfaces
• Interfacial erosion: a three-dimensional numerical model
• Trapping of Planar Brownian Motion: Full First Passage Time Distributions by Kinetic Monte Carlo, Asymptotic, and Boundary Integral Methods
• Numerical modelling of interfacial soil erosion with viscous incompressible flows
• Quadrature error estimates for kernels with logarithmic singularity
• Boundary integral methods for particle diffusion in complex geometries: shielding, confinement, and escape

This page was built for publication: A boundary-integral framework to simulate viscous erosion of a porous medium

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2002225)