LAUSR

Abstract Soil erosion is a critical process that is being studied in soil science, hydraulic engineering, and geotechnical engineering. Among many societal and environmental impacts, soil erosion is a major cause for the failures of bridges. The erodibility of soil is determined by its physical and geochemical properties and is also affected by surrounding biological activities. In most of the current models for soil erosion, erodibility of non-cohesive soil is characterized by its median grain size (D50), density, and porosity. The contribution to erodibility of the irregular shape of soil grains, which plays an important role in the mechanical and hydraulic properties of coarse-grained soils, is generally ignored. In this paper, a coupled computational fluid dynamics and discrete element method model is developed to analyze the influence of the shape of sand grain on soil erodibility. A numerical model for the drag force on spherical and non-spherical particles is verified by using the results from physical free settling experiments. Erosion of sand grains of different shapes is simulated in a virtual erosion function apparatus, a laboratory device used to measure soil erodibility. The simulation results indicate that the grain shape has major effects on erodibility. Spherical particles do not show a critical velocity because of their low rolling resistance, but a critical velocity does exist for angular particles owing to grain interlocking. The erosion rate is proportional to the flow velocity for both spherical and non-spherical particles. The simulation result for angular particle erosion is fairly consistent with the experimental observations, implying that grain shape is an important factor affecting the erodibility of non-cohesive soils.