Effects of chain resolution on the configurational and rheological predictions of dilute polymer solutions in flow fields with hydrodynamic interactions

Abstract

In a recent study, the resolution of a polymer chain model was shown to significantly affect rheological predictions from Brownian dynamics (BD) simulations [Kumar and Dalal, “Effects of chain resolution on the configurational and rheological predictions from Brownian dynamics simulations of an isolated polymer chain in flow,” J. Non-Newtonian Fluid Mech. 315, 105017 (2023)], even in the absence of hydrodynamic interactions (HI) and excluded volume. In this study, we investigate the effects of chain resolution in the presence of HI. Toward this, we perform BD simulations of a long polymer chain, with the discretization level varying from a single Kuhn step (bead–rod model) to several tens of Kuhn-steps (bead–spring model). The chain models were subjected to flow fields of uniaxial extension (purely stretching) and steady shear (equal rates of stretching and rotation). Broadly, our results indicate an amplification of the differences observed between the differently resolved bead–rod and bead–spring models, in the presence of HI. Interestingly, all rheological predictions qualitatively fall in two groups for extensional flow, with the predictions from the bead–spring model with HI being close to those of the bead–rod model without HI. This indicates significantly reduced sensitivity of coarser bead–spring models to HI, relative to the one resolved to a single Kuhn step. However, in shear flow, the bead–spring rheological predictions fall between those of the bead–rod model with and without HI, forming a third group. This is linked to the presence of stretched and coiled states in the ensemble for shear flow. HI effects are large for the coiled states and weak for the stretched states, thereby yielding predictions that are intermediate between those for no HI and dominant HI. Thus, quite surprisingly, the quality of predictions of the bead–spring models is strongly affected by the physics of the flow field, irrespective of the parameterization.