Accurately assessing the erodibility of geomaterials is of great significance for the design of earthen structures and the prevention of the associated failure induced by seepage force. Recently, the un-resolved Computational Fluid Dynamics–Discrete Element Method (CFD–DEM) has been widely used to investigate internal erosion. However, due to the use of wall boundary and the fact that the fixed CFD domain cannot be changed wit...