Mischen in viskoelastischen Flüssigkeiten mittels elastischer Turbulenz

Grundlegende Informationen

• Papier-ID: 2409.06391
• Titel: Mixing in viscoelastic fluids using elastic turbulence
• Autoren: Reinier van Buel, Holger Stark (Technische Universität Berlin)
• Klassifizierung: physics.flu-dyn
• Veröffentlichungsdatum: September 2024
• Papierlink: https://arxiv.org/abs/2409.06391

Zusammenfassung

Diese Arbeit untersucht die Auswirkungen elastischer Turbulenz auf die Mischung von Skalarkonzentrationsfeldern in viskoelastischen Flüssigkeiten in einer zweidimensionalen Taylor-Couette-Geometrie mittels numerischer Lösungen des Oldroyd-B-Modells. Der Strömungszustand wird durch einen Ordnungsparameter der Sekundärströmung bestimmt, der die Region elastischer Turbulenz anzeigt. Beim Übergang vom turbulenten Zustand und anschließender Reduktion der Weissenberg-Zahl erscheint unterhalb von Wi_c eine schwach chaotische Strömung. Die Konvektion sowohl im turbulenten als auch im schwach chaotischen Strömungszustand induziert Mischung, veranschaulicht durch die zeitliche Entwicklung der Standardabweichung der Lösungskonzentration relativ zur gleichmäßigen Verteilung. Insbesondere ist die Mischung im elastisch turbulenten Zustand stark und wird durch die Mischungsrate, die Mischungszeit und die Mischungseffizienz quantifiziert. Alle drei Größen folgen Skalierungsgesetzen. Wichtig ist, dass der Ordnungsparameter stark mit der Mischungsrate korreliert und daher auch ein guter Indikator für die Mischung in der Flüssigkeit ist.

Forschungshintergrund und Motivation

Problemhintergrund

1. Mischungsherausforderungen bei niedrigen Reynolds-Zahlen: In Newtonschen Flüssigkeiten ist die Strömung bei niedrigen Reynolds-Zahlen laminar, und der Mischungsprozess ist langsam und wird hauptsächlich durch Diffusion kontrolliert. Dies ist eine wichtige Einschränkung in Mikrofluidik-Geräten.
2. Entdeckung elastischer Turbulenz: Das Hinzufügen hochmolekularer Polymere zu Flüssigkeiten macht diese viskoelastisch und ermöglicht elastische Turbulenzphänomene, selbst wenn Trägheitskräfte vernachlässigbar sind.
3. Bedeutung von Mischung im Mikromaßstab: In Lab-on-a-Chip-Geräten ist effiziente Flüssigkeitsmischung entscheidend, besonders im Mikrometerbereich.

Forschungsmotivation

• Erkundung des Potenzials elastischer Turbulenz zur Verbesserung der Mischungseffizienz
• Verständnis der physikalischen Natur von Mischungsmechanismen in viskoelastischen Flüssigkeiten
• Bereitstellung theoretischer Richtlinien für Mikrofluidik-Anwendungen

Kernbeiträge

1. Bestimmung der Übergangsnatur: Identifizierung des Übergangs von laminarer Strömung zu elastischer Turbulenz als subkritischen Übergang und Entdeckung des Zustands schwach chaotischer Strömung
2. Etablierung von Mischungscharakterisierungsmethoden: Einführung von drei quantitativen Indikatoren – Mischungsrate, Mischungszeit und Mischungseffizienz – die alle Skalierungsgesetzen folgen
3. Entdeckung starker Korrelation zwischen Ordnungsparameter und Mischung: Nachweis starker Korrelation zwischen dem Ordnungsparameter der Sekundärströmung und der Mischungsrate, der als guter Indikator für den Mischungsgrad in Flüssigkeiten dienen kann
4. Verifizierung des Batchelor-Mischungsmechanismus: Beobachtung exponentiellen Zerfalls im elastisch turbulenten Zustand, entsprechend der Batchelor-Mischungsregion, konsistent mit experimentellen Ergebnissen

Methodische Details

Geometrische Anordnung und Kontrollgleichungen

Die Studie verwendet eine zweidimensionale Taylor-Couette-Geometrie mit zwei konzentrischen Zylindern:

• Außenzylinder-Radius: $r_o = 10\mu m$
• Innenzylinder-Radius: $r_i = 2.5\mu m$
• Außenzylinder rotiert gegen den Uhrzeigersinn mit Winkelgeschwindigkeit $\Omega = 2\pi s^{-1}$

Die Flüssigkeitsbewegung wird durch die verallgemeinerte Navier-Stokes-Gleichung beschrieben: $\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}$

Oldroyd-B-Konstitutivmodell

Der Polymerspannungstensor $\boldsymbol{\tau}$ folgt dem Oldroyd-B-Modell: $\boldsymbol{\tau} + \lambda \stackrel{\nabla}{\boldsymbol{\tau}} = \eta_p[\nabla \otimes \mathbf{u} + (\nabla \otimes \mathbf{u})^T]$

wobei $\stackrel{\nabla}{\boldsymbol{\tau}}$ die obere konvektive Ableitung ist.

Mischungsmodellierung

Die Entwicklung des Skalarkonzentrationsfeldes $c(\mathbf{r}, t)$ wird durch die Advektions-Diffusions-Gleichung kontrolliert: $\frac{\partial c}{\partial t} + (\mathbf{u} \cdot \nabla)c = D\nabla^2 c$

In dieser Studie wird $D = 0$ gesetzt, wodurch molekulare Diffusion vernachlässigt wird.

Wichtige dimensionslose Parameter

• Reynolds-Zahl: $Re = \rho\Omega r_o^2/\eta_s = 0.628 \times 10^{-3}$
• Weissenberg-Zahl: $Wi = \lambda\Omega$
• Viskositätsverhältnis: $\beta = \eta_p/\eta_s = 1.5$

Experimentelle Anordnung

Numerische Methode

• Verwendung des OpenFOAM®-Finite-Volumen-Lösers mit offener Quelle
• RheoTool-Speziallöser für viskoelastische Strömungen
• Netz: Speichennetz mit 100 radialen und 120 azimutalen Elementen
• Zeitschritt: $\delta t = 10^{-5}s$

Definition des Ordnungsparameters

Definition der Sekundärströmungsintensität als Ordnungsparameter: $\Phi = \sigma = \sqrt{\langle[\mathbf{u}(\mathbf{r},t) - \mathbf{u}_0(\mathbf{r})]^2\rangle}/u_0^{\max}$

Mischungscharakterisierungsindikatoren

Normalisierte Standardabweichung: $\Delta c = \sqrt{\frac{\langle[c(\mathbf{r},t) - c_\infty]^2\rangle}{\langle[c(\mathbf{r},0) - c_\infty]^2\rangle}}$

wobei $c_\infty$ die Konzentration im vollständig gemischten Zustand ist.

Experimentelle Ergebnisse

Strömungszustandsübergang

1. Subkritischer Übergang: Der Übergang von laminarer Strömung zu elastischer Turbulenz zeigt subkritische Charakteristiken mit Hysterese
2. Kritische Weissenberg-Zahl:
   • Aufwärtszweig: $Wi_c = 9.95$
   • Abwärtszweig: $Wi_c^I = 9.26$
3. Schwach chaotischer Zustand: Unterhalb von $Wi_c^I$ wird schwach chaotische Strömung beobachtet

Mischungscharakteristiken-Analyse

Schwach chaotischer Strömungszustand ($5.03 < Wi \leq 8.17$)

• Standardabweichung zeigt Potenzgesetz-Zerfall mit Exponent nahe $-1/6$
• Leistungsspektrum zeigt Potenzgesetz-Skalierung mit Exponent $-5$ oder größer

Elastisch turbulenter Zustand ($Wi \geq 9.42$)

• Anfänglicher Potenzgesetz-Zerfall (Exponent etwa $1/5$), gefolgt von exponentiellem Zerfall
• Exponentielle Zerfallsform: $\Delta c \sim \exp(-\alpha_{mix}t)$
• Leistungsspektrum zeigt Batchelor-Spektrum-Charakteristiken ($k^{-1}$-Skalierung)

Skalierungsgesetze

1. Mischungszeit: $t_{mix} \sim (Wi - Wi_c^I)^{-0.52}$
2. Region-II-Konzentration: $c_{II}/c_\infty \sim (Wi - Wi_c^I)^{0.28}$
3. Skalierungsbeziehung: $c_{II} \sim t_{mix}^{-1/2}$

Korrelation zwischen Ordnungsparameter und Mischungsrate

Die experimentellen Ergebnisse zeigen starke Korrelation zwischen dem Ordnungsparameter $\Phi_d$ und der Mischungsrate $\alpha_{mix}$, was darauf hindeutet, dass erhöhte Geschwindigkeitsschwankungen die Mischungseffizienz in der Flüssigkeit direkt verbessern.

Verwandte Arbeiten

Forschung zu elastischer Turbulenz

• Erste Entdeckung elastischer Turbulenz in Polymerlösungen durch Groisman & Steinberg
• Lineare und nichtlineare Stabilitätsanalyse in Taylor-Couette-Geometrie
• Beobachtung spiraliger und bandartiger Strömungsmuster in dreidimensionalen numerischen Simulationen

Forschung zur Mischungsverstärkung

• Streckströmungsmischung in Vier-Walzen-Kalander
• Elastische Instabilität in Kreuzschlitz-Geometrie
• Konvektive Wärmeverstärkung in Serpentinen-Kanälen

Schlussfolgerungen und Diskussion

Hauptschlussfolgerungen

1. Übergangsnatur: Bestätigung der subkritischen Natur des elastischen Turbulenzübergangs in Taylor-Couette-Strömung
2. Mischungsmechanismus: Bestätigung, dass elastische Turbulenz effiziente Mischung bei extrem niedrigen Reynolds-Zahlen ermöglicht
3. Batchelor-Region: Beobachtung von Batchelor-Mischungscharakteristiken im elastisch turbulenten Zustand
4. Vorhersage-Indikator: Der Ordnungsparameter kann als zuverlässiger Vorhersage-Indikator für Mischungseffizienz dienen

Einschränkungen

1. Zweidimensionale Beschränkung: Die Studie ist auf zweidimensionale Geometrie beschränkt; dreidimensionale Effekte könnten unterschiedliche Ergebnisse liefern
2. Netzabhängigkeit: Das Hysterese-Verhalten verschwindet bei feineren Netzen, was auf Auswirkungen numerischer Genauigkeit hindeutet
3. Empfindlichkeit gegenüber Anfangsbedingungen: Mischungscharakteristiken hängen von der Wahl der anfänglichen Konzentrationsverteilung ab

Zukünftige Richtungen

1. Einfluss elastischer Wellen: Untersuchung der Auswirkungen elastischer Wellen auf Mischungseffizienz, Mischungsrate und Mischungszeit
2. Dreidimensionale Erweiterung: Verifizierung der entdeckten Skalierungsgesetze in dreidimensionaler Geometrie
3. Praktische Anwendungen: Erkundung praktischer Anwendungspotenziale in Mikrofluidik-Geräten

Tiefgreifende Bewertung

Stärken

1. Systematische Forschung: Umfassende Analyse der Mischungscharakteristiken während des Übergangsprozesses von laminarer Strömung zu elastischer Turbulenz
2. Quantitative Charakterisierung: Einführung mehrerer quantitativer Indikatoren zur Charakterisierung des Mischungsgrades und der Effizienz
3. Theoretische Verifizierung: Numerische Ergebnisse stimmen mit Batchelor-Mischungstheorie und experimentellen Beobachtungen überein
4. Praktischer Wert: Bietet wichtige theoretische Richtlinien für Mikrofluidik-Anwendungen

Mängel

1. Numerische Genauigkeitsprobleme: Die Netzabhängigkeit des Hysterese-Phänomens deutet auf mögliche Einschränkungen der numerischen Methode hin
2. Mechanismusverständnis: Die physikalische Mechanik des Potenzgesetz-Zerfalls im schwach chaotischen Zustand ist nicht ausreichend erklärt
3. Parameterbereich: Der untersuchte Parameterbereich ist relativ begrenzt und es fehlt eine umfassendere Erkundung des Parameterraums

Auswirkungen

1. Theoretischer Beitrag: Bietet wichtige Erkenntnisse zum Verständnis von Mischungsmechanismen in viskoelastischen Flüssigkeiten
2. Anwendungsaussichten: Von großer Bedeutung für das Design von Mikrofluidik-Geräten
3. Methodologischer Wert: Die etablierten Charakterisierungsmethoden können auf andere Strömungsgeometrien verallgemeinert werden

Anwendungsszenarien

• Mischungsoptimierung in Mikrofluidik-Geräten
• Mischungskontrolle in Polymerverarbeitungsprozessen
• Mischung im Mikromaßstab in biomedizinischen Anwendungen

Literaturverzeichnis

Das Papier zitiert 57 relevante Referenzen, die elastische Turbulenz, Stabilität viskoelastischer Strömungen, Mischungstheorie und andere Aspekte abdecken und eine solide theoretische Grundlage für die Forschung bieten.