Mélange dans les fluides viscoélastiques utilisant la turbulence élastique

Informations de base

• ID de l'article: 2409.06391
• Titre: Mixing in viscoelastic fluids using elastic turbulence
• Auteurs: Reinier van Buel, Holger Stark (Technische Universität Berlin)
• Classification: physics.flu-dyn
• Date de publication: Septembre 2024
• Lien de l'article: https://arxiv.org/abs/2409.06391

Résumé

Cet article étudie l'effet de la turbulence élastique sur le mélange du champ de concentration scalaire dans les fluides viscoélastiques au sein d'une géométrie bidimensionnelle de Couette-Taylor, en utilisant des solutions numériques du modèle d'Oldroyd-B. L'état d'écoulement est déterminé par un paramètre d'ordre de flux secondaire, indiquant la région de turbulence élastique. Lors du passage d'un état turbulent et de la réduction subséquente du nombre de Weissenberg, un flux faiblement chaotique apparaît en dessous de Wic. L'advection dans les états de turbulence élastique et de flux faiblement chaotique induit le mélange, illustré par l'évolution temporelle de l'écart-type de la concentration de soluté par rapport à une distribution uniforme. En particulier, le mélange est intense dans l'état de turbulence élastique, quantifié par le taux de mélange, le temps de mélange et l'efficacité de mélange. Ces trois quantités obéissent toutes à des lois d'échelle. Importantly, le paramètre d'ordre est fortement corrélé au taux de mélange, et constitue donc un bon indicateur du mélange au sein du fluide.

Contexte et motivation de la recherche

Contexte du problème

1. Défis du mélange à faible nombre de Reynolds: Dans les fluides newtoniens, les conditions à faible nombre de Reynolds entraînent un écoulement laminaire, et le processus de mélange est lent, contrôlé principalement par la diffusion. Ceci représente une limitation importante dans les dispositifs microfluidiques.
2. Découverte de la turbulence élastique: L'ajout de polymères de haut poids moléculaire au fluide le rend viscoélastique, permettant l'apparition de phénomènes de turbulence élastique même lorsque l'inertie est négligeable.
3. Importance du mélange à l'échelle microscopique: Dans les dispositifs de laboratoire sur puce (lab-on-a-chip), un mélange efficace des fluides est crucial, particulièrement à l'échelle micrométrique.

Motivation de la recherche

• Explorer le potentiel de la turbulence élastique pour améliorer l'efficacité de mélange
• Comprendre la nature physique des mécanismes de mélange dans les fluides viscoélastiques
• Fournir des orientations théoriques pour les applications microfluidiques

Contributions principales

1. Identification de la nature de la transition: Identification de la transition du régime laminaire à la turbulence élastique comme une transition sous-critique, avec découverte de l'existence d'un état de flux faiblement chaotique
2. Établissement de méthodes de caractérisation du mélange: Introduction de trois indicateurs quantitatifs—taux de mélange, temps de mélange et efficacité de mélange—qui obéissent tous à des lois d'échelle
3. Découverte d'une forte corrélation entre le paramètre d'ordre et le mélange: Démonstration que le paramètre d'ordre de flux secondaire est fortement corrélé au taux de mélange et peut servir de bon indicateur du degré de mélange du fluide
4. Vérification du mécanisme de mélange de Batchelor: Observation d'une décroissance exponentielle dans l'état de turbulence élastique, correspondant à la région de mélange de Batchelor, en accord avec les résultats expérimentaux

Détails méthodologiques

Configuration géométrique et équations de gouvernance

L'étude utilise une géométrie bidimensionnelle de Couette-Taylor contenant deux cylindres concentriques:

• Rayon du cylindre externe: $r_o = 10\mu m$
• Rayon du cylindre interne: $r_i = 2.5\mu m$
• Le cylindre externe tourne dans le sens antihoraire à une vitesse angulaire $\Omega = 2\pi s^{-1}$

Le mouvement du fluide est décrit par l'équation de Navier-Stokes généralisée: $\rho\left(\frac{\partial \mathbf{u}}{\partial t} + \mathbf{u} \cdot \nabla \mathbf{u}\right) = -\nabla p + \eta_s \nabla^2 \mathbf{u} + \nabla \cdot \boldsymbol{\tau}$

Modèle constitutif d'Oldroyd-B

Le tenseur de contrainte polymère $\boldsymbol{\tau}$ utilise le modèle d'Oldroyd-B: $\boldsymbol{\tau} + \lambda \stackrel{\nabla}{\boldsymbol{\tau}} = \eta_p[\nabla \otimes \mathbf{u} + (\nabla \otimes \mathbf{u})^T]$

où $\stackrel{\nabla}{\boldsymbol{\tau}}$ est la dérivée convectée supérieure.

Modélisation du mélange

L'évolution du champ de concentration scalaire $c(\mathbf{r}, t)$ est gouvernée par l'équation d'advection-diffusion: $\frac{\partial c}{\partial t} + (\mathbf{u} \cdot \nabla)c = D\nabla^2 c$

Dans cette étude, $D = 0$, négligeant la diffusion moléculaire.

Paramètres adimensionnels clés

• Nombre de Reynolds: $Re = \rho\Omega r_o^2/\eta_s = 0.628 \times 10^{-3}$
• Nombre de Weissenberg: $Wi = \lambda\Omega$
• Rapport de viscosité: $\beta = \eta_p/\eta_s = 1.5$

Configuration expérimentale

Méthode numérique

• Utilisation du solveur de volume fini open-source OpenFOAM®
• Emploi du solveur spécialisé RheoTool pour traiter les écoulements viscoélastiques
• Maillage: grille radiale de 100 éléments, grille en rayon de 120 éléments
• Pas de temps: $\delta t = 10^{-5}s$

Définition du paramètre d'ordre

Définition de l'intensité du flux secondaire comme paramètre d'ordre: $\Phi = \sigma = \sqrt{\langle[\mathbf{u}(\mathbf{r},t) - \mathbf{u}_0(\mathbf{r})]^2\rangle}/u_0^{\max}$

Indicateurs de caractérisation du mélange

Écart-type normalisé: $\Delta c = \sqrt{\frac{\langle[c(\mathbf{r},t) - c_\infty]^2\rangle}{\langle[c(\mathbf{r},0) - c_\infty]^2\rangle}}$

où $c_\infty$ est la concentration à l'état de mélange complet.

Résultats expérimentaux

Transition d'état d'écoulement

1. Transition sous-critique: Découverte que la transition du régime laminaire à la turbulence élastique présente une nature sous-critique, avec présence d'hystérésis
2. Nombre de Weissenberg critique:
   • Branche ascendante: $Wi_c = 9.95$
   • Branche descendante: $Wi_c^I = 9.26$
3. État faiblement chaotique: Observation d'un état d'écoulement faiblement chaotique pour $Wi_c^I$ en dessous

Analyse des caractéristiques de mélange

État de flux faiblement chaotique ($5.03 < Wi \leq 8.17$)

• Décroissance en loi de puissance de l'écart-type, avec exposant proche de $-1/6$
• Spectre de puissance montrant une loi d'échelle avec exposant de $-5$ ou supérieur

État de turbulence élastique ($Wi \geq 9.42$)

• Décroissance initiale en loi de puissance (exposant environ $1/5$), suivie d'une décroissance exponentielle
• Forme de décroissance exponentielle: $\Delta c \sim \exp(-\alpha_{mix}t)$
• Spectre de puissance présentant les caractéristiques du spectre de Batchelor (échelle $k^{-1}$)

Relations de loi d'échelle

1. Temps de mélange: $t_{mix} \sim (Wi - Wi_c^I)^{-0.52}$
2. Concentration de la région II: $c_{II}/c_\infty \sim (Wi - Wi_c^I)^{0.28}$
3. Relation d'échelle: $c_{II} \sim t_{mix}^{-1/2}$

Corrélation entre le paramètre d'ordre et le taux de mélange

Les résultats expérimentaux montrent une forte corrélation entre le paramètre d'ordre $\Phi_d$ et le taux de mélange $\alpha_{mix}$, indiquant que l'augmentation des fluctuations de vitesse améliore directement l'efficacité de mélange du fluide.

Travaux connexes

Recherche sur la turbulence élastique

• Première découverte de la turbulence élastique dans les solutions polymères par Groisman & Steinberg
• Analyses de stabilité linéaire et non-linéaire dans la géométrie de Couette-Taylor
• Observation de motifs d'écoulement hélicoïdaux et en bandes dans les simulations numériques tridimensionnelles

Recherche sur l'amélioration du mélange

• Mélange par écoulement d'extension dans les laminoirs à quatre rouleaux
• Instabilité élastique dans la géométrie de fente croisée
• Amélioration du transfert de chaleur par convection dans les canaux sinueux

Conclusions et discussion

Conclusions principales

1. Nature de la transition: Confirmation de la nature sous-critique de la transition de turbulence élastique dans l'écoulement de Couette-Taylor
2. Mécanisme de mélange: Confirmation que la turbulence élastique peut réaliser un mélange efficace à des nombres de Reynolds extrêmement faibles
3. Région de Batchelor: Observation des caractéristiques de mélange de Batchelor dans l'état de turbulence élastique
4. Indicateur prédictif: Le paramètre d'ordre peut servir d'indicateur fiable de l'efficacité de mélange

Limitations

1. Limitation bidimensionnelle: L'étude est limitée à la géométrie bidimensionnelle; les effets tridimensionnels pourraient produire des résultats différents
2. Dépendance au maillage: Le comportement d'hystérésis disparaît avec un maillage plus fin, indiquant l'influence de la précision numérique
3. Sensibilité aux conditions initiales: Les caractéristiques de mélange dépendent du choix de la distribution de concentration initiale

Directions futures

1. Effets des ondes élastiques: Étude de l'influence des ondes élastiques sur l'efficacité de mélange, le taux de mélange et le temps de mélange
2. Extension tridimensionnelle: Vérification des lois d'échelle découvertes dans les géométries tridimensionnelles
3. Applications pratiques: Exploration du potentiel d'application pratique dans les dispositifs microfluidiques

Évaluation approfondie

Points forts

1. Étude systématique: Analyse complète des caractéristiques de mélange lors de la transition du régime laminaire à la turbulence élastique
2. Caractérisation quantitative: Introduction de multiples indicateurs quantitatifs pour caractériser le degré et l'efficacité de mélange
3. Vérification théorique: Les résultats numériques sont en accord avec la théorie de mélange de Batchelor et les observations expérimentales
4. Valeur pratique: Fournit des orientations théoriques importantes pour les applications microfluidiques

Insuffisances

1. Problèmes de précision numérique: La dépendance au maillage du phénomène d'hystérésis suggère des limitations potentielles de la méthode numérique
2. Compréhension des mécanismes: L'explication du mécanisme physique de la décroissance en loi de puissance dans l'état faiblement chaotique est insuffisante
3. Plage de paramètres: La plage de paramètres étudiée est relativement limitée, manquant d'exploration d'un espace de paramètres plus large

Impact

1. Contribution théorique: Fournit des perspectives importantes pour comprendre les mécanismes de mélange dans les fluides viscoélastiques
2. Perspectives d'application: Possède une importance significative pour la conception de dispositifs microfluidiques
3. Valeur méthodologique: Les méthodes de caractérisation établies peuvent être généralisées à d'autres géométries d'écoulement

Scénarios applicables

• Optimisation du mélange dans les dispositifs microfluidiques
• Contrôle du mélange dans les processus de transformation des polymères
• Mélange à l'échelle microscopique dans les applications biomédicales

Références bibliographiques

L'article cite 57 références connexes, couvrant plusieurs domaines tels que la turbulence élastique, la stabilité des écoulements viscoélastiques et la théorie du mélange, fournissant une base théorique solide pour la recherche.