Helicity subgrid-scale models and their numerical validation

Dit onderzoek toont aan dat het integreren van SGS-heliciteitseffecten in het Smagorinsky-model, naast de wervelviscositeit, de nauwkeurigheid van grootschalige turbulentiesimulaties verbetert, zoals geverifieerd via directe numerieke simulaties.

De kern

De Heliciteit-SGS-modellen en hun validatie: Een uitleg voor iedereen

Stel je voor dat je probeert het weer te voorspellen, de stroming van een rivier te begrijpen of hoe een tornado ontstaat. De natuur is hierbij extreem chaotisch. Water en lucht bewegen in een wirwar van draaikolken, van enorme stormsystemen tot microscopisch kleine werveltjes. Dit noemen we turbulentie.

Het probleem voor wetenschappers is dat computers niet krachtig genoeg zijn om elk van die werveltjes, van de grootte van een berg tot die van een zandkorrel, tegelijkertijd te berekenen. Het zou te veel tijd en energie kosten.

Om dit op te lossen, gebruiken ze een slimme truc: Large-Eddy Simulations (LES).

• De Grote Dingen: De computer rekent de grote stromingen (de "eddies" of wervels) die je kunt zien, wel uit.
• De Kleine Dingen: De heel kleine werveltjes die de computer niet kan zien, worden niet berekend, maar geschat met een wiskundig model. Dit noemen we het Subgrid-Scale (SGS) model.

Het oude probleem: De "Smagorinsky"

Voor decennia was het meest gebruikte model om die kleine werveltjes te schatten het Smagorinsky-model. Je kunt dit vergelijken met een simpele schatting: "Als de stroming hard beweegt, moet er veel wrijving zijn, dus we gaan veel energie verliezen."

Maar dit model heeft een groot nadeel: het is te pessimistisch. Het denkt dat er te veel wrijving is. Het "dempt" de stroming te snel, alsof je een auto in de modder rijdt terwijl hij eigenlijk over glad ijs glijdt. In de echte wereld zie je vaak sterke, draaiende structuren (zoals in een tornado of langs de wanden van een pijp) die dit simpele model niet goed kan nabootsen. Wetenschappers moesten daarom de instellingen van het model handmatig aanpassen voor elke situatie, wat niet ideaal is.

De nieuwe oplossing: Heliciteit (De "Schroef")

In deze paper stellen de auteurs een nieuwe, betere manier voor. Ze kijken niet alleen naar hoe snel de stroming beweegt, maar ook naar de vorm van die beweging.

Ze introduceren een concept dat heliciteit heet.

• Analogie: Stel je voor dat je een touw vasthoudt. Als je het recht trekt, is er alleen spanning (energie). Maar als je het in een schroefdraad of een dop draait, heb je een draaiende beweging. Die draaiing is heliciteit.
• In de natuur zie je dit vaak: tornado's, draaikolken in de oceaan, of zelfs de wind die om een gebouw draait. Deze dingen hebben een sterke "schroefvorm".

Het oude Smagorinsky-model keek alleen naar de snelheid (hoe hard de wind waait). Het nieuwe Heliciteit-SGS-model kijkt ook naar de draaiing (hoe sterk de wind schroeft).

Wat hebben ze gedaan? (De "Test")

De auteurs (Yokoi, Mininni, en anderen) hebben dit nieuwe model getest. Ze deden dit niet in een echte rivier, maar in een superkrachtige computer-simulatie die ze Direct Numerical Simulation (DNS) noemen.

• Dit is als een "god-simulatie": hierin wordt alles berekend, van de grootste storm tot het kleinste werveltje. Het is de "waarheid".
• Vervolgens hebben ze hun nieuwe model gebruikt om te proberen deze "waarheid" na te bootsen, maar dan alsof ze alleen de grote dingen zagen (net als in een echte simulatie).

Ze hebben twee situaties getest:

1. Zonder draaiing: Een willekeurige, chaotige stroming.
2. Met draaiing: Een stroming die wordt beïnvloed door rotatie (zoals op de Aarde, waar de Corioliskracht werkt).

De Resultaten: Waarom is dit beter?

De uitkomsten waren verrassend goed:

1. De "Schroef" maakt het verschil: In situaties waar de stroming draait (zoals bij rotatie of in een pijp), zag het oude Smagorinsky-model de stroming verkeerd. Het dacht dat de energie verdween, terwijl hij eigenlijk bleef bestaan in die draaiende structuren.
2. Het nieuwe model pakt het op: Door de "schroefvorm" (heliciteit) mee te nemen, kon het nieuwe model precies voorspellen hoe die draaiende wervels zich gedragen. Het model werd veel nauwkeuriger, vooral bij de complexe, schuine krachten in de stroming.
3. Een universele oplossing? Het mooie is dat dit nieuwe model misschien wel werkt voor alle situaties zonder dat je de instellingen hoeft aan te passen. Het oude model had voor elke situatie een andere "knop" nodig. Dit nieuwe model zou met één instelling voor alles kunnen werken, omdat het de fysica van de draaiing begrijpt.

De Conclusie in het Kort

Deze paper zegt eigenlijk: "We hebben al die tijd te veel gekeken naar hoe hard de wind waait, en te weinig naar hoe hij draait."

Door de draaiing (heliciteit) toe te voegen aan de wiskundige modellen die we gebruiken om stromingen te simuleren, krijgen we veel betere voorspellingen. Het is alsof je van een simpele platte kaart overstapt naar een 3D-kaart met diepte en vorm. Dit helpt ons niet alleen om betere weermodellen te maken, maar ook om beter te begrijpen hoe tornado's ontstaan, hoe sterren energie transporteren, en hoe we efficiëntere vliegtuigen of schepen kunnen ontwerpen.

Kortom: Meer draaiing in de wiskunde betekent minder fouten in de simulatie.

Technische samenvatting

Titel: Heliciteit SGS-modellen en hun validatie

Auteurs: N. Yokoi et al.
Onderwerp: Turbulente stromingen, Large-Eddy Simulatie (LES), Subgrid-scale (SGS) modellering.

---

1. Het Probleem

Directe Numerieke Simulaties (DNS) van turbulente stromingen bij hoge Reynolds-getallen zijn vaak onmogelijk vanwege de enorme rekenkracht die nodig is om alle schalen op te lossen. Large-Eddy Simulatie (LES) is een krachtig alternatief waarbij grote schalen expliciet worden opgelost en kleine schalen (Subgrid-scale of SGS) worden gemodelleerd.

De meest gebruikte SGS-model is het Smagorinsky-model. Dit model heeft echter fundamentele beperkingen:

• Over-dissipatie: Het model dissipeert vaak te veel energie, vooral in gebieden met grote rek (strain).
• Gebrek aan universaliteit: De Smagorinsky-constante ($C_S$) moet handmatig worden aangepast afhankelijk van de stroming (bijv. $0.18$ voor isotrope turbulentie, maar slechts $0.10$ voor wandstromingen).
• Beperkte fysica: Het Smagorinsky-model baseert zich uitsluitend op de turbulentie-intensiteit (energie) en de rek-tensor. Het negeert de structurele of geometrische informatie van de turbulentie, zoals de correlatie tussen snelheid en vorticiteit (wervelbeweging).

In wandstromingen en andere complexe stromingen ontstaan sterke, langsgestroomde wervelstructuren (streaks). De sterkte van deze structuren varieert per stromingstype, wat de noodzaak voor het aanpassen van $C_S$ verklaart. De auteurs stellen dat heliciteit (de lokale correlatie tussen snelheid en vorticiteit, $u \cdot \omega$) de sleutel is om deze structurele informatie te vangen en de over-dissipatie van het Smagorinsky-model te verminderen.

2. Methodologie

De auteurs ontwikkelen en valideren nieuwe heliciteit SGS-modellen (H-SGS) door middel van vergelijkingen met Directe Numerieke Simulaties (DNS).

• Theoretische Basis:

  • Gebaseerd op de Twee-Schaal Directe Interactie Benadering (TSDIA) voor inhomogene turbulentie.
  • De Reynolds-spanning (en analoog de SGS-spanning) wordt niet alleen uitgedrukt in termen van de rek-tensor (symmetrisch deel), maar ook in termen van de absolute vorticiteit (antisymmetrisch deel) gekoppeld aan de gradient van de turbulentie-heliciteit.
  • De nieuwe SGS-spanningsterm ($\tau_{SGS}$) bevat een extra term: $\eta_S [\nabla H_S \cdot \omega^* + \dots]$, waarbij $H_S$ de SGS-heliciteit is en $\omega^*$ de absolute vorticiteit (relatieve vorticiteit + systeemrotatie).

• Modelvarianten:
  De auteurs presenteren drie niveaus van modellen, afhankelijk van hoe SGS-variabelen worden berekend:

  1. Twee-vergelijking model: Oplossing van transportvergelijkingen voor SGS-energie ($K_S$) en SGS-heliciteit ($H_S$).
  2. Eén-vergelijking model: Oplossing van alleen de $H_S$-vergelijking, met $K_S$ afgeleid uit lokale evenwichtshypothese.
  3. Nul-vergelijking model (Zero-equation): Geen transportvergelijkingen. Alle SGS-variabelen ($K_S, H_S, \nu_S, \eta_S$) worden algebraïsch uitgedrukt in termen van de filterbreedte ($\Delta$) en de opgeloste (GS) veldgrootheden. Dit is het meest vergelijkbaar met het Smagorinsky-model.

• Numerieke Validatie:

  • Gebruik van de GHOST-code (een pseudo-spectrale code) voor DNS in een driedimensionaal periodiek domein ($1024^3$ roosterpunten).
  • Configuratie: Twee scenario's: zonder rotatie (vergelijkbaar met isotrope turbulentie met grote schaal inhomogeniteit) en met rotatie ($\omega_F$).
  • Afdwinging (Forcing): Mechanische forcing wordt toegepast om een grote-schaal inhomogene heliciteit te creëren (een hyperbolische tangens profiel in de x-richting).
  • Filtering: Een Gauss-filter wordt gebruikt om de DNS-data te scheiden in opgeloste (GS) en niet-opgeloste (SGS) componenten.
  • Validatiemethode: "Scatter plots" worden gebruikt om de correlatie te meten tussen de SGS-spanningen berekend uit de DNS-data en de voorspellingen van de modellen (zowel Smagorinsky als H-SGS).

3. Belangrijkste Bijdragen

1. Incorporatie van Heliciteit in SGS: De auteurs introduceren een nieuwe term in de SGS-spanning die de gradient van de SGS-heliciteit koppelt aan de absolute vorticiteit. Dit introduceert structurele informatie in het model.
2. Algebraïsche Evaluatie van SGS-Heliciteit: Voor het nul-vergelijking model wordt een methode ontwikkeld om de SGS-heliciteit ($H_S$) direct te schatten op basis van de forcing en de lokale rek, zonder extra transportvergelijkingen te hoeven oplossen.
3. Systematische Validatie: Een uitgebreide vergelijking tussen DNS-data en modelvoorspellingen voor zowel diagonale als niet-diagonale componenten van de SGS-spanningstensor, onder verschillende filterbreedten en rotatiecondities.

4. Resultaten

De numerieke resultaten tonen duidelijke verbeteringen van het H-SGS-model ten opzichte van het standaard Smagorinsky-model:

• Diagonale Componenten: Voor de diagonale componenten van de SGS-spanning (gerelateerd aan druk en isotrope druk) is er weinig verschil tussen de modellen. Dit is te verwachten omdat de heliciteitsterm in de huidige opstelling voornamelijk invloed heeft op de niet-diagonale (shear) componenten.
• Niet-diagonale Componenten (Off-diagonal):
  • Het Smagorinsky-model toont geen significante correlatie met de DNS-data voor de niet-diagonale spanningen ($\tau_{13}$), ongeacht de filterbreedte of rotatie.
  • Het H-SGS-model toont een sterke verbetering in correlatie. Vooral bij kleinere filterbreedten ($k_C=14$) en in aanwezigheid van sterke heliciteit-gradienten, bereikt de correlatie een hoge waarde ($r \approx 0.90$).
  • De verbetering is het grootst in gebieden met grote SGS-heliciteit-gradienten, wat bevestigt dat de nieuwe term essentieel is voor het modelleren van de momentumtransport in deze gebieden.
• Invloed van Rotatie: De verbetering door het H-SGS-model is ook aanwezig in rotatie-stromingen, wat aangeeft dat het model robuust is voor zowel niet-rotatie als rotatie-gebaseerde turbulentie.
• SGS-variabelen: De tijdschalen, energie en dissipatie van SGS kunnen goed worden gemodelleerd met algebraïsche relaties gebaseerd op de filterbreedte en rek, zelfs met rotatie.

5. Betekenis en Conclusie

De studie concludeert dat het inhomogene heliciteitseffect een cruciale rol speelt in het momentumtransport van turbulente stromingen, vooral in aanwezigheid van absolute vorticiteit (rotatie of wandstromingen).

• Oplossing voor Over-dissipatie: Door de heliciteitsterm toe te voegen, kan het model de over-dissipatie van het Smagorinsky-model verminderen. De heliciteitsterm werkt tegen de eddy-viscositeitsterm en helpt de menging van impuls te onderdrukken waar nodig.
• Universaliteit: Het model biedt de belofte om de Smagorinsky-constante ($C_S$) universeel te maken (bijv. vasthouden aan $0.18$ voor alle stromingen), in plaats van deze handmatig aan te passen voor wandstromingen. De heliciteitsterm compenseert dan het verschil in dissipatiebehoeften.
• Toekomstperspectief: Hoewel de huidige validatie in periodieke dozen succesvol was, wordt aanbevolen om het model toe te passen op complexere stromingen zoals kanaalstromingen (channel flows) en menglagen, waar sterke langsgestroomde wervelstructuren voorkomen. Als het model daar werkt, zou het een doorbraak kunnen betekenen voor het modelleren van wandturbulentie zonder ad-hoc aanpassingen.

Kortom, het artikel levert een sterke fysieke onderbouwing en numerieke validatie voor het integreren van heliciteit in subgrid-modellen, wat leidt tot nauwkeurigere voorspellingen van turbulente stromingen.