Leveraging unstructured grids for direct numerical simulations of wall turbulence

Dit artikel introduceert de {\eta}-grid, een ongestructureerd raamwerk voor het genereren van roosters voor directe numerieke simulaties van wandturbulentie dat roostergroottes schaalt met de lokale Kolmogorov-schaal, waarbij nauwkeurigheid wordt bereikt die vergelijkbaar is met conventionele Cartesische roosters, terwijl de rekenkosten aanzienlijk worden verlaagd, met name bij hoge Reynoldsgetallen en over complexe riblet-geometrieën.

De kern

Stel je voor dat je probeert te simuleren hoe lucht over een auto stroomt of hoe water langs de romp van een schip beweegt. Om dit nauwkeurig op een computer te doen, gebruiken wetenschappers een techniek die Direct Numerical Simulation (DNS) heet. Denk aan DNS als het maken van een gigantische, 3D digitale microscoop die het fluïdum (lucht of water) opdeelt in miljoenen kleine, onzichtbare kubussen (een rooster). De computer berekent vervolgens hoe elke enkele kubus beweegt en interageert met zijn buren.

Het probleem is dat stroming dicht bij een oppervlak (zoals de zijkant van een schip) ongelooflijk chaotisch en gedetailleerd is. Om een duidelijk beeld te krijgen, heb je een enorm aantal van deze kleine kubussen direct naast het oppervlak nodig. Echter, naarmate je verder van het oppervlak verwijderd bent, kalmeert het chaos en heb je niet meer zulke kleine kubussen nodig.

De Oude Weg: De "Stijve Bakstenen Muur"

Traditioneel gebruikten wetenschappers een Cartesiaans rooster. Stel je voor dat je een muur bouwt van identieke, stijve bakstenen.

• Het Probleem: Om de kleine details dicht bij het oppervlak te zien, moet je de bakstenen onderaan heel klein maken. Maar omdat deze bakstenen stijf zijn en in een rechte lijn verbonden, word je gedwongen om diezelfde kleine bakstenen te gebruiken tot aan de top van je muur, zelfs waar de details niet belangrijk zijn.
• Het Resultaat: Je eindigt met een muur van miljarden kleine bakstenen, waarvan de meeste onnodig zijn. Dit maakt de computersimulatie ongelooflijk traag en duur, alsof je probeert elk zandkorreltje op een strand te tellen om alleen maar het getij te meten.

De Nieuwe Oplossing: Het "Slimme, Rekbaar Net"

Dit artikel introduceert een nieuwe methode die het $\eta$-rooster (eta-rooster) wordt genoemd. In plaats van stijve bakstenen, stel je je een slim, rekbaar visnet voor.

• Hoe het werkt: De auteurs hebben een systeem ontworpen waarbij de grootte van de gaten in het net automatisch verandert op basis van hoeveel detail er nodig is.
  • Dicht bij het oppervlak (De "Binnenste Laag"): Het net heeft zeer kleine, strakke gaten om de kleine, chaotische wervelingen van het fluïdum te vangen.
  • Verder weg (De "Buitenste Laag"): Naarmate het fluïdum kalmeert, rekt het net automatisch uit, waardoor de gaten veel groter worden.
• Het Geheime Ingrediënt: De grootte van deze gaten is gebaseerd op iets dat de Kolmogorov-schaal (aangeduid als $\eta$) wordt genoemd. Denk hierbij aan de "kleinste mogelijke werveling" die op een bepaald moment in het fluïdum kan bestaan. Het nieuwe rooster zegt simpelweg: "Maak de gatgrootte net groot genoeg om de kleinste werveling op deze specifieke hoogte te vangen, en niet groter."

Waarom Dit Een Grote Zaken Is

De auteurs hebben dit "slimme net" getest op twee verschillende soorten computercodes (de ene lijkt op een spectrale elementmethode, de andere op een finite volume-methode) en vergeleken het met de oude "stijve baksteen"-methode.

1. Nauwkeurigheid: De resultaten waren bijna identiek. Het "slimme net" ving de fysica net zo goed op als de "stijve bakstenen", met minder dan 1% verschil in belangrijke metingen zoals wrijving en snelheid.
2. Enorme Besparingen: Hier gebeurt de magie.
   • Voor gladde oppervlakken (zoals een vlakke muur) verminderde het nieuwe rooster het aantal benodigde "bakstenen" (roosterpunten) met ongeveer 90% bij hoge snelheden.
   • Voor ruwe oppervlakken (zoals een muur met kleine groeven die "riblets" worden genoemd en zijn ontworpen om weerstand te verminderen) waren de besparingen nog dramatischer: tot 97% minder roosterpunten.

De Analogie van de "Gegroefde Muur"

Om het riblet-gedeelte te begrijpen, stel je je een muur voor die bedekt is met kleine, parallelle groeven (zoals de textuur van een golfbal of de huid van een haai).

• De Oude Weg: Om dit te simuleren, moest de stijve baksteenmethode de bakstenen overal klein houden omdat de groeven het rooster dwongen om tot ver omhoog fijn te zijn. Het was alsof je probeerde elke enkele draad in een trui te tellen, zelfs de delen ver weg van de stof.
• De Nieuwe Weg: Het "slimme net" weet dat zodra je een paar centimeter boven die kleine groeven bent, de stroming weer glad wordt. Het rekt de gaten direct boven de groeven uit en negeert de kleine details die niet langer belangrijk zijn.

De Conclusie

De auteurs hebben een raamwerk gecreëerd dat fungeert als een slimme zoomlens voor stromingssimulaties. Het richt zijn rekenkracht precies daar waar het nodig is (dicht bij de muur) en ontspant waar het dat niet is.

• Voor Gladde Muren: Het schaalt de inspanning veel langzamer naarmate de simulatie groter wordt.
• Voor Ruwe Muren: Het schaalt nog beter, waardoor simulaties van complexe, weerstandsverminderende oppervlakken haalbaar worden op computers die ze eerder niet aankonden.

Kortom, ze hebben een manier gevonden om hetzelfde hoogwaardige werk te verrichten met een fractie van de rekenkracht, waardoor een taak die een supercomputer misschien een maand zou kosten, in een paar dagen kan worden gedaan.

Technische samenvatting

Technische Samenvatting: Benutting van Ongestructureerde Roosters voor Directe Numerieke Simulaties van Wandturbulentie

Probleemstelling
Directe Numerieke Simulatie (DNS) van wandturbulentie is rekenkundig duur, voornamelijk vanwege de vereiste om de Kolmogorov-lengteschaal ($\eta$) over het gehele domein op te lossen. Conventionele Cartesische roosteroplossers, die wijdverspreid worden gebruikt voor gladde en onregelmatige oppervlakken (bijvoorbeeld via Immersed Boundary Methods), hanteren doorgaans een vaste of langzaam variërende streamwise ($\Delta x^+$) en spanwise ($\Delta z^+$) roosterresolutie over het hele domein. Deze aanpak leidt tot aanzienlijke overresolutie in de buitenste gebieden waar de Kolmogorov-schaal toeneemt, wat een rekenkundige knelpunt creëert. Dit is bijzonder acuut voor simulaties die complexe oppervlakken zoals ribbels (riblets) omvatten, waarbij het rooster extreem fijn moet zijn nabij het oppervlak om micro-groeven op te lossen, maar deze fijne resolutie in Cartesische kaders onnodig wordt gehandhaafd over de volledige dikte van de grenslaag.

Methodologie
De auteurs stellen een nieuw raamwerk voor het genereren van ongestructureerde roosters voor, genaamd het $\eta$-rooster. Het kernprincipe is om de wandnormale ($\Delta y^+$) en spanwise ($\Delta z^+$) roostergroottes evenredig in te stellen met de lokale Kolmogorov-schaal ($\eta^+$), waardoor roostervergroving mogelijk is in de buitenste stromingsgebieden waar de turbulentieschalen groter zijn.

1. Roostervormulering:

   • Binnenste laag: Een dunne laag ($\sim 50$ viskeuze eenheden) aangrenzend aan de wand maakt gebruik van op viskeuze schaal gebaseerde roostergroottes die vergelijkbaar zijn met conventionele DNS ($0.3 \lesssim \Delta y^+ \lesssim 4$, $\Delta z^+ \simeq 5$). Voor onregelmatige oppervlakken (bijvoorbeeld ribbels) lost het rooster de kleinste oppervlaktegolflengte ($\ell^+$) op met $\ell^+/30 \lesssim \Delta y^+, \Delta z^+ \lesssim 4$.
   • Buitenste gebied: Boven de binnenste laag nemen de roostergroottes toe evenredig met de lokale Kolmogorov-schaal ($\Delta y^+ \simeq \Delta z^+ \simeq 2\eta^+$).
   • Schaalwetten: De auteurs verfijnen semi-empirische fits voor $\eta^+$ in zowel het logaritmische als het wake-gebied voor turbulente kanaalstroming en turbulente grenslaagstromingen met nul-drukgradiënt (ZPG TBL). Zij stellen stuksgewijze machtswet-fits voor $\eta^+$ voor die overgaan bij $y^+ \approx \delta^+/2$.
   • Ribbel-aanpassing: Voor ribbelgeometrieën introduceert het raamwerk een "ribbelsublaag" ($\delta_r^+$) gevuld met vierkante cellen om secundaire stromingen en Kelvin-Helmholtz-rollen op te lossen, die soepel overgaat in de op $\eta$ gebaseerde schaling hierboven.

2. Implementatie:

   • Het raamwerk maakt gebruik van het open-source pakket Gmsh om multi-blok ongestructureerde roosters te genereren.
   • De roosters worden getest met twee verschillende oplossers: SOD2D (Spectrale Elementen Methode) en OpenFOAM (Finite Volume Methode).
   • Simulaties bestrijken turbulente kanaalstroming en ZPG TBL's over gladde wanden en diverse ribbelgeometrieën (driehoekige groeven) tot een wrijvings-Reynoldsgetal $\delta_0^+ = 1000$.

Belangrijkste Bijdragen

• Ontwikkeling van het Raamwerk: Een robuuste, multi-blok strategie voor het genereren van ongestructureerde roosters die de roosterafstand dynamisch aanpast aan de lokale Kolmogorov-schaal, toepasbaar op zowel gladde als complexe ruwe oppervlakken.
• Gefinancierde Schaalwetten: Bijgewerkte semi-empirische fits voor de Kolmogorov-schaal in het wake-gebied van TBL's en kanaalstromingen, waardoor de nauwkeurigheid van de roostergroottebepaling in de buitenste laag verbetert.
• Oplosser-onafhankelijkheid: Aantonen dat het $\eta$-rooster-raamwerk consistente, hoge-trouw resultaten oplevert over verschillende numerieke discretisatieschema's (SEM en FVM).
• Specifieke Formulering voor Ribbels: Een op maat gemaakt roosterstrategie voor ribbels die de fysica van de sublaag (secundaire stromingen, KH-rollen) oplost, terwijl snelle vergroving in de buitenste stroming mogelijk is, waarmee de beperkingen van Cartesische/IBM-aanpakken worden overwonnen.

Resultaten

• Nauwkeurigheid: Het $\eta$-rooster produceert resultaten met minder dan 1% verschil in vergelijking met conventionele Cartesische roosters en referentie-DNS/experimentele data. Metrieken omvatten wrijvingscoëfficiënten ($C_f$), gemiddelde snelheidsprofielen ($U^+$), turbulente spanningen en spectrogrammen.
• Roosterefficiëntie:
  • Gladde Wand: Het aantal roosterpunten ($N_\eta$) schaalt als $\propto \delta_0^{+2.5}$, vergeleken met $\propto \delta_0^{+3.0}$ voor Cartesische roosters met hyperbolische tangens-verlenging. Bij $\delta_0^+ = 6000$ vereist het $\eta$-rooster slechts $\sim 10\%$ van het aantal roosterpunten van een standaard Cartesisch rooster ($N_\eta/N_{Tanh} \simeq 0.1$).
  • Ribbels: Het winst in efficiëntie is nog uitgesprokener. Voor wrijvingsreducerende driehoekige ribbels schaalt $N_\eta$ als $\propto \delta_0^{+2.0}$, terwijl Cartesische roosters schalen als $\propto \delta_0^{+3.0}$. Bij $\delta_0^+ = 6000$ vereist het $\eta$-rooster slechts $\sim 3\%$ van het aantal roosterpunten van een Cartesisch rooster ($N_\eta/N_{Tanh} \simeq 0.03$).
• Ribbelfysica: Het raamwerk slaagt erin complexe stromingsmechanismen over ribbels vast te leggen, inclusief de virtuele oorsprong, secundaire stromingen en Kelvin-Helmholtz-rollen, met roosterconvergentie bereikt met gematigde roostereparameters ($C_y=2.0, C_z=2.5, y_{in}^+=50$).
• TBL-validatie: DNS van ZPG TBL's over ribbels toonde goede overeenstemming met experimentele data (Baron & Quadrio; Choi & Orchard) wat betreft percentages wrijvingsreductie en snelheidsprofielen, ondanks inherente discrepanties tussen experiment en DNS in wake-gebieden.

Betekenis
Het artikel beweert dat het $\eta$-rooster-raamwerk een "enorme roosterbewaring" biedt, waardoor DNS van wandturbulentie bij hoge Reynoldsgetallen over complexe oppervlakken rekenkundig haalbaar wordt. Door de roosterresolutie af te stemmen op de lokale fysieke schalen (Kolmogorov-schaal) in plaats van vaste resoluties te handhaven, verlaagt de methode de rekenkosten (roosterpunten en tijdstapbeperkingen) drastisch zonder nauwkeurigheid op te offeren. De auteurs benadrukken dat deze aanpak het betaalbaar maakt om DNS van TBL's over complexe oppervlakken uit te voeren bij Reynoldsgetallen die overeenkomen met experimentele omstandigheden, een taak die eerder werd gehinderd door de verbodsdreigende kosten van op Cartesische methoden gebaseerde benaderingen. Het werk valideert de doeltreffendheid van ongestructureerde roosters voor hoge-trouw turbulentiesimulaties, en gaat voorbij de beperkingen van gestructureerde Cartesische roosters.