“…The explicit numerical scheme, which can provide a precise description of both particle–particle and particle–wall interactions, is used in DEM to investigate individual particle behavior. – Considering computational efficiency and successful experience according to the previous literature, the typical Hertz-Mindlin particle contact model is adopted in this work. Based on Newton’s laws of motion, the translational motion and rotational motion of particle i are, respectively, determined by m i normald v i normald t = m i bold-italicg + ∑ j = 1 n i false( F boldn , i j + F boldt , i j false) I i normald ω i normald t = ∑ j = 1 n i false( T normalt , i j + T normalr , i j false) where v i , ω i , and I i are the translational and angular velocities and the moment of inertia of particle i , respectively; g is the gravitational acceleration and t is time; m i is the mass of particle I ; n i is the total number of particles in contact with particle i ; F n, ij , F t, ij , T t, ij , and T r, ij are the normal contact force, the tangential contact force, the tangential torque, and the rolling torque, respectively.…”