Abstract:Considering the cell subdivision method and the modified velocityVerlet algorithm, the implementation of LeesEdwards boundary conditions for DPD method was investigated. The simulated results show that the profiles of velocities, densities, temperatures, pressures and stresses agree well with the expected. Furthermore, when dissipative coefficient γ was increased to 100, the velocity profile in the system was homogeneous, indicating that the LE boundary conditions still valid for higher dissipation rate. Viewpoints different from Ref [21] are proposed.