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

Sintesi Tecnica: Studio Numerico del Sistema di Boussinesq Kaup-Broer-Kuperschmidt 2D

Enunciato del Problema
Questo lavoro indaga la versione ben posta del sistema bidimensionale (2D) Kaup-Broer-Kuperschmidt (KBK), un caso specifico dei sistemi di Boussinesq 'abcd'. Mentre il sistema KBK unidimensionale (1D) è completamente integrabile e possiede soluzioni d'onda solitaria note con stabilità orbitale stabilita, le proprietà qualitative della sua generalizzazione 2D rimangono in gran parte inesplorate. Il sistema 2D non è integrabile e mancano risultati rigorosi riguardanti il comportamento a lungo termine, l'esistenza globale o l'esplosione in tempo finito. L'obiettivo principale è determinare numericamente se esistono strutture localizzate stabili (come solitoni 2D o solitoni lineari) nel sistema KBK 2D e caratterizzare l'evoluzione di dati iniziali localizzati generali.

Metodologia
Gli autori impiegano un approccio numerico ad alta precisione che combina metodi spettrali nello spazio con schemi di differenziazione esponenziale nel tempo (ETD).

• Formulazione: Il sistema è riformulato utilizzando un potenziale di velocità $V$ (dove $v = \nabla V$) e trasformato in un sistema di equazioni di Schrödinger non lineari (NLS) frazionarie nello spazio di Fourier. Questa formulazione facilita l'uso di integratori esponenziali, in particolare lo schema ETD del quarto ordine di Cox e Matthews, necessario per gestire la rigidità del sistema derivante dalla dispersione di ordine superiore.
• Soluzioni Statiche: Per costruire soluzioni statiche localizzate, gli autori risolvono prima l'equazione differenziale ordinaria (ODE) radiale associata in coordinate polari utilizzando una discretizzazione spettrale di Chebyshev e un metodo Newton-Krylov con rilassamento. Questa soluzione radiale 1D viene quindi interpolata su una griglia di Fourier 2D per servire come ipotesi iniziale per risolvere il completo sistema non lineare 2D.
• Evoluzione Temporale: L'evoluzione temporale è simulata su un dominio doppiamente periodico (2-toro) per approssimare funzioni a decadimento rapido. Gli autori monitorano la conservazione relativa dell'energia ($\Delta E$) per garantire l'accuratezza numerica, mirando a errori significativamente inferiori a $10^{-3}$.
• Analisi dell'Esplosione: Per investigare potenziali singolarità, gli autori utilizzano una tecnica di riscalamento dinamico. Ciò comporta il riscalamento delle variabili spaziali e del tempo per tracciare il comportamento delle soluzioni vicino a un potenziale tempo di esplosione $t^*$. Il fattore di scala $L(t)$ e gli esponenti delle norme $L^\infty$ e $L^2$ sono adattati a leggi di potenza della forma $(t^* - t)^\alpha$ per determinare la natura dell'esplosione (algebrica vs esponenziale).

Contributi e Risultati Chiave

1. Esistenza di una Soluzione Statica: Gli autori costruiscono con successo una soluzione statica localizzata e radialmente simmetrica per il sistema KBK 2D. Questa soluzione corrisponde a uno stato fondamentale di un'equazione NLS vettoriale.
2. Instabilità delle Soluzioni Statiche: Le perturbazioni numeriche di questa soluzione statica rivelano una dicotomia nel comportamento:
   • Dispersione: Le perturbazioni con una massa (o norma $L^2$) leggermente inferiore a quella della soluzione statica portano alla dispersione dell'onda all'infinito.
   • Esplosione in Tempo Finito: Le perturbazioni con una massa leggermente superiore risultano in un rapido processo di auto-focalizzazione che porta all'esplosione in tempo finito. I dati numerici suggeriscono che l'esplosione è autosimile, con il fattore di scala $L$ che si comporta come $L \propto (t^* - t)^{1/2}$, e le norme che crescono come $\|\eta\| \sim (t^* - t)^\alpha$ con $\alpha \approx -1$. Questo comportamento è analogo al regime di focalizzazione nelle equazioni NLS $L^2$-supercritiche.
3. Instabilità Trasversale dei Solitoni Lineari: Lo studio esamina i solitoni lineari 1D (soluzioni indipendenti dalla coordinata trasversale $y$) estesi trivialmente a 2D. Queste strutture risultano fortemente instabili rispetto alle perturbazioni trasversali. Sia le perturbazioni positive che quelle negative portano all'esplosione in tempo finito, caratterizzata dalla formazione di picchi acuti e da un riscalamento autosimile simile al caso della soluzione statica 2D. Si nota che l'esplosione si verifica anche quando la massa perturbata è inferiore a quella della soluzione statica costruita, un comportamento distinto dall'equazione NLS $L^2$-critica.
4. Dati Localizzati Generali: Le simulazioni con dati iniziali gaussiani generali (dove non viene utilizzata alcuna soluzione statica specifica come base) mostrano che il sistema non supporta strutture localizzate stabili. A seconda dei parametri, i dati si disperdono o subiscono un'esplosione.

Significato e Affermazioni
Il documento conclude che il sistema KBK 2D non supporta strutture localizzate stabili. Gli autori affermano che:

• Le soluzioni per dati iniziali lisci e localizzati possono esplodere in tempo finito tramite un meccanismo autosimile, rispecchiando il comportamento delle equazioni NLS supercritiche.
• Le onde solitarie lineari sono fortemente instabili e soggette a esplosione sotto perturbazioni trasversali.
• Non esistono strutture stabili localizzate in due dimensioni per questo sistema.

Gli autori dichiarano esplicitamente che i criteri precisi affinché i dati esplodano e i profili esatti delle soluzioni di esplosione ($U^\infty$ e $V^\infty$) rimangono sconosciuti. Il lavoro funge da caratterizzazione numerica della dinamica del sistema, evidenziando l'assenza di stabilità in 2D nonostante le proprietà di integrabilità e stabilità del suo corrispettivo 1D.