We propose a simple and effective boundary model in a nonequilibrium molecular dynamics(NEMD) simulation to study the out-of-equilibrium dynamics of polymer fluids. The present boundary model can effectively weaken the depletion effect and the slip effect near the boundary, and remove the unwanted heat instantly. The validity of the boundary model is checked by investigating the flow behavior of dilute polymer solution driven by an external force. Reasonable density distributions of both polymer and solvent particles, velocity profiles of the solvent and temperature profiles of the system are obtained. Furthermore, the studied polymer chain shows a cross-streaming migration towards center of the tube,which is consistent with that predicted in previous literatures. These numerical results give powerful evidences for the validity of the present boundary model. Besides, the boundary model can also be used in other flows in addition to the Poiseuille flow.
A single polymer chain in solvent confined in a slit formed by two parallel plates is studied by using molecular dynamics simulation method. The square radii of gyration and diffusion behaviors of polymers are greatly affected by the distance between the two plates, but they do not follow the same way. The chain size decays drastically with increasing h (h is the distance between two plates), until a basin occurs, and a universal h/ 0 dependence for polymer chains with different degrees of polymerization can be obtained. While, for the chain's diffusion coefficient, it decays monotonously and there is no such basin-like behavior. Furthermore, we studied the radial distribution function of confined polymer chains to explain the reason why there is a difference for the decay behaviors between dynamic properties and static properties. Besides, we also give the degree of confinement dependence of the static scaling exponent for a single polymer chain. Our work provides an efficient way to estimate the dynamics and static properties of confined polymer chains, and also helps us to understand the behavior of polymer chains under confinement.
