Beyond binary scission: a generalized three-species cascade breakage model for wormlike micellar solutions

Résumé Technique : Un modèle de rupture en cascade à trois espèces généralisé pour les solutions de micelles filamenteuses

Énoncé du problème
Les solutions de micelles filamenteuses (WLM) présentent des comportements rhéologiques complexes, notamment la viscoélasticité, le amincissement par cisaillement et le bandement de cisaillement (shear banding), pilotés par la rupture et la recombinaison continues des réseaux de micelles auto-assemblées. Les cadres théoriques existants, particulièrement les modèles à deux espèces largement utilisés (par exemple, le modèle de Vasquez-Cook-McKinley ou VCM), reposent sur une représentation binaire grossière de la population micellaire (distinguant typiquement uniquement les longues chaînes enchevêtrées des fragments courts). Bien que réussis pour capturer certains traits macroscopiques, ces modèles possèdent une limitation fondamentale : ils ne parviennent pas à résoudre explicitement la rupture en cascade multi-échelle et les régimes de relaxation continus inhérents aux réseaux micellaires réels. Cette approximation binaire filtre les états structurels intermédiaires, conduisant à une incapacité à décrire précisément les spectres de relaxation larges, particulièrement le redressement à haute fréquence des modules dynamiques observé dans les expériences, et limite l'interprétation physique des réponses transitoires et des mécanismes de bandement de cisaillement.

Méthodologie
Pour remédier à ces limitations, les auteurs développent un modèle de rupture en cascade à trois espèces qui introduit un état micellaire intermédiaire dans le cadre cinétique. Le modèle suppose que le réseau micellaire est composé de trois états structurels discrets :

1. Espèce A : Un réseau de gel (primaire, hautement enchevêtré).
2. Espèce B : Des chaînes longues (état intermédiaire).
3. Espèce C : Des chaînes courtes (fragments désenchevêtrés).

La formulation mathématique est basée sur une approche de théorie cinétique utilisant des distributions continues dans l'espace des phases. Les auteurs dérivent des équations de Smoluchowski couplées pour les fonctions de densité de probabilité des trois espèces, incorporant la convection d'étirement, la diffusion et les termes de réaction. Les principales caractéristiques méthodologiques incluent :

• Cinétique hiérarchique : Le modèle définit un chemin de cascade où le réseau de gel (A) se rompt en chaînes longues (B), qui se rompent ensuite en chaînes courtes (C). Des voies de recombinaison sont également définies pour assurer la conservation de la masse.
• Scission dépendante de l'écoulement : Les taux de rupture du réseau de gel et des chaînes longues sont couplés au flux macroscopique local via un mécanisme de convection de brins à extension partielle (PEC). Cela crée un mécanisme de rupture dépendant de la contrainte où la déformation accélère la scission.
• Équations de comportement : Le modèle dérive des équations d'évolution adimensionnelles pour les densités numériques et les tenseurs de conformation, couplées à l'équation de quantité de mouvement de Cauchy macroscopique. La contrainte totale est la superposition des contraintes élastiques de toutes les trois espèces.
• Implémentation numérique : Les auteurs analysent le modèle sous des conditions d'écoulement homogène (pour étudier la viscoélasticité linéaire et la réponse de cisaillement stationnaire) et d'écoulement de Couette cylindrique inhomogène (pour étudier le bandement de cisaillement, le démarrage transitoire et la relaxation).

Contributions clés et résultats

1. Réponse viscoélastique à trois modes :
   Sous écoulement homogène, le modèle produit analytiquement une structure de Maxwell à trois modes discrets. Contra�iquement aux modèles à deux espèces, ce cadre capture avec succès le redressement à haute fréquence des modules de stockage ($G'$) et de perte ($G''$) dynamiques. Les auteurs démontrent que le mode de relaxation supplémentaire fourni par l'espèce intermédiaire (B) et les chaînes courtes (C) améliore significativement l'ajustement aux courbes de maîtrise expérimentales des solutions CPyCl/NaSal, particulièrement dans le régime de haute fréquence.

2. Courbes d'écoulement non monotones et bandement de cisaillement :
   Le modèle prédit une courbe constitutive non monotone présentant une région de pente strictement négative, un trait caractéristique des instabilités constitutives. Dans les simulations d'écoulement de Couette inhomogène, cette instabilité conduit à la formation spontanée de bandes de cisaillement coexistantes (régions de taux de cisaillement élevé et faible). Le modèle à trois espèces révèle que l'espèce intermédiaire (B) forme une couche de transition distincte à l'interface des bandes, reliant la bande de faible cisaillement dominée par le réseau de gel et la bande de haut cisaillement dominée par les chaînes courtes. Cela fournit une description microstructurale plus résolue de l'interface de bande que les modèles binaires le permettent.

3. Dynamique transitoire et hystérésis :
   Le modèle capture des comportements transitoires complexes :

• Dépassement de contrainte (Stress overshoot) : Lors du démarrage, la compétition entre le chargement élastique et la rupture en cascade accélérée conduit à des dépassements de contrainte prononcés.
• Hystérésis d'écoulement : Le modèle présente une hystérésis, où l'état stationaire final dépend de l'historique de l'écoulement (démarrage à partir du repos vs descente par pali de cisaillement élevé).
• Relaxation multi-étapes : Lors de la cessation du cisaillement, le modèle prédit un processus de relaxation multi-étapes. Au lieu d'une décroissance exponentielle unique, la contrainte se relaxe par étapes séquentielles : libération rapide des chaînes courtes (C), ajustement intermédiaire par les chaînes longues (B), et récupération lente du réseau de gel (A). Cela résout le « temps de récupération prolongé » souvent observé dans les expériences mais difficile à capturer avec des modèles réduits.

4. Cohérence mécanique continue :
   Le modèle satisfait la relation de Lodge-Meissner ($N_1/\sigma_{r\theta} = \gamma_0$) lors de grands paliers de déformation imposée, confirmant sa cohérence tensorielle et sa validité physique sous déformation finie.

Signification et affirmations
L'article affirme que le cadre proposé à trois espèces comble le fossé entre les modèles de réseau phénoménologiques discrets et les théories cinétiques continues (telles que la théorie originale de reptation-réaction de Cates). En résolvant explicitement la cinétique de cascade régissant les états structurels intermédiaires, le modèle fournit un cadre physiquement interprétable pour le bandement de cisaillement des micelles filamenteuses.

Les auteurs soutiennent que cette approche capture la connexion intrinsèque entre :

• L'évolution microstructurale en cascade (A $\to$ B $\to$ C).
• La dynamique de relaxation large (trois échelles de temps distinctes).
• La localisation de l'écoulement macroscopique (bandement de cisaillement).

La signification réside dans la capacité du modèle à maintenir une tractabilité computationnelle macroscopique tout en récupérant les modes de relaxation intermédiaires manquants. Cela offre une base physique plus complète pour interpréter la richesse des dynamiques transitoires et stationnaires des fluides micellaires complexes, pouvant potentiellement aider à l'interprétation de futures expériences microfluidiques et d'applications industrielles impliquant des WLM. L'article ne propose pas de nouveaux protocoles expérimentaux mais offre plutôt un outil théorique pour mieux interpréter les phénomènes rhéologiques existants.