Toward an Improved Description of Constraint Release and Contour Length Fluctuations in Tube Models for Entangled Polymer Melts Guided by Atomistic Simulations

Predicting High-Density Polyethylene Melt Rheology Using a Multimode Tube Model Derived Using Non-Equilibrium Thermodynamics

Based on the Generalized bracket, or Beris-Edwards, formalism of non-equilibrium thermodynamics, we recently proposed a new differential constitutive model for the rheological study of entangled polymer melts and solutions. It amended the shortcomings of a previous model that was too strict regarding the values of the convective constraint release parameter for the model not to violate the second law of thermodynamics, and it has been shown capable of predicting a transient stress undershoot (following the overshoot) at high shear rates. In this study, we wish to further examine this model's capability to predict the rheological response of industrial polymer systems by extending it to its multiple-mode version. The comparison with industrial rheological data (High-Density Polyethylene resins), which was based on comparison with experimental data available in (a) Small Amplitude Oscillatory shear, (b) start-up shear, and (c) start-up uniaxial elongation, was noted to be good.

Thermodynamically consistent incorporation of entanglement spatial fluctuations in the slip-link model

We evaluate the thermodynamic consistency of the anisotropic mobile slip-link model for entangled flexible polymers. The level of description is that of a single chain, whose interactions with other chains are coarse grained to discrete entanglements. The dynamics of the model consist of the motion of entanglements through space and of the chain through the entanglements, as well as the creation and destruction of entanglements, which are implemented in a mean-field way. Entanglements are modeled as discrete slip links, whose spatial positions are confined by quadratic potentials. The confinement potentials move with the macroscopic velocity field, hence the entanglements fluctuate around purely affine motion. We allow for anisotropy of these fluctuations, described by a set of shape tensors. By casting the model in the form of the general equation for the nonequilibrium reversible-irreversible coupling from nonequilibrium thermodynamics, we show that (i) since the confinement potentials contribute to the chain free energy, they must also contribute to the stress tensor, (ii) these stress contributions are of two kinds: one related to the virtual springs connecting the slip links to the centers of the confinement potentials and the other related to the shape tensors, and (iii) these two kinds of stress contributions cancel each other if the confinement potentials become anisotropic in flow, according to a lower-convected evolution of the confinement strength or, equivalently, an upper-convected evolution of the shape tensors of the entanglement spatial fluctuations. In previous publications, we have shown that this cancellation is necessary for the model to obey the stress-optical rule and the Green-Kubo relation, and simultaneously to agree with plateau modulus predictions of multichain models and simulations.

Determine Mesh Size through Monomer Mean-Square Displacement

A dynamic method to determine the main parameter of the tube theory through monomer mean-square displacement is discussed in this paper. The tube step length can be measured from the intersection of the slope-12 line and the slope-14 line in log-log plot, and the tube diameter can be obtained by recording the time at which g1 data start to leave the slope-12 regime. According to recent simulation data, the ratio of the tube step length to the tube diameter was found to be about 2 for different entangled polymer systems. Since measuring the tube diameter does not require g1 data to reach the slope-14 regime, this could be the best way to find the entanglement length from microscopic consideration.

Comparative Analysis of Different Tube Models for Linear Rheology of Monodisperse Linear Entangled Polymers

The aim of the present paper is to analyse the differences between tube-based models which are widely used for predicting the linear viscoelasticity of monodisperse linear polymers, in comparison to a large set of experimental data. The following models are examined: Milner–McLeish, Likhtman–McLeish, the Hierarchical model proposed by the group of Larson, the BoB model of Das and Read, and the TMA model proposed by the group of van Ruymbeke. This comparison allows us to highlight and discuss important questions related to the relaxation of entangled polymers, such as the importance of the contour-length fluctuations (CLF) process and how it affects the reptation mechanism, or the contribution of the constraint release (CR) process on the motion of the chains. In particular, it allows us to point out important approximations, inherent in some models, which result in an overestimation of the effect of CLF on the reptation time. On the contrary, by validating the TMA model against experimental data, we show that this effect is underestimated in TMA. Therefore, in order to obtain accurate predictions, a novel modification to the TMA model is proposed. Our current work is a continuation of earlier research (Shchetnikava et al., 2014), where a similar analysis is performed on well-defined star polymers.

Assessment of the Tumbling-Snake Model against Linear and Nonlinear Rheological Data of Bidisperse Polymer Blends

We have recently solved the tumbling-snake model for concentrated polymer solutions and entangled melts in the academic case of a monodisperse sample. Here, we extend these studies and provide the stationary solutions of the tumbling-snake model both analytically, for small shear rates, and via Brownian dynamics simulations, for a bidisperse sample over a wide range of shear rates and model parameters. We further show that the tumbling-snake model bears the necessary capacity to compare well with available linear and non-linear rheological data for bidisperse systems. This capacity is added to the already documented ability of the model to accurately predict the shear rheology of monodisperse systems.

From intermediate anisotropic to isotropic friction at large strain rates to account for viscosity thickening in polymer solutions

The steady-state extensional viscosity of dense polymeric liquids in elongational flows is known to be peculiar in the sense that for entangled polymer melts it monotonically decreases-whereas for concentrated polymer solutions it increases-with increasing strain rate beyond the inverse Rouse time. To shed light on this issue, we solve the kinetic theory model for concentrated polymer solutions and entangled melts proposed by Curtiss and Bird, also known as the tumbling-snake model, supplemented by a variable link tension coefficient that we relate to the uniaxial nematic order parameter of the polymer. As a result, the friction tensor is increasingly becoming isotropic at large strain rates as the polymer concentration decreases, and the model is seen to capture the experimentally observed behavior. Additional refinements may supplement the present model to capture very strong flows. We furthermore derive analytic expressions for small rates and the linear viscoelastic behavior. This work builds upon our earlier work on the use of the tumbling-snake model under shear and demonstrates its capacity to improve our microscopic understanding of the rheology of entangled polymer melts and concentrated polymer solutions.

Tumbling-Snake Model for Polymeric Liquids Subjected to Biaxial Elongational Flows with a Focus on Planar Elongation

We have recently solved the tumbling-snake model for concentrated polymer solutions and entangled melts in the presence of both steady-state and transient shear and uniaxial elongational flows, supplemented by a variable link tension coefficient. Here, we provide the transient and stationary solutions of the tumbling-snake model under biaxial elongation both analytically, for small and large elongation rates, and via Brownian dynamics simulations, for the case of planar elongational flow over a wide range of rates, times, and the model parameters. We show that both the steady-state and transient first planar viscosity predictions are similar to their uniaxial counterparts, in accord with recent experimental data. The second planar viscosity seems to behave in all aspects similarly to the shear viscosity, if shear rate is replaced by elongation rate.