On multiple stable states in Taylor-Couette flow with realistic end-wall boundary conditions

Dit onderzoek toont aan dat realistische no-slip randvoorwaarden aan de eindwanden in Taylor-Couette-stroming leiden tot het ontstaan van meerdere stabiele toestanden met hysteresis, waarbij de stroming bij toenemende Reynolds-getal een specifieke reeks structurele overgangen doorloopt die worden gekwantificeerd via een uitgebreid hoekmoment-fluxmodel.

De kern

Titel: De dans van de vloeistof: Waarom één draaisnelheid niet altijd hetzelfde resultaat geeft

Stel je voor dat je een grote, doorzichtige cilinder hebt, zoals een gigantische blikje soep. In het midden zit een kleinere cilinder die ronddraait, terwijl de buitenste cilinder stil blijft staan. De ruimte ertussen is gevuld met een vloeistof, bijvoorbeeld water of olie. Dit systeem heet een Taylor-Couette-stroming. Het klinkt als een ingewikkeld natuurkundig proefje, maar het is eigenlijk heel simpel: hoe gedraagt vloeistof zich als je het laat draaien?

In dit onderzoek kijken wetenschappers van het Max Planck Instituut naar wat er gebeurt als je dit systeem echt bouwt, met deksels boven en onder, in plaats van te doen alsof het oneindig lang is.

Hier is wat ze hebben ontdekt, vertaald naar alledaags taal:

1. Het probleem met de "oneindige" cilinder

Vroeger maakten wetenschappers het zich makkelijk. Ze dachten: "Laten we doen alsof deze cilinder oneindig lang is, zonder boven- en onderkant." Dat is makkelijk om te simuleren op de computer, maar in het echte leven bestaat dat niet. Elke cilinder heeft een bodem en een deksel.

Deze onderzoekers zeiden: "Nee, we doen het echt." Ze bouwden een simulatie met echte, stilstaande deksels boven en onder.

• De analogie: Stel je voor dat je een dansvloer hebt. Als je doet alsof de vloer oneindig groot is, dansen de mensen in een perfect patroon. Maar als je echte muren en een plafond hebt, botsen de dansers tegen de muren. Die muren veranderen de hele dans!

2. De verrassing: Meerdere stabiele danspassen

Het meest fascinerende wat ze vonden, is dat je dezelfde snelheid kunt hebben, maar toch twee totaal verschillende patronen kunt krijgen.

• De analogie: Denk aan een stapel kaarten. Als je de kaarten netjes opstapelt, blijft ze staan. Maar als je ze een beetje schuift, kunnen ze in een heel ander patroon vallen, en ook daar blijven staan.
• In hun experiment bleek dat als je de binnenste cilinder met een bepaalde snelheid laat draaien, de vloeistof in één geval kan kiezen voor een patroon met 18 draaikolken (ronde wervels), en in een ander geval voor 24 draaikolken. Beide patronen zijn stabiel en blijven zo, zolang je de snelheid niet verandert.
• Het hangt er dus van af hoe je begint. Als je de vloeistof rustig laat beginnen, krijg je patroon A. Als je de vloeistof een klein beetje "schudt" (een verstoring toevoegt) voordat je begint, krijg je patroon B.

Dit noemen ze hysteresis: de geschiedenis van het systeem bepaalt de toekomst. Het is alsof je een deur hebt die twee kanten op kan openen; welke kant het opgaat, hangt af van waar je hem vandaan duwt.

3. De rol van de deksels (De "Eindmuur" effecten)

Waarom gebeurt dit? Omdat de deksels boven en onder de vloeistof vasthouden.

• De analogie: Stel je voor dat je een emmer water ronddraait. De vloeistof aan de bodem kan niet meedraaien omdat de bodem stil staat. Dit zorgt voor wrijving en zorgt dat er kleine wervels ontstaan bij de bodem. Deze wervels sturen de hele stroming in de emmer.
• De onderzoekers ontdekten dat deze "wrijving aan de bodem" de vloeistof dwingt om asymmetrische patronen aan te nemen. Soms zijn deze patronen zelfs efficiënter in het transporteren van energie dan de perfecte, symmetrische patronen die je in theorie zou verwachten.

4. Van rustig naar chaotisch (en weer terug)

Ze lieten de binnenste cilinder steeds sneller draaien (van heel rustig tot heel snel).

1. Rustig: De vloeistof maakt mooie, ronde ringen (Taylor-vortexen).
2. Iets sneller: Deze ringen beginnen te wiebelen en golven (Wavy Vortex Flow).
3. Snel: Het wordt chaotisch en turbulent, alsof er een storm in de emmer waait.
4. Zeer snel: En hier komt het verrassende deel: op een gegeven moment wordt het turbulentie weer geordend! De vloeistof vindt weer grote, stabiele patronen, zelfs als het heel snel draait.

Waarom is dit belangrijk?

Dit klinkt misschien als pure theorie, maar het heeft grote gevolgen voor de praktijk:

• Industrie: Veel machines (zoals mengvaten in de chemische industrie of lagers in turbines) werken op dit principe. Als je weet dat je door de startmethode een ander patroon kunt kiezen, kun je machines ontwerpen die efficiënter mengen of minder energie verbruiken.
• Voorspellen: Het laat zien dat complexe systemen (zoals weerpatronen of stroming in pijpleidingen) niet altijd één voorspelbaar resultaat hebben. Soms hangt het resultaat af van kleine details in het begin.

Conclusie

Deze wetenschappers hebben laten zien dat als je een vloeistof laat draaien in een echte cilinder met deksels, de wereld veel interessanter is dan we dachten. Er is niet één "goede" manier om te stromen; er zijn meerdere manieren, en welke je krijgt, hangt af van hoe je begint en hoe de wanden eruitzien. Het is een prachtige herinnering aan dat in de natuur, en vooral in vloeistoffen, de geschiedenis telt.

Technische samenvatting

Titel: Over meervoudige stabiele toestanden in Taylor–Couette-stroming met realistische eindwand-grensvoorwaarden

Auteurs: M. Kriening, Z. Yao, M. S. Emran, J. Song, A. Teimurazov en O. Shishkina
Affiliatie: Max Planck Institute for Dynamics and Self-Organization, Göttingen, Duitsland

---

1. Probleemstelling

Taylor–Couette (TC) stroming, de stroming van een vloeistof tussen twee concentrische cilinders, is een canoniek model voor het bestuderen van hydrodynamische stabiliteit en turbulentie. Traditionele numerieke studies en veel theoretische modellen gebruiken vaak periodieke randvoorwaarden in de axiale richting om de eindwanden te simuleren. Dit benadert een oneindig lange cilinder.

In de praktijk echter hebben TC-systemen fysieke boven- en onderplaten (eindwanden) waar no-slip randvoorwaarden gelden (de vloeistof heeft nul snelheid ten opzichte van de wand). Deze realistische randvoorwaarden introduceren:

• Ekman-laag effecten: Axiale transport van hoekmomentum door de eindwanden.
• Asymmetrie: De stroming wordt axiaal niet-uniform, zelfs in de gemiddelde stroming.
• Meervoudige toestanden: Het bestaan van verschillende, langlevende stromingspatronen (met een verschillend aantal rollen) bij exact dezelfde besturingsparameters (Reynoldsgetal), wat leidt tot hysterese.

De kernvraag is hoe deze realistische eindwanden de dynamica, stabiliteit en het transport van hoekmomentum beïnvloeden ten opzichte van de ideale, periodieke benadering.

2. Methodologie

De auteurs voerden Directe Numerieke Simulaties (DNS) uit met de in-house code goldfish.

• Systeemconfiguratie:
  • Binnenste cilinder roterend, buitenste cilinder stationair.
  • Boven- en onderplaten stationair (geen rotatie).
  • Twee geometrieën onderzocht: Aspectratio $\Gamma = 11$ en $\Gamma = 30$ (met $\Gamma = H/d$).
  • Reynoldsgetal bereik: $90 \le Re \le 7500$.
• Randvoorwaarden en Singulariteiten:
  • Er werden no-slip voorwaarden toegepast op alle oppervlakken.
  • Singulariteitsprobleem: Op de hoek waar de roterende binnenste cilinder de stationaire eindplaat raakt, ontstaat een snelheidsdiscontinuïteit. Om numerieke instabiliteit te voorkomen, werd een gladmakende functie (sigmoid) gebruikt om de hoeksnelheid aan de binnenwand over een kleine axiale afstand ($\epsilon$) te dempen naar nul.
  • De auteurs analyseerden systematisch de invloed van deze gladmakende breedte en kalibreerden de resultaten om overeen te komen met experimentele data (door een vaste verschuiving in $Re$ en een multiplicatieve correctie in de flux toe te passen).
• Initiële Voorwaarden:
  • Om meervoudige toestanden te onderzoeken, werden simulaties gestart met verschillende initieel opgelegde rollen ($n$) en met toegevoegde witte ruis om overgangen tussen toestanden te triggeren.
• Theoretische Uitbreiding:
  • Het klassieke model voor hoekmomentumflux (Eckhardt-Grossmann-Lohse, 2007), dat geldt voor periodieke randen, werd uitgebreid. Er werd een axiale transportterm toegevoegd om de invloed van de Ekman-lagen bij de eindwanden te accounten. Dit leidt tot een totale flux $J_\omega$ die radiaal constant is, in tegenstelling tot de radiale flux alleen.

3. Belangrijkste Bijdragen

1. Theoretische Correctie: Afleiding van een uitgebreide formule voor de totale hoekmomentumflux die axiaal transport door eindwanden incorporeert. Dit verklaart waarom de radiale flux niet constant is in realistische systemen, maar de totale flux wel.
2. Validatie van DNS: Het aantonen dat DNS met realistische randvoorwaarden uitstekende kwantitatieve overeenkomst toont met experimentele data (Martínez-Arias et al., 2014; Ramesh et al., 2019) na correctie voor de numerieke gladmakende effecten.
3. Kartering van Meervoudige Stabiele Toestanden: Een systematische verkenning van de parameter ruimte die aantoont dat er bij identieke Reynoldsgetallen meerdere stabiele toestanden kunnen coëxisteren, afhankelijk van de initiële condities (geschiedenisafhankelijkheid).

4. Resultaten

• Invloed van Randvoorwaarden:

  • Realistische eindwanden veranderen de stromingsdynamica aanzienlijk ten opzichte van periodieke simulaties. Ze stabiliseren asymetrische rollenconfiguraties en verhogen het hoekmomentumtransport in vergelijking met symmetrische toestanden.
  • De grenslaagverhouding ($\delta_o/\delta_i$) wordt beïnvloed door de axiale transporttermen, wat leidt tot een betere overeenkomst met theorie dan de periodieke benadering.

• Meervoudige Stabiele Toestanden en Hysterese:

  • Er werden meervoudige stabiele toestanden gevonden bij dezelfde $Re$. Het systeem kan vastlopen in een staat met een specifiek aantal rollen ($n$), afhankelijk van hoe het systeem is gestart.
  • Hysterese-lussen: Bij het variëren van $Re$ vertoont het systeem hysterese; de overgang van de ene naar de andere toestand gebeurt bij een andere $Re$ dan de terugkeer.
  • Fasenruimte: De toegankelijke fasruimte (het aantal mogelijke stabiele rollen) varieert met $Re$:
    • Bij lage $Re$($<900$): Veel coëxisterende toestanden.
    • Bij middelbare $Re$($900-2000$): De fasruimte krimpt; het systeem convergeert naar één specifieke toestand (vaak door het samensmelten van rollen).
    • Bij hoge $Re$($>2500$): De toegankelijke fasruimte groeit weer aan ("regrowth"). Turbulentie onderdrukt secundaire instabiliteiten, waardoor grote schaal coherente structuren opnieuw stabiel worden.

• Overgangsdynamica:

  • De overgang van Taylor-vortex flow (TVF) naar wavy vortex flow (WVF) en uiteindelijk naar turbulente toestanden verloopt via een cascade van overgangen.
  • Rol-samensmelting: Instabiele configuraties evolueren langzaam (over 200-400 rotaties) door het samensmelten van aangrenzende rollen. Dit gebeurt bij voorkeur in het midden van de cilinder, terwijl de rollen dicht bij de eindwanden door de no-slip condities worden gestabiliseerd.
  • Energiebudget: Analyse van modale energie toont aan dat tijdens overgangen de niet-axiale symmetrische modi tijdelijk energie opbouwen, terwijl de basis-axiale modus ($m=0$) energie verliest voordat de nieuwe configuratie zich stabiliseert.

• Transport Efficiëntie:

  • De efficiëntie van het hoekmomentumtransport hangt af van het aantal rollen. Toestanden met meer rollen (kleinere aspectratio van de rollen) vertonen vaak een hogere transportefficiëntie in het onderzochte regime. Dit verschilt van 2D Rayleigh-Bénard convectie, waar de schalingsexponent onafhankelijk is van de rolconfiguratie.

5. Betekenis en Conclusie

Dit onderzoek benadrukt dat realistische randvoorwaarden essentieel zijn voor een correct begrip van Taylor–Couette stroming, vooral voor industriële toepassingen en het interpreteren van experimenten.

• Fundamenteel: Het toont aan dat TC-stroming, zelfs met deterministische vergelijkingen en vaste parameters, geen unieke turbulente toestand heeft. Het systeem is geschiedenisafhankelijk (multistabiliteit).
• Praktisch: Het bestaan van meerdere stabiele toestanden met verschillende transportefficiënties biedt kansen voor stroomcontrole. Door de initiële condities of perturbaties te manipuleren, kan men mogelijk de toestand selecteren die de beste menging of warmteoverdracht biedt.
• Toekomst: De bevindingen leggen de basis voor verder onderzoek naar de mechanismen die overgangen tussen deze toestanden sturen en de relevantie voor andere complexe, aan wanden gebonden turbulente systemen.

Samenvattend corrigeert deze studie het beeld van Taylor–Couette stroming door de kritische rol van eindwanden en de complexiteit van meervoudige stabiele toestanden bloot te leggen, wat leidt tot een nauwkeurigere theoretische beschrijving en betere overeenkomst met experimenten.