Inviscid scaling in the Kuramoto-Sivashinsky equation from functional renormalization group and direct numerical simulations

Technische Samenvatting: Invisscale Scaling in de Kuramoto-Sivashinsky-vergelijking

Probleemstelling
De Kuramoto-Sivashinsky (KS)-vergelijking beschrijft de dynamica van een één-dimensionaal scalair veld $h(t, x)$ dat wordt beheerst door negatieve viscositeit ($\nu < 0$), positieve vierde-orde dissipatie ($\tau > 0$) en een niet-lineaire term. Hoewel het gedrag op grote schaal (infrarood) van de KS-vergelijking algemeen wordt geaccepteerd als behorend tot de Kardar-Parisi-Zhang (KPZ)-universaliteitsklasse met een dynamische exponent $z = 3/2$, zijn de schalingsregimes op intermediaire schaal minder goed begrepen. Eerdere numerieke studies rapporteerden vaak diffussieve schaling ($z=2$) geassocieerd met de Edwards-Wilkinson-vergelijking, een resultaat dat wordt toegeschreven aan eindige-grootte-effecten in deterministische simulaties waarbij de effectieve ruisamplitude aanvankelijk nul is. Het specifieke gedrag van het systeem naarmate de effectieve viscositeit evolueert van zijn microscopische negatieve waarde naar een macroscopische positieve waarde, en het potentiële bestaan van een distinct schalingsregime op intermediaire schalen, waren niet volledig gekarakteriseerd.

Methodologie
De auteurs hanteren twee complementaire benaderingen om de schalingsregimes van de KS-vergelijking te onderzoeken:

1. Functionele Renormalisatiegroep (FRG):

   • De studie maakt gebruik van het niet-perturbatieve FRG-formalisme, specifiek de Next-to-Leading-Order (NLO)-benadering die is ontwikkeld voor de KPZ-vergelijking.
   • Om FRG toe te passen op de deterministische KS-vergelijking, wordt een stochastische ruis term geïntroduceerd (volgend op eerdere werken [19, 20]) om de gemiddelde over beginvoorwaarden te vervangen door een gemiddelde over ruisrealisaties. De limiet van een verdwijnende ruisamplitude wordt opgemerkt als onderwerp voor toekomstig werk, maar de ruis blijkt geen invloed te hebben op intermediaire schalen.
   • De methode volgt de evolutie van schaalafhankelijke functies voor effectieve viscositeit ($\nu_\kappa$) en ruisamplitude ($D_\kappa$) naarmate de impulschaal $\kappa$ stroomt van het ultraviolette (microscopische) naar het infrarode (macroscopische) gebied.
   • De stroomvergelijkingen worden numeriek opgelost met een specifiek schema dat gekoppelde roosters voor dimensieloze en dimensie-bevattende grootheden omvat, waardoor de beschrijving van het volledige bereik aan golfgetallen en frequenties mogelijk wordt.

2. Directe Numerieke Simulaties (DNS):

   • De 1D KS-vergelijking wordt numeriek geïntegreerd met behulp van een pseudo-spectrale methode gekoppeld aan een Exponential Time Differencing vierde-orde Runge-Kutta (ETDRK4)-schema.
   • Simulaties worden uitgevoerd op een periodiek domein met grootte $L = 2^{17}$ met $2^{18}$ roosterpunten, waarbij GPU-versnelling (NVIDIA Tesla P100) wordt gebruikt om de rekenlast te hanteren.
   • De studie genereert 102 onafhankelijke realisaties beginnend bij willekeurige Gaussische begindata. De ruimtetijd-correlatiefunctie $C(t, p)$ wordt berekend door te middelen over de stationaire toestand en vervolgens over verschillende realisaties om statistische convergentie te waarborgen.

Belangrijkste Bijdragen en Resultaten

• Ontdekking van een Intermediair Schalingsregime: De belangrijkste bevinding is het bestaan van een robuust schalingsregime op intermediaire schalen, gekenmerkt door een dynamische exponent $z = 1$. Dit regime verschijnt tussen de schaling op grote schaal (KPZ, $z = 3/2$) en het niet-universele gedrag op kleine schaal (in de buurt van de meest onstabiele modus).
• Mechanisme van Invisscale Scaling: Het $z=1$-regime is intrinsiek voor de KS-dynamica. Het ontstaat specifiek wanneer de effectieve viscositeit $\nu_\kappa$, evoluerend vanuit de microscopische negatieve waarde ($\nu < 0$), door nul gaat voordat de macroscopische positieve KPZ-waarde wordt bereikt. Op de schaal waar $\nu_\kappa \approx 0$, worden de correlaties van het systeem beheerst door een effectieve vergelijking met verdwijnende viscositeit.
• Identificatie van de Universaliteitsklasse: De auteurs identificeren dit $z=1$-regime met de Inviscid Burgers (IB)-universaliteitsklasse. Dit komt overeen met het nulpunt-viscositeit vaste punt van de KPZ-vergelijking, distinct van het standaard KPZ-vaste punt.
  • FRG-bewijs: De stroom van de effectieve viscositeit kruist nul, en de correlatiefuncties in het intermediaire bereik komen overeen met de schalingsfuncties die zijn voorspeld voor het IB-vaste punt. De analyse toont aan dat tijd-omkeer-symmetrie, die op microscopisch KS-niveau is verbroken, op grote afstanden weer optreedt, consistent met het KPZ-vaste punt, maar dat het intermediaire regime wordt gedomineerd door het IB-vaste punt.
  • DNS-bewijs: Numerieke data voor de twee-punts correlatiefunctie $C(t, p)$ en decorrelatietijden $\tau_\alpha(p)$ vertonen duidelijk een gebied waar de schalingsvariabele $pt$ is (implicerend $z=1$). De data collapseert op de exacte asymptotische schalingsfunctie die is afgeleid voor het IB-vaste punt (een Gaussische vorm voor kleine tijden), distinct van de Prah¨ofer-Spohn schalingsfunctie van het KPZ-vaste punt.
• Robuustheid: De $z=1$-schaling wordt waargenomen ongeacht het uiteindelijke infrarode regime (of het KPZ of Edwards-Wilkinson is) en vereist geen fijnafstelling van parameters. Het beheerst een uitgebreid bereik van intermediaire impulsen tot aan de meest onstabiele modus.

Betekenis en Claims
Het artikel claimt het eerste bewijs en de karakterisering te leveren van dit "tot nu toe over het hoofd geziene" schalingsregime in de KS-vergelijking. Door FRG en DNS te combineren, tonen de auteurs aan dat het verdwijnen van de effectieve viscositeit een universele $z=1$-schaling afdrukt op correlaties op intermediaire schalen.

De auteurs stellen dat dit gedrag behoort tot de Inviscid Burgers universaliteitsklasse, wat overeenkomt met het nulpunt-viscositeit vaste punt van de KPZ-vergelijking. Zij benadrukken dat dit regime generiek is voor KS-dynamica, natuurlijk voortvloeiend uit de overgang van negatieve naar positieve effectieve viscositeit, en niet het gevolg is van artefacten van systeemgrootte of ruisintroductie. Het werk suggereert dat eerdere Renormalisatiegroep-analyses van de benadering van het infrarode gedrag opnieuw moeten worden bezien in het licht van dit intermediaire IB-regime. De bevindingen worden gepresenteerd als een directe karakterisering van de KS-vergelijking zelf, ondersteund door zowel theoretische stroomanalyse als hoog-precisie numerieke simulaties.