Numerical study of Lagrangian velocity structure functions using acceleration statistics and a spatial-temporal perspective

Dit artikel gebruikt directe numerieke simulaties van isotrope turbulentie om te tonen dat de waargenomen gedragingen van de tweede-orde Lagrangiaanse snelheidsstructuurfunctie sterk worden beïnvloed door zowel het beperkte tijdsbereik als de deeltjesverplaatsingen, waarbij de schaalconstante $C_0$ mogelijk gevoelig is voor Lagrangiaanse intermitentie maar asymptotisch constant kan blijven bij nog hogere Reynoldsgetallen.

De kern

De dans van de turbulentie: Waarom waterdeeltjes zich soms raar gedragen

Stel je voor dat je in een drukke, chaotische menigte loopt. Je probeert een vriend te vinden, maar iedereen duwt en duwt. Soms word je snel vooruit geduwd, soms sta je stil, en soms word je teruggeblazen. Dit is wat er gebeurt in een vloeistof die turbulent is, zoals een woelige rivier of een stormachtige luchtstroom.

Wetenschappers proberen al decennia lang de regels van dit chaos te begrijpen. Een beroemde theorie uit de jaren '40 (van Kolmogorov) zegt dat als je heel goed kijkt naar hoe snel een deeltje versnelt of vertraagt over een korte tijd, er een heel simpel, voorspelbaar patroon zou moeten zijn. Het zou moeten lijken op een rechte lijn op een grafiek.

Maar in de praktijk is het niet zo simpel. De grafieken zijn vaak krom, hebben pieken en gedragen zich niet zoals de theorie voorspelt.

In dit nieuwe onderzoek kijken twee wetenschappers van de Georgia Tech (Rohini Uma-Vaideswaran en P. K. Yeung) naar waarom dit zo is. Ze hebben supercomputers gebruikt om een virtuele, wervelende vloeistof te simuleren. Hier is hun ontdekking, vertaald in begrijpelijke taal:

1. Het probleem: De "piek" in plaats van de "vlakte"

De theorie voorspelt een vlakte (een plateau): een periode waarin het gedrag van de deeltjes stabiel en voorspelbaar is.
De computersimulaties tonen echter vaak een piek: het gedrag verandert snel, bereikt een hoogtepunt en daalt dan weer. Het is alsof je verwacht dat een auto met constante snelheid rijdt, maar in plaats daarvan ziet je dat hij constant accelereert en remt.

2. De eerste verklaring: De versnellingskracht (De "Remmen en Gas")

De auteurs kijken naar de versnelling van de deeltjes. Stel je voor dat een deeltje een auto is.

• De snelheid van de auto is de snelheid.
• Het gaspedaal en de rem zijn de versnelling.

Ze ontdekten dat de "piek" in de grafiek direct samenhangt met hoe sterk de versnelling fluctueert. Bij hoge turbulentie (veel chaos) is de versnelling extreem onvoorspelbaar. Het deeltje wordt soms enorm hard weggeblazen en dan weer abrupt gestopt. Deze extreme schokken zorgen ervoor dat het "stabiele patroon" dat we zoeken, nooit lang genoeg duurt om een echte rechte lijn te vormen. Het is alsof je probeert een rechte lijn te tekenen terwijl je hand constant schokt.

3. De tweede verklaring: De dubbele beweging (De "Vervoerders")

Dit is het meest interessante deel. Een deeltje in een stroming verandert van snelheid om twee redenen tegelijk:

1. De lokale verandering: De stroming op die specifieke plek verandert van snelheid (alsof de wind plotseling harder waait op je positie).
2. De verplaatsing: Het deeltje beweegt naar een nieuwe plek waar de stroming anders is (alsof je van een rustige straat naar een drukke kruising loopt).

In de natuurkunde noemen ze dit lokaal en convectief.
Het geheim is dat deze twee krachten vaak tegen elkaar werken.

• Stel je voor dat je op een roterende carrousel staat. Als je naar voren loopt (verplaatsing), voel je een duw. Maar als de carrousel op dat moment vertraagt (lokale verandering), voel je een duw in de tegenovergestelde richting.
• In de turbulentie heffen deze twee krachten elkaar bijna op, maar niet helemaal.

De auteurs ontdekten dat deze "bijna-annihilatie" de oorzaak is van de rare pieken. Omdat ze elkaar bijna opheffen, is het resultaat heel klein, maar omdat ze niet perfect opheffen, blijven er kleine, extreme restjes over. Deze restjes zorgen voor de intermittentie (het plotselinge, extreme gedrag) dat we zien.

4. De reis van het deeltje: Hoe ver komt het?

Een ander belangrijk punt is hoe ver een deeltje reist in een heel korte tijd.

• Bij lage turbulentie (rustig weer) blijft een deeltje dicht bij zijn startpunt.
• Bij hoge turbulentie (storm) kan een deeltje in een fractie van een seconde al een enorme afstand afleggen, ver buiten de "rustige zone".

De auteurs tonen aan dat deeltjes zo snel verplaatsen dat ze binnen een oogwenk de "inertiale zone" bereiken (een zone waar de wetten van de grote stroming gelden). Maar omdat ze zo snel verplaatsen, verlaten ze die zone ook weer heel snel.
De analogie: Het is alsof je probeert een foto te maken van een raket die net is opgestart. Je wilt de fase zien waarin het stabiel vliegt, maar de raket is zo snel dat hij die fase al voorbij is voordat je de camera hebt ingeschakeld. De "stabiele zone" is zo kort, dat we hem in onze metingen vaak missen of alleen een piek zien.

Conclusie: Waarom is dit belangrijk?

Deze studie legt uit waarom we in de natuurkunde nog steeds worstelen met het voorspellen van turbulentie.

1. Tijd is te kort: De tijd waarin de "wetten" gelden, is zo kort dat het moeilijk te meten is.
2. Beweging is te snel: Deeltjes verplaatsen zich zo snel dat ze constant van omgeving wisselen, waardoor de lokale regels verwarren.
3. Krachten heffen elkaar bijna op: De twee soorten beweging (lokaal en verplaatsing) vechten tegen elkaar, wat zorgt voor extreme, onvoorspelbare schokken.

Kortom: Turbulentie is niet alleen een kwestie van "hoe snel gaat het?", maar ook van "hoe ver is het gegaan en hoe snel is de omgeving veranderd?". De wetenschappers concluderen dat we waarschijnlijk nog veel langere en krachtigere computers nodig zullen hebben om het perfecte, stabiele patroon te vinden dat de oude theorie voorspelde. Tot die tijd blijft turbulentie een fascinerend, chaotisch dansje.

Technische samenvatting

Titel

Numerisch onderzoek naar Lagrangiaanse snelheidsstructuurfuncties met behulp van versnellingsstatistieken en een ruimtetijdperspectief.

1. Het Probleem

In de theorie van turbulente stromingen, specifiek binnen het kader van de Kolmogorov-theorie (K41), wordt voorspeld dat de tweede-orde Lagrangiaanse snelheidsstructuurfunctie, $D_2^L(\tau) = \langle (\Delta_\tau u^+)^2 \rangle$, lineair schaalt met de tijdvertraging $\tau$ in het trajectgebied (inertial range) bij voldoende hoge Reynolds-getallen:
$$\frac{1}{3}\langle (\Delta_\tau u^+)^2 \rangle = C_0 \langle \epsilon \rangle \tau$$
Hierbij is $C_0$ de universele Lagrangiaanse Kolmogorov-constante.

Echter, in praktische Directe Numerieke Simulaties (DNS) en experimenten vertonen de data geen duidelijk plateau voor $C_0$. In plaats daarvan zien onderzoekers een gladde lokale piek in het bereik van 5 tot 10 Kolmogorov-tijdschalen ($\tau_\eta$), gevolgd door een snelle afname. De oorzaken hiervan zijn niet volledig begrepen. Uitdagingen omvatten:

• De beperkte duur van simulaties bij hoge Reynolds-getallen, wat statistische onzekerheid introduceert.
• De sterke intermitentie (extreme waarden) van Lagrangiaanse versnellingen.
• De complexiteit van het onderscheid tussen lokale tijdsveranderingen en convectieve veranderingen door deeltjesverplaatsing.

2. Methodologie

De auteurs hebben gebruik gemaakt van Directe Numerieke Simulaties (DNS) van gedwongen isotrope turbulentie.

• Reynolds-getallen: Het bereik van Taylor-schaal Reynolds-getallen ($R_\lambda$) loopt van 140 tot 1300. Dit is een significant bereik, waarbij de hoogste waarden drie keer zo hoog zijn als in veel eerdere literatuur.
• Oplossing: Simulaties zijn uitgevoerd op roosters tot $12288^3$ punten op de supercomputer Frontera (TACC), met tot $10^9$ deeltjes voor statistische nauwkeurigheid.
• Twee complementaire benaderingen:
  1. Versnellingsstatistieken: De snelheidsverandering wordt geanalyseerd als de dubbele integraal van de versnellingsautocorrelatie. Dit omzeilt de noodzaak van zeer lange tijdsreeksen voor de structuurfunctie zelf, omdat versnelling kortere tijdschalen heeft.
  2. Ruimtetijdperspectief: De Lagrangiaanse snelheidsverandering $\Delta_\tau u^+$ wordt ontbonden in twee componenten (volgens vergelijking 5):
     • Lokaal ($v_L$): De verandering van het snelheidsveld op een vaste positie in de tijd ($\partial u / \partial t$).
     • Convectief ($v_C$): De verandering door de verplaatsing van het deeltje naar een nieuwe locatie ($u \cdot \nabla u$).
       De auteurs analyseren hoe deze twee componenten elkaar wederzijds opheffen (canceleren) en hoe de deeltjesverplaatsing $\ell(\tau)$ de statistieken beïnvloedt.

3. Belangrijkste Bijdragen en Resultaten

A. Rol van Versnellingsautocorrelatie en $C_0$

• De piekwaarde van de genormaliseerde structuurfunctie ($C_0^*$) is sterk gerelateerd aan de genormaliseerde versnellingsvariantie ($a_0$).
• De data toont aan dat $a_0$ toeneemt met het Reynolds-getal (ongeveer als $R_\lambda^{0.25}$) als gevolg van intermitentie.
• $C_0^*$ neemt ook toe met $R_\lambda$, maar langzamer ($R_\lambda^{0.2}$).
• Conclusie over $C_0$: Het is onwaarschijnlijk dat $C_0$ een asymptotische constante bereikt binnen de bereikbare Reynolds-getallen van DNS. Als dit wel gebeurt, zou dit gebeuren bij Reynolds-getallen die ver buiten het bereik van toekomstige simulaties liggen, en zou de waarde waarschijnlijk hoger zijn dan eerdere schattingen (rond 8 in plaats van 7).
• Tijdsduur: Een simulatietijd van ongeveer 10 $\tau_\eta$ is voldoende om betrouwbare resultaten voor $C_0^*$ te verkrijgen, ondanks de beperkte duur ten opzichte van de integrale tijdschaal $T_L$.

B. Ruimtetijdperspectief en Wederzijdse Opheffing

• Er is een sterke, maar onvolledige, wederzijdse opheffing tussen de convectieve ($v_C$) en lokale ($v_L$) bijdragen. Dit bevestigt het "random-sweeping"-hypothese van Tennekes.
• De correlatie tussen $v_C$ en $v_L$ begint dicht bij -1,0 (volledige tegenstelling) bij kleine tijdsvertragingen, maar zwakt af naarmate $\tau$ toeneemt.
• Deze onvolledige opheffing is de fysieke oorzaak van de sterke intermitentie in de totale Lagrangiaanse snelheidsverandering. Extreme waarden ontstaan wanneer grote waarden van $v_C$ en $v_L$ niet perfect elkaar opheffen of zelfs dezelfde teken hebben.

C. De Rol van Deeltjesverplaatsing

• De verplaatsing van een deeltje $\ell(\tau)$ volgt een Chi-verdeling (orde 3).
• Bij hoge Reynolds-getallen kunnen deeltjes zich in zeer korte tijd (enkele $\tau_\eta$) verplaatsen over afstanden die veel groter zijn dan de Kolmogorov-schaal $\eta$, en zelfs het trajectgebied (inertial range) binnendringen.
• Schaling: Als men kijkt naar de convectieve bijdrage $v_C$ voor deeltjes die een specifieke afstand hebben afgelegd (binnen het trajectgebied), vertoont deze een duidelijke schaling van $\tau^{2/3}$ (overeenkomend met Euleriaanse ruimtelijke schaling).
• Beperking: Omdat de totale structuurfunctie een mix is van deeltjes met verschillende verplaatsingen, en omdat de verplaatsing snel de trajectgebiedgrenzen overschrijdt, wordt het schalingsbereik voor de totale structuurfunctie snel "weggespoeld". Dit verklaart waarom het plateau van $C_0$ zo kort is en snel verdwijnt.

4. Significatie en Conclusies

Dit werk biedt nieuwe inzichten in waarom de klassieke Lagrangiaanse Kolmogorov-schaling ($D_2^L \propto \tau$) in simulaties zo moeilijk te observeren is als een stabiel plateau:

1. Beperkte tijdschaal: Het bereik van tijdschalen in DNS is smaller dan dat van lengteschalen, wat het vinden van een breed inertieel gebied bemoeilijkt.
2. Deeltjesverplaatsing: De snelle verplaatsing van deeltjes naar het trajectgebied en daarbuiten zorgt ervoor dat de convectieve bijdrage snel de Euleriaanse schaling ($\tau^{2/3}$) vertoont, maar dat de totale Lagrangiaanse grootheid door de complexe interactie en onvolledige opheffing met de lokale term een ander gedrag vertoont.
3. Intermitentie: De onvolledige opheffing tussen convectieve en lokale termen is de drijvende kracht achter de extreme waarden in Lagrangiaanse statistieken.

De studie concludeert dat zowel de beperkte tijdschalen als de dynamiek van deeltjesverplaatsing cruciale rollen spelen in het gedrag van de tweede-orde Lagrangiaanse snelheidsstructuurfunctie. Dit onderstreept de noodzaak van een ruimtetijdperspectief om Lagrangiaanse turbulentie volledig te begrijpen, en suggereert dat de universele constante $C_0$ mogelijk niet binnen het bereik van huidige of nabije toekomstige simulaties kan worden vastgesteld.