Wavemaker and endogeneity of gravitationally stretched weakly viscoelastic jets

Cet article présente un modèle unifié en une dimension pour les jets viscoélastiques étirés par la gravité, démontrant que l'élasticité modifie les conditions critiques d'instabilité et identifiant, grâce à l'analyse de la sensibilité structurelle, la région proche de la buse comme le lieu dominant de la réceptivité du jet.

L'essentiel

🌊 Le Jet de Liquide : Une Danse entre la Gravité et l'Élasticité

Imaginez que vous versez du miel ou de l'eau depuis un robinet. Parfois, le liquide tombe en un filet lisse et régulier (c'est le jet). Parfois, il se brise immédiatement en gouttes qui tombent une par une (c'est le goutte-à-goutte).

Les scientifiques s'intéressent à ce moment précis où le liquide passe d'un filet stable à une série de gouttes. C'est ce qu'on appelle la transition "jet-goutte".

Dans cet article, les chercheurs étudient un cas particulier : un jet de liquide qui s'étire sous l'effet de la gravité (comme un fil de miel qui s'allonge en tombant) et qui contient un peu de polymère (comme de la colle ou du plastique dissous), ce qui le rend légèrement élastique.

🧶 L'Analogie du Fil de Pâte à Modeler

Pour comprendre la différence entre un liquide normal (comme l'eau) et un liquide viscoélastique (comme une solution de polymères), utilisons une analogie :

1. L'Eau (Newtonienne) : Imaginez un fil de pâte à modeler très humide et gluant. Si vous l'étirez, il s'amincit uniformément jusqu'à ce qu'il casse net. C'est prévisible. La science avait déjà très bien compris ce comportement.
2. Le Liquide Élastique (Viscoélastique) : Maintenant, imaginez un fil de pâte à modeler qui contient des élastiques cachés à l'intérieur. Si vous l'étirez, il résiste un peu plus, comme s'il voulait revenir en arrière. Cette "mémoire" de l'élasticité change la façon dont le fil s'amincit et se brise.

🔍 Ce que les chercheurs ont fait

L'équipe a créé un modèle mathématique très précis pour prédire exactement quand ce fil élastique va se briser en gouttes. Ils ont combiné deux idées :

• La forme exacte du fil (qui n'est pas un simple cylindre parfait, mais qui se courbe et s'amincit).
• La tension interne causée par les élastiques (les polymères).

Ils ont ensuite utilisé des outils mathématiques avancés (qu'ils appellent "analyse de stabilité globale") pour voir comment de petites perturbations (un tout petit choc au départ) se propagent le long du fil.

🎯 Les Découvertes Clés (Traduites en langage simple)

Voici les trois grandes conclusions de l'étude, expliquées avec des métaphores :

1. L'Élasticité est un "Frein"

Quand on ajoute des polymères (élasticité) au liquide, le jet devient plus stable.

• L'image : C'est comme si vous ajoutiez un petit frein à un vélo qui descend une pente. Le vélo (le jet) va plus lentement et met plus de temps à décoller ou à se briser.
• Le résultat : Pour qu'un jet élastique se brise en gouttes, il faut que le débit soit encore plus faible que pour l'eau pure. L'élasticité retarde le moment de la rupture.

2. Le "Cœur" de l'instabilité s'étire

C'est la découverte la plus intéressante.

• Pour l'eau : Le "déclencheur" de la rupture (là où le problème commence) est très localisé, juste à la sortie du robinet. C'est comme un point de tension précis.
• Pour le liquide élastique : À cause de la "mémoire" des polymères, le déclencheur s'étire sur une plus grande longueur du fil.
• L'image : Imaginez que pour l'eau, c'est un seul muscle qui se contracte pour casser le fil. Pour le liquide élastique, c'est tout le bras qui se contracte doucement sur plusieurs mètres avant que le fil ne casse. L'instabilité est répartie sur une plus grande zone.

3. La Fréquence change

Le rythme auquel les gouttes se forment ralentit aussi.

• L'image : Si l'eau fait des gouttes comme un métronome rapide (tictac-tictac), le liquide élastique, à cause de son élasticité, fait un métronome plus lent (tac... tac... tac). Les polymères ajoutent un délai, comme un élastique qui met du temps à se détendre.

💡 Pourquoi est-ce important ?

Comprendre cela, ce n'est pas juste de la théorie. Cela aide dans la vie réelle pour :

• La médecine : Créer des gouttes de médicaments très précises (micro-gouttes) pour des inhalateurs ou des injections.
• L'industrie : Fabriquer des fibres ultra-fines pour les vêtements ou les filtres.
• L'impression 3D : Contrôler comment l'encre ou le plastique coule pour éviter les défauts.

En résumé

Les chercheurs ont montré que si vous ajoutez un peu d'élasticité à un filet de liquide qui tombe, il devient plus difficile à briser, il se brise plus lentement, et le "point de rupture" ne se situe plus juste à la sortie du tuyau, mais s'étale sur une plus grande distance.

C'est comme si la gravité tirait sur le fil, mais que les polymères, comme de petits élastiques internes, résistaient et redistribuaient la tension, changeant toute la dynamique de la chute.

Résumé technique

1. Problématique et Contexte

La génération contrôlée de jets liquides étirés par la gravité est fondamentale pour de nombreuses applications technologiques (micro-gouttes, filaments, électrofilage). Pour les fluides newtoniens, la transition entre un régime d'écoulement stable (jetting) et un régime oscillatoire (dripping) a été bien caractérisée par des équations de jets minces unidimensionnels (1D) incluant la courbure interfaciale exacte.

Cependant, les fluides viscoélastiques (solutions polymériques) présentent des dynamiques de minceur et de rupture distinctes, notamment des régimes « perles sur un fil » (beads-on-a-string) et des retards à la rupture dus aux contraintes élastiques. Bien que des modèles locaux (analyse de stabilité spatio-temporelle) aient été utilisés pour étudier ces jets, ils échouent souvent à prédire avec précision les seuils critiques et les fréquences d'oscillation dans les régimes fortement étirés où les conditions aux limites et la non-parallélisme de l'écoulement de base jouent un rôle décisif.

Objectif de l'article : Développer un cadre unifié de stabilité globale pour les jets viscoélastiques étirés par la gravité, en combinant une description 1D complète de la courbure capillaire avec une loi de comportement viscoélastique (Giesekus/Oldroyd-B), et analyser les mécanismes d'instabilité via des diagnostics de sensibilité structurelle.

2. Méthodologie

A. Modélisation Mathématique

Les auteurs formulent un modèle 1D pour un jet axisymétrique incompressible, caractérisé par son rayon $r(z,t)$ et sa vitesse moyenne $u(z,t)$.

• Équations de conservation : Masses et quantités de mouvement (équations de type slender-jet) incluant la gravité, la tension de surface (avec courbure complète) et les contraintes visqueuses.
• Constitution viscoélastique : Le comportement du polymère est décrit par le modèle de Giesekus. Pour simplifier l'analyse, le paramètre de mobilité $\alpha$ est fixé à zéro, ce qui ramène le modèle à la loi d'Oldroyd-B. Les contraintes polymériques ($\tau_{zz}, \tau_{rr}$) sont résolues via des équations d'évolution convectées.
• Conditions aux limites :
  • Entrée ($z=0$) : Rayon fixe, débit imposé, et contraintes polymériques nulles (hypothèse de relaxation complète à la sortie du réservoir).
  • Sortie ($z=L$) : Conditions de jet droit (courbure nulle).

B. Non-dimensionnalisation

Le système est mis sous forme adimensionnelle en utilisant les échelles gravito-capillaires ($\ell_\sigma, U_\sigma$). Les paramètres clés sont :

• Nombre de Weber ($We$) : Rapport inertie/tension de surface.
• Nombre de Bond ($Bo$) : Gravité/tension de surface.
• Nombre de Kapitza ($\Gamma$) : Viscosité du solvant.
• Nombre de Deborah ($De$) : Temps de relaxation du polymère / temps caractéristique gravitaire.
• Ratio de viscosité ($\beta$) : Viscosité du solvant / viscosité totale.

C. Analyse de Stabilité Globale

Au lieu d'une analyse de stabilité locale (parallèle), les auteurs adoptent une approche globale :

1. Décomposition modale : Perturbation linéaire de l'état stationnaire sous la forme $\hat{q}(z)e^{\omega t}$.
2. Problème aux valeurs propres : Résolution du système linéarisé $(\mathbf{A} - \omega \mathbf{B})\hat{q} = 0$ par collocation spectrale (Chebyshev). Cela permet d'obtenir des paires propres globales (taux de croissance $\omega_r$ et fréquence $\omega_i$).
3. Diagnostics de sensibilité (Direct/Adjoint) :
   • Fonction d'onde (Wavemaker) : Calculée à partir du produit des modes directs et adjoints ($S = \|\mathbf{Q}\hat{q}^\dagger\| \|\hat{q}\|$). Elle identifie la région spatiale où l'instabilité est la plus sensible à une perturbation locale (cœur de la boucle de rétroaction).
   • Endogénéité ($E$) : Décomposition spatiale de la contribution de l'opérateur non perturbé à la valeur propre. Elle permet de quantifier les contributions locales (capillaires, cinématiques, élastiques) à la croissance et à la fréquence.

3. Résultats Clés

A. Validation Newtonienne

Le modèle reproduit avec précision les résultats de référence pour les fluides newtoniens (Rubio-Rubio et al., 2013), notamment les spectres marginaux, les profils d'écoulement de base et la dépendance du seuil critique vis-à-vis des nombres $Bo$ et $\Gamma$. Cela confirme la validité de l'implémentation numérique et de l'approche globale.

B. Impact de l'Élasticité sur la Stabilité

Pour des nombres de Deborah modérés ($De = O(1)$), typiques des solutions polymériques diluées :

• Déplacement du seuil critique : L'augmentation de l'élasticité ($De$) déplace la frontière jetting-dripping vers des nombres de Weber critiques ($We_c$) plus faibles. Autrement dit, l'élasticité stabilise le jet, permettant un écoulement stable à des débits plus faibles que pour un fluide newtonien équivalent.
• Fréquence d'oscillation : La fréquence critique d'oscillation ($\omega_i$) diminue avec l'augmentation de $De$. L'élasticité introduit un temps de mémoire qui ralentit la boucle de rétroaction globale.
• Dépendance à la viscosité : L'effet stabilisateur est plus marqué pour les jets à viscosité modérée à élevée ($\Gamma \ge 5.83$). Pour les jets très peu visqueux ($\Gamma = 0.84$), l'effet de l'élasticité est atténué car l'instabilité est contrôlée par un équilibre inertio-capillaire rapide localisé à l'entrée.

C. Mécanismes d'Instabilité et Sensibilité Structurelle

L'analyse par fonctions d'onde et endogénéité révèle des changements profonds dans la nature de l'instabilité :

• Localisation de la réceptivité : Dans les deux cas (newtonien et viscoélastique), le mode adjoint (réceptivité) reste fortement localisé près de la buse (entrée). Cela signifie que les perturbations à l'entrée sont les plus efficaces pour exciter l'instabilité.
• Élargissement de la zone sensible : Contrairement au cas newtonien où la zone de sensibilité (wavemaker) est très localisée, l'élasticité élargit significativement la région sensible en aval. Les contraintes polymériques sont advectées et se relaxent sur une distance finie, couplant ainsi des points axiaux distants.
• Boucle de rétroaction couplée : L'endogénéité montre que l'état marginal résulte d'un équilibre entre :
  1. Une contribution capillaire/cinématique (dominante en amont).
  2. Une contribution élastique ($E_{\tau_{zz}}$) qui s'ajoute et compense partiellement la première.
     La fréquence d'oscillation reste déterminée par le cœur proche de la buse, mais la croissance/decay est contrôlée par une boucle de rétroaction spatialement distribuée incluant la tension polymérique.

4. Contributions Majeures

1. Modélisation unifiée : Développement d'un modèle 1D global intégrant la courbure exacte et la rhéologie viscoélastique (Giesekus/Oldroyd-B) pour des jets étirés par gravité.
2. Validation quantitative : Confirmation que l'élasticité modifie systématiquement les seuils critiques et les fréquences d'instabilité, validant l'approche globale par rapport aux analyses locales.
3. Insight mécanistique : Démonstration que l'effet principal de la viscoélasticité n'est pas seulement un déplacement de la valeur propre, mais une redistribution spatiale de la sensibilité. L'instabilité passe d'un mécanisme purement newtonien localisé à une boucle de rétroaction capillaire-élastique étendue.
4. Outils diagnostics : Application réussie des diagnostics direct/adjoint (wavemaker, endogénéité) pour décomposer les contributions physiques (tension polymérique vs tension de surface) à l'instabilité globale.

5. Signification et Perspectives

Ce travail fournit un cadre prédictif robuste pour la conception de procédés impliquant des jets de polymères (impression 3D, fabrication de fibres, microfluidique). Il montre que pour contrôler la taille des gouttes ou la stabilité des jets, il faut considérer non seulement les paramètres rhéologiques globaux, mais aussi la manière dont l'élasticité étend la zone de couplage dynamique le long du jet.

Les auteurs suggèrent plusieurs extensions futures :

• Intégration de conditions d'entrée plus réalistes (effets de gonflement à la sortie du die).
• Analyse de la réponse forcée et de l'amplification du bruit (analyse resolvent).
• Incorporation d'effets supplémentaires comme la relaxation du profil de vitesse à la sortie ou les surfactants.

En résumé, l'article démontre que la viscoélasticité transforme la nature de l'instabilité des jets étirés en créant une boucle de rétroaction globale étendue, où la tension polymérique joue un rôle stabilisateur crucial et modifie la dynamique de sélection de la fréquence.