Fully Turbulent Wakes at Low Reynolds Numbers: the Case of the Thin Flat Plate

Dit artikel toont door middel van directe numerieke simulatie en experimentele vergelijking aan dat de zogflow achter een dun tweedimensionaal plat plaat volledig turbulent wordt bij een relatief laag Reynoldsgetal van 400, waarbij het statistische en spectrale kenmerken vertoont die niet te onderscheiden zijn van turbulente zogstromen bij hogere Reynoldsgetallen, een transitiepad dat significant verschilt van dat van canonieke circulaire of vierkante cilinders.

De kern

Stel je voor dat je een dun, plat stuk karton (zoals een speelkaart) vasthoudt in een sterke wind. Terwijl de wind tegen het kaartje blaast, ontstaat er een rommelig, kolkend spoor van lucht achter het object, een zogenaamde "wake" (vaan). Een lange tijd geloofden wetenschappers dat de wind heel hard moest blazen, of dat het object een specifieke vorm zoals een ronde buis of een vierkant blok moest hebben, voordat deze wake echt chaotisch en "turbulent" zou worden.

Dit artikel vertelt een ander verhaal. De onderzoekers ontdekten dat als je een dunne, platte plaat gebruikt, de lucht erachter veel bij een veel lagere windsnelheid chaotisch en turbulent wordt dan men verwacht. Sterker nog, dit gebeurt bij een snelheid waarbij de lucht bij andere vormen nog relatief kalm en geordend is.

Hier is een overzicht van hun bevindingen met eenvoudige analogieën:

1. De verrassing van de "Turbulentiedrempel"

Denk aan turbulentie als een drukke dansvloer.

• Het oude geloof (Ronde cilinders): Als je een ronde paal in de wind hebt, begint de lucht erachter als een kalme, ritmische dans (heen en weer wiegend). Het kost veel energie (hoge snelheid) voordat de dansers tegen elkaar opbotsen, wild ronddraaien en een chaotische bende veroorzaken (turbulentie). Deze overgang vindt geleidelijk plaats over een breed scala aan snelheden.
• De nieuwe ontdekking (Dunne platte plaat): De onderzoekers ontdekten dat voor een dunne platte plaat de "dansvloer" bijna onmiddellijk van kalm naar een wilde moshpit gaat. Zelfs bij een relatief lage windsnelheid (Reynoldsgetal 400) is de lucht achter de plaat al volledig chaotisch. Het gaat niet door de langzame, ritmische fasen die ronde palen wel doorlopen. Het springt direct naar het feestje.

2. Hoe ze het bewezen hebben

Om er zeker van te zijn dat ze het niet alleen maar verbeeldden, gedroegen de onderzoekers zich als detectives die een plaats delict vergelijken.

• De Simulatie (Het virtuele lab): Ze gebruikten supercomputers om de wind die tegen de plaat botst bij lage snelheden te simuleren (Re 150 en Re 400).
• De Test in de echte wereld (De windtunnel): Ze keken ook naar echte experimenten waarbij de wind veel harder blies (Re 12.500 en Re 19.700).
• De Match: Toen ze de simulatie bij lage snelheid (Re 400) vergeleken met de experimenten met hoge snelheid in de echte wereld, kwamen de patronen perfect overeen. De "vingerafdrukken" van de turbulentie — hoe de lucht bewoog, hoeveel energie het had en hoe het kolkte — waren identiek.
• De controlegroep: Toen ze naar de simulatie bij een nog lagere snelheid keken (Re 150), waren de patronen totaal anders. Het bevond zich nog steeds in de "kalme" fase, nog niet chaotisch. Dit bewees dat de overgang naar chaos ergens tussen 150 en 400 plaatsvindt.

3. De "Vingerafdruk" van turbulentie

Hoe weet je of een stroming echt turbulent is? Het artikel zoekt naar specifieke "tekens van leven" in de data:

• Het Energiespectrum (Het geluid van de wind): In een kalme stroming is de energie geconcentreerd in een paar specifieke tonen (zoals een fluit die een enkele toon speelt). In een turbulente stroming klinkt het als witte ruis of statische elektriciteit, met energie verspreid over een enorm bereik van frequenties. De onderzoekers vonden dat bij Re 400 het "geluid" van de wind achter de plaat al vol zat met deze chaotische ruis, net als in de experimenten met hoge snelheid.
• De "Intermittentie" (De incidentele schreeuw): In een echt turbulente stroming draait de lucht niet alleen maar rustig rond; het heeft plotselinge, intense uitbarstingen van snelheid en rotatie. De onderzoekers vonden deze "schreeuwen" in de data bij Re 400, maar deze waren afwezig bij Re 150.

4. Waarom is dit anders?

Het artikel suggereert dat de reden voor deze plotselinge sprong de vorm van het object is.

• Ronde/Vierkante objecten: Wanneer de wind tegen een rond of vierkant object botst, werkt het achterste deel van het object als een schild dat de luchtstroom erachter stabiliseert. Het kost veel energie om die stabiliteit te doorbreken.
• De dunne plaat: Omdat de plaat zo dun is, is er geen "achterkant" om de lucht af te schermen. De drukfluctuaties (het duwen en trekken van de lucht) zijn vanaf het begin direct verbonden met de kolkende wervelingen. Het is als het proberen te balanceren van een potlood op zijn punt versus het balanceren van een bowlingbal; het potlood (de dunne plaat) is inherent onstabiel en kantelt veel sneller over in chaos.

De Kern van het Verhaal

Dit artikel verandert ons begrip van hoe lucht rond platte objecten stroomt. Het bewijst dat dunne platte platen bij verrassend lage snelheden volledig turbulente wakes creëren, veel lager dan ronde of vierkante objecten. De overgang naar chaos is voor deze vormen geen langzaam, geleidelijk proces; het is een plotselinge, fundamentele verschuiving die heel vroeg in het snelheidsbereik plaatsvindt.

De onderzoekers hebben niet besproken hoe dit van toepassing is op de bouw van bruggen, het ontwerpen van auto's of medische hulpmiddelen. Ze hebben zich strikt gericht op het bewijzen dat dit fenomeen bestaat en hoe de fysica van de luchtstroom verschilt van wat we voorheen dachten.

Technische samenvatting

Technische Samenvatting: Volledig turbulente wakedatastromen bij lage Reynoldsgetallen: de casus van de dunne platte plaat

Probleemstelling
De numerieke simulatie van de nominaal tweedimensionale (2D) stroming langs een dunne platte plaat die loodrecht op een uniforme stroming staat, bleek historisch gezien uitdagend, waarbij CFD-studies tegenstrijdige resultaten opleverden met betrekking tot de wakedynamica bij lage Reynoldsgetallen ($Re < 1000$). Terwijl experimentele gegevens suggereren dat integrale parameters (zoals weerstandscoëfficiënten en Strouhal-getallen) zwak afhankelijk zijn van $Re$ in het 3D-regime, hebben simulaties zeer verschillende wakedynamica gerapporteerd. Een centrale vraag blijft of de wake van een dunne platte plaat direct na het ontstaan van 3D-instabiliteiten overgaat naar een volledig turbulente staat, vergelijkbaar met canonieke cilindrische waken, of dat deze een afzonderlijk, meerfasig transitiepad volgt. De auteurs onderzoeken of de wake langs een dunne platte plaat al volledig turbulent is bij een relatief laag Reynoldsgetal van $Re = 400$, een regime waarin canonieke circulaire en vierkante cilindrische waken typisch een sterke $Re$-afhankelijkheid en niet-turbulente fluctuatiekenmerken vertonen.

Methodologie
De studie maakt gebruik van een veelzijdige aanpak die Direct Numerical Simulations (DNS) combineert met experimentele validatie:

1. Numerieke Simulaties: Er werd gebruikgemaakt van twee onafhankelijke DNS-codes om robuustheid te waarborgen:

   • Spectral/hp Element Method (SEM): Geïmplementeerd in Nektar++, gebruikt voor de primaire hoog-nauwkeurige DNS bij $Re = 400$ (DNS400) en $Re = 150$ (DNS150).
   • Finite Volume Method (FVM): Geïmplementeerd in OpenFOAM, gebruikt voor mesh-convergentiestudies en kruisvalidatie (DNS400V).
   • Het computationele domein werd afgekapt en onderworpen aan no-slip condities op de plaat, vrije instroomcondities, en specifieke uitlaatcondities voorgesteld door Dong et al. Periodiciteit werd afgedwongen in de spanwise richting.
   • Mesh-convergentie werd rigoureus getest, waarbij werd gewaarborgd dat de resolutie van de Kolmogorov-lengteschaal ($\eta$) met $\Delta x/\eta$-waarden vergelijkbaar met of beter dan bestaande literatuur voor soortgelijke geometrieën werd bereikt.

2. Experimentele Data: Twee experimentele datasets werden verkregen met behulp van tijd-opgeloste stereoscopische Particle Image Velocimetry (PIV) in een windtunnel aan de University of Calgary:

   • EX12K: $Re = 12.500$
   • EX20K: $Re = 19.700$
   • Deze datasets vertegenwoordigen volledig turbulente regimes en dienen als benchmark voor de lage-$Re$ simulaties.

3. Analyse: De studie vergelijkt globale integrale grootheden, gemiddelde snelheidvelden, Reynolds-spanningen, turbulente kinetische energie (TKE) budgetten, energiespectra en statistische eigenschappen van snelheidgradiënten (specifiek niet-Gaussianiteit) tussen de verschillende $Re$-gevallen.

Belangrijkste Bijdragen en Resultaten

• Vroegtijdig Ontstaan van Turbulentie: De belangrijkste bevinding is dat de wake van een dunne platte plaat volledig turbulent wordt bij $Re = 400$. Dit ligt aanzienlijk lager dan de transitie waargenomen bij canonieke circulaire ($Re \approx 10^5$) of vierkante ($Re \approx 6000$) cilindrische waken.
• Kwantitatieve Overeenkomst met Hoog-$Re$ Experimenten: De DNS-resultaten bij $Re = 400$ (DNS400) vertonen een opmerkelijke overeenkomst met de experimentele data bij $Re = 12.500$ en $19.700$ (EX12K, EX20K) met betrekking tot:
  • Globale Grootheden: De gemiddelde weerstandscoëfficiënt ($C_D$) en het Strouhal-getal ($St$) zijn traag variërende functies van $Re$ voor $Re \ge 400$.
  • Stromingsstructuur: Gemiddelde stroomlijnen, snelheidsprofielen en Reynolds-spanningverdelingen ($u'^2, v'^2, w'^2$) zijn statistisch ononderscheidbaar tussen $Re=400$ en de hoog-$Re$ experimenten binnen de meetonzekerheid.
  • Energiebudgetten: De ruimtelijke verdeling van TKE-productie en advectietermen in de $k$-transportvergelijking is consistent over de $Re=400$ en de hoog-$Re$ gevallen.
• Turbulente Signaturen bij Laag $Re$: Het $Re=400$ geval vertoont definitieve kenmerken van turbulentie die afwezig zijn in het $Re=150$ geval:
  • Energiespectra: De fluctuerende snelheidcomponenten vertonen breedbandige spectra met een inert bereik-schaling die ongeveer voldoet aan $f^{-5/3}$ (Kolmogorovs wet), wat duidt op een energiespectrum. In contrast hiermee vertonen de $Re=150$ spectra energieoverdrachten tussen slechts enkele modi.
  • Niet-Gaussianiteit: Waarschijnlijkheidsdichtheidsfuncties (PDF's) van fluctuerende rek ($\Sigma$) en rotatiesnelheden ($\Omega$) vertonen zware staarten, wat wijst op intermittentie en zeldzame, hoog-intensieve gebeurtenissen die karakteristiek zijn voor turbulentie. Het $Re=150$ geval vertoont Gaussianische verdelingen.
• Afwijkend Transitiepad: De transitie bij $Re=150$ omvat 3D-instabiliteiten maar mist de coherente, volledig turbulente structuur die gezien wordt bij $Re=400$. De $Re=150$ wake kenmerkt zich door zwakkere terugstroming, een langzamere wakereductie en een niet-monotone evolutie van de wakediameter, wat fundamenteel verschilt van het gedrag van canonieke cilinders waar 3D-instabiliteiten (Modi A, B, C) over een breed $Re$-bereik evolueren voordat volledige turbulentie optreedt.

Betekenis en Claims
Het artikel betoogt dat het pad naar de transitie naar turbulentie voor een dunne platte plaat fundamenteel anders is dan dat van canonieke cilinders. Terwijl cilindrische waken een langdurige sequentie van instabiliteiten (Modus A, B, C) ondergaan over een breed Reynoldsgetal voordat ze een volledig turbulente staat bereiken, transiteert de wake van een dunne platte plaat bijna onmiddellijk na het ontstaan van 3D-instabiliteiten (tussen $Re=150$ en $Re=400$).

De auteurs suggereren dat dit verschil toegeschreven kan worden aan het gebrek aan een stroomopwaartse achterzijde (afterbody) in de geometrie. Voor cilinders stabiliseert de achterzijde de stroming en ontkoppelt deze drukfluctuaties van het scheidingspunt. Voor dunne platte platen zijn drukfluctuaties echter direct gekoppeld aan de vortexvormingsregio, wat potentieel een eerdere en snellere transitie naar turbulentie drijft. Deze bevindingen helpen de discrepanties in de eerdere literatuur over de gevoeligheid van dunne plaat-wakes voor aspectratio's en Reynoldsgetallen te verzoenen, en stellen vast dat het "volledig turbulente" regime voor deze geometrie bij een veel lager $Re$ begint dan voorheen aangenomen voor canonieke bluff bodies.