On the deformation of a shear thinning viscoelastic drop in a steady electric field

Cette étude numérique, utilisant le solveur Basilisk, examine la déformation et la rupture d'une goutte viscoélastique de type LPTT dans un champ électrique uniforme, révélant que l'élasticité modifie significativement la dynamique de déformation et les régimes de rupture selon les régions de conductivité et de permittivité, avec des comportements non monotones observés dans certaines zones du diagramme de phase.

L'essentiel

🧪 Quand une goutte de plastique rencontre l'électricité : Une histoire de formes et de résistance

Imaginez que vous tenez une goutte d'eau (ou plutôt, une goutte de "plastique liquide" un peu élastique) entre vos mains. Maintenant, imaginez que vous appliquez une très forte force électrique autour d'elle. Que va-t-il se passer ?

Cette étude, réalisée par des chercheurs de l'Institut indien des sciences, répond à cette question. Ils ont simulé le comportement de gouttes viscoélastiques (des fluides qui sont à la fois liquides et élastiques, comme du miel très épais ou des solutions de polymères) soumises à un champ électrique constant.

Voici les grandes lignes de leur découverte, expliquées avec des métaphores du quotidien.

1. Le décor : La goutte et le champ électrique

Imaginez une petite goutte de gelée suspendue dans un bain d'huile.

• Le champ électrique agit comme un vent invisible et puissant qui pousse la goutte.
• La tension de surface est comme une peau élastique qui veut garder la goutte ronde et compacte.
• La viscoélasticité est la "mémoire" de la goutte : elle veut se déformer sous la poussée, mais elle a aussi envie de revenir à sa forme initiale, comme un ressort.

Les chercheurs ont utilisé un super-ordinateur (le logiciel Basilisk) pour voir comment ces gouttes se comportent selon la force du "vent électrique" et la "fermeté" de la goutte.

2. Les trois scénarios possibles (Le jeu des formes)

Selon les propriétés électriques de la goutte par rapport à son environnement, trois destins différents s'offrent à elle :

• Le scénario "Élastique calme" (Régions PR- et OB+) :
  Imaginez une goutte très résistante. Même si le vent électrique souffle fort, la goutte s'étire un peu, comme un élastique qu'on tire doucement, puis elle se stabilise. Elle ne casse pas. Dans ce cas, le fait que la goutte soit "élastique" ne change pas grand-chose par rapport à une goutte d'eau normale.

• Le scénario "L'explosion en plusieurs morceaux" (Région PR+A) :
  Ici, la goutte s'étire beaucoup. Si le vent électrique est trop fort, au lieu de se casser net, elle commence à former des "lobes" (comme des petits ballons attachés les uns aux autres).

  • La surprise : Plus la goutte est élastique (plus elle a de "ressort"), plus elle résiste ! C'est comme si vous aviez un élastique très épais : il faut beaucoup plus de force pour le casser. Les chercheurs ont vu que plus la goutte est élastique, plus elle peut supporter un vent électrique violent avant de se briser.

• Le scénario "Le cône pointu" (Région PR+B) :
  C'est le plus spectaculaire. La goutte s'étire et finit par former des pointes très fines à ses extrémités, comme un stylo ou un cône de glace.

  • Le paradoxe : Ici, l'élasticité joue un rôle bizarre. À un niveau moyen d'élasticité, la goutte résiste bien. Mais si elle devient trop élastique, elle finit par se déformer plus facilement et former ces pointes dangereuses. C'est comme si un ressort trop tendu finissait par céder d'un coup.

• Le scénario "Le disque plat" (Région OB-) :
  Au lieu de s'allonger, la goutte s'écrase et devient plate, comme une crêpe ou un disque de hockey.

  • Le comportement en "montagne russe" : Pour cette forme, l'élasticité a un effet curieux. Au début, quand on augmente l'élasticité, la goutte s'écrase plus vite. Mais si on continue d'augmenter l'élasticité, elle redevient plus résistante et s'écrase moins. C'est une relation non-linéaire : ni trop, ni trop peu, il faut trouver le bon équilibre.

3. La grande différence : Le modèle "Oldroyd-B" vs Le modèle "LPTT"

Les chercheurs ont comparé deux façons de modéliser ces gouttes élastiques :

1. Le modèle "Oldroyd-B" (L'ancien modèle) : Il imagine que les molécules de la goutte peuvent s'étirer à l'infini. C'est un peu comme imaginer un élastique magique qui ne casse jamais, même si on le tire jusqu'à ce qu'il devienne infiniment fin. Cela donne des résultats parfois irréalistes.
2. Le modèle "LPTT" (Le nouveau modèle de cette étude) : Il prend en compte la réalité : les molécules ont une longueur maximale. On ne peut pas les étirer à l'infini. C'est comme un vrai élastique qui a une limite avant de se rompre.

Le verdict :
Le nouveau modèle (LPTT) montre que les gouttes sont souvent plus stables et plus résistantes que ce que l'ancien modèle prédisait, surtout quand le champ électrique est très fort. L'ancien modèle sous-estimait la capacité de la goutte à résister à la rupture.

4. Pourquoi est-ce important ? (À quoi ça sert ?)

Comprendre comment ces gouttes se déforment et se cassent est crucial pour notre quotidien et l'industrie :

• Impression 3D et encres : Pour imprimer des gouttes de plastique parfaites sur une page.
• Médicaments et biologie : Pour manipuler des gouttelettes contenant des médicaments ou des cellules sans les abîmer.
• Météo : Pour mieux comprendre comment les gouttes de pluie se brisent dans les orages (ce qui crée la foudre !).

En résumé

Cette étude nous apprend que la "mémoire élastique" d'une goutte est un super-pouvoir complexe. Selon la force du champ électrique et la nature de la goutte, l'élasticité peut soit la protéger (en retardant sa rupture), soit la rendre plus fragile (en favorisant la formation de pointes).

C'est comme si on apprenait à piloter un bateau dans une tempête : parfois, avoir un moteur puissant (l'élasticité) vous aide à résister aux vagues, mais parfois, si vous forcez trop, vous risquez de vous retourner. Les chercheurs ont maintenant la carte pour savoir exactement quand et comment cela se produit.

Résumé technique

Voici un résumé technique détaillé de l'article de recherche intitulé « On the deformation of a shear thinning viscoelastic drop in a steady electric field » (Sur la déformation d'une goutte viscoélastique à amincissement de cisaillement dans un champ électrique stationnaire).

1. Problématique et Contexte

L'étude se concentre sur la dynamique de déformation et de rupture d'une goutte viscoélastique soumise à un champ électrique uniforme stationnaire. Ce phénomène est crucial pour diverses applications industrielles et naturelles, notamment l'impression jet d'encre, la microfluidique, le déshuilage et l'électro-spraying.

Bien que le comportement des gouttes newtoniennes et des gouttes viscoélastiques modélisées par le modèle d'Oldroyd-B ait été largement étudié, ces modèles présentent des limitations physiques :

• Le modèle d'Oldroyd-B suppose une extensibilité infinie des chaînes polymères, ce qui conduit à une divergence non physique de la viscosité en extension et ne capture pas l'amincissement de cisaillement (shear-thinning).
• De nombreux fluides réels (solutions polymères, fluides biologiques) présentent un comportement d'amincissement de cisaillement et une extensibilité finie.

L'objectif de cette recherche est donc d'investiguer la déformation d'une goutte décrite par le modèle linéaire de Phan-Thien–Tanner (LPTT), qui intègre à la fois l'amincissement de cisaillement et l'extensibilité finie des polymères, et de comparer ses résultats avec ceux du modèle d'Oldroyd-B.

2. Méthodologie

• Modélisation Mathématique :

  • La goutte (phase interne) est modélisée comme un fluide viscoélastique incompressible régi par l'équation constitutive LPTT.
  • Le fluide environnant (phase externe) est newtonien.
  • Les équations de conservation de la masse, de la quantité de mouvement et de l'électricité (lois de Gauss et Faraday) sont résolues en couplant les contraintes Maxwell (forces électriques) avec la dynamique des fluides.
  • Le modèle LPTT introduit un paramètre d'extensibilité ($\varepsilon$) et un temps de relaxation ($\lambda$), permettant de limiter la croissance des contraintes à haut taux de déformation.

• Paramètres Adimensionnels :
  L'étude explore l'espace des paramètres défini par :

  • Le nombre capillaire électrique ($Ca_E$) : rapport entre les forces électriques et la tension de surface.
  • Le nombre de Deborah ($De$) : rapport entre le temps de relaxation du polymère et l'échelle de temps de l'écoulement.
  • Les rapports de conductivité ($\sigma_r$) et de permittivité ($\epsilon_r$).

• Approche Numérique :

  • Les simulations sont réalisées avec le solveur open-source Basilisk, utilisant une méthode de volume de fluide (VOF) géométrique pour capturer l'interface.
  • L'équation constitutive viscoélastique est résolue via l'approche du tenseur de conformation logarithmique (log-conformation tensor) pour assurer la stabilité numérique.
  • Une étude de convergence de maillage et d'indépendance du domaine a été effectuée pour valider les résultats.

• Stratégie d'Étude :
  Les simulations couvrent six régions caractéristiques de l'espace des phases $(\sigma_r, \epsilon_r)$ (notées $PR^+_A, PR^+_B, PR^-_A, PR^-_B, OB^+, OB^-$), identifiées par des analyses asymptotiques précédentes.

3. Résultats Clés et Contributions

Les résultats sont analysés en fonction des régions de l'espace des phases :

A. Régions $PR^-_A, PR^-_B$ et $OB^+$

Dans ces régions, où les coefficients de déformation d'ordre 1 et 2 ont des signes opposés, les gouttes LPTT montrent un comportement très proche de celui des gouttes newtoniennes. La viscoélasticité a un effet négligeable sur la dynamique de déformation et de rupture.

B. Région $PR^+_A$ (Déformation prolate et rupture)

• Comportement : En dessous d'un $Ca_E$ critique, la goutte forme un sphéroïde prolate stable. Au-delà, elle subit une rupture ou forme des structures multi-lobe.
• Effet de l'élasticité ($De$) : L'augmentation de $De$ s'oppose à la déformation, augmentant ainsi le $Ca_E$ critique nécessaire à la rupture.
• Comparaison Oldroyd-B vs LPTT :
  • Les deux modèles montrent une réduction de la déformation avec l'augmentation de $De$.
  • Cependant, les gouttes LPTT présentent une rupture plus retardée à haut $De$ et forment des structures multi-lobe stables uniquement à faible $De$, tandis que les gouttes Oldroyd-B peuvent maintenir des formes multi-lobe stables sur une plage plus large de $De$.

C. Région $PR^+_B$ (Formation de pointes)

• Comportement : La goutte forme des sphéroïdes stables jusqu'à un seuil critique, puis développe des extrémités pointues (pointed ends).
• Dépendance non monotone : Contrairement au modèle d'Oldroyd-B (où la déformation diminue monotone avec $De$), la déformation LPTT présente un comportement non monotone : elle augmente initialement avec $De$ (due au durcissement par déformation ou strain-hardening) avant de diminuer à des $De$ élevés (due à la saturation des contraintes et à l'extensibilité finie).
• Résistance à la déformation : Les gouttes LPTT montrent une plus grande déformation globale que les gouttes Oldroyd-B pour les $De$ non nuls, indiquant une résistance moindre dans cette région spécifique due à la réponse non linéaire du modèle LPTT.

D. Région $OB^-$ (Déformation oblate et rupture)

• Comportement : La goutte s'aplatit (oblate) jusqu'à un seuil critique, puis se rompt.
• Stabilité : Les gouttes LPTT sont plus stables que les gouttes Oldroyd-B. Elles résistent mieux à la rupture à haut $De$.
• Dynamique de rupture : Pour les gouttes Oldroyd-B, l'augmentation de $De$ favorise la déformation et accélère la rupture. Pour les gouttes LPTT, la déformation augmente initialement puis diminue à haut $De$, rendant la goutte plus résistante à la rupture. Le temps de rupture des gouttes LPTT augmente généralement avec $De$, contrairement aux gouttes Oldroyd-B qui se rompent plus rapidement.

4. Signification et Conclusion

Cette étude met en évidence l'importance cruciale de l'extensibilité finie et de l'amincissement de cisaillement dans la prédiction précise de la dynamique des gouttes viscoélastiques sous champ électrique.

• Limites du modèle Oldroyd-B : Le modèle Oldroyd-B surestime souvent la déformation et sous-estime la stabilité des gouttes à haut nombre de Deborah, en raison de son incapacité à capturer la saturation des contraintes.
• Apport du modèle LPTT : Le modèle LPTT offre une description plus réaliste, montrant que la viscoélasticité peut soit stabiliser, soit déstabiliser une goutte selon la région de l'espace des phases et le niveau d'élasticité.
• Implications : Ces résultats sont essentiels pour le contrôle précis des gouttes dans les systèmes électrohydrodynamiques (EHD), permettant d'optimiser des procédés comme l'impression jet d'encre ou la manipulation de fluides complexes où la rupture prématurée ou la déformation excessive doivent être évitées.

En résumé, l'interplay complexe entre la viscoélasticité (via le modèle LPTT) et les contraintes électriques détermine non seulement la morphologie finale des gouttes, mais aussi leur stabilité et leur seuil de rupture, offrant des perspectives nouvelles pour la manipulation contrôlée de fluides complexes.