Numerical study of the 2D Kaup-Broer-Kuperschmidt Boussinesq system

Technische Samenvatting: Numerieke Studie van het 2D Kaup-Broer-Kuperschmidt Boussinesq Systeem

Probleemstelling
Dit werk onderzoekt de goedgestelde versie van het tweedimensionale (2D) Kaup-Broer-Kuperschmidt (KBK) systeem, een specifiek geval van de Boussinesq 'abcd' systemen. Waar het 1D KBK-systeem volledig integreerbaar is en bekende solitaire golfoplossingen bezit met gevestigde orbitale stabiliteit, blijven de kwalitatieve eigenschappen van de 2D-generalisatie grotendeels onontgonnen. Het 2D-systeem is niet-integreerbaar, en er ontbreken rigoureuze resultaten betreffende gedrag op lange termijn, globale existentie of eindtijd-blow-up. Het primaire doel is om numeriek te bepalen of stabiele, gelokaliseerde structuren (zoals 2D-solitons of lijnsolitons) bestaan in het 2D KBK-systeem en om de evolutie van algemene gelokaliseerde beginvoorwaarden te karakteriseren.

Methodologie
De auteurs maken gebruik van een hoogprecisie numerieke aanpak die spectrale methoden in de ruimte combineert met exponentiële tijdsdifferentiatie (ETD) schema's in de tijd.

• Formulering: Het systeem wordt herschreven met behulp van een snelheidspotentiaal $V$ (waarbij $v = \nabla V$) en getransformeerd naar een systeem van fractionele niet-lineaire Schrödinger (NLS) vergelijkingen in de Fourier-ruimte. Deze formulering faciliteert het gebruik van exponentiële integratoren, specifiek het vierde-orde ETD-schema van Cox en Matthews, wat noodzakelijk is om de stijfheid van het systeem te hanteren die voortvloeit uit dispersie van hoge orde.
• Statische Oplossingen: Om gelokaliseerde statische oplossingen te construeren, lossen de auteurs eerst de bijbehorende radiale gewone differentiaalvergelijking (ODE) in poolcoördinaten op met behulp van een Chebyshev-spectrale discretisatie en een Newton-Krylov-methode met relaxatie. Deze 1D-radiële oplossing wordt vervolgens geïnterpoleerd op een 2D Fourier-rooster om te dienen als een beginvoorstel voor het oplossen van het volledige 2D niet-lineaire systeem.
• Tijds-evolutie: De tijds-evolutie wordt gesimuleerd op een dubbel periodiek domein (2-torus) om snel afnemende functies te benaderen. De auteurs bewaken het behoud van relatieve energie ($\Delta E$) om numerieke nauwkeurigheid te waarborgen, met als doel fouten die aanzienlijk onder de $10^{-3}$ liggen.
• Blow-up Analyse: Om potentiële singulariteiten te onderzoeken, maken de auteurs gebruik van een dynamische herschalingstechniek. Dit houdt in dat ruimtelijke variabelen en tijd worden herschaald om het gedrag van oplossingen in de buurt van een potentiële blow-up-tijd $t^*$ te volgen. De schalingsfactor $L(t)$ en de exponenten van de $L^\infty$- en $L^2$-normen worden aangepast aan machtswetten van de vorm $(t^* - t)^\alpha$ om de aard van de blow-up te bepalen (algebraïsch versus exponentieel).

Belangrijkste Bijdragen en Resultaten

1. Existentie van een Statische Oplossing: De auteurs construeren succesvol een radiaal symmetrische, gelokaliseerde statische oplossing voor het 2D KBK-systeem. Deze oplossing komt overeen met een grondtoestand van een vector NLS-vergelijking.
2. Instabiliteit van Statische Oplossingen: Numerieke perturbaties van deze statische oplossing onthullen een dichotomie in gedrag:
   • Dispersie: Perturbaties met een massa (of $L^2$-norm) iets kleiner dan de statische oplossing leiden tot dispersie van de golf naar oneindig.
   • Eindtijd-blow-up: Perturbaties met een iets grotere massa resulteren in een snel zelf-focuserend proces dat leidt tot eindtijd-blow-up. De numerieke data suggereert dat de blow-up zelf-similaar is, waarbij de schalingsfactor $L$ zich gedraagt als $L \propto (t^* - t)^{1/2}$, en de normen groeien als $\|\eta\| \sim (t^* - t)^\alpha$ met $\alpha \approx -1$. Dit gedrag is analoog aan het focuserend regime in $L^2$-supercritische NLS-vergelijkingen.
3. Transversale Instabiliteit van Lijnsolitons: De studie onderzoekt 1D-lijnsolitons (oplossingen die onafhankelijk zijn van de transversale coördinaat $y$) die triviaal tot 2D zijn uitgebreid. Deze structuren blijken sterk instabiel te zijn tegen transversale perturbaties. Zowel positieve als negatieve perturbaties leiden tot eindtijd-blow-up, gekenmerkt door de vorming van scherpe pieken en zelf-similaire schaling vergelijkbaar met het geval van de 2D-statische oplossing. Opmerkelijk is dat blow-up optreedt zelfs wanneer de geperturbeerde massa kleiner is dan die van de geconstrueerde statische oplossing, een gedrag dat verschilt van de $L^2$-kritieke NLS-vergelijking.
4. Algemene Gelokaliseerde Data: Simulaties met algemene Gaussische beginvoorwaarden (waarbij geen specifieke statische oplossing als basis wordt gebruikt) tonen aan dat het systeem geen stabiele gelokaliseerde structuren ondersteunt. Afhankelijk van de parameters, ofwel dispergeert de data of ondergaat deze blow-up.

Betekenis en Beweringen
Het artikel concludeert dat het 2D KBK-systeem geen stabiele gelokaliseerde structuren ondersteunt. De auteurs stellen het volgende:

• Oplossingen voor gelokaliseerde gladde beginvoorwaarden kunnen in eindige tijd blow-up ondergaan via een zelf-similaar mechanisme, wat het gedrag van supercritische NLS-vergelijkingen weerspiegelt.
• Lijn-solitaire golven zijn sterk instabiel en vatbaar voor blow-up onder transversale perturbaties.
• Er zijn geen stabiele structuren die in twee dimensies gelokaliseerd zijn voor dit systeem.

De auteurs stellen expliciet dat de precieze criteria voor data om te blow-upen en de exacte profielen van de blow-up-oplossingen ($U^\infty$ en $V^\infty$) onbekend blijven. Het werk dient als een numerieke karakterisering van de dynamiek van het systeem, waarbij de afwezigheid van stabiliteit in 2D wordt benadrukt, ondanks de integreerbaarheid en stabiliteitseigenschappen van het 1D-tegenhanger.