Laminar-to-Turbulent Transition of Yield-Stress Fluids in Pipe and Channel Flows

Dit artikel presenteert directe numerieke simulaties van de overgang van laminaire naar turbulente stroming voor Herschel-Bulkley-vloeistoffen in pijpen en kanalen, waarbij wordt aangetoond dat turbulentie optreedt wanneer lokale Reynolds-spanningen de vloeigrens overschrijden en een overgangsregime wordt geïdentificeerd dat overeenkomt met experimentele trends.

De kern

De Grote Strijd: Wanneer Vloeistoffen van Rust naar Chaos Overschakelen

Stel je voor dat je een grote bak honing of een dikke soep hebt. Als je die rustig laat staan, gedraagt het zich als een stevige gel: het beweegt niet, het houdt zijn vorm. Maar zodra je er hard genoeg op duwt (of er een lepel doorheen roert), wordt het plotseling vloeibaar en stroomt het. Dit is wat wetenschappers een Herschel-Bulkley-vloeistof noemen. Denk aan tandpasta, verf, of zelfs modder.

In dit onderzoek kijken we naar wat er gebeurt als je deze "dikke" vloeistoffen door een pijp of een kanaal stuurt. De vraag is: Wanneer stopt het met stromen als een rustige, gladde stroom en begint het te borrelen en te draaien als een wilde, turbulente stroom?

Hier is wat de onderzoekers hebben ontdekt, vertaald naar alledaagse taal:

1. De "Kern" die Weigert te Bewegen

In een rustige stroom (we noemen dit laminaire stroming) vormt deze vloeistof een kern of een "plug" in het midden van de pijp.

• De Analogie: Denk aan een ijsblokje in een glas water. Het water stroomt eromheen, maar het ijsblokje zelf beweegt niet mee; het glijdt gewoon als een blok. Bij deze vloeistoffen is het midden van de pijp zo'n "ijsblokje". Het is te dik en te stijf om te vervormen, dus het glijdt als een star blokje door de pijp, terwijl alleen de vloeistof vlak tegen de wanden in de war raakt.

2. De Overgang: De "Krakende" Kern

Naarmate je harder pompt (de snelheid omhoog gaat), gebeurt er iets fascinerends. De druk wordt zo groot dat de "ijsblokken" in het midden beginnen te barsten.

• Het Moment van Chaos: Er komt een punt waarop de kracht van de stroming groot genoeg is om de stijfheid van de vloeistof te breken. De kern begint te smelten en te vervormen.
• De Turbulentie: Zodra de kern breekt, ontstaan er kleine draaikolken (wirrels) langs de wanden. Deze draaikolken groeien en versmelten, totdat de hele pijp vol zit met een chaotische dans van vloeistof. Dit is de turbulente overgang.

3. Wat hebben de computersimulaties laten zien?

De onderzoekers gebruikten superkrachtige computers om dit proces in detail te bekijken (zonder echte vloeistof, puur wiskunde). Ze keken naar twee situaties:

1. Ronde pijpen (zoals waterleidingen).
2. Rechthoekige kanalen (zoals een afvoerbak).

Ze ontdekten drie duidelijke fasen, afhankelijk van hoe hard ze pompen:

• Fase 1: De Rustige Kern (Laminaar). Alles is stil. Het midden is een star blokje. Geen chaos.
• Fase 2: De Onzekere Overgang. De druk is net hoog genoeg om de kern te laten wankelen. Hier zie je soms rust, soms plotseling een explosie van draaikolken. Het is een onstabiele fase.
• Fase 3: De Wilde Dans (Turbulent). De kern is volledig weggesmolten. De hele vloeistof is nu in beweging, met veel draaikolken en wervelingen.

4. De "Magische Drempel"

Een van de belangrijkste bevindingen is dat er een magische drempel is.

• Zolang de wervelkracht (de energie van de stroming) niet sterk genoeg is om de "stijfheid" van de vloeistof te overwinnen, blijft het midden een star blokje.
• Pas op het moment dat de kracht van de stroming groter is dan de stijfheid van de vloeistof, breekt de kern en begint de echte turbulentie.

5. Waarom is dit belangrijk?

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

• Industrie: Bedrijven die verf, chocopasta, of afvalwater (slib) transporteren, moeten weten hoeveel pompen ze nodig hebben. Als ze niet weten dat de vloeistof eerst een "stijf blok" vormt, kunnen ze hun pompen te klein kiezen.
• Energie: Als je weet precies wanneer de stroming verandert, kun je energie besparen. Je wilt niet harder pompen dan nodig is, maar je wilt ook niet dat de vloeistof vastloopt.
• Medisch: Zelfs bloed (dat ook een beetje als deze vloeistof werkt) gedraagt zich zo. Begrijpen hoe het stroomt helpt bij het ontwerpen van betere kunstmatige hartkleppen of injectiespuiten.

Samenvatting

Kortom, deze paper vertelt het verhaal van een vloeistof die eerst als een stijf blokje door een pijp glijdt, en dan, als je er hard genoeg op duwt, plotseling "smelt" en verandert in een wilde, draaiende storm. De onderzoekers hebben de exacte momenten berekend waarop dit gebeurt, wat helpt om alles van verfvervoer tot bloedcirculatie beter te begrijpen en te optimaliseren.

Technische samenvatting

Hier is een gedetailleerde technische samenvatting van het artikel "Laminar-to-Turbulent Transition of Yield-Stress Fluids in Pipe and Channel Flows", geschreven in het Nederlands.

Technische Samenvatting: Overgang van Laminair naar Turbulent bij Vloeistoffen met Vloeigrens

Probleemstelling
Niet-Newtonse vloeistoffen met een vloeigrens (yield-stress), zoals die beschreven worden door het Herschel-Bulkley (HB) model, komen veel voor in industriële toepassingen (bijv. Carbopol-gels, slurry's, cosmetica). Een fundamentele uitdaging in de stromingsmechanica is het begrijpen van de overgang van laminair naar turbulent gedrag in deze vloeistoffen. In tegenstelling tot Newtonse vloeistoffen, waar de overgang goed gedocumenteerd is, is het gedrag van HB-vloeistoffen complexer vanwege de co-existentie van "geyieldde" (vloeibare) en "ongeyieldde" (plug-achtige) gebieden.
Bestaande literatuur mist een unificerend beeld van de volledige overgangsreeks (van plug-vorming tot volledige turbulentie) voor zowel pijpstroming als kanaalstroming. Er is behoefte aan direct numeriek onderzoek (DNS) om de mechanismen van plug-afbraak en de rol van vloeigrens in de turbulentie-ontwikkeling te ontrafelen.

Methodologie
De auteurs voerden Direct Numerical Simulations (DNS) uit voor Herschel-Bulkley vloeistoffen in twee configuraties:

1. Pijpstroming: Een cirkelvormige doorsnede met een lengte van $4\pi D$.
2. Kanaalstroming: Een rechthoekig kanaal met een stroomrichtingslengte van 0,2 m.

• Rheologisch Model: De vloeistof werd gemodelleerd als een Herschel-Bulkley materiaal met vloeigrens ($\tau_y$), consistentie-index ($k$) en power-law index ($n$). Om numerieke stijfheid bij zeer lage schuifsnelheden te voorkomen, werd een bi-viscositeits-regularisatie toegepast. Hierbij wordt de viscositeit afgesneden onder een bepaalde cut-off schuifsnelheid ($\dot{\gamma}_0$), waardoor de ongeyieldde zones een zeer hoge maar eindige viscositeit hebben.
• Numeriek Oplossingskader: De incompressibele momentumvergelijkingen werden opgelost met een schuif-afhankelijke viscositeit. Voor stabiliteit werd de niet-lineaire term in een skew-symmetrische vorm geïmplementeerd. De simulaties gebruikten gestructureerde roosters (O-grid voor pijpen, kartesiaans voor kanalen) met een resolutie van ongeveer 4,2 miljoen cellen voor het kanaal.
• Reynoldsgetallen: Omdat de viscositeit variabel is, werd een generaliseerd Reynoldsgetal ($Re_G$) gebruikt, gebaseerd op de wandviscositeit. De simulaties bestreken een breed bereik van $Re_G$ van 378 tot 5300.
• Validatie: De numerieke resultaten voor pijpstroming werden gevalideerd tegen experimentele LDV-data van Guzel et al. (2009) voor 0,1% Carbopol-oplossingen.

Belangrijkste Bijdragen

1. Eerste Unificerende DNS: Dit werk biedt de eerste DNS die de volledige overgang van laminair naar turbulent in HB-vloeistoffen oplost voor zowel pijp- als kanaalconfiguraties.
2. Mechanisme van Overgang: Het onderzoek identificeert dat de overgang wordt bepaald door de lokale concurrentie tussen turbulente momentumtransport en de vloeigrens. Turbulentie kan alleen ontstaan en in stand blijven wanneer de lokale Reynolds-spanningen de vloeigrens overschrijden, waardoor de plug-achtige kern vervormt en afbreekt.
3. Kwantificering van Regimes: De studie definieert duidelijke grenzen voor de stromingsregimes op basis van $Re_G$, wat overeenkomt met trends waargenomen bij Carbopol-vloeistoffen.

Resultaten

• Pijpstroming:

  • Laminaire regime ($Re_G < 1735$): De stroming wordt gedomineerd door een sterke, onbeweeglijke plug in de kern. Turbulentie is verwaarloosbaar.
  • **Transitieregio ($1735 < Re_G < 2920$):** Er treedt een scherpe stijging op in turbulentie-intensiteit en RMS-snelheidsfluctuaties ($u'_{rms}$). De plug begint te breken, en er ontstaan lokale instabiliteiten. Er wordt een asymmetrie waargenomen in de turbulentie-intensiteit, wat overeenkomt met eerdere experimentele bevindingen.
  • Turbulent regime ($Re_G > 2920$): De plug is volledig afgebroken. De stroming wordt gedomineerd door wandturbulentie, hoewel de kern minder sterk is dan bij Newtonse vloeistoffen.
  • De simulaties reproduceerden nauwkeurig de experimentele trends voor turbulentie-intensiteit en vorticale structuren (gevisualiseerd via het Q-criterium).

• Kanaalstroming:

  • De simulaties bevestigden dezelfde fysica: vorming van een centrale plug, groei van wand-nabije turbulentie (zichtbaar via Q-criterium structuren) en een progressieve toename van de turbulentie-intensiteit en kinetische energie met toenemende $Re_G$.
  • De snelheidsprofielen evolueerden van een afgevlakte "plug-vorm" bij lage $Re_G$ naar een voller profiel met steilere wandgradiënten bij hoge $Re_G$, kenmerkend voor turbulente wandstroming.

Betekenis en Conclusie
Dit onderzoek levert een cruciale bijdrage aan het fundamentele begrip van niet-Newtonse stromingen. Het bewijst dat de vloeigrens een actieve regulator is in de turbulentie-ontwikkeling: zolang de vloeigrens de lokale schuifspanning kan weerstaan, wordt turbulentie onderdrukt.
De resultaten bieden een robuust numeriek raamwerk dat experimentele observaties valideert en nieuwe inzichten geeft in de schaalafhankelijkheid van plug-afbraak. Dit is essentieel voor het ontwerp van pijpleidingen en industriële processen waarbij dergelijke vloeistoffen worden getransporteerd. Toekomstig werk omvat experimentele validatie met PIV (Particle Image Velocimetry) in een 3,6 m kanaal om de voorspellingen verder te verifiëren en geometrie-afhankelijke drempelwaarden te kwantificeren.