SUMMARYThis study focuses on non-coaxial flow behavior of cohesionless soil undergoing cyclic rotational shear, with a special interest in the effects of particle-scale characteristics. To this end, we perform a series of 2D discrete element simulations with various particle shapes, inter-particle coefficient of friction, initial density, and stress ratios. The validity and efficacy of the numerical model is established by systematically comparing numerical simulation results with existing laboratory testing results. Such comparison shows that the numerical simulations are capable of capturing mechanical behavior observed in laboratory testing under rotational shear. We further demonstrate and quantify a strong yet simple relationship between the deviatoric part of the normalized strain increment and the non-coaxial angle, denoted by Δε ð Þ R q and ψ, respectively. This quantitative correlation between ψ and Δε ð Þ R q is independent of applied stress ratio, initial and current void ratio, and the number of cycles applied, but dependent on the principal stress orientation and particle-scale characteristics. At the same Δε ð Þ R q , specimens with higher inter-particle friction angle or smaller particle aspect ratio show greater non-coaxial angles. A simple model ψ ¼ 45 ∘ m Δε ð Þ R q is able to fit this ψ-Δε ð Þ R q relationship well, which provides a useful relationship that can be exploited in developing constitutive models for rotational shearing.