Continuum modeling of Soft Glassy Materials under shear

Cet article présente une introduction pédagogique à un modèle continu basé sur une approche de fluidité spatialement résolue, qui rend compte des effets non locaux et capture quantitativement les phénomènes complexes tels que le dépassement de contrainte et le cisaillement hétérogène lors de l'écoulement de matériaux vitreux mous.

L'essentiel

🍯 Les Matériaux "Glaceux" : Quand le Solide devient Liquide

Imaginez une foule très dense de gens dans une salle de concert. Tant que la musique est douce, tout le monde reste immobile, serré les uns contre les autres : c'est un solide. Mais si la musique commence à devenir très forte (une contrainte mécanique), les gens commencent à bouger, à se pousser, et soudain, la foule se met à couler comme un liquide.

C'est exactement ce qui se passe avec les Matériaux Glaceux Mous (ou Soft Glassy Materials). Ce sont des substances comme le dentifrice, la mayonnaise, les gels de carbopol (utilisés dans les cosmétiques) ou même certaines encres d'impression 3D. À l'arrêt, ils sont solides et tiennent leur forme. Mais si on les pousse trop fort ou trop vite, ils se transforment en liquide.

Les scientifiques Roberto Benzi et ses collègues ont écrit un article pour expliquer comment cette transformation se produit, surtout quand on les étire ou qu'on les cisaille (qu'on les fait glisser entre deux plaques).

🎢 Le "Saut de Stress" (L'effet de ressort)

Quand on commence à faire bouger ces matériaux, il se passe quelque chose de curieux, appelé le "stress overshoot" (ou dépassement de contrainte).

• L'analogie du ressort : Imaginez que vous tirez sur un élastique très épais. Au début, il résiste fort (c'est le comportement élastique). La tension monte, monte, monte... jusqu'à atteindre un pic maximum. À ce moment précis, l'élastique "casse" intérieurement et se détend soudainement pour couler plus facilement.
• Dans le matériau : Le matériau résiste d'abord, atteint un pic de tension, puis se relâche pour commencer à couler. Les chercheurs veulent comprendre : Plus on tire vite, plus ce pic est haut ? Et combien de temps faut-il pour que tout le matériau se mette à couler ?

🌊 Le Modèle de la "Fluidité" : Une Carte de Trafic

Pour prédire ce comportement, les auteurs ont créé un modèle mathématique qu'ils appellent le modèle de "fluidité".

Au lieu de regarder chaque petite particule individuellement (ce qui serait trop compliqué, comme compter chaque grain de sable sur une plage), ils regardent le matériau comme un tout, mais avec une carte spéciale.

• La fluidité (f) : C'est comme un indicateur de "déblocage".
  • Si f = 0, c'est un embouteillage total (le matériau est solide, rien ne bouge).
  • Si f > 0, c'est que la circulation reprend (le matériau est devenu liquide).

L'idée clé : La contagion (Coopérativité)
Dans un vrai embouteillage, si une voiture avance, elle libère de l'espace pour la suivante. De même, dans ces matériaux, si une petite zone commence à couler, elle aide ses voisines à couler aussi. Les chercheurs ont ajouté une "distance de contagion" (appelée longueur de coopérativité) dans leur équation. Cela signifie qu'un mouvement local influence ses voisins, créant des vagues de fluidité.

🚧 Les Bandes de Cisaillement : Une Autoroute à deux vitesses

Lorsqu'on commence à faire couler le matériau, il ne devient pas liquide partout en même temps. C'est comme si, sur une autoroute, une seule voie se débloquait soudainement alors que les autres restent bloquées.

• La bande de cisaillement (Shear Band) : Une fine couche de matériau près de la paroi qui bouge (le liquide), tandis que le reste du matériau reste figé (le solide).
• L'évolution : Cette bande de liquide grandit petit à petit, comme une tache d'huile, jusqu'à ce qu'elle finisse par envahir tout le matériau. C'est ce qu'on appelle la fluidisation.

Le modèle des auteurs montre très bien comment cette "tache de liquide" grandit et comment la vitesse de croissance dépend de la force appliquée.

🛞 Le Secret du Glissement : L'effet "Elasto-Hydrodynamique"

Il y a un détail crucial que les chercheurs ont ajouté pour rendre leur modèle encore plus réaliste : le frottement contre les parois.

• L'analogie du patin à glace : Quand vous glissez sur la glace, il y a une fine couche d'eau entre le patin et la glace qui vous permet de glisser. Dans les matériaux mous (comme les microgels), quand on les écrase contre une paroi lisse, il se forme une minuscule couche de liquide entre les particules et la paroi.
• L'effet EHD : Cette couche lubrifie le mouvement. Si la paroi est lisse, le matériau glisse beaucoup plus facilement que si la paroi est rugueuse.
• La découverte : En incluant cet effet dans leur équation, les chercheurs ont pu expliquer pourquoi, sur des surfaces lisses, le pic de tension (le "saut") se comporte différemment. C'est comme si la "règle du jeu" changeait selon la texture du sol sur lequel on glisse.

🎯 Pourquoi est-ce important ?

Ce travail est une sorte de "manuel de conduite" pour les ingénieurs qui travaillent avec ces matériaux bizarres.

1. Prédiction : Le modèle permet de prédire exactement combien de temps il faudra pour qu'un matériau (comme une colle ou une peinture) se mette à couler, selon la vitesse à laquelle on l'applique.
2. Conception : Cela aide à concevoir de meilleurs procédés industriels, comme l'impression 3D ou le revêtement de surfaces, en évitant que le matériau ne se bloque ou ne se brise de manière imprévisible.
3. Comprendre la rupture : Le modèle montre aussi que selon les conditions aux bords (lisses ou rugueuses), le matériau peut se comporter comme un métal ductile (qui se déforme doucement) ou comme du verre (qui casse brutalement).

En résumé

Les auteurs ont créé une carte mathématique intelligente qui décrit comment un matériau solide et collant se transforme en liquide sous l'effet de la force. Ils ont découvert que :

• La transformation se propage comme une vague de déblocage.
• La vitesse de cette transformation dépend de la force appliquée.
• La texture des parois (lisses ou rugueuses) joue un rôle de lubrifiant invisible qui change toute la donne.

C'est une avancée majeure pour comprendre et maîtriser ces matériaux omniprésents dans notre vie quotidienne, de la crème glacée aux pneus de voiture !

Résumé technique

1. Problématique et Contexte

Les Matériaux Mous Vitreux (SGM) (ex : microgels, mousses, émulsions) sont des assemblages denses de particules colloïdales qui présentent des propriétés solides à l'équilibre mais peuvent se fluidifier sous l'effet d'une contrainte externe suffisante. Ce phénomène de transition solide-liquide induite par le cisaillement est complexe et se manifeste par :

• Des réponses de contrainte non monotones, notamment un pic de contrainte (stress overshoot) suivi d'une relaxation.
• Des écoulements spatialement hétérogènes, tels que le cisaillement en bandes (shear-banding) ou des ruptures fragiles.
• Une dépendance forte aux conditions aux limites et à l'histoire mécanique.

L'objectif de l'article est de proposer une introduction pédagogique et une extension d'un modèle continu basé sur l'approche de la fluidité (fluidity approach) pour rationaliser ces phénomènes, en particulier la transition de fluage et les effets non locaux.

2. Méthodologie : Le Modèle de Fluidité Non Local

Les auteurs développent un modèle continu en deux dimensions (géométrie de cisaillement simple) où l'état du matériau est décrit par un paramètre d'ordre local : la fluidité $f(y,t)$, qui représente le taux d'événements plastiques locaux.

Les composantes clés du modèle sont :

• Équation de dynamique de la fluidité : L'évolution de $f$ est régie par une équation de type gradient de fonctionnelle libre, incluant un terme de diffusion non local (échelle de coopérativité $\xi$) et des termes non linéaires :
  $$\frac{\partial f}{\partial \tilde{t}} = f \left( \xi^2 \Delta f + m f - f^{3/2} \right)$$
  où $\tilde{t}$ est le temps adimensionné et $m$ dépend de la contrainte $\Sigma$. Le choix de la mobilité $\kappa[f] \propto f$ permet la coexistence d'une région fluidifiée ($f>0$) et d'une région solide ($f=0$).
• Équation d'évolution de la contrainte : Couplée à la fluidité via un modèle de Maxwell généralisé, reliant la contrainte totale à une partie élastique et une partie plastique (déformation plastique proportionnelle à la moyenne spatiale de la fluidité $\langle f \rangle$).
• Conditions aux limites : Le modèle explore l'impact crucial des conditions aux limites, notamment la "fluidité de paroi" ($f_w$) imposée par le mouvement de la paroi, qui contrôle la nucléation de la bande de cisaillement.

3. Contributions et Résultats Clés

A. Reproduction de la Phénoménologie du Cisaillement

Le modèle reproduit quantitativement les observations expérimentales sur les microgels de Carbopol :

• Pic de contrainte : Une croissance linéaire initiale (réponse élastique) suivie d'un maximum $\Sigma_M$ et d'une relaxation.
• Cisaillement en bandes transitoires : Formation d'une bande fluidifiée de largeur $\ell_b(t)$ qui grandit depuis la paroi mobile, coexistant avec le matériau solide, jusqu'à l'homogénéisation complète à un temps de fluidisation $\tau_f$.
• Profil de vitesse : Transition d'un profil hétérogène (bande de cisaillement) vers un profil homogène (loi de Herschel-Bulkley) à l'état stationnaire.

B. Lois d'Échelle du Pic de Contrainte ($\Sigma_M$)

L'analyse théorique permet de dériver des lois d'échelle reliant le pic de contrainte à la vitesse de cisaillement imposée $\dot{\Gamma}$ :
$$\Sigma_M - 1 \sim B \left( \frac{\dot{\Gamma}}{\xi \tau^{1/2}} \right)^{\beta(n)} + C \left( \frac{\dot{\Gamma}}{\xi \tau} \right)^{\alpha(n)}$$

• Deux régimes sont identifiés : un régime diffusif (petites vitesses, exposant $\beta = 2n/3$) et un régime asymptotique (grandes vitesses, exposant $\alpha = 4n/(9-n)$).
• Ces prédictions sont en excellent accord avec les données expérimentales sur les microgels, validant le rôle de l'échelle de coopérativité $\xi$.

C. Intégration des Interactions Élasto-Hydrodynamiques (EHD)

Une contribution majeure de l'article est l'intégration des interactions EHD, cruciales pour les particules déformables (microgels) près des parois.

• Le modèle est modifié pour inclure une déformation supplémentaire $\Gamma_{EHD}$ proportionnelle à $\Sigma^2$.
• Résultat : Pour des vitesses de cisaillement faibles, les effets EHD dominent et modifient la loi d'échelle du pic de contrainte en $\Sigma_M \sim \dot{\Gamma}^{1/2}$, indépendamment de l'exposant de rhéologie Herschel-Bulkley ($n$).
• Cela explique les écarts observés expérimentalement entre des surfaces rugueuses (comportement HB classique) et lisses (comportement dominé par le glissement/EHD).

D. Temps de Fluidisation ($\tau_f$)

Le modèle prédit la durée de la phase transitoire de cisaillement en bandes :

• En cisaillement imposé : $\tau_f \sim \dot{\Gamma}^{-9/4}$.
• En contrainte imposée : $\tau_f \sim (\Sigma - 1)^{-9/4n}$.
• Le rapport des exposants de mise à l'échelle entre les deux modes de contrôle est égal à l'exposant $n$ de la loi de Herschel-Bulkley, ce qui est confirmé par l'expérience.

4. Importance et Implications

Ce travail offre un cadre théorique unifié et polyvalent pour comprendre la rhéologie complexe des SGMs :

1. Validation du modèle non local : Il confirme que les effets non locaux (coopérativité) sont essentiels pour décrire la transition de phase dynamique et la formation de bandes de cisaillement.
2. Rôle des conditions aux limites : L'article met en lumière que le type de condition aux limites (fluidité fixée vs gradient de fluidité fixé) dicte la nature de la rupture (ductile vs fragile), expliquant la variabilité des comportements observés en laboratoire.
3. Prédiction des effets EHD : L'incorporation des interactions élasto-hydrodynamiques permet de résoudre des incohérences entre les modèles classiques et les données expérimentales sur les surfaces lisses, reliant la microstructure (déformation des particules) à la macro-rhéologie.
4. Perspectives : Le modèle ouvre la voie à l'étude de protocoles complexes (rampe de vitesse, hystérésis) et à une meilleure compréhension des mécanismes microscopiques aux interfaces.

En résumé, cette étude établit que le modèle de fluidité, enrichi par des effets non locaux et des interactions EHD, constitue un outil puissant pour prédire quantitativement la dynamique de fluage et les transitions de phase dans les matériaux mous vitreux.