Entanglements in quiescent and sheared polymer melts

Effect of shear flow on the transverse thermal conductivity of polymer melts

The effect of shear flows on the thermal conductivity of polymer melts is investigated using a reversed nonequilibrium molecular-dynamics (RNEMD) method. We extended the original RNEMD method to simultaneously produce spatial gradients of temperature and flow velocity in a single direction. This method enables accurate measurement of the thermal conductivity in the direction transverse to shear flow. The Weissenberg number defined with the shear rate and the relaxation time of the polymer conformation can uniformly differentiate the occurrence of shear rate dependence of the thermal conductivity across different chain lengths. The stress-thermal rule (STR) (i.e., the linear relationship between anisotropic parts of the stress tensor and the thermal conductivity tensor) holds for entangled polymer melts even under shear flows but not for unentangled polymer melts. Furthermore, once entanglements form in polymer chains, the stress-thermal coefficient in the STR remains independent of the polymer chain length. These observations align with the theoretical foundation of the STR, which focuses on energy transmission along the network structure of entangled polymer chains [B. van den Brule, Rheol. Acta 28, 257 (1989)0035-451110.1007/BF01329335]. However, under driven shear flows, the stress-thermal coefficient is notably smaller than that measured in the literature for a quasiquiescent state without external forces. Although the mechanism of the STR in shear flows has yet to be fully elucidated, our study confirmed the validity of the STR in shear flows. This allows us to use the STR as a constitutive equation for computational thermofluid dynamics of polymer melts, thus having broad engineering applications.

Molecular Dynamics Simulation of Entangled Melts at High Rates: Identifying Entanglement Lockup Mechanism Leading to True Strain Hardening

In the present work, molecular dynamics simulations are carried out based on the bead-spring model to indicate how the entanglement lockup manifests in the late stage of fast Rouse-Weissnberg number (Wi_R >>1) uniaxial melt stretching of entangled polymer melts. At high strains, distinct features show up to reveal the emergence of an increasingly tightened entanglement network. Chain tension can build up, peaking at the middle of the chain, to a level for chain scission, through accumulated interchain interactions, as if there is a tug-of-war ongoing for each load-bearing chain. Thanks to the interchain uncrossability, network junctions form by the pairing of two or more hairpins. It is hypothesized that the interchain entanglement at junctions can lockup through prevailing twist-like interchain couplings as long as Wi_R > 9. In this limit, a significant fraction of chains act like cyclic chains to form a network held by interchain uncrossability, and appreciable chain tension emerges.

Fingerprinting Brownian Motions of Polymers under Flow

We present a quantitative approach to the self-dynamics of polymers under steady flow by employing a set of complementary reference frames and extending the spherical harmonic expansion technique to dynamic density correlations. Application of this method to nonequilibrium molecular dynamics simulations of polymer melts reveals a number of universal features. For both unentangled and entangled melts, the center-of-mass motions in the flow frame are described by superdiffusive, anisotropic Gaussian distributions, whereas the isotropic component of monomer self-dynamics in the center-of-mass frame is strongly suppressed. Spatial correlation analysis shows that the heterogeneity of monomer self-dynamics increases significantly under flow.

Decoding polymer self-dynamics using a two-step approach

The self-correlation function and corresponding self-intermediate scattering function in Fourier space are important quantities for describing the molecular motions of liquids. This work draws attention to a largely overlooked issue concerning the analysis of these space-time density-density correlation functions of polymers. We show that the interpretation of non-Gaussian behavior of polymers is generally complicated by intrachain averaging of distinct self-dynamics of different segments. By the very nature of the mathematics involved, the averaging process not only conceals critical dynamical information, but also contributes to the observed non-Gaussian dynamics. To fully expose this issue and provide a thorough benchmark of polymer self-dynamics, we perform analyses of coarse-grained molecular dynamics simulations of linear and ring polymer melts as well as several theoretical models using a "two-step" approach, where interchain and intrachain averagings of segmental self-dynamics are separated. While past investigations primarily focused on the average behavior, our results indicate that a more nuanced approach to polymer self-dynamics is clearly required.

Spatial correlations of entangled polymer dynamics

The spatial correlations of entangled polymer dynamics are examined by molecular dynamics simulations and neutron spin-echo spectroscopy. Due to the soft nature of topological constraints, the initial spatial decays of intermediate scattering functions of entangled chains are, to the first approximation, surprisingly similar to those of an unentangled system in the functional forms. However, entanglements reveal themselves as a long tail in the reciprocal-space correlations, implying a weak but persistent dynamic localization in real space. Comparison with a number of existing theoretical models of entangled polymers suggests that they cannot fully describe the spatial correlations revealed by simulations and experiments. In particular, the strict one-dimensional diffusion idea of the original tube model is shown to be flawed. The dynamic spatial correlation analysis demonstrated in this work provides a useful tool for interrogating the dynamics of entangled polymers. Lastly, the failure of the investigated models to even qualitatively predict the spatial correlations of collective single-chain density fluctuations points to a possible critical role of incompressibility in polymer melt dynamics.

Synchronized Molecular-Dynamics Simulation of the Thermal Lubrication of an Entangled Polymeric Liquid

The thermal lubrication of an entangled polymeric liquid in wall-driven shear flows between parallel plates is investigated by using a multiscale hybrid method, coupling molecular dynamics and hydrodynamics (i.e., the synchronized molecular dynamics method). The temperature of the polymeric liquid rapidly increases due to viscous heating once the drive force exceeds a certain threshold value, and the rheological properties drastically change at around the critical drive force. In the weak viscous-heating regime, the conformation of polymer chains is dominated by the flow field so that the polymers are more elongated as the drive force increases. However, in the large viscous-heating regime, the conformation dynamics is dominated by the thermal agitation of polymer chains so that the conformation of polymers recovers more uniform and random structures as the drive force increases, even though the local shear flows are further enhanced. Remarkably, this counter-intuitive transitional behavior gives an interesting re-entrant transition in the stress–optical relation, where the linear stress–optical relation approximately holds even though each of the macroscopic quantities behaves nonlinearly. Furthermore, the shear thickening behavior is also observed in the large viscous-heating regime—this was not observed in a series of previous studies on an unentangled polymer fluid. This qualitative difference of the thermo-rheological property between the entangled and unentangled polymer fluids gives completely different velocity profiles in the thermal lubrication system.

Simulating the Flow of Entangled Polymers

To optimize automation for polymer processing, attempts have been made to simulate the flow of entangled polymers. In industry, fluid dynamics simulations with phenomenological constitutive equations have been practically established. However, to account for molecular characteristics, a method to obtain the constitutive relationship from the molecular structure is required. Molecular dynamics simulations with atomic description are not practical for this purpose; accordingly, coarse-grained models with reduced degrees of freedom have been developed. Although the modeling of entanglement is still a challenge, mesoscopic models with a priori settings to reproduce entangled polymer dynamics, such as tube models, have achieved remarkable success. To use the mesoscopic models as staging posts between atomistic and fluid dynamics simulations, studies have been undertaken to establish links from the coarse-grained model to the atomistic and macroscopic simulations. Consequently, integrated simulations from materials chemistry to predict the macroscopic flow in polymer processing are forthcoming.

Computational study of the propagation of the longitudinal velocity in a polymer melt contained within a cylinder using a scale-bridging method

The "constitutive equation"-free scale-bridging method connecting nonequilibrium molecular dynamics and continuum fluid mechanics, that had hitherto been applied only to a parallel-plates geometry, is extended to study the flow of a polymer melt in a cylindrical pipe subject to a velocity in the direction parallel to the cylinder's axis. The system, initially at rest, is given a velocity at the cylinder's surface, and the evolution of the velocity profile within the fluid is studied, along with the time taken for the velocity to propagate toward the cylinder's axis. The said time of propagation is found to increase with the boundary velocity-a fact in contrast with the case of a Newtonian fluid for which the time of propagation is expected to be independent of the boundary velocity. For a fixed value of the boundary velocity, the propagation time is found to increase with the cylinder radius according to a power law with an exponent that is smaller than the corresponding exponent for a Newtonian fluid. For the lower values of the boundary velocity and the lower values of the radius studied, a velocity overshoot is observed at the cylinder's axis-a manifestation of elastic behavior of the fluid.

Dynamics of relaxation of entangled polymers in shear flow

The application of shear flow to entangled polymer melts can strongly modify its rheological and physicochemical behaviors, giving rise to an acceleration of several chemical processes such as diffusion-controlled reactions. In the present work, we investigate the modification of conformational and diffusive properties of an entangled polymer in shear flow by numerical methods. The flow affects both the conformational and diffusive properties of the system, giving rise to a quasinematic ordering of the macromolecules which take prolate spheroid shape with the main axis aligned to the shear direction. The shear flow is found to accelerate the overall diffusion of the chains in all directions at times longer than the polymer relaxation time. The polymer chains display a quite peculiar displacement behavior in direction parallel to the flow. At the same conditions, the linear relation between the diffusion constant in direction perpendicular to the flow and the inverse of the relaxation time, usually adopted in equilibrium regimes, is shown to hold even in the presence of flow.