Direct numerical simulation of wave propagation in saturated random granular packings using coupled LBM-DEM

¹ Faculty of Engineering Technology, University of Twente, P.O. Box 217, 7500 AE Enschede, the Netherlands
² Forschungszentrum Jülich, IEK-11, Fürther Straße 248, 90429 Nürnberg, Germany
³ Department of Applied Physics, TU/e, P.O. Box 513, 5600 MB Eindhoven, the Netherlands

^* e-mail: h.cheng@utwente.nl

Published online: 7 June 2021

Abstract

Poroelasticity theory predicts wave velocities in a saturated porous medium through a coupling between the bulk deformation of the solid skeleton and porous fluid flow. The challenge emerges below the characteristic wavelengths at which hydrodynamic interactions between grains and pore fluid become important. We investigate the pressure and volume fraction dependence of compressional- and shear-wave velocities in fluid-saturated, random, isotropic, frictional granular packings. The lattice Boltzmann method (LBM) and discrete element method (DEM) are two-way coupled to capture the particle-pore fluid interactions; an acoustic source is implemented to insert a traveling wave from the fluid reservoir to the saturated medium. We extract wave velocities from the acoustic branches in the wavenumber-frequency space, for a range of confining pressures and volume fractions. For random isotropic granular media the pressure-wave velocity data collapse on a single curve when scaled properly by the volume fraction.