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

Technische Zusammenfassung: Numerische Studie des 2D Kaup-Broer-Kuperschmidt-Boussinesq-Systems

Problemstellung
Diese Arbeit untersucht die wohlgestellte Version des zweidimensionalen (2D) Kaup-Broer-Kuperschmidt (KBK)-Systems, einen spezifischen Fall der Boussinesq-„abcd"-Systeme. Während das 1D-KBK-System vollständig integrabel ist und bekannte Solitärwellenlösungen mit etablierter orbitale Stabilität besitzt, bleiben die qualitativen Eigenschaften seiner 2D-Verallgemeinerung weitgehend unerforscht. Das 2D-System ist nicht integrabel, und es liegen keine rigorosen Ergebnisse bezüglich des Langzeitverhaltens, der globalen Existenz oder des Endzeit-Blow-ups vor. Das Hauptziel besteht darin, numerisch zu bestimmen, ob stabile, lokalisierte Strukturen (wie 2D-Solitonen oder Liniensolitonen) im 2D-KBK-System existieren, und die Entwicklung allgemeiner lokalisierter Anfangsdaten zu charakterisieren.

Methodik
Die Autoren wenden einen hochpräzisen numerischen Ansatz an, der Spektralmethoden im Raum mit exponentiellen Zeitschrittverfahren (ETD) in der Zeit kombiniert.

• Formulierung: Das System wird unter Verwendung eines Geschwindigkeitspotenzials $V$ (wobei $v = \nabla V$) umformuliert und in ein System fraktionaler nichtlinearer Schrödinger-Gleichungen (NLS) im Fourier-Raum transformiert. Diese Formulierung ermöglicht die Verwendung exponentieller Integratoren, speziell des ETD-Schemas vierter Ordnung von Cox und Matthews, das notwendig ist, um die Steifigkeit des Systems zu handhaben, die durch hochordentliche Dispersion entsteht.
• Statische Lösungen: Um lokalisierte statische Lösungen zu konstruieren, lösen die Autoren zunächst die zugehörige radiale gewöhnliche Differentialgleichung (ODE) in Polarkoordinaten unter Verwendung einer Chebyshev-Spektraldiskretisierung und einer Newton-Krylov-Methode mit Relaxation. Diese 1D-radiale Lösung wird dann auf ein 2D-Fourier-Gitter interpoliert, um als Anfangsschätzung für die Lösung des vollen 2D-nichtlinearen Systems zu dienen.
• Zeitentwicklung: Die Zeitentwicklung wird auf einem doppelt periodischen Gebiet (2-Torus) simuliert, um schnell abklingende Funktionen zu approximieren. Die Autoren überwachen die relative Energieerhaltung ($\Delta E$), um die numerische Genauigkeit sicherzustellen, mit dem Ziel, Fehler deutlich unter $10^{-3}$ zu halten.
• Blow-up-Analyse: Um potenzielle Singularitäten zu untersuchen, nutzen die Autoren eine dynamische Reskalierungstechnik. Dies beinhaltet das Reskalieren von Raumvariablen und Zeit, um das Verhalten von Lösungen in der Nähe einer potenziellen Blow-up-Zeit $t^*$ zu verfolgen. Der Skalierungsfaktor $L(t)$ und die Exponenten der $L^\infty$- und $L^2$-Normen werden an Potenzgesetze der Form $(t^* - t)^\alpha$ angepasst, um die Natur des Blow-ups (algebraisch vs. exponentiell) zu bestimmen.

Hauptbeiträge und Ergebnisse

1. Existenz einer statischen Lösung: Die Autoren konstruieren erfolgreich eine rotationssymmetrische, lokalisierte statische Lösung des 2D-KBK-Systems. Diese Lösung entspricht einem Grundzustand einer vektoriellen NLS-Gleichung.
2. Instabilität statischer Lösungen: Numerische Störungen dieser statischen Lösung offenbaren eine Dichotomie im Verhalten:
   • Dispersion: Störungen mit einer Masse (oder $L^2$-Norm), die geringfügig kleiner als die der statischen Lösung ist, führen zur Dispersion der Welle ins Unendliche.
   • Endzeit-Blow-up: Störungen mit einer geringfügig größeren Masse resultieren in einem schnellen Selbstfokussierungsprozess, der zu einem Endzeit-Blow-up führt. Die numerischen Daten deuten darauf hin, dass der Blow-up selbstähnlich ist, wobei der Skalierungsfaktor $L$ sich wie $L \propto (t^* - t)^{1/2}$ verhält und die Normen wie $\|\eta\| \sim (t^* - t)^\alpha$ mit $\alpha \approx -1$ wachsen. Dieses Verhalten ist analog zum Fokussierungsregime in $L^2$-superkritischen NLS-Gleichungen.
3. Transversale Instabilität von Liniensolitonen: Die Studie untersucht 1D-Liniensolitonen (Lösungen, die unabhängig von der transversalen Koordinate $y$ sind), die trivial auf 2D erweitert werden. Diese Strukturen erweisen sich als stark instabil gegenüber transversalen Störungen. Sowohl positive als auch negative Störungen führen zu einem Endzeit-Blow-up, der durch die Bildung scharfer Spitzen und selbstähnliche Skalierung gekennzeichnet ist, ähnlich wie im Fall der 2D-statischen Lösung. Bemerkenswerterweise tritt der Blow-up sogar dann auf, wenn die gestörte Masse kleiner ist als die der konstruierten statischen Lösung, ein Verhalten, das sich von der $L^2$-kritischen NLS-Gleichung unterscheidet.
4. Allgemeine lokalisierte Daten: Simulationen mit allgemeinen Gaußschen Anfangsdaten (bei denen keine spezifische statische Lösung als Basis verwendet wird) zeigen, dass das System keine stabilen lokalisierten Strukturen unterstützt. Je nach Parametern dispergiert die Daten entweder oder erleidet einen Blow-up.

Bedeutung und Behauptungen
Die Arbeit kommt zu dem Schluss, dass das 2D-KBK-System keine stabilen lokalisierten Strukturen unterstützt. Die Autoren behaupten:

• Lösungen für lokalisierte glatte Anfangsdaten können über einen selbstähnlichen Mechanismus in endlicher Zeit einen Blow-up erleiden, was das Verhalten superkritischer NLS-Gleichungen widerspiegelt.
• Liniensolitärwellen sind stark instabil und anfällig für Blow-up unter transversalen Störungen.
• Es gibt keine stabilen Strukturen, die für dieses System in zwei Dimensionen lokalisiert sind.

Die Autoren stellen ausdrücklich fest, dass die genauen Kriterien dafür, wann Daten einen Blow-up erleiden, und die exakten Profile der Blow-up-Lösungen ($U^\infty$ und $V^\infty$) unbekannt bleiben. Die Arbeit dient als numerische Charakterisierung der Dynamik des Systems und hebt das Fehlen von Stabilität in 2D hervor, trotz der Integrabilität und Stabilitätseigenschaften seines 1D-Gegenstücks.