A Solid-Based Approach for Modeling Simple Yield-Stress Fluids: Rheological Transitions, Overshoot and Relaxation

Cette étude propose un modèle constitutif basé sur un solide viscoélastique de type Zener qui prédit avec succès les comportements rhéologiques complexes des fluides à seuil, notamment le dépassement de contrainte lors du cisaillement, en l'attribuant à un mécanisme homogène d'amplification du tenseur des contraintes par les différences de contraintes normales.

L'essentiel

🧱 Le Secret des Fluides qui "Décident" de Couler : Une Nouvelle Recette Mathématique

Imaginez que vous avez un pot de moutarde ou de peinture épaisse. Si vous appuyez doucement avec votre doigt, rien ne bouge : c'est dur comme un solide. Mais si vous appuyez fort ou si vous secouez le pot, ça coule comme un liquide.

C'est ce qu'on appelle un fluide à seuil (ou yield-stress fluid). Ces matériaux sont partout : dans les boues des rivières, les pâtes dentifrices, et même dans les encres pour imprimer des batteries.

Le problème, c'est que les scientifiques avaient du mal à prédire exactement comment ces matériaux se comportent quand on les fait bouger vite, surtout quand ils montrent des comportements bizarres comme un pic de résistance avant de couler (comme si le matériau disait "Non, je résiste !" avant de céder).

C'est là que cette étude intervient. Les auteurs ont créé une nouvelle "recette mathématique" (une équation) pour prédire ce comportement avec une précision incroyable.

---

🏗️ L'Analogie du "Gâteau et du Sirop"

Pour comprendre leur modèle, imaginez que le matériau (comme le Carbopol, une sorte de gel utilisé dans les expériences) est composé de deux choses mélangées :

1. Le "Gâteau" (la partie solide) : Ce sont les micro-gouttelettes de gel qui sont serrées les unes contre les autres, comme des billes dans un bocal. Elles forment une structure solide.
2. Le "Sirop" (le liquide entre les billes) : C'est l'eau qui coule entre les billes.

L'ancienne façon de voir les choses :
Les anciens modèles traitaient ces matériaux comme des liquides compliqués. Ils disaient : "Si on pousse, ça coule". Mais ils ne savaient pas bien expliquer pourquoi, quand on arrête de pousser, le matériau garde une certaine forme (comme un solide) au lieu de tout s'effondrer.

La nouvelle idée des auteurs :
Ils ont dit : "Attendez, ce matériau se comporte d'abord comme un solide élastique (comme un ressort), et seulement si on pousse trop fort, il commence à couler."

Ils ont construit leur modèle mathématique comme un assemblage de pièces mécaniques :

• Un ressort (le Gel) : Il représente la structure solide. Il peut se déformer et revenir en place.
• Un amortisseur (le Sirop) : Il représente l'eau qui frotte et dissipe l'énergie.
• Une règle de "cassure" : Quand la force est trop forte, le ressort ne revient plus tout à fait en place (c'est la partie plastique).

---

🎢 Les Trois Phénomènes Magiques Expliqués

Leur nouvelle équation réussit à prédire trois choses que les anciens modèles rataient souvent :

1. Le "Pic de Stress" (Stress Overshoot) 📈

Quand vous commencez à mélanger une peinture épaisse, vous devez pousser très fort au début. La résistance monte, atteint un pic (le moment où vous sentez que ça va céder), puis baisse un peu pour se stabiliser.

• L'explication du papier : Ce pic n'est pas dû à un changement de structure interne complexe ou à de la "fatigue" du matériau. C'est simplement une conséquence mathématique de la façon dont les forces s'organisent dans les trois dimensions de l'espace. C'est comme si le matériau se tendait comme un élastique avant de se détendre.

2. Le Repos qui ne finit jamais (Relaxation) 🛑

Si vous étirez ce gel et que vous l'arrêtez net, il ne revient pas à zéro. Il garde une tension résiduelle.

• L'explication : Le modèle montre que le matériau a une "mémoire" solide. Même au repos, il garde une certaine rigidité, comme un ressort qui reste légèrement tendu. De plus, plus vous l'avez étiré vite avant de l'arrêter, plus cette tension résiduelle est faible (comme un élastique qu'on a trop étiré).

3. Le Passage de Solide à Liquide (Fluage) 🧊➡️💧

Si vous posez un poids léger sur le gel, il reste immobile (comportement solide). Si vous posez un poids lourd, il coule lentement mais sûrement (comportement liquide).

• L'explication : Le modèle trouve exactement le seuil de poids où le matériau bascule. En dessous, il est solide. Au-dessus, il devient fluide. C'est une frontière très nette que leur équation décrit parfaitement.

---

🌍 Pourquoi est-ce important ?

Avant, pour simuler ces matériaux sur ordinateur (par exemple pour concevoir une usine de batteries ou prévoir une inondation de boue), il fallait faire des approximations qui ne marchaient pas toujours.

Grâce à cette nouvelle équation :

• Les ingénieurs peuvent mieux concevoir des encres pour batteries (pour qu'elles s'écoulent bien dans les machines sans boucher les tuyaux).
• On peut mieux comprendre les glissements de terrain ou les coulées de boue.
• On a une vision plus juste de la nature : ces matériaux sont des solides qui savent se transformer en liquides, et non pas des liquides bizarres.

🎯 En résumé

Les auteurs ont dit : "Arrêtons de traiter ces matériaux comme des liquides compliqués. Traitez-les comme des solides élastiques qui ont un mécanisme de sécurité pour se transformer en liquide quand on les pousse trop fort."

Leur modèle est comme un GPS mathématique très précis qui nous dit exactement comment ces matériaux réagiront, que ce soit en les étirant doucement ou en les secouant violemment. C'est une avancée majeure pour comprendre la physique de la matière molle !

Résumé technique

1. Problématique et Contexte

Les fluides à seuil (yield-stress fluids), tels que les dispersions de Carbopol ou les boues d'électrodes pour batteries, sont omniprésents dans l'industrie et la nature. Bien que leur comportement en écoulement stationnaire soit souvent décrit par le modèle d'Herschel-Bulkley, leur réponse transitoire (démarrage de cisaillement, relaxation de contrainte, fluage) reste difficile à modéliser avec précision.

Les défis principaux identifiés par les auteurs sont :

• L'absence de modèles unifiés : Peu de modèles parviennent à prédire simultanément les comportements stationnaires et transitoires.
• Le problème du « dépassement de contrainte » (Stress Overshoot) : Lors du démarrage d'un écoulement de cisaillement, la contrainte atteint souvent un pic avant de se stabiliser. Les modèles existants attribuent souvent ce phénomène à l'écrouissage isotrope ou à l'hétérogénéité microstructurale (thixotropie), ce qui n'est pas toujours le cas pour des fluides simples comme le Carbopol (faible thixotropie).
• La contrainte relaxée non nulle : Contrairement aux fluides viscoélastiques classiques qui relaxent totalement la contrainte à zéro, les fluides à seuil comme le Carbopol présentent une contrainte résiduelle non nulle après l'arrêt de l'écoulement, un comportement typique des solides.
• La transition solide-fluide : En fluage, ces matériaux passent d'un comportement solide (saturation de la déformation) à un comportement fluide (déformation continue) selon l'amplitude de la contrainte appliquée.

L'objectif de l'étude est de développer une équation constitutive capable de reproduire ces quatre caractéristiques clés sans recourir à des hypothèses complexes d'hétérogénéité spatiale ou de thixotropie forte.

2. Méthodologie

Les auteurs proposent une approche basée sur la mécanique des solides viscoélastiques, inspirée par la microstructure des dispersions de Carbopol (assemblage jammed de microgels dans un solvant).

A. Modèle 1D (Déformation simple)

Le modèle initial est construit sur une analogie mécanique :

• Élément Zener (Phase gel) : Un ressort linéaire (déformation élastique) en série avec un élément de Kelvin-Voigt (ressort + amortisseur, représentant la déformation plastique et le réarrangement des microgels).
• Amortisseur linéaire parallèle (Phase solvant) : Représente la dissipation visqueuse du solvant interstitiel, ajoutant une contrainte proportionnelle au taux de déformation total.
• Décomposition de la déformation : Additive ($\gamma = \gamma_e + \gamma_p$) pour le modèle 1D.
• Loi d'écoulement et évolution de la contrainte de retour (Back stress) : Utilisation d'une loi de type Armstrong-Frederick pour l'évolution de la contrainte de retour ($\tau_B$), liée à la déformation plastique.
• Modèles de viscosité : Deux modèles sont testés pour décrire la viscosité dépendante de la contrainte : le modèle d'Eyring et le modèle de Carreau-Yasuda (converti en fonction de la contrainte).

B. Extension 3D (Grandes déformations)

Pour respecter l'invariance du cadre matériel et traiter des déformations finies, le modèle est étendu en 3D :

• Décomposition multiplicative (Kröner-Lee) : Le gradient de déformation $F$ est décomposé en $F = F_e \cdot F_p$.
• Tenseurs : Utilisation de tenseurs de Cauchy pour les contraintes et de tenseurs de déformation de gauche (B) pour la partie élastique.
• Invariance : Les équations d'évolution sont formulées pour être invariantes sous les rotations rigides.
• Loi d'écoulement tensorielle : La vitesse de déformation plastique est liée à la différence de contrainte déviatorique entre la phase gel et la contrainte de retour.

3. Contributions Clés

1. Approche « Solide » unifiée : Le modèle traite le fluide à seuil comme un solide viscoélastique avec une dissipation solvant parallèle, expliquant naturellement la contrainte relaxée non nulle.
2. Origine du dépassement de contrainte (Overshoot) : L'étude démontre que le dépassement de contrainte en cisaillement de démarrage peut émerger d'un mécanisme homogène dans un modèle continu 3D. Il ne nécessite ni écrouissage isotrope ni hétérogénéité spatiale. Le phénomène est dû à l'interaction entre les contraintes normales (qui se développent même en cisaillement simple en 3D) et la loi de viscosité, accélérant temporairement la réponse plastique.
3. Modélisation de la transition de fluage : Le modèle prédit correctement le seuil de transition entre un comportement de fluage solide (saturation) et fluide (écoulement continu).
4. Dépendance de la contrainte relaxée : Contrairement au modèle 1D, le modèle 3D prédit que la contrainte relaxée finale dépend du taux de cisaillement pré-shearing, en accord avec les observations expérimentales récentes.

4. Résultats

Les simulations numériques et les analyses analytiques ont validé le modèle contre les critères suivants :

• Démarrage de cisaillement (Start-up shear) :

  • Le modèle 1D ne montre pas de dépassement de contrainte (convergence monotone).
  • Le modèle 3D reproduit qualitativement et quantitativement le dépassement de contrainte, dont l'amplitude augmente avec le taux de cisaillement, en accord avec les données expérimentales sur le Carbopol.
  • L'analyse montre que le dépassement résulte d'un moment où le terme de relaxation (plastique) dépasse temporairement le terme de chargement élastique, favorisé par les contraintes normales dans la formulation 3D.

• Comportement stationnaire :

  • À faible taux de cisaillement, la contrainte stationnaire suit la loi d'Herschel-Bulkley.
  • À haut taux de cisaillement, le modèle capture la déviation vers un deuxième plateau newtonien (dû à la viscosité du solvant), ce que le modèle d'Herschel-Bulkley simple ne fait pas.

• Relaxation de contrainte :

  • Après l'arrêt de l'écoulement, la contrainte ne tombe pas à zéro mais converge vers une contrainte relaxée non nulle.
  • Dans le modèle 3D, cette contrainte relaxée diminue lorsque le taux de cisaillement pré-shearing augmente, reproduisant les résultats de Vinutha et al. (2024).

• Essai de fluage (Creep) :

  • Le modèle prédit une transition nette : pour des contraintes appliquées inférieures au seuil de yield (estimé par le fit d'Herschel-Bulkley), la déformation sature (comportement solide). Au-delà, la déformation augmente indéfiniment avec un taux de déformation non nul (comportement fluide).

5. Signification et Perspectives

Cette étude apporte une avancée significative dans la rhéologie des fluides complexes :

• Changement de paradigme : Elle valide l'hypothèse que les fluides à seuil simples peuvent être modélisés comme des solides viscoélastiques, offrant une alternative aux modèles fluides basés sur la thixotropie.
• Simplicité et Robustesse : Le modèle parvient à capturer des phénomènes complexes (overshoot, relaxation non nulle) avec un nombre limité de paramètres physiques, sans avoir besoin de supposer des mécanismes microstructuraux hétérogènes.
• Applications industrielles : Ce cadre constitutif est particulièrement pertinent pour la simulation numérique (éléments finis) de procédés impliquant des fluides à seuil, tels que l'impression 3D de encres (Direct Ink Writing) ou la fabrication de batteries (pâtes d'électrodes), où les écoulements sont souvent complexes et transitoires.

En conclusion, l'équation constitutive proposée offre un cadre théorique robuste et généralisable pour prédire le comportement dynamique des fluides à seuil, comblant le fossé entre les observations expérimentales et la modélisation mécanique.