Brownian dynamics simulations of shear-thickening in dilute polymer solutions

A versatile model describing the shear thickening behaviour of dilute polymer solutions in high shear flows is presented. The polymer macromolecules are modelled as Hookean elastic dumbbells which deform affinely during flow. In addition, the dumbbells feel a retractive anisotropic hydrodynamic drag and an isotropic Brownian force. Furthermore, it is assumed that high shear rate increases the probability of molecules forming associations and this is described through expressions for the frequencies of association and dissociation, without explicitly accounting for finite extensibility, hydrodynamic interaction or excluded volume effects. Thus, a reversible kinetic process is incorporated into the model, which results in two diffusion equations for the associated and dissociated dumbbells. Numerical simulations predict shear thickening for specific range of parameters, which are physically meaningful and related to molecular characteristics of the polymer. A comparison against experimental data reported in the literature revealed very promising results, thus confirming the ability of this model to predict shear thickening under a wide range of conditions, for various polymer models.Nomenclature A A factor in the frequency of association - B Frequency of dissociation - B₀ Reference frequency of dissociation - c Concentration of polymer solution - c_i Concentration of singlets (i = 1) and doublets (i = 2) in the solution - c^* The overlap concentration - D_t Translation coefficient of molecule - F_i(Q) Spring force for a singlet (i = 1) and for a doublet (i = 2) - F Frequency of association - F₀ Reference frequency of association - H_i Dumbbell spring constant for a singlet (i = 1) and for a doublet (i = 2) - k Boltzman's constant - k_H Huggins constant - MW Molecular weight - MW_c Critical molecular weight for formation of entanglements - n Number density of molecules in the polymer solution - n₀ Number density of dumbbells at equilibrium - n_i Number density of singlets (i = 1) and doublets (i = 2) - Q Vector defining the size and orientation of a dumbbell - t Time - T Absolute temperature - x Degree of multimerization - W Interaction energy between the two components of a doubletGreek letters a Dimensionless anisotropy parameter - Shear rate - _i Friction coefficient of singlets (i = 1) and doublets (i = 2) - _i Intrinsic viscosity of singlets (i = 1) and doublets (i = 2) - _red Reduced viscosity of solution - _sp Specific viscosity - Viscosity of the polymer solution of concentration c - _s Viscosity of the solvent - (t) White noise - K^T Velocity gradient tensor - _Hi Time constant of a singlet (i = 1) and a doublet (i = 2) - ₁ Length scale of singlets (standard deviation of singlet lengths at equilibrium) - ₂ Length scale of doublets - T_p Stress tensor - T_xy Shear Stress (xy element of T_p) - T_pi Contributions to the stress tensor of singlets (i = 1) and doublets (i = 2) - ₀ Equilibrium configuration distribution function of Q - _i Configuration distribution function of singlets (i = 1) and doublets (i = 2)