Mixing of miscible liquids: Dimensionless scaling for intermediate-to-large density differences in a stirred tank

Diese Studie nutzt numerische Simulationen eines Rührbehälters mit mischbaren Flüssigkeiten, um zu zeigen, dass die Mischzeit zwar positiv mit der Richardson-Zahl korreliert, eine abgeleitete exponentielle Skalierungsgesetzgebung auf Basis von Leistung, Froude- und Richardson-Zahl jedoch alle Daten für mittlere bis große Dichteunterschiede erfolgreich auf eine einzige Masterkurve zusammenführt.

Das Wesentliche

Stellen Sie sich vor, Sie versuchen, einen riesigen Kessel mit zwei verschiedenen Flüssigkeiten zu mischen: einen schweren, dicken Sirup am Boden und einen leichteren, dünneren Saft oben. Sie werfen ein riesiges, sich drehendes Rührwerk (ein Laufrad) hinein, um sie zu vermischen.

In der realen Welt ist dies eine gängige Aufgabe in Fabriken, die alles von Medikamenten bis hin zur Abwasserbehandlung herstellen. Aber hier liegt der Haken: Da die Flüssigkeiten unterschiedliche Gewichte (Dichten) haben, möchte die schwere unten bleiben und die leichte oben schwimmen. Dies erzeugt einen „Kampf" zwischen dem sich drehenden Rührwerk, das sie mischen will, und der Schwerkraft, die versucht, sie getrennt zu halten.

Dieser Artikel ist wie eine Detektivgeschichte, in der Wissenschaftler leistungsstarke Computersimulationen nutzten, um genau herauszufinden, wie lange es dauert, bis diese beiden Flüssigkeiten perfekt gemischt sind, und wie man diese Zeit vorhersagen kann, ohne jedes Mal einen physischen Tank zu bauen und teure Tests durchzuführen.

Das Setup: Eine digitale Testküche

Die Forscher bauten eine virtuelle Version eines standardmäßigen industriellen Mischbehälters.

• Der Tank: Ein großer Zylinder mit Wänden und vier vertikalen Flossen (Leitblechen), um zu verhindern, dass die Flüssigkeit nur wie ein träger Fluss in einer Kreiselbewegung rotiert.
• Das Rührwerk: Eine sich drehende Schaufel in der Mitte.
• Die Flüssigkeiten: Sie simulierten eine 50/50-Mischung aus einer schweren und einer leichteren Flüssigkeit. Sie verwendeten keine echten Chemikalien; sie behandelten sie einfach als „schwere" und „leichte" Fluide mit der gleichen Zähigkeit (Viskosität).
• Die Methode: Anstatt Standard-Mathematikgleichungen zu verwenden, nutzten sie einen cleveren Trick namens Gitter-Boltzmann-Methode. Stellen Sie sich dies vor wie die Simulation der Flüssigkeit nicht als kontinuierliche Masse, sondern als Milliarden winziger, unsichtbarer Billardkugeln, die herumhüpfen und kollidieren. Dies ermöglichte ihnen, genau zu sehen, wie sich die Turbulenz (das chaotische Wirbeln) verhielt.

Die große Frage: Wie schnell können wir mischen?

Das Hauptziel war es, eine „magische Formel" zu finden, um die Mischzeit vorherzusagen.

• Die Variablen: Sie veränderten zwei Hauptfaktoren:
  1. Wie schnell sich das Rührwerk dreht (Reynolds-Zahl): Schnelleres Drehen bedeutet normalerweise mehr Turbulenz und schnellere Mischung.
  2. Wie unterschiedlich die Gewichte sind (Richardson-Zahl): Wenn die Flüssigkeiten fast das gleiche Gewicht haben, mischen sie sich leicht. Wenn eines viel schwerer ist, bekämpft die Schwerkraft die Mischung und erzeugt Schichten, die schwer zu durchbrechen sind.

Die Entdeckung: Die Schlacht „Schwerkraft gegen Drehung"

Die Forscher fanden einige interessante Muster:

1. Wenn die Schwerkraft keine Rolle spielt (gleiches Gewicht):
   Wenn die beiden Flüssigkeiten exakt das gleiche Gewicht haben, ist die Mischzeit überraschend konsistent. Unabhängig davon, wie schnell Sie das Rührwerk drehen (innerhalb eines bestimmten Bereichs), bleibt die „dimensionslose Mischzeit" (eine ausgefallene Art zu sagen, „wie viele Umdrehungen des Rührwerks es braucht") konstant bei etwa 20 Umdrehungen. Es ist wie ein Naturgesetz: Sobald das Wasser ausreichend wirbelt, führt schnelleres Drehen nicht dazu, dass es sich in Bezug auf Rührwerksumdrehungen schneller mischt.

2. Wenn die Schwerkraft sich zur Wehr setzt (unterschiedliche Gewichte):
   Sobald die Flüssigkeiten unterschiedliche Gewichte haben, möchte die schwere unten bleiben. Je größer der Gewichtsunterschied, desto schwieriger ist die Mischung.

   • Der Trend: Je unterschiedlicher die Gewichte sind, desto länger dauert es, bis sie gemischt sind.
   • Die überraschende Wendung: Wenn Sie den „Gewichtsunterschied" konstant halten und das Rührwerk einfach schneller drehen, geht die Mischzeit nicht immer zurück. Manchmal führt schnelleres Drehen tatsächlich dazu, dass es länger dauert, einen bestimmten Mischgrad zu erreichen.
   • Warum? Stellen Sie sich vor, die schwere Flüssigkeit ist wie ein dicker Deckel. Wenn Sie das Rührwerk zu schnell drehen, erzeugen Sie viel Energie, aber die schwere Flüssigkeit bildet eine stabile „Kappe", die die leichtere Flüssigkeit nicht durchdringen kann. Die Energie wird verschwendet, indem die obere Schicht gewirbelt wird, während die untere Schicht abgedichtet bleibt. Es ist, als würde man versuchen, einen Topf Suppe zu rühren, in dem sich das schwere Gemüse zu einem festen Block am Boden gesetzt hat; wenn man den Löffel schneller dreht, spritzt man nur die Brühe oben herum, ohne den Gemüseblock zu zerbrechen.

Die Lösung: Eine neue „Master-Kurve"

Die größte Leistung des Teams bestand darin, eine einzelne, einfache Formel zu erstellen, die all diese Faktoren kombiniert. Sie stellten fest, dass, wenn man die Mischzeit durch die Linse von drei spezifischen Zahlen betrachtet (Leistung, Froude-Zahl und Richardson-Zahl), all ihre unübersichtlichen Datenpunkte auf eine glatte, exponentielle Kurve zusammenfallen.

Stellen Sie es sich so vor: Früher mussten Ingenieure raten oder Hunderte von Tests durchführen, um zu sehen, wie sich eine neue Flüssigkeit mischen würde. Jetzt haben sie ein „Rezept". Wenn Sie ihnen den Gewichtsunterschied und die geplante Drehgeschwindigkeit nennen, sagt diese Formel die Mischzeit mit hoher Genauigkeit voraus.

Das Fazit

Der Artikel kommt zu dem Schluss, dass für diese spezifischen industriellen Tanks:

• Turbulenz ist der Schlüssel: Sobald die Flüssigkeit vollständig wirbelt, ist das Mischverhalten vorhersagbar.
• Die Schwerkraft ist der Boss: Wenn die Flüssigkeiten unterschiedliche Dichten haben, erzeugt die Schwerkraft eine „Schichtung" (Stratifizierung), die der Mischung widersteht.
• Schneller ist nicht immer besser: In Systemen mit großen Dichteunterschieden garantiert das bloße Hochdrehen der Motorgeschwindigkeit nicht eine schnellere Mischung; manchmal erzeugt es nur eine stabilere Trennung.

Die Autoren stellen diese neue Formel bereit, um Ingenieuren zu helfen, bessere Mischprozesse zu entwickeln, ohne zuerst teure Prototypen bauen zu müssen. Sie planen, diese Formel in Zukunft an verschiedenen Tankformen und Rührwerkstypen zu testen, aber vorläufig funktioniert sie perfekt für die standardmäßigen Tanks, die sie simuliert haben.

Technische Zusammenfassung

Technische Zusammenfassung: Dimensionslose Skalierung für mittlere bis große Dichteunterschiede in einem gerührten Behälter

Problemstellung
Das Mischen mischbarer Flüssigkeiten ist ein kritischer verfahrenstechnischer Schritt in industriellen Prozessen, von der pharmazeutischen Produktion bis zur Abwasserbehandlung. Während turbulente Mischung in einphasigen, homogenen Systemen gut charakterisiert ist, führt das Vorhandensein von Dichteunterschieden zwischen Fluiden zu einer durch Auftrieb verursachten Schichtung. Dieses Phänomen erzeugt vertikal geschichtete Bereiche, in denen Gravitationskräfte der konvektiven Strömung entgegenwirken, was die vertikale Durchmischung potenziell unterdrückt und die Zeit verlängert, die zur Erreichung von Homogenität erforderlich ist. Obwohl empirische Korrelationen für einphasiges Tracermischen und spezifische geschichtete Fälle existieren, bleibt ein allgemein gültiges Rahmenwerk, das Gravitations-, Trägheits-, Viskositäts- und Auftriebseffekte über verschiedene Strömungsregime und Geometrien hinweg konsistent berücksichtigt, weiterhin schwer fassbar. Diese Studie adressiert die Herausforderung, Mischzeiten in baffleierten Rührbehältern vorherzusagen, die zwei mischbare Flüssigkeiten mit mittleren bis großen Dichteunterschieden enthalten.

Methodik
Die Autoren setzten Computational Fluid Dynamics (CFD) unter Verwendung der Gitter-Boltzmann-Methode (LBM) ein, um turbulente Mischung in einem baffleierten Rührbehälter zu simulieren. Das System bestand aus einem Volumenverhältnis von 50/50 zweier mischbarer Flüssigkeiten mit identischer dynamischer Viskosität, jedoch unterschiedlichen Dichten ($\rho_l$ und $\rho_h$), wobei die leichtere Flüssigkeit anfänglich über der schwereren Flüssigkeit geschichtet war.

• Numerischer Ansatz: Die Simulationen lösten die inkompressiblen Navier-Stokes-Gleichungen unter Verwendung eines D3Q19-Gitterschemas. Um hohe Reynolds-Zahlen ($Re > 10^5$) zu handhaben, wurde ein Smagorinsky-Lilly-Großwirbel-Simulationsmodell (S-L LES) implementiert, um Subgitter-Turbulenz zu modellieren.
• Skalartransport: Ein skalarer Transportmodell auf Basis des Massenflusses wurde eingeführt, um die Konzentration der leichteren Flüssigkeit zu verfolgen. Dieses Modell koppelt das Skalarfeld über eine konzentrationsabhängige gravitative Volumenkraft an die Impulsgleichung, was die Simulation von Auftriebseffekten ermöglicht, die aus lokalen Dichteschwankungen resultieren.
• Simulationsraum: Eine systematische Parameterstudie wurde durchgeführt, indem die Rotationsgeschwindigkeit des Rührers ($N$) und die Dichte der leichten Phase ($\rho_l$) variiert wurden. Dies ergab einen breiten Bereich dimensionsloser Zahlen: Reynolds-Zahlen ($Re$) von 50.000 bis 125.000 und Richardson-Zahlen ($Ri$) von 0 bis 14.
• Validierung: Die numerische Genauigkeit wurde durch eine Grid Convergence Index (GCI)-Analyse verifiziert, die die numerische Unsicherheit auf dem Produktionsgitter auf etwa 2,39 % schätzte. Zusätzlich wurden die Simulationsergebnisse gegen experimentelle Tracerdaten aus der Literatur validiert, wobei eine gute Übereinstimmung bei der Sondenantwort und den Definitionen der Mischzeit festgestellt wurde.

Hauptergebnisse
Die Studie quantifizierte die Mischzeit ($t_m$), definiert als die Zeit, die erforderlich ist, um 90 % Homogenität zu erreichen, und drückte sie dimensionslos aus ($t^*_m = N \cdot t_m$).

1. Homogener Basisfall ($Ri = 0$): In Abwesenheit von Dichteunterschieden kollabierte die dimensionslose Mischzeit über den gesamten untersuchten Reynolds-Zahlenbereich auf einen konstanten Wert von $t^*_m \approx 20$. Dies bestätigt, dass die Strömung vollständig turbulent ist und die Trägheitsmischung dominiert, was mit etablierter Literatur für baffleierte Behälter übereinstimmt.
2. Einfluss der Schichtung ($Ri > 0$): Mit zunehmender Richardson-Zahl, was auf größere Dichteunterschiede hindeutet, nahm die dimensionslose Mischzeit systematisch zu. Die Daten folgten über alle Reynolds-Zahlen hinweg einem klaren exponentiellen Trend.
3. Vorgeschlagene Skalierungsgesetz: Die Autoren leiteten eine Masterkurve ab, die alle Simulationsdaten auf eine einzige exponentielle Skalierungsbeziehung kollabiert:
   $$t^*_m = 20 \cdot e^{0.58 \cdot N_p \cdot Fr \cdot Ri}$$
   Dabei ist $N_p$ die Leistungszahl, $Fr$ die Froude-Zahl und $Ri$ die Richardson-Zahl. Der Vorfaktor 20 repräsentiert die Trägheitsbasis, während der exponentielle Term den stabilisierenden Effekt des Auftriebs erfasst.
4. Nicht-monotones Verhalten: Bei der Analyse der dimensionsbehafteten Mischzeiten ($t_m$) bei konstantem $Ri$ beobachtete die Studie eine nicht-monotone Abhängigkeit von $Re$ für höhere Richardson-Zahlen ($Ri \geq 8$). Während eine Erhöhung von $Re$ (und somit der Rührerdrehzahl) die Mischzeit zunächst durch verstärkte Turbulenz reduziert, führen weitere Erhöhungen zu einer stärkeren Dichteschichtung (aufgrund der Kopplung der Parameter im Simulationsaufbau), was den vertikalen Transport unterdrückt. Dieses Zusammenspiel führt zu einer minimalen Mischzeit bei mittleren Reynolds-Zahlen.
5. Strömungsstruktur: Die Analyse des Geschwindigkeitsfeldes ergab, dass bei hohen $Ri$ die leichtere Phase vom Energieeintrag des Rührers „abgeschnitten" werden kann, wobei vertikale Strukturen primär in der dichteren Flüssigkeit auftreten. Die Analyse der Leistungsdichtespektren (PSD) zeigte, dass, obwohl der Auftrieb die räumliche Organisation der Strömung beeinflusst, die vom Rührer angetriebene Dynamik den Energiekaskadenprozess weiterhin dominiert.

Bedeutung und Behauptungen
Die Arbeit behauptet, eine kompakte, prädiktive Korrelation für turbulente Mischzeiten in anfänglich dichte-geschichteten Systemen im vollständig turbulenten Regime bereitzustellen. Das abgeleitete Skalierungsgesetz, basierend auf dem kombinierten Einfluss von Leistungszahl, Froude-Zahl und Richardson-Zahl, bietet Ingenieuren und Anwendungsexperten einen einfachen analytischen Ausdruck.

Die Autoren stellen fest, dass ihnen nach Kenntnisstand keine vergleichbare Korrelation für geschichtete Systeme im vollständig turbulenten Regime zuvor berichtet wurde. Die Bedeutung dieser Arbeit liegt in ihrer Fähigkeit, die Auslegung von Mischprozessen zu straffen, indem die Abhängigkeit von teuren realen Tests reduziert wird. Das vorgeschlagene Modell ermöglicht die Schätzung von Mischzeiten in Systemen, in denen Auftriebseffekte signifikant sind, und erfasst das komplexe Wettspiel zwischen verstärkter Turbulenz und verstärkter Schichtung. Die Autoren schlagen bescheiden vor, dass zukünftige Arbeiten den Test dieser Korrelation an verschiedenen Behältergeometrien und Rührertypen, einschließlich schräger Rührer, umfassen werden, um die Robustheit der vorgeschlagenen Skalierung zu bewerten.