Shear thinning behavior of linear polymer melts under shear flow via nonequilibrium molecular dynamics

The properties of both untangled and entangled linear polymer melts under shear flow are studied by nonequilibrium molecular dynamics simulations. The results reveal that the dependence of shear viscosity η on shear rate γ, expressed by n ~ γ(-n), exhibits three distinct regimes. The first is the well-known Newtonian regime, namely, η independent of shear rate at small shear rates γ < τ0(-1) (where τ0 is the longest polymer relaxation time at equilibrium). In the non-Newtonian regime (γ > τ0(-1)) the shear dependence of viscosity exhibits a crossover at a critical shear rate γc dividing this regime into two different regimes, shear thinning regime I (ST-I) and II (ST-II), respectively. In the ST-I regime (τ0(-1) < γ < γc), the exponent n increases with increasing chain length N, while in the ST-II regime (γ > γc) a universal power law n ~ γ(-0.37) is found for considered chain lengths. Furthermore, the longer the polymer chain is, the smaller the shear viscosity for a given shear rate in the ST-II regime. The simulation also shows that a characteristic chain length, below which γc will be equal to τ0(-1), lies in the interval 30 < N < 50. For all considered chain lengths in the ST-II regime, we also find that the first and second normal stress differences N1 and N2 follow power laws of N1 ~ γ(2/3) and N2 ~ γ(0.82), respectively; the orientation resistance parameter mG follows the relation mG ~ γ(0.75) and the tumbling frequency ftb follows ftb ~ γ(0.75). These results imply that the effects of entanglement on the shear dependences of these properties may be negligible in the ST-II regime. These findings may shed some light on the nature of shear thinning in flexible linear polymer melts.