Molecular dynamics simulation of entangled polymers in shear flow

Large scale molecular dynamics simulation is conducted for a system of entangled polymers in shear flow. The polymer consists of 100, 200 and 400 beads, which is 3-10 times larger than the number of beads between the entanglement points. The simulation reproduces many characteristic features of the rheological properties of real polymeric liquids. The steady state viscosity eta((gamma) over dot) plotted against the shear rate (gamma) over dot, approaches a power law curve eta(gamma) over dot similar to (gamma) over dot(-n) which is independent of the molecular weight with the exponent n similar or equal to 1. The second normal stress coefficient Psi(2)((gamma) over dot) is negative and its ratio to the first normal stress coefficiently Psi(1)((gamma) over dot), - Psi(2)/Psi(1) approaches to zero with the increase of the shear rate. The bond orientation is measured as a function of the position of the bond along the chain, and the profile is consistent with the recent theory of Mead et al. for the convective constraint release. (C) 2000 Elsevier Science Ltd. All rights reserved.