Shear Banding in Simulations of Polymer Melts

This paper compares numerical simulations of sheared polymer melts with a theoretical model coupling Rolie-Poly tube dynamics and Convected Constraint Release, finding semi-quantitative agreement that predicts shear banding occurs when the entanglement number exceeds a stiffness-dependent critical threshold, while highlighting limitations of tube models in strongly flowing, partially-disentangled liquids.

The Gist

The Big Picture: Spaghetti in a Blender

Imagine you have a giant bowl of cooked spaghetti. If you stir it gently, the noodles move together smoothly. But if you start spinning a spoon really fast, something weird happens. Instead of the whole bowl spinning at the same speed, the spaghetti might suddenly split into two distinct zones:

1. A fast zone where the noodles are whipping around wildly.
2. A slow zone where the noodles are barely moving.

This phenomenon is called Shear Banding. It's like the fluid is saying, "I can't handle this speed uniformly, so I'm going to split into a fast lane and a slow lane."

Scientists have long debated whether this happens in real-world plastics (polymer melts) or just in specific solutions. This paper uses computer simulations to figure out exactly when and why this happens, and whether our current math models can predict it.

---

The Cast of Characters

To understand the study, we need to meet the "actors" in this drama:

1. The Polymers (The Spaghetti): Long, chain-like molecules. In a melt (like melted plastic), they are all tangled up with each other, like a bowl of spaghetti that has been left to sit for hours.
2. The Entanglements (The Knots): The points where the chains get stuck on each other. The more tangled they are, the harder it is to move them.
3. The Shear (The Stirring): The force applied to make the fluid flow.
4. The "Tube" (The Hallway): In physics, we imagine each spaghetti noodle is trapped in a hallway made of its neighbors. It can wiggle forward and backward (reptation), but it can't easily jump out sideways.
5. CCR (The "Get Out of Jail Free" Card): This stands for Convected Constraint Release. Imagine that when you stir the spaghetti fast, the "hallways" (tubes) get destroyed or rearranged because the noodles are being pulled so hard. This frees the noodles up to move faster. The paper studies a specific knob, called $\beta$, that controls how fast these hallways get destroyed.

---

The Experiment: A Digital Lab

The researchers didn't use a real bowl of spaghetti. They built a virtual world inside a supercomputer (using software called LAMMPS).

• The Setup: They created thousands of digital chains (beads connected by springs).
• The Variables: They changed two main things:
  • Stiffness: How bendy are the chains? (Some are floppy like cooked noodles; others are stiff like uncooked spaghetti).
  • Length/Entanglement: How many chains are there, and how tangled are they?
• The Action: They "stirred" the virtual fluid at different speeds and watched what happened.

They used two different methods to stir:

1. The "Strict Teacher" (SLLOD): This forces the fluid to move at a perfectly uniform speed. It's like a teacher forcing every student to walk at the exact same pace. This method prevents banding, so they used it to measure how the chains react to stress.
2. The "Free-Range" (DPD): This lets the fluid move however it wants. If the fluid wants to split into fast and slow lanes, it's allowed to do so. This is where they actually saw the banding.

---

The Main Discovery: The "Tangled" Threshold

The researchers compared their computer results with a mathematical model called the DO Model (named after Dolata and Olmsted). This model tries to predict when the fluid will split into bands.

The Prediction:
The model says: "If your spaghetti is tangled enough (a high number of knots, called $Z$), and the 'hallways' don't get destroyed too quickly (a low $\beta$), then banding will happen."

The Result:
The computer simulations agreed with the model!

• Short/Few Tangles: The fluid flowed smoothly. No banding.
• Long/Many Tangles: The fluid split into fast and slow lanes. Banding occurred.

It was a "semi-quantitative" match, meaning the numbers weren't perfect, but the general trend was spot on. The model successfully predicted which digital polymers would break into bands and which wouldn't.

---

The Twist: Why Real Plastics Might Not Band

Here is the most interesting part. The model predicts that banding should happen in many real-world plastics. But in real life, we rarely see it in solid plastics (melts). Why?

The paper suggests the culprit is Stiffness.

• The Analogy: Imagine the difference between a floppy wet noodle and a stiff, dry spaghetti stick.
• The Finding: Stiffer chains (like those found in many real plastics) have a higher $\beta$. This means their "hallways" get destroyed very quickly when you stir them.
• The Consequence: When the hallways get destroyed fast, the chains can relax and move easily without getting stuck in a "traffic jam." This smooths out the flow and prevents the splitting into bands.

The authors speculate that real-world polymers (like polyethylene) are stiff enough that they "escape" the banding condition. However, very flexible polymers (like silicone rubber, or PDMS) might actually show this banding behavior if we look closely enough.

---

The "History" Problem: It Depends on How You Start

The paper also found that banding is history-dependent.

• The Analogy: Think of a crowd of people trying to leave a stadium.
  • If you start the exit slowly and gradually speed up, the crowd might suddenly surge and split into a fast lane and a slow lane.
  • But if you start the crowd running at full speed immediately, they might just keep running uniformly without splitting.

In the simulations, if they started the fluid at a high speed and then slowed it down, banding sometimes didn't happen, even when the model said it should. This suggests that in real experiments, the way you start the machine might hide the banding effect, making it hard to see.

---

Summary: What Does This Mean for Us?

1. We Can Predict It: We now have a better mathematical tool (the DO model) to predict when complex fluids will act weirdly and split into layers.
2. It's About Tangles: The more tangled the molecules, the more likely they are to band.
3. Stiffness Saves the Day: Stiffer molecules are less likely to band because they can "unstick" themselves faster.
4. Real vs. Virtual: While computer simulations show clear banding, real-world experiments are messy. Things like wall slip (fluid sliding off the container walls) and edge fractures (fluid breaking at the edges) often hide the banding, making it hard to prove it exists in real life.

The Takeaway: This paper bridges the gap between fancy math and computer simulations, helping us understand why some plastics flow smoothly while others might suddenly "jam" and split into lanes. It suggests that if you want to find a plastic that does this, look for something very flexible and very long, not something stiff and short.

Technical Summary

1. Problem Statement

The paper addresses the long-standing discrepancy between theoretical predictions and experimental observations regarding shear banding in entangled polymer melts.

• Theoretical Context: Classical tube models (e.g., Doi-Edwards) predict that under high shear rates, the relationship between shear stress ($\sigma_{xy}$) and shear rate ($\dot{\gamma}$) becomes non-monotonic. This non-monotonicity leads to flow instabilities known as shear banding, where the fluid separates into regions of high and low shear rates.
• The Conflict: While shear banding is frequently observed in polymer solutions and in molecular dynamics (MD) simulations of melts, experimental evidence for steady-state shear banding in polymer melts is scarce or controversial. Experiments often show a slight stress increase rather than a plateau, and phenomena are often confounded by wall slip and edge fracture.
• The Gap: There is a need to reconcile the clear simulation evidence of banding with the lack of experimental observation, specifically by understanding the role of Convected Constraint Release (CCR) and how polymer stiffness affects the critical conditions for banding.

2. Methodology

The authors employed a multi-faceted approach combining continuum modeling with large-scale molecular dynamics simulations.

A. Theoretical Framework (The DO Model)

• They utilized the Dolata-Olmsted (DO) model, a thermodynamically consistent extension of the Rolie-Poly (RP) model.
• Key Mechanism: The model couples the evolution of the tube conformation tensor ($A$) with the reduced number of entanglements ($\nu$). It incorporates CCR, where flow-induced stretching accelerates the release of entanglements.
• Parameters: The model is governed by:
  • $Z$: Equilibrium entanglement number.
  • $\beta$: The CCR parameter (strength of flow-induced disentanglement).
  • $\alpha$: Relaxation anisotropy parameter (controls normal stress differences).
  • $\lambda_{max}$: Maximum stretch ratio.
• Prediction: The model predicts a "phase diagram" where shear banding occurs only if $Z$ exceeds a critical value $Z_c$, which depends on $\beta$ and $\alpha$. Specifically, higher $\beta$ (stronger CCR) is predicted to inhibit banding by eliminating the non-monotonicity in the stress-strain curve.

B. Molecular Dynamics Simulations

• System: Semi-flexible Kremer-Grest (KG) bead-spring polymers simulated using LAMMPS.
• Variables:
  • Chain lengths ($N_{mon}$) from 200 to 2000 beads.
  • Bending stiffness ($k_\theta$) varied (0, 1.5, 2.0) to tune the entanglement length and stiffness.
  • Shear Protocols: Two methods were used:
    1. SLLOD: Imposes a uniform velocity profile (suppresses banding, used for measuring constitutive properties).
    2. Fix Deform (DPD Thermostat): Allows the velocity profile to evolve naturally, permitting shear banding.
• Analysis:
  • Entanglement Measurement: Used the $Z1+$ code to calculate the number of primitive path kinks ($Z_k$) to determine the disentanglement ratio $\nu = Z_k/Z_{k,eq}$.
  • Parameter Fitting: The CCR parameter $\beta$ was determined by fitting the simulation data of $\nu$ vs. Rouse Weissenberg number ($Wi_R$) to the DO model predictions.
  • Banding Observation: Velocity profiles were analyzed to detect the formation of distinct fast and slow bands.

3. Key Contributions

1. Quantification of $\beta$: The study successfully mapped the CCR parameter $\beta$ to polymer stiffness (quantified by the ratio of Kuhn length to packing length, $b_K/p$). They found $\beta$ increases with stiffness, following a power law $\beta \propto (b_K/p)^{1.6}$.
2. Validation of the DO Model: The simulations confirmed that the DO model correctly predicts the existence of a critical entanglement number ($Z_c$) above which banding occurs.
3. Identification of Discrepancies: The authors identified specific areas where the simplified DO model diverges from simulation reality, particularly regarding the dependence of $Z_c$ on $\beta$ and the breakdown of tube models at high shear rates.
4. History Dependence: The study highlighted the metastable nature of shear banding in simulations, showing that shear history (starting from rest vs. ramping down from high shear) significantly impacts whether bands form.

4. Key Results

• Agreement on Banding Existence: Simulations show clear steady-state shear banding for sufficiently entangled chains ($Z > Z_c$). The region of banding in the $(Z, \beta)$ parameter space agrees semi-quantitatively with the DO model predictions.
• The $\beta$ vs. $Z_c$ Discrepancy:
  • Model Prediction: Higher $\beta$ (stronger CCR) should increase the critical entanglement number $Z_c$, making banding harder to achieve.
  • Simulation Observation: The simulations showed a slight trend where $Z_c$ decreases as $\beta$ increases (stiffer chains with higher $\beta$ banded at lower $Z$). This contradicts the model's mechanism that CCR suppresses banding.
• Shear Rate Dependence: Banding was observed in an intermediate shear rate window. At very high shear rates, the system became uniform again (likely due to complete disentanglement/tube breakdown), and at very low rates, it was uniform.
• Normal Stress Ratio ($-N_2/N_1$): The simulations and experimental data for polymer melts suggest that the normal stress ratio is higher in melts than in solutions, implying a higher relaxation anisotropy $\alpha$ (likely $\alpha > 0.5$). However, the linear regime limit is difficult to measure, leaving $\alpha$ somewhat uncertain.
• Finite Size Effects: The authors noted that finite simulation box sizes can obscure banding, particularly when the width of the high-shear band becomes comparable to the molecular size (few monomers).

5. Significance and Implications

• Explaining the "Missing" Experimental Banding: The results suggest that physical polymers with specific microscopic structures (like United-Atom Polyethylene with dihedral angles) may have very high $\beta$ values. If the trend holds that high $\beta$ raises $Z_c$ significantly, then many experimental polymer melts may simply have $Z < Z_c$ under typical conditions, explaining why banding is rarely seen in experiments compared to simulations.
• Model Limitations: The study highlights that while the DO model captures the qualitative physics, it fails quantitatively at high shear rates where the "tube" concept breaks down due to significant disentanglement. The assumption of a single relaxation mode and the specific form of the CCR term may be insufficient for strongly non-linear flows.
• Future Directions: The work suggests that to observe shear banding in physical melts, one should look for flexible polymers with low stiffness (low $\beta$) and high entanglement density ($Z$). It also underscores the need for more precise measurements of the normal stress ratio in the linear regime to fix the $\alpha$ parameter in theoretical models.

In summary, the paper provides a rigorous bridge between microscopic simulations and continuum rheology, confirming that shear banding is a real phenomenon in polymer melts under specific conditions, while simultaneously revealing the limitations of current tube-based models in describing the complex interplay of stiffness, entanglement, and flow history.