This review article focuses on the modeling of complex granular flows employing the discrete element method (DEM) approach. The specific topic discussed is the application of DEM models for the study of the flow behavior of nonspherical, flexible, or cohesive particles, including particle breakage. The major sources of particle cohesion—liquid induced, electrostatics, van der Waals forces—and their implementation into DEM simulations are covered. These aspects of particle flow are of great importance in practical applications and hence are the significant foci of research at the forefront of current DEM modeling efforts. For example, DEM simulations of nonspherical grains can provide particle stress information needed to develop constitutive models for continuum-based simulations of large-scale industrial processes.

The effect of particle aspect ratio and surface geometry on granular flows is assessed by performing numerical simulations of rod-like particles in simple shear flows using the discrete element method (DEM). The effect of particle surface geometry is explored by adopting two types of particles: glued-spheres particles and true cylindrical particles. The particle aspect ratio varies from one to six. Compared to frictionless spherical particles, smaller stresses are obtained for the glued-spheres and cylindrical particle systems in dilute and moderately dense flows due to the loss of translational energy, which is partially converted to rotational energy, for the non-spherical particles. For dilute granular flows of non-spherical particles, stresses are primarily affected by the particle aspect ratio rather than the surface geometry. As the particle aspect ratio increases, the effective particle projected area in the plane perpendicular to the flow direction increases, so that the probability of the occurrence of the particle collisions increases, leading to a reduction in particle velocity fluctuation and therefore a decrease in the stresses. Hence, a simple modification is made to the kinetic theory for granular flows to describe the stress tensors for dilute flows of non-spherical particles by incorporating a normalized effective particle projected area to account for the effect of particle collision probability. For dense granular flows, the stresses depend on both the particle aspect ratio and the surface geometry. Sharp stress increases at high solid volume fractions are observed for the glued-spheres particles with large aspect ratios due to the bumpy surfaces, which impede the flow. However, smaller stresses are obtained for the true cylindrical particles with large aspect ratios at high solid volume fractions. This trend is attributed to the combined effects of the smooth particle surfaces and the particle alignments such that the major/long axes of particles are aligned in the flow direction. In addition, the apparent friction coefficient, defined as the ratio of shear to normal stresses, is found to decrease as the particle aspect ratio increases and/or the particle surface becomes smoother at high solid volume fractions.