Shear Banding in Simulations of Polymer Melts

Les simulations numériques de polymères en écoulement de cisaillement montrent un accord semi-quantitatif avec un modèle couplant la dynamique des tubes et la libération de contraintes convectées (CCR), révélant que le cisaillement hétérogène (shear banding) survient lorsque le nombre d'enchevêtrements dépasse une valeur critique dépendante de la rigidité de flexion.

L'essentiel

🧶 Le Grand Défi des Spaghetti Énervés

Imaginez un bol rempli de spaghettis cuits, très longs et très emmêlés. C'est ce qu'on appelle un polymère fondu (comme du plastique fondu ou du caoutchouc). Quand vous essayez de les mélanger doucement avec une cuillère, tout se déplace ensemble. C'est fluide et régulier.

Mais que se passe-t-il si vous fouettez ce mélange à toute vitesse ? C'est là que la magie (et le problème) opère.

🚦 Le Phénomène de "Banding" (Les Bandes de Trafic)

Les chercheurs de cet article ont découvert quelque chose de fascinant : si vous fouettez assez vite, le mélange ne reste pas uniforme. Il se sépare en deux zones distinctes :

1. Une zone où les spaghettis glissent très vite (comme une autoroute libre).
2. Une zone où ils bougent très lentement (comme un embouteillage).

C'est ce qu'on appelle le "shear banding" (ou bandage par cisaillement). C'est comme si, sur une autoroute, une moitié des voitures roulait à 200 km/h et l'autre moitié à 5 km/h, sans que les deux ne se mélangent.

🔍 Le Problème : Pourquoi ne le voit-on pas toujours ?

Dans la vraie vie (en laboratoire avec de vrais plastiques), il est très difficile de voir ces bandes. Les expériences sont souvent gâchées par des effets de bord (comme le plastique qui glisse sur les parois du récipient) ou par des fractures.

C'est pourquoi les auteurs ont utilisé des simulations informatiques. Ils ont créé des "spaghettis virtuels" dans un ordinateur pour observer ce qui se passe sans les problèmes des expériences réelles. Ils ont vu clairement ces bandes apparaître.

🧠 La Théorie : Le Modèle "Tube" et le "Détachement"

Pour prédire quand ces bandes vont apparaître, les scientifiques utilisent une théorie appelée le modèle du tube.

• L'idée : Imaginez que chaque spaghetti est coincé dans un tube imaginaire fait par ses voisins. Il ne peut bouger que dans ce tube, comme un ver dans un tuyau.
• Le problème : Quand on fouette trop fort, les spaghettis s'étirent et se détachent de leurs voisins. Ils sortent de leur "tube".

Les chercheurs ont comparé leurs simulations à un modèle mathématique (le modèle DO) qui prend en compte ce détachement. Le modèle dit : "Si vous avez assez de spaghettis emmêlés (un nombre élevé de nœuds, noté Z) et que vous les fouettez à la bonne vitesse, des bandes vont se former."

🎯 La Découverte Principale : La Rigidité est la Clé

C'est ici que ça devient intéressant. Le modèle prédit que plus les polymères sont rigides (raides comme des baguettes) et plus ils sont emmêlés, plus ils devraient former de bandes.

Mais les simulations ont révélé une surprise :

• Les polymères flexibles (comme des spaghettis mous) forment des bandes assez facilement.
• Les polymères rigides (comme des baguettes en bois) sont plus difficiles à faire bander, même s'ils sont très emmêlés.

Pourquoi ? Parce que les polymères rigides ont une propriété spéciale appelée CCR (Relâchement Contraint Convectif). C'est un peu comme si, quand on les étire, ils se "détachent" si vite de leurs voisins que le système reste stable et ne se sépare pas en bandes. C'est comme si le spaghetti rigide trouvait un moyen de glisser hors de son tube avant que le chaos ne s'installe.

🌍 Pourquoi est-ce important pour nous ?

Cette recherche est cruciale pour l'industrie du plastique.

• Si vous fabriquez des bouteilles ou des fibres, vous ne voulez pas que le plastique se sépare en bandes pendant la fabrication, car cela créerait des défauts et des faiblesses.
• En comprenant exactement quels types de plastiques (flexibles ou rigides) et à quelles vitesses de production ces bandes risquent de se former, les ingénieurs peuvent ajuster leurs machines pour éviter le désastre.

En Résumé

Les chercheurs ont utilisé un ordinateur pour simuler des millions de spaghettis virtuels. Ils ont confirmé qu'ils peuvent se séparer en zones rapides et lentes (des bandes) quand on les fouette trop fort. Ils ont aussi découvert que la rigidité de ces spaghettis joue un rôle de "frein" : plus ils sont raides, plus ils résistent à cette séparation, ce qui explique pourquoi on ne voit pas toujours ce phénomène dans les plastiques du quotidien.

C'est une victoire pour la science : on passe de "ça semble bizarre" à "voici exactement pourquoi ça arrive et comment l'éviter".

Résumé technique

1. Problématique et Contexte

Le phénomène de banding de cisaillement (shear banding) dans les liquides polymères fortement enchevêtrés est une instabilité d'écoulement où le fluide se sépare en bandes de vitesses de cisaillement différentes sous une contrainte uniforme. Bien que ce phénomène soit bien documenté expérimentalement pour certaines solutions polymères, son existence dans les mélanges polymères purs (melt) reste débattue et difficile à observer en raison d'artefacts expérimentaux tels que le glissement aux parois et la fracture des bords.

Les modèles théoriques classiques, comme le modèle de Doi-Edwards (DE), prédisent une relation non monotone entre la contrainte de cisaillement et le taux de cisaillement, ce qui favorise le banding. Cependant, l'introduction du Relâchement de Contrainte Convecté (CCR - Convected Constraint Release) dans des modèles avancés (comme le modèle Rolie-Poly étendu) tend à supprimer cette non-monotonie pour des polymères très enchevêtrés, rendant le banding improbable.

L'objectif de cet article est de comparer les prédictions d'un modèle théorique récent (le modèle DO, couplant la dynamique des tubes à la dynamique d'enchevêtrement via le CCR) avec des simulations de dynamique moléculaire non équilibrée (NEMD) de mélanges polymères, afin de déterminer les conditions physiques sous lesquelles le banding se produit réellement dans les mélanges.

2. Méthodologie

A. Modélisation Théorique (Modèle DO)
Les auteurs utilisent le modèle développé par Dolata et Olmsted (DO), qui couple l'évolution du tenseur de conformation du tube local ($A$) au nombre réduit de nœuds d'enchevêtrement ($\nu$).

• Le modèle intègre le mécanisme CCR, contrôlé par un paramètre $\beta$.
• Il inclut l'anisotropie de relaxation ($\alpha$) et la finitude de l'extensibilité des chaînes.
• La prédiction du banding dépend de la non-monotonie de la courbe constitutive (contrainte vs taux de cisaillement). Le modèle prédit que le banding se produit si le nombre d'enchevêtrements à l'équilibre $Z$ dépasse une valeur critique $Z_c(\beta, \alpha)$.

B. Simulations Numériques
Les auteurs ont réalisé des simulations de dynamique moléculaire (MD) non équilibrée de polymères semi-flexibles de type Kremer-Grest (KG) (modèle bille-ressort) en utilisant le logiciel LAMMPS.

• Paramètres : Ils ont varié la rigidité de flexion ($k_\theta$) et la longueur des chaînes ($N_{mon}$) pour modifier le nombre d'enchevêtrements $Z$.
• Protocoles de cisaillement : Deux méthodes ont été utilisées :
  1. SLLOD : Impose un profil de vitesse uniforme (empêche le banding, utile pour mesurer les propriétés rhéologiques globales et le paramètre $\beta$).
  2. Déformation de boîte (fix deform) avec thermostat DPD : Permet la formation de profils de vitesse non uniformes, donc l'observation du banding.
• Analyse : Le nombre d'enchevêtrements $Z_k$ a été mesuré en temps réel via l'analyse du chemin primitif (code Z1+). Le paramètre $\beta$ a été ajusté en comparant la décroissance de $Z_k$ avec le nombre de Weissenberg de Rouse ($Wi_R$) aux prédictions du modèle DO.

3. Contributions Clés

1. Détermination du paramètre CCR ($\beta$) : L'article établit une relation quantitative entre la rigidité de la chaîne ($k_\theta$) et le paramètre $\beta$. Il confirme que $\beta$ augmente avec la rigidité de la chaîne (rapport longueur de Kuhn / longueur d'empilement $b_K/p$), bien que la loi de puissance trouvée diffère légèrement de travaux antérieurs.
2. Cartographie du banding : Les auteurs construisent un diagramme de phase reliant $Z$, $\beta$ et $\alpha$ pour prédire la zone de banding, et la comparent directement aux résultats des simulations.
3. Analyse de l'histoire de cisaillement : L'étude met en évidence la forte dépendance de l'état stationnaire vis-à-vis de l'histoire de déformation (métastabilité), expliquant pourquoi certaines simulations ne montrent pas de banding même dans des conditions théoriquement favorables.

4. Résultats Principaux

• Accord Semi-Quantitatif : Les simulations montrent un accord semi-quantitatif avec le modèle DO. Pour des valeurs de $\beta$ faibles à modérées (typiques des chaînes KG flexibles et semi-flexibles, $\beta \approx 0.03 - 0.3$), le modèle prédit correctement que le banding se produit pour des chaînes avec un nombre d'enchevêtrements $Z \gtrsim 10$.
• Déviation sur la Rigidité : Une divergence notable apparaît concernant l'effet de la rigidité. Le modèle prédit que l'augmentation de $\beta$ (due à une plus grande rigidité) devrait augmenter le seuil critique $Z_c$, rendant le banding plus difficile à observer. Cependant, les simulations montrent une tendance inverse ou neutre : les systèmes plus rigides ($k_\theta=2.0$) présentent du banding à des $Z$ plus faibles que prévu par le modèle, ou le seuil $Z_c$ ne suit pas la tendance théorique attendue.
• Rôle de l'Histoire et de la Taille Finie :
  • Le banding est fortement dépendant de l'histoire de cisaillement (métastabilité). Un système démarrant d'un état de repos peut bander, tandis qu'un système passant d'un taux de cisaillement élevé (non-banding) à un taux plus faible peut rester homogène.
  • Les effets de taille finie dans les simulations (boîtes de taille limitée) peuvent masquer le banding, en particulier lorsque la largeur de la bande rapide devient inférieure à quelques monomères.
• Contraintes Normales : L'analyse des rapports de contraintes normales ($-N_2/N_1$) suggère que le paramètre d'anisotropie $\alpha$ pourrait être universel pour différents types de polymères, mais les données à faible taux de cisaillement sont insuffisantes pour le fixer avec précision.

5. Signification et Conclusion

Cet article apporte des éclairages cruciaux sur la nature du banding de cisaillement dans les mélanges polymères :

• Validation des modèles : Le modèle DO, bien que simplifié (version mono-mode), capture correctement la physique globale du banding dans les mélanges simulés, validant l'approche couplant dynamique des tubes et CCR.
• Explication de l'absence de banding expérimental : Les auteurs suggèrent que les polymères physiques réels (comme le polyéthylène) possèdent des mécanismes microscopiques (rotations dièdres) qui pourraient conduire à des valeurs de $\beta$ beaucoup plus élevées que celles des modèles de chaînes simples. Un $\beta$ élevé augmenterait le seuil critique $Z_c$ au-delà des nombres d'enchevêtrements typiques des polymères expérimentaux, expliquant pourquoi le banding est rarement observé dans les mélanges réels, contrairement aux solutions ou aux modèles de simulation simplifiés.
• Limites des modèles de tube : Les écarts observés à haut taux de cisaillement soulignent les limites des modèles de tube (comme DO ou Rolie-Poly) lorsque les enchevêtrements sont fortement détruits (désenchevêtrement partiel) et que des modes dynamiques plus rapides entrent en jeu.

En résumé, cette étude démontre que le banding dans les mélanges polymères est un phénomène réel et simulable, mais qu'il est extrêmement sensible aux paramètres microscopiques spécifiques au polymère (rigidité, mécanismes de relâchement) et aux conditions historiques de l'écoulement, ce qui complique son observation expérimentale.