In this paper the Discrete Element Method (DEM) is coupled with the Lattice-Boltzmann Method (LBM) to model the undrained condition of dense granular media that display significant dilation under highly confined loading. DEM-only models are commonly used to simulate the micromechanics of an undrained specimen by applying displacements at the domain boundaries so that the specimen volume remains constant. While this approach works well for uniform strain conditions found in laboratory tests, it doesn’t realistically represent non-uniform strain conditions that exist in the majority of real geotechnical problems. The LBM offers a more realistic approach to simulate the undrained condition since the fluid can locally conserve the system volume. To investigate the ability of the DEM-LBM model to effectively represent the undrained constraint while conserving volume and accurately calculating the stress path of the system, a two dimensional biaxial test is simulated using the coupled DEM-LBM model, and the results are compared with those attained from a DEM-only constant volume simulation. The compressibility of the LBM fluid was found to play an important role in the model response. The compressibility of the fluid is expressed as an apparent Skempton’s pore pressure parameter B. The biaxial test, both with and without fluid, demonstrated particle-scale instabilities associated with shear band development. The results show that the DEM-LBM model offers a promising technique for a variety of geomechanical problems that involve particle-fluid mixtures undergoing large deformation under shear loading.
