Multi-domain spectral approach for Zakharov-Kuznetsov equations in 3D with cylindrical symmetry

Probleemstelling

Het artikel richt zich op het numeriek bestuderen van de drie-dimensionale (3D) gegeneraliseerde Zakharov-Kuznetsov (gZK) vergelijking:
$$u_t + (\Delta u + u^p)_x = 0$$
waarbij $u = u(t, x, y, z)$ een reële waarde heeft en $p > 1$ een niet-geheel getal of een geheel getal is. De auteurs focussen specifiek op de $L^2$-kritieke case in 3D, waarbij de macht van de niet-lineariteit $p = 7/3$ is.

De uitdagingen in dit probleem zijn:

1. Hoge dimensie en symmetrie: De vergelijking beschrijft golven in een magnetisch plasma met voorkeur voor de $x$-as. Door cilindrische symmetrie rond de $x$-as ($u = u(x, \rho)$ met $\rho = \sqrt{y^2+z^2}$) kan het probleem worden gereduceerd, maar de term $1/\rho$ in de Laplaciaan introduceert een singulariteit op de as ($\rho=0$).
2. Fractionele niet-lineariteit: De macht $p=7/3$ is een breuk. Dit veroorzaakt numerieke instabiliteiten, vooral wanneer de oplossing negatieve waarden aanneemt of wanneer complexe getallen ontstaan door Fourier-transformaties (bij het berekenen van $u^{7/3}$).
3. Dynamisch gedrag: In de kritieke regime is de stabiliteit van solitonoplossingen (grondtoestanden) een open vraag. Er moet worden bepaald of oplossingen globaal bestaan (dispersie) of in eindige tijd "blow-up" vertonen (instabiliteit), afhankelijk van de initiële massa in vergelijking met die van de grondtoestand.

Methodologie

De auteurs presenteren een nieuw numeriek raamwerk dat een multi-domein spectrale aanpak combineert met een robuuste tijdsintegratie.

1. Ruimtelijke Discretisatie:

   • In de $x$-richting: Er wordt gebruikgemaakt van een Fourier-spectrale methode op een torus (periodieke randvoorwaarden). Dit is efficiënt omdat de oplossing als lokaal gelokaliseerd wordt beschouwd.
   • In de $\rho$-richting (radiaal): Vanwege de singulariteit bij $\rho=0$ en de variabele coëfficiënten, wordt een Chebyshev-collocatiemethode gebruikt.
   • Multi-domein Strategie: Het domein wordt opgesplitst in twee gebieden om de singulariteit en de resolutie te optimaliseren:
     • Domein I ($\rho < \rho_0$): Hier wordt de variabele $s = \rho^2$ gebruikt. Dit verwijdert de singulariteit ($1/\rho$) en maakt de vergelijking minder singulier, wat cruciaal is voor de nauwkeurigheid op de as. De oplossing is hier analytisch in $\rho^2$.
     • Domein II ($\rho_0 < \rho < \rho_1$): Hier wordt de standaard variabele $\rho$ gebruikt.
   • Koppeling: De oplossingen in beide domeinen worden gekoppeld via $C^1$-continuïteitsvoorwaarden (continuïteit van de functie en de eerste afgeleide) op het interfacepunt $\rho_0$. Dit wordt geïmplementeerd via de $\tau$-methode.

2. Tijdsintegratie:

   • Vanwege de stijfheid van het systeem (door de derde afgeleide in $x$ en de Chebyshev-matrices) wordt een impliciete Runge-Kutta methode van de vierde orde (IRK4), ook wel het 2-staps Gauss-methode (Hammer-Hollingsworth), gebruikt.
   • Het resulterende niet-lineaire systeem wordt opgelost met een vereenvoudigde Newton-iteratie.

3. Omgaan met Fractionele Machten:

   • Om problemen met complexe getallen bij het berekenen van $u^{7/3}$ te voorkomen (vooral door kleine imaginaire delen na inverse FFT), wordt de oplossing na elke inverse transformatie geforceerd tot reële waarden.
   • De niet-lineariteit wordt berekend als $|u|^{7/3}\text{sign}(u)$ om de tak van de wortel op de negatieve reële as correct te behandelen.

Belangrijkste Bijdragen

• Novel Numeriek Kader: De ontwikkeling van een efficiënt algoritme dat specifiek is ontworpen voor 3D dispersieve vergelijkingen met cilindrische symmetrie en fractionele niet-lineariteiten.
• Omgaan met Singulariteiten: De succesvolle implementatie van een domein-decompositie met een veranderde variabele ($s=\rho^2$) om de singulariteit op de symmetrie-as te elimineren zonder verlies van nauwkeurigheid.
• Behandeling van $L^2$-kritieke gevallen: Het vermogen om de $p=7/3$ case nauwkeurig te simuleren, wat eerder als zeer uitdagend werd beschouwd vanwege numerieke instabiliteiten.

Resultaten

De auteurs hebben de methode getest op de constructie van de grondtoestand (soliton) en de dynamiek van perturbaties:

1. Grondtoestand Constructie:

   • De numeriek gevonden grondtoestand $Q$ heeft een massa $\|Q\|_{L^2} \approx 7.98$ (kwadraat $\approx 63.78$) en energie $E \approx 0$, wat consistent is met theoretische verwachtingen voor de kritieke case.
   • De spectrale coëfficiënten dalen exponentieel, wat aantoont dat de ruimtelijke resolutie zeer hoog is.

2. Stabiliteit en Blow-up:

   • Subkritieke massa ($\lambda < 1$): Voor initiële data $u_0 = \lambda Q$ met $\lambda = 0.99$ (kleinere massa dan de grondtoestand), disperseren de oplossingen en verdwijnen ze naar oneindig. De $L^\infty$-norm neemt monotoon af.
   • Kritieke/Superkritieke massa ($\lambda > 1$):
     • Bij kleine perturbaties ($\lambda = 1.01$) op een eindig domein lijkt de oplossing stabiel te zijn door reflecties van straling, maar dit is een artefact van het eindige domein.
     • Bij grotere perturbaties ($\lambda = 1.1$) treedt er duidelijk blow-up in eindige tijd op. De $L^\infty$-norm groeit exponentieel.
   • Gaussische Data: Voor Gaussische initiële data werd de drempelwaarde bevestigd. Data met massa onder de drempel ($\lambda=5$) disperseren, terwijl data boven de drempel ($\lambda=6.5$) in eindige tijd instabiel worden en blow-up vertonen.

3. Blow-up Profiel:

   • De simulaties tonen aan dat bij blow-up de energie zich concentreert, terwijl er ook straling wordt uitgestraald in een kegelvormige structuur rond de negatieve $x$-as met een opening hoek van ongeveer $60^\circ$, wat overeenkomt met eerdere theoretische bevindingen.

Significantie

Deze studie is van groot belang voor de plasmafysica en de wiskundige analyse van niet-lineaire dispersieve vergelijkingen:

• Voorspelling van Instabiliteit: Het bevestigt dat voor de 3D kritieke ZK-vergelijking de massa van de grondtoestand de scherpe drempel vormt tussen globale bestaan en eindige-tijd instabiliteit (blow-up).
• Plasma-toepassingen: De resultaten bieden inzicht in het gedrag van golven in gemagnetiseerde plasma's (zoals in laboratoria of de ruimte), waar de balans tussen dispersie en niet-lineariteit kan leiden tot energie-localisatie en instabiliteit.
• Methodologische Vooruitgang: De voorgestelde multi-domein spectrale methode biedt een robuust instrument voor het bestuderen van andere complexe, niet-lineaire problemen in hogere dimensies met symmetrie, waar traditionele methoden vaak falen door singulariteiten of fractionele machten.

Samenvattend biedt dit artikel een krachtige numerieke oplossing voor een langdurig probleem in de analyse van 3D niet-lineaire golven en levert het overtuigend bewijs voor de stabiliteitsdrempel in de kritieke regime.