Multi-scale modelling of dilute viscoelastic liquids: Atomistic to mesoscale mapping of polymer solutions

Abstract

A framework predicting the rheological (storage and loss moduli, first normal stress coefficient, and relaxation time) and transport (viscosity, diffusivity) properties of non-Newtonian dilute polymer solutions at mesoscales (e.g. ∼ns to μs) from the atomistic-scale molecular behaviour is presented. More specifically, the rheological behaviour differences of OCP and PMA polymer solutions in PAO-2 oil are simulated using both atomistic molecular dynamics (MD) and many-body dissipative particle dynamics (mDPD) within a temperature range of 313–373 K. The simulation methodology described is able to distinguish itself from the standard DPD model by accurately reproducing the shear-thinning with high sensitivity, as seen in the atomistic MD simulations at high shear rates (e.g. 108−1013 s−1). It is shown that the model is well-suited to compute properties such as first normal stress differences and relaxation times that are difficult to estimate at atomistic scales due to the low signal-to-noise ratio. Moreover, the Schmidt numbers (>103) are predicted with high accuracy when compared with the values from atomistic-scale simulations. The proposed model is able to predict relaxation times of dilute polymer mixtures that are difficult to be obtained using state-of-the-art rheometers. Finally, it is found that the terminal relaxation time of PMA polymer chain does not vary monotonically as a function of temperature, unlike in the case of OCP; this is significant for describing viscoelastic behaviour at macroscales where satisfactory constitutive equations are not available.