The first normal stress difference of non-Brownian hard-sphere suspensions in the oscillatory shear flow near the liquid and crystal coexistence region

We carry out a numerical study to investigate the dynamics of non-Brownian hard-sphere suspensions near the liquid and crystal coexistence region in small to large amplitude oscillatory shear flow. The first normal stress difference (N1) and related rheological functions are carefully analyzed, focusing on the strain stiffening phenomenon, which occurs in the large strain amplitude region. Under oscillatory shear, we observe several unique behaviors of N1. A negative nonzero mean value of N1 (N1,0) is observed for the applied strain amplitudes. The change of the sign, from negative to positive, at the maximum value of N1 (N1,max) is observed at a specific point, which is not consistent with the critical strain amplitude (γ0,c) at which the modulus begins to deviate from linear viscoelasticity. The behavior of N1 in the oscillatory shear flow is different from that of N1 in steady shear f…