Abstract:The efficient on-lattice Monte Carlo (MC) method is extended to simulate, for the first time, the simple shear flow of multiple macromolecular chains with self-avoiding walk (SAW) on a molecular level by introducing a pseudopotential to describe the flow field. This pseudopotential makes sense only for the potential difference associated with each local microrelaxational movement of a bead in the chain strictly defined by the four-site lattice model and bond fluctuation approach. The free-draining bead-spring model is thus investigated at low and high shear rates, and the resultant shear stress and first normal stress difference are obtained by statistics according to the sampled configurational distributions under flow. As the first paper of a corresponding series, the pseudopotential is checked in detail and confirmed by the simulation outputs for both a single SAW chain and multiple SAW chains in two dimensions. The simulated velocity profile is found to greatly satisfy the requirement of the simple shear flow unless the shear rate is unreasonably high. The material functions (apparent viscosity and first normal stress coefficient) show the shear-rate dependence at high shear rate, which is, in turn, explained by the large chain deformation on the microscopic origin. Both Newtonian regime and non-Newtonian regime are found in our on-lattice MC simulation. The nonlinear rheological behaviors at high shear rates revealed in this MC simulation agree with the experimental observations in literature for most polymers. On the other hand, the asymptotic behavior about the chain-length dependence of zero-shear viscosity can be reasonably explained by the present linear viscoelastic theory. Consequently, this paper puts forward a novel and unified simulation approach to study the chain conformation and chain dynamics under shear flow and the nonlinear viscoelasticity of polymeric fluids for both dilute and concentrated solutions or polymer melts.

Monte Carlo Simulation of Sat in Three Dimensional Lattice

Worm-type Monte Carlo Simulation of the Ashkin-Teller Model on the Triangular Lattice

Monte Carlo Simulation of SATs in 2D

Critical Exponents of Sat in Two Dimensional Square Lattice

Universal Critical Properties of the Eulerian Bond-Cubic Model

Simultaneous Analysis of Several Models in the Three-Dimensional Ising Universality Class.

XY model on a self-avoiding walk

Monte Carlo Study of the Site-Percolation Model in Two and Three Dimensions.

Monte Carlo Simulation of Self-Avoiding Lattice Chains Subject to Oscillatory Shear Flow in Polymer Solution

The Trail Problem on Finitely Ramified Fractals

The three-dimensional randomly dilute Ising model: Monte Carlo results

High-Precision Monte Carlo Study Of Several Models In The Three-Dimensional U(1) Universality Class

A Monte Carlo study for the universality in short-time critical dynamics

Monte Carlo Study Of The Triangular Lattice Gas With First- And Second-Neighbor Exclusions

Monte Carlo simulation of electrostatic interactions in inhomogeneous dielectric media: Correct sampling for the local lattice simulation algorithm

Monte Carlo Simulation of Self-Avoiding Lattice Chains Subject to Simple Shear Flow. I. Model and Simulation Algorithm

Simplicial Lattice Study of the 2d Ising CFT

Monte Carlo simulation of self-avoiding lattice chains subject to simple shear flow Part II. Three-dimensional results and comparison with experiments

Logarithmic Finite-Size Scaling of the Self-Avoiding Walk at Four Dimensions.

Cluster Monte Carlo Simulation of the Transverse Ising Model

Critical structure factor in Ising systems