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

Deze studie maakt gebruik van numerieke simulaties van een geroerde tank met mengbare vloeistoffen om aan te tonen dat, hoewel de mengtijd positief correleert met het Richardson-getal, een afgeleide exponentiële schalingswet gebaseerd op het vermogen, het Froude-getal en het Richardson-getal alle gegevens voor intermediaire tot grote dichtheidsverschillen succesvol op één enkele mastercurve reduceert.

De kern

Stel je voor dat je probeert een gigantische tank met twee verschillende vloeistoffen te mengen: een zware, dikke stroop onderaan en een lichtere, dunner sap bovenop. Je laat een gigantisch roterend roerblad (een impeller) vallen om ze samen te mengen.

In de echte wereld is dit een veelvoorkomende taak in fabrieken die van alles produceren, van medicijnen tot afvalwaterbehandeling. Maar hier zit de adder onder het gras: omdat de vloeistoffen verschillende gewichten (dichtheden) hebben, wil de zware onderaan blijven en de lichte bovenaan drijven. Dit creëert een "strijd" tussen het roterende roerblad dat ze probeert te mengen en de zwaartekracht die probeert ze gescheiden te houden.

Dit artikel is als een detectiveverhaal waarin wetenschappers krachtige computersimulaties gebruikten om precies uit te vinden hoe lang het duurt voordat deze twee vloeistoffen perfect gemengd zijn, en hoe die tijd te voorspellen zonder elke keer een fysieke tank te bouwen en dure tests uit te voeren.

De Opzet: Een Digitale Testkeuken

De onderzoekers bouwden een virtuele versie van een standaard industriële mengtank.

• De Tank: Het is een grote cilinder met wanden en vier verticale vinnen (geleiders) om te voorkomen dat de vloeistof gewoon in een cirkel draait als een lui riviertje.
• Het Roerblad: Een roterend mes in het midden.
• De Vloeistoffen: Ze simuleerden een 50/50 mengsel van een zware vloeistof en een lichtere vloeistof. Ze gebruikten geen echte chemicaliën; ze behandelden ze gewoon als "zware" en "lichte" vloeistoffen met dezelfde dikte (viscositeit).
• De Methode: In plaats van standaard wiskundige vergelijkingen te gebruiken, gebruikten ze een slimme truc genaamd de Lattice Boltzmann-methode. Denk hierbij aan het simuleren van de vloeistof niet als een continue klomp, maar als biljoenen tiny, onzichtbare biljartballen die rondhollen en botsen. Dit stelde hen in staat om precies te zien hoe de turbulentie (het chaotische wervelen) zich gedroeg.

De Grote Vraag: Hoe Snel Kunnen We Maken?

Het hoofddoel was een "magische formule" te vinden om de mengtijd te voorspellen.

• De Variabelen: Ze veranderden twee hoofdzaakken:
  1. Hoe snel het roerblad draait (Reynoldsgetal): Sneller draaien betekent meestal meer turbulentie en snellere menging.
  2. Hoe verschillend de gewichten zijn (Richardsongetal): Als de vloeistoffen bijna hetzelfde gewicht hebben, mengen ze gemakkelijk. Als één veel zwaarder is, vecht de zwaartekracht tegen het mengen, waardoor lagen ontstaan die moeilijk te breken zijn.

De Ontdekking: De "Zwaartekracht versus Spin"-Strijd

De onderzoekers vonden enkele interessante patronen:

1. Wanneer Zwaartekracht Niet Uitmaakt (Zelfde Gewichten):
   Als de twee vloeistoffen exact hetzelfde gewicht hebben, is de mengtijd verrassend consistent. Het maakt niet uit hoe snel je het roerblad draait (binnen een bepaald bereik), de "dimensieloze mengtijd" (een chique manier om te zeggen "hoeveel omwentelingen van het roerblad het kost") blijft constant op ongeveer 20 omwentelingen. Het is als een natuurwet: zodra het water genoeg kolt, maakt sneller draaien het niet sneller gemengd in termen van roerblaad-omwentelingen.

2. Wanneer Zwaartekracht Terugvecht (Verschillende Gewichten):
   Zodra de vloeistoffen verschillende gewichten hebben, wil de zware onderaan blijven. Hoe groter het gewichtsverschil, hoe moeilijker het is om te mengen.

   • De Trend: Hoe groter het gewichtsverschil, hoe langer het duurt om te mengen.
   • De Verrassende Wending: Als je het "gewichtsverschil" constant houdt en gewoon het roerblad sneller draait, gaat de mengtijd niet altijd omlaag. Soms duurt het juist langer om een specifiek punt van menging te bereiken als je sneller draait.
   • Waarom? Stel je voor dat de zware vloeistof als een dikke deken is. Als je het roerblad te snel draait, creëer je veel energie, maar vormt de zware vloeistof een stabiele "kap" die de lichtere vloeistof niet kan penetreren. De energie wordt verspild aan het wervelen van de bovenste laag terwijl de onderste laag afgesloten blijft. Het is alsof je een pan soep probeert te roeren waarbij de zware groenten tot een stevig blok onderaan zijn gezakt; het sneller draaien van de lepel plakt alleen de bouillon bovenop zonder het groenteblok te breken.

De Oplossing: Een Nieuwe "Meestercurve"

De grootste prestatie van het team was het creëren van één enkele, simpele formule die al deze factoren combineert. Ze ontdekten dat als je de mengtijd bekijkt door de lens van drie specifieke getallen (Vermogen, Froudgetal en Richardsongetal), al hun rommelige datapunten samenvallen tot één gladde, exponentiële curve.

Denk hierbij aan dit: Vroeger moesten ingenieurs raden of honderden tests uitvoeren om te zien hoe een nieuwe vloeistof zou mengen. Nu hebben ze een "recept". Als je hen het gewichtsverschil vertelt en hoe snel je van plan bent te draaien, voorspelt deze formule de mengtijd met hoge nauwkeurigheid.

De Conclusie

Het artikel concludeert dat voor deze specifieke industriële tanks:

• Turbulentie is cruciaal: Zodra de vloeistof volledig kolkt, is het menggedrag voorspelbaar.
• Zwaartekracht is de baas: Als de vloeistoffen verschillende dichtheden hebben, creëert de zwaartekracht een "stratificatie" (laagvorming) die het mengen weerstaat.
• Sneller is niet altijd beter: In systemen met grote dichtheidsverschillen garandeert het simpelweg harder draaien van de motor niet snellere menging; soms creëert het juist een stabielere scheiding.

De auteurs leveren deze nieuwe formule om ingenieurs te helpen betere mengprocessen te ontwerpen zonder eerst dure prototypes te hoeven bouwen. Ze plannen om deze formule in de toekomst te testen op verschillende tankvormen en roerbladtypes, maar voor nu werkt het perfect voor de standaardtanks die ze simuleerden.

Technische samenvatting

Technische Samenvatting: Dimensieloze Schaling voor Intermediaire tot Grote Dichtheidsverschillen in een Geroerde Tank

Probleemstelling
Het mengen van mengbare vloeistoffen is een kritieke eenheidsoperatie in industriële processen, variërend van farmaceutische productie tot afvalwaterbehandeling. Hoewel turbulente menging in éénfasige, homogene systemen goed gekarakteriseerd is, introduceert de aanwezigheid van dichtheidsverschillen tussen vloeistoffen door opwaartse krachten gedreven stratificatie. Dit fenomeen creëert verticaal gelaagde gebieden waar zwaartekrachtsterken convectief transport tegenwerken, wat het verticale mengen potentieel onderdrukt en de tijd verlengt die nodig is om homogeniteit te bereiken. Hoewel empirische correlaties bestaan voor éénfasige tracer-menging en specifieke gestratificeerde gevallen, ontbreekt een algemeen geldig raamwerk dat consistent rekening houdt met zwaartekracht-, traagheids-, viskeuze en opwaartse krachten-effecten over verschillende stromingsregimes en geometrieën. Deze studie adresseert de uitdaging om mengtijden te voorspellen in afgeleide geroerde tanks die twee mengbare vloeistoffen bevatten met intermediaire tot grote dichtheidsverschillen.

Methodologie
De auteurs gebruikten Computational Fluid Dynamics (CFD) met behulp van de Lattice Boltzmann-methode (LBM) om turbulente menging in een afgeleide geroerde tank te simuleren. Het systeem bestond uit een volumeverhouding van 50/50 van twee mengbare vloeistoffen met identieke dynamische viscositeiten maar variërende dichtheden ($\rho_l$ en $\rho_h$), waarbij de lichtere vloeistof aanvankelijk gestratificeerd boven de zwaardere vloeistof zat.

• Numerieke Aanpak: De simulaties losten de incompressibele Navier-Stokes-vergelijkingen op met behulp van een D3Q19 rooster-schema. Om hoge Reynolds-getallen ($Re > 10^5$) te hanteren, werd een Smagorinsky-Lilly Large Eddy Simulation (S-L LES) model geïmplementeerd om sub-rooster turbulentie te modelleren.
• Scalar Transport: Er werd een op massavlakte gebaseerd scalar transportmodel aangenomen om de concentratie van de lichtere vloeistof te volgen. Dit model koppelde het scalarveld aan de impulsvergelijking via een concentratie-afhankelijke zwaartekrachtskracht, waardoor de simulatie van opwaartse krachteffecten die voortvloeien uit lokale dichtheidsvariaties mogelijk werd.
• Simulatieruimte: Een systematische parameterstudie werd uitgevoerd door de rotatiesnelheid van de roerder ($N$) en de dichtheid van de lichte fase ($\rho_l$) te variëren. Dit resulteerde in een breed scala aan dimensieloze getallen: Reynolds-getallen ($Re$) van 50.000 tot 125.000 en Richardson-getallen ($Ri$) van 0 tot 14.
• Validatie: De numerieke nauwkeurigheid werd geverifieerd via een Grid Convergence Index (GCI) analyse, die numerieke onzekerheid schatte op ongeveer 2,39% op het productie-rooster. Bovendien werden simulatieresultaten gevalideerd tegen experimentele tracerdata uit de literatuur, waarbij goede overeenkomst werd getoond in proefrespons en definities van mengtijden.

Belangrijkste Resultaten
De studie kwantificeerde de mengtijd ($t_m$), gedefinieerd als de tijd nodig om 90% homogeniteit te bereiken, en drukte deze uit in dimensieloze vorm ($t^*_m = N \cdot t_m$).

1. Homogene Baseline ($Ri = 0$): Bij afwezigheid van dichtheidsverschillen klapte de dimensieloze mengtijd samen tot een constante waarde van $t^*_m \approx 20$ over het hele onderzochte Reynolds-getalbereik. Dit bevestigt dat de stroming volledig turbulent is en dat traagheidsmengen overheerst, in overeenstemming met gevestigde literatuur voor afgeleide tanks.
2. Effect van Stratificatie ($Ri > 0$): Naarmate het Richardson-getal toenam, wat grotere dichtheidsverschillen aangeeft, nam de dimensieloze mengtijd systematisch toe. De data volgden een duidelijke exponentiële trend over alle Reynolds-getallen.
3. Voorgestelde Schalingswet: De auteurs leidden een mastercurve af die alle simulatiedata samenvoegt tot één enkele exponentiële schalingsrelatie:
   $$t^*_m = 20 \cdot e^{0.58 \cdot N_p \cdot Fr \cdot Ri}$$
   Waar $N_p$ het vermogensgetal is, $Fr$ het Froude-getal, en $Ri$ het Richardson-getal. De pre-factor van 20 vertegenwoordigt de traagheidsbaseline, terwijl de exponentiële term het stabiliserende effect van opwaartse krachten vastlegt.
4. Niet-monotoon Gedrag: Bij analyse van dimensionale mengtijden ($t_m$) bij constante $Ri$, observeerde de studie een niet-monotone afhankelijkheid van $Re$ voor hogere Richardson-getallen ($Ri \geq 8$). Hoewel het verhogen van $Re$ (en dus de roersnelheid) aanvankelijk de mengtijd verkort door turbulentie te versterken, leiden verdere verhogingen tot sterkere dichtheidsstratificatie (door de koppeling van parameters in de simulatieopstelling), wat verticaal transport onderdrukt. Deze wisselwerking resulteert in een minimale mengtijd bij intermediaire Reynolds-getallen.
5. Stromingsstructuur: Analyse van het snelheidsveld onthulde dat bij hoge $Ri$ de lichtere fase kan worden "afgesloten" van de vermogensinvoer van de roerder, waarbij verticale structuren voornamelijk in de dichte vloeistof verschijnen. Power Spectral Density (PSD) analyse gaf aan dat hoewel opwaartse krachten de ruimtelijke organisatie van de stroming beïnvloeden, de door de roerder gedreven dynamiek nog steeds de energiecascade domineert.

Betekenis en Claims
Het artikel claimt een compacte, voorspellende correlatie te bieden voor turbulente mengtijden in aanvankelijk door dichtheid gestratificeerde systemen binnen het volledig turbulente regime. De afgeleide schalingswet, gebaseerd op de gecombineerde invloed van het Vermogens-, Froude- en Richardson-getal, biedt een eenvoudige analytische uitdrukking voor ingenieurs en toepassingsdeskundigen.

De auteurs merken op dat, voor zover hen bekend, geen vergelijkbare correlatie voor gestratificeerde systemen in het volledig turbulente regime eerder is gerapporteerd. De betekenis van dit werk ligt in het vermogen om het ontwerp van mengprocessen te stroomlijnen door de afhankelijkheid van dure real-world testen te verminderen. Het voorgestelde model maakt het mogelijk mengtijden te schatten in systemen waar opwaartse krachteffecten significant zijn, en legt de complexe concurrentie tussen versterkte turbulentie en versterkte stratificatie vast. De auteurs suggereren bescheiden dat toekomstig werk zal bestaan uit het testen van deze correlatie op verschillende tankgeometrieën en roerertypes, inclusief schuine roerders, om de robuustheid van de voorgestelde schaling te evalueren.