Towards a full-scale version of Yakhot's model of strong turbulence

Dit artikel presenteert een uitgebreide, parameterloze versie van Yakhot's model voor sterke turbulentie die structurele functies van de tweede en derde orde nauwkeurig beschrijft van de dissipatieve schaal tot de systeemgrootte, gebaseerd op een empirische relatie voor even-orde structuurfuncties en Kolmogorovs vier-vijfde wet.

De kern

Stel je voor dat je naar een stromende rivier kijkt. Op grote schaal zie je misschien grote, rustige draaikolven. Maar als je heel dichtbij kijkt, zie je dat het water niet rustig stroomt, maar een chaotische dans van kleine wervelingen uitvoert. Dit is turbulentie. Het is een van de lastigste problemen in de natuurkunde: hoe beschrijf je die chaos precies?

Deze paper, geschreven door Christoph Renner, is als het ware een nieuwe, betere recept voor het koken van die turbulentie. Hij bouwt voort op een bestaand recept (het model van Yakhot) en voegt daar de ontbrekende ingrediënten aan toe om het te laten werken, van de grootste stromingen tot de allerkleinste druppels.

Hier is een uitleg in gewone taal, met een paar creatieve vergelijkingen:

1. Het Probleem: De "Gaten" in het Recept

Vroeger hadden wetenschappers een model (Yakhot's model) dat goed werkte voor de "middelmatige" schaal van de turbulentie (de inertiale range).

• Het grote gat: Het model wist niet hoe het zich moest gedragen op de grootste schaal (waar de stroming begint) en op de kleinste schaal (waar wrijving het water tot stilstand brengt).
• De oplossing van Renner: Hij heeft twee nieuwe stukjes toegevoegd om het recept compleet te maken:
  1. Een stukje voor de grote schaal (al eerder bedacht door hemzelf).
  2. Een nieuw stukje voor de kleine schaal (de focus van dit artikel).

2. De Nieuwe Ontdekking: Een Geheime Link

Renner keek naar meetgegevens van heliumgas (een heel koud, vloeibaar gas dat als een ideale vloeistof werkt). Hij zocht naar een patroon tussen verschillende soorten "stootjes" in de stroming.

• De Analogie: Stel je voor dat je de snelheid van het water meet op twee punten. Het verschil tussen die snelheden noemen we een "stootje".
  • Soms meet je het gemiddelde van die stootjes (evene machten).
  • Soms meet je het gemiddelde van de richting (oneven machten).
• De Doorbraak: Renner ontdekte een verrassende regel: De manier waarop de "stootjes" veranderen op de allerkleinste schaal, hangt direct samen met de "stootjes" van de volgende grootteschaal.
  • Het is alsof je ziet dat de trilling van een kleine golfje (de dissipatie) precies voorspelt hoe het volgende, iets grotere golfje zich gedraagt.
  • Hij noemt dit een "empirische relatie". In het Nederlands: "We hebben gekeken, gemeten en een regel gevonden die werkt."

3. De "Overgangszone": De Drempel

Een belangrijk nieuw element in zijn model is een nieuwe lengte.

• De Metafoor: Denk aan een auto die van een snelweg (de grote, chaotische stroming) afrijdt en de oprijlaan oprijdt naar een parkeerplaats (waar de wrijving de auto tot stilstand brengt).
• Er is een punt waarop de auto van "snelweg-modus" naar "parkeer-modus" schakelt.
• Renner introduceert een nieuwe maatstaf (noem het $\rho$) die precies dat punt aangeeft.
  • Boven deze maatstaf: De stroming volgt de oude, bekende wetten van chaos.
  • Onder deze maatstaf: De wrijving (viscositeit) begint te winnen en de stroming wordt rustiger.
• Het mooie is: Deze maatstaf hangt af van de Reynoldsgetal (een maat voor hoe "snel" en "chaotisch" de stroming is). Hoe sneller de stroming, hoe kleiner dit puntje is.

4. Het Resultaat: Een Volledig Recept

Door deze nieuwe regels te combineren met de oude wetten (zoals de beroemde "vier-vijfde wet" van Kolmogorov), heeft Renner formules opgesteld die:

1. Geen losse parameters hebben: Je hoeft niet te gokken met instellingen; het model wordt volledig bepaald door de eigenschappen van de stroming zelf.
2. Overal werken: Ze beschrijven het gedrag van de allerkleinste werveltjes (waar de energie verdwijnt als warmte) tot aan de grootste stromingen.
3. De data volgen: Als je de formules vergelijkt met de echte meetgegevens van het heliumgas, passen ze perfect. Het is alsof je een kaart hebt die precies de route van de rivier volgt, van bron tot monding.

5. Wat is er nog niet perfect? (De "Gaten" in de theorie)

Hoewel het model voor de twee belangrijkste soorten metingen (de tweede en derde orde) perfect werkt, is er nog een klein probleem:

• Het werkt geweldig voor de "evene" metingen (zoals de gemiddelde snelheidsverschillen).
• Voor de "oneven" metingen (waar de richting belangrijk is) is de link nog niet helemaal gevonden. Het is alsof het recept voor de soep perfect is, maar je weet nog niet precies hoe je de kruiden moet mengen voor de salade.
• Renner suggereert dat we misschien niet alleen naar de stroming in één richting hoeven te kijken, maar dat we de interactie tussen verschillende richtingen beter moeten begrijpen.

Samenvatting

Christoph Renner heeft een puzzelstukje gevonden dat ontbrak in de theorie van turbulente stroming. Hij heeft een brug geslagen tussen de chaos van de grote stroming en de rust van de kleine wervels.

• Vroeger: We hadden modellen die alleen op de snelweg werkten of alleen in de parkeergarage.
• Nu: We hebben één model dat de hele rit beschrijft, van het begin tot het einde, zonder dat we hoeven te gokken met instellingen.

Het is een stap dichter bij het volledig begrijpen van één van de meest complexe verschijnselen in de natuur: hoe vloeistoffen en gassen zich gedragen als ze in volle chaos verkeren.

Technische samenvatting

Probleemstelling

Turbulentie in vloeistoffen blijft een van de meest uitdagende problemen in de natuurkunde. Bestaande modellen, zoals die van Yakhot (1998), beschrijven succesvol de overgang van kleine schalen naar grote schalen (inertiebereik) en voorspellen correct het gedrag van oneven-orde structuurfuncties (zoals het afnemen naar nul bij grote schalen). Echter, deze modellen hebben twee belangrijke tekortkomingen:

1. Ze beschrijven niet de convergentie van even-orde structuurfuncties naar hun constante waarden op grote schalen.
2. Ze missen volledig de dissipatieve schalen (de kleinste schalen waar viskeuze krachten dominant worden). Yakhot's originele model en de recente uitbreiding voor grote schalen [9] zijn niet geldig in het dissipatiebereik, waar de schaalwetten fundamenteel anders zijn dan in het inertiebereik.

Het doel van dit artikel is om Yakhot's model uit te breiden tot een volledig schaalmodel (van dissipatieve schalen tot het systeemgrootte) voor even en oneven orde structuurfuncties, zonder gebruik te maken van vrije aanpasbare parameters.

Methodologie

De auteur baseert de analyse op experimentele data afkomstig van een cryogene, axiale heliumgasjet. De methode omvat de volgende stappen:

1. Dimensionalisering en Uitgangspunten:
   De vergelijkingen worden uitgedrukt in dimensieloze grootheden ($r = l/L$, $u = v/\sigma$, $S_n$), waarbij $L$ de systeemgrootte is en $\sigma$ de RMS van de snelheidsfluctuaties. Er wordt gebruikgemaakt van de eerder ontwikkelde grote-schaal uitbreiding van Yakhot's model, die extra termen ($D$ en $C$) introduceert om de convergentie op grote schalen te beschrijven.

2. Analyse van Experimentele Data:
   De auteur onderzoekt de relatie tussen de ruimtelijke afgeleiden van even-orde structuurfuncties ($\partial_r S_n$) en de volgende oneven-orde structuurfunctie ($S_{n+1}$).

   • Er wordt een hypothese getoetst (H2) dat de term $d(r)$ (die verantwoordelijk is voor symmetrie en grote schalen) geen bijdrage levert aan de dynamiek op zeer kleine schalen.
   • Door de residuen ($R_n$) van de bestaande vergelijkingen te analyseren, wordt een nieuwe empirische relatie gevonden voor even-orde functies.

3. Afleiding van een Nieuwe Schaalwet:
   Voor even-orde structuurfuncties ($n=2, 4$) wordt ontdekt dat de gecompenseerde residuen een machtswet volgen:
   $$-n \frac{R_n(r)}{S_{n+1}(r)} = \frac{\tau}{r^2}$$
   Hierbij is $\tau \approx 0.026$ een constante. Deze relatie geldt van het dissipatiebereik tot ver in het inertiebereik.

4. Sluiten van de Vergelijkingen:

   • Voor $n=2$ en $n=3$ wordt het stelsel gesloten met behulp van Kolmogorov's 4/5-wet (die een exacte relatie geeft tussen $S_3$ en de afgeleide van $S_2$).
   • Een nieuwe karakteristieke lengteschaal $\rho$ wordt geïntroduceerd, die de overgang markeert tussen het dissipatie- en inertiebereik.
   • De functie $d(r)$ wordt gekozen als een sigmoid-functie die overgaat van 0 (bij kleine schalen) naar 1 (bij grote schalen).

Belangrijkste Bijdragen

1. Empirische Relatie voor Even Orde: De ontdekking van de machtswet-relatie (vergelijking 22) tussen de afgeleide van even-orde structuurfuncties en de volgende oneven-orde structuurfunctie. Dit vult een cruciale ontbrekende schakel in Yakhot's model voor kleine schalen in.
2. Volledig Schaalmodel: De afleiding van gesloten, analytische uitdrukkingen voor de tweede en derde orde structuurfuncties ($S_2$ en $S_3$) die geldig zijn van de kleinste dissipatieve schalen ($r \to 0$) tot de systeemgrootte ($r \to 1$).
3. Geen Vrije Parameters: Het model bevat geen vrij aanpasbare parameters. De nieuwe schaalparameter $\rho$ (of $\tau$) wordt uitsluitend bepaald door de Reynolds-getal ($Re$), de dimensieloze dissipatiesnelheid ($\epsilon$) en de parameters van het grote-schaal model.
4. Karakteristieke Lengteschaal: De introductie van $\rho$ als de overgangsschaal van dissipatie naar inertie. De auteur toont aan dat $\rho$ een vergelijkbare afhankelijkheid van het Reynolds-getal heeft als de Kolmogorov-microschaal ($\eta$), maar een orde van grootte groter is, wat overeenkomt met het begin van het dissipatiebereik.

Resultaten

• Overeenkomst met Data: Het afgeleide model voor $S_2(r)$ en $S_3(r)$ toont uitstekende overeenkomst met de experimentele data over het volledige bereik van schalen.
  • In het dissipatiebereik ($r \ll \rho$) volgt het model correct de verwachte schaalwetten: $S_2 \propto r^2$ en $S_3 \propto r^3$.
  • In het inertiebereik gaat het model over in de bekende machtswetten met exponenten $\zeta_n$.
  • Op grote schalen convergeert het model correct naar de waarden voorspeld door het uitgebreide Yakhot-model.
• Validatie van Hypothesen: De aanname dat de term $d(r)$ verwaarloosbaar is op kleine schalen wordt bevestigd door de consistentie van het model met de kleine-schaal limieten en de experimentele data.
• Beperkingen: De gevonden relatie (22) geldt strikt genomen alleen voor even orde. Voor oneven orde ($n=3, 5$) vertonen de residuen een complexer gedrag dat niet door een eenvoudige machtswet wordt beschreven, waardoor het model voor willekeurige $n$ nog niet volledig gesloten is. Echter, voor de praktische toepassingen van $S_2$ en $S_3$ (via de 4/5-wet) is dit geen belemmering.
• Reynolds-getal Afhankelijkheid: Voor lagere Reynolds-getallen (beneden $Re \approx 2 \cdot 10^4$) beginnen afwijkingen in de machtswet (22) op te treden op grotere schalen, wat suggereert dat het model het beste werkt bij hoge Reynolds-getallen.

Betekenis

Dit werk vormt een aanzienlijke stap in de theoretische beschrijving van volledig ontwikkelde turbulentie. Door Yakhot's model uit te breiden naar het dissipatiebereik, biedt het voor het eerst een gesloten, analytisch model dat de volledige dynamiek van snelheidsincrementen beschrijft, van de kleinste viskeuze schalen tot de grootste systeemschalen, zonder gebruik te maken van empirische aanpassingen voor elke schaal.

De resultaten onderstrepen dat de overgang van dissipatie naar inertie kan worden gemodelleerd met een enkele karakteristieke lengteschaal $\rho$ die direct gerelateerd is aan het Reynolds-getal. Dit biedt een krachtig instrument voor het begrijpen van de structuur van turbulentie en kan nuttig zijn voor het verbeteren van turbulentiemodellen in CFD-simulaties (Large Eddy Simulation), waar de overgang tussen schalen vaak een uitdaging vormt.