Escalade transitoire de tige dans un fluide d'Oldroyd-B

Informations fondamentales

• ID de l'article: 2501.01016
• Titre: Escalade transitoire de tige dans un fluide d'Oldroyd-B
• Auteurs: Tachin Ruangkriengsin, Rodolfo Brandão, Katie Wu, Jonghyun Hwang, Evgeniy Boyko, Howard A. Stone
• Classification: physics.flu-dyn (mécanique des fluides)
• Date de publication: 2 janvier 2025 (soumis à arXiv)
• Journal: Journal of Fluid Mechanics (en cours de soumission)
• Lien de l'article: https://arxiv.org/abs/2501.01016

Résumé

Cet article étudie l'effet Weissenberg dans les fluides non-newtoniens (phénomène d'escalade de tige), c'est-à-dire l'ascension de l'interface du fluide le long d'une tige en rotation. Les auteurs adoptent le modèle d'Oldroyd-B et, par analyse asymptotique et variables adimensionnelles, étudient la hauteur d'interface dépendante du temps. L'étude révèle que pour les petites déformations, les équations de gouvernance indiquent l'existence d'une couche limite temporelle, où l'interface monte rapidement sur une courte échelle de temps, puis s'approche progressivement de l'état stationnaire. En introduisant une variable de temps étirée, les auteurs déduisent le champ de vitesse transitoire et le profil d'interface correspondant, et réintroduisent les petits effets d'inertie finis pour étudier leur interaction avec la viscoélasticité.

Contexte et motivation de la recherche

Définition du problème

Le phénomène d'escalade de tige (effet Weissenberg) est l'un des comportements les plus caractéristiques des fluides non-newtoniens. Lorsqu'une tige tourne dans un fluide viscoélastique, la surface du fluide grimpe le long de la tige. Ce phénomène manque encore d'une compréhension théorique approfondie.

Importance

1. Signification théorique: Ce phénomène est un concept fondamental pour comprendre le comportement des fluides complexes, avec une importance majeure dans l'enseignement et la recherche
2. Valeur pratique: Applications importantes dans les mesures rhéologiques, permettant de déterminer les paramètres matériaux tels que le coefficient de contrainte normale et le temps de relaxation
3. Applications d'ingénierie: Fournit des orientations pour les applications industrielles des solutions polymères et autres fluides complexes

Limitations des méthodes existantes

1. Limitations théoriques: Les recherches antérieures se concentraient principalement sur le comportement stationnaire, utilisant des modèles de fluides du second ordre
2. Problèmes de modèle: Les modèles de fluides du second ordre ne sont valides que pour les petits nombres de Deborah ou de Weissenberg, et sont hautement instables
3. Dépendance temporelle: Absence d'analyse théorique de la dynamique transitoire

Motivation de la recherche

Les auteurs choisissent le modèle d'Oldroyd-B pour éviter les problèmes des modèles de fluides du second ordre dans les écoulements non-stationnaires, et étudient pour la première fois systématiquement la dépendance temporelle du phénomène d'escalade de tige.

Contributions principales

1. Première analyse théorique: Fournit la première analyse théorique transitoire complète du phénomène d'escalade de tige dans un fluide d'Oldroyd-B
2. Théorie de la couche limite temporelle: Découvre et analyse la structure de la couche limite temporelle, révélant le processus en deux étapes de montée rapide et d'approche lente vers l'état stationnaire
3. Cadre d'analyse adimensionnelle: Établit une méthode d'analyse adimensionnelle systématique basée sur le nombre de Weissenberg
4. Étude des effets d'inertie: Première étude de l'interaction entre les petits effets d'inertie finis et la viscoélasticité
5. Critère d'escalade: Propose un critère théorique pour déterminer les conditions d'occurrence de l'escalade de tige: $Re < 4\beta_p Wi$

Détails méthodologiques

Définition de la tâche

Étudier l'évolution temporelle de la hauteur de la surface libre du fluide $h(r,t)$ lorsqu'une tige cylindrique infiniment longue (rayon $a$) tourne à une vitesse angulaire $\Omega$ dans une solution polymère viscoélastique.

Modèle mathématique

Équations de gouvernance

Adoption du modèle constitutif d'Oldroyd-B:

• Équation d'évolution du tenseur de conformation: $\frac{\partial\mathbf{A}}{\partial t} + \mathbf{u} \cdot \nabla\mathbf{A} - \mathbf{A} \cdot \nabla\mathbf{u} - (\nabla\mathbf{u})^T \cdot \mathbf{A} = -\frac{1}{\lambda}(\mathbf{A} - \mathbf{I})$
• Tenseur de contrainte: $\boldsymbol{\sigma} = -p\mathbf{I} + 2\mu_s\mathbf{E} + G(\mathbf{A} - \mathbf{I})$
• Équation de quantité de mouvement: $\nabla \cdot \boldsymbol{\sigma} = \rho g\mathbf{e}_z$

Adimensionnalisation

Introduction des paramètres adimensionnels clés:

• Nombre de Weissenberg: $Wi = \lambda\Omega$
• Paramètre de gravité: $G = \frac{\rho ga}{(\mu_s + \mu_p)\Omega}$
• Nombre de Reynolds: $Re = \frac{\rho a^2\Omega}{\mu_s + \mu_p}$

Limitation caractéristique

Les auteurs identifient la limitation caractéristique clé: $G = Wi \cdot \mathcal{G}$, où $\mathcal{G} = O(1)$, ce qui simplifie l'expansion asymptotique à deux paramètres en une expansion à un seul paramètre.

Points d'innovation technique

1. Séparation des échelles de temps

Découverte de deux échelles de temps:

• Temps externe: $T = O(1)$, correspondant au comportement stationnaire
• Temps interne: $\tau = T/Wi$, correspondant à la couche limite transitoire

2. Méthode de perturbation de domaine

Pour les petites déformations d'interface, expansion en série de Taylor de la condition aux limites de $Z = H(R,T)$ à $Z = 0$, simplifiant le traitement des conditions aux limites de surface libre.

3. Transformation de Weber

Lors du traitement des effets d'inertie, application de la transformation de Weber pour traiter l'opérateur différentiel de Bessel sur un domaine semi-infini, convertissant l'équation aux dérivées partielles en équation différentielle ordinaire.

Configuration expérimentale

Cadre d'analyse théorique

Cet article est une étude purement théorique utilisant la méthode d'analyse asymptotique:

1. Petit paramètre: $Wi \ll 1$ (viscoélasticité faible)
2. Plage de gravité: $G \gg Wi^{-1}$ (effets gravitationnels forts)
3. Hypothèse de petite déformation: Déformation d'interface relativement faible

Stratégie d'analyse

1. Section 4: Analyse externe - profil d'interface stationnaire
2. Section 5: Analyse interne - profil d'interface transitoire
3. Section 6: Analyse incluant les petits effets d'inertie

Résultats expérimentaux

Résultats principaux

1. Profil d'interface stationnaire

$H_{steady}(R) = \frac{2\beta_p Wi^2}{\mathcal{G}R^4}$

Ce résultat est cohérent avec celui obtenu par Joseph & Fosdick (1973) en utilisant un modèle de fluide du second ordre.

2. Profil d'interface transitoire

$H_{transient}(R,\tau) = \frac{2\beta_p(1-(\tau+1)e^{-\tau})Wi^2}{\mathcal{G}R^4}$

où $\tau = T/Wi$ est la variable de temps étirée.

3. Critère d'escalade

Lorsque les effets d'inertie sont inclus, la condition d'occurrence de l'escalade de tige est: $4El\beta_p > 1$ où $El = Wi/Re$ est le nombre d'élasticité.

Caractéristiques de l'évolution temporelle

1. Phase de montée rapide: L'interface monte rapidement dans un intervalle de temps $\tau \approx 6$
2. Approche stationnaire: La solution transitoire converge vers la solution stationnaire aux temps longs
3. Comportement non-monotone: Avec l'inertie incluse, la hauteur d'interface peut d'abord augmenter puis diminuer

Découvertes clés

1. Contrainte azimutale: La contrainte azimutale produite par l'étirement du polymère est le mécanisme principal entraînant l'escalade
2. Compétition inertielle: La force centrifuge produite par l'inertie entre en compétition avec la contrainte azimutale de la viscoélasticité
3. Structure de couche limite: L'existence de la couche limite temporelle explique le comportement en deux étapes de transitoire rapide et de stationnaire lent

Travaux connexes

Développement historique

1. Serrin (1959): Première tentative théorique, utilisant un modèle de fluide du second ordre
2. Giesekus (1961): Équations constitutives plus générales négligeant les effets d'inertie
3. Joseph & Fosdick (1973): Cadre systématique de la méthode de perturbation de domaine
4. Oldroyd (1950, 1951): Applications précoces du modèle d'Oldroyd-B

Avantages de cet article

1. Dépendance temporelle: Première analyse transitoire complète
2. Stabilité du modèle: Évite les problèmes d'instabilité des modèles de fluides du second ordre
3. Cadre adimensionnel: Analyse systématique des paramètres adimensionnels
4. Effets d'inertie: Première étude de l'interaction entre inertie et viscoélasticité

Conclusions et discussion

Conclusions principales

1. Structure de couche limite: Le phénomène d'escalade de tige présente une couche limite temporelle distincte, avec une montée rapide de l'interface suivie d'une approche lente vers l'état stationnaire
2. Séparation d'échelles: L'échelle de temps transitoire est $O(Wi)$, l'échelle de temps stationnaire est $O(1)$
3. Mécanisme d'escalade: La contrainte azimutale produite par l'étirement du polymère est le mécanisme fondamental entraînant l'escalade
4. Influence de l'inertie: Les petits effets d'inertie peuvent influencer significativement le comportement d'escalade, avec un nombre d'élasticité critique

Limitations

1. Hypothèse de petite déformation: Applicable uniquement aux petites déformations d'interface
2. Viscoélasticité faible: Limitée à la plage $Wi \ll 1$
3. Gravité forte: Nécessite la condition $G \gg Wi^{-1}$
4. Modèle simplifié: Néglige les effets de tension superficielle et d'angle de contact

Directions futures

1. Nombre de Weissenberg élevé: Étudier le comportement pour $Wi \gg 1$
2. Modèles constitutifs complexes: Extension aux modèles FENE-P, Giesekus, etc.
3. Déformations finies: Traiter les cas de déformation d'interface importante
4. Validation expérimentale: Comparaison avec les données expérimentales

Évaluation approfondie

Points forts

1. Rigueur théorique: Dérivations mathématiques rigoureuses, application appropriée de la méthode d'analyse asymptotique
2. Perspicacité physique: Révèle la structure de la couche limite temporelle et le mécanisme d'escalade
3. Innovation méthodologique: Combinaison ingénieuse de la transformation de Weber et de la méthode de perturbation de domaine
4. Valeur pratique: Fournit une base théorique pour les mesures rhéologiques

Insuffisances

1. Plage de paramètres: Plage de paramètres applicables relativement limitée
2. Comparaison expérimentale: Absence de comparaison directe avec les données expérimentales
3. Vérification numérique: Certains résultats ne sont donnés que par calcul numérique, manquent d'expressions analytiques

Impact

1. Contribution académique: Comble le vide dans la théorie transitoire de l'escalade de tige
2. Valeur méthodologique: Les méthodes d'analyse fournies peuvent être généralisées à d'autres problèmes de fluides complexes
3. Support d'ingénierie: Fournit un support théorique pour la conception de rhéomètres et la caractérisation des matériaux

Scénarios applicables

1. Solutions polymères diluées: Fluides viscoélastiques à viscosité de cisaillement constante tels que les fluides de Boger
2. Faible vitesse de rotation: Conditions expérimentales avec petit nombre de Weissenberg
3. Mesure rhéologique: Détermination du coefficient de contrainte normale et du temps de relaxation

Références

1. Joseph, D.D. & Fosdick, R.L. (1973). The free surface on a liquid between cylinders rotating at different speeds. Arch. Rat. Mech. Anal., 49, 321-380.
2. Oldroyd, J.G. (1950). On the formulation of rheological equations of state. Proc. R. Soc. A, 200, 523-541.
3. Bird, R.B., Armstrong, R.C. & Hassager, O. (1987). Dynamics of Polymeric Liquids, Vol. 1, 2nd edn. John Wiley and Sons.

Évaluation globale: Cet article est un travail théorique de haute qualité en mécanique des fluides, qui analyse pour la première fois systématiquement le comportement transitoire du phénomène d'escalade de tige dans un fluide d'Oldroyd-B. L'article présente des dérivations mathématiques rigoureuses, des perspectives physiques profondes et des méthodes innovantes, apportant des contributions importantes à la théorie des fluides complexes et aux applications rhéologiques. Bien qu'il existe des limitations dans la plage de paramètres, son cadre théorique et ses méthodes d'analyse possèdent une excellente valeur de généralisation.