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

Résumé technique : Étude numérique du système de Boussinesq Kaup-Broer-Kuperschmidt en deux dimensions

Énoncé du problème
Ce travail examine la version bien posée du système bidimensionnel (2D) Kaup-Broer-Kuperschmidt (KBK), un cas particulier des systèmes de Boussinesq « abcd ». Alors que le système KBK unidimensionnel (1D) est complètement intégrable et possède des solutions d'ondes solitaires connues avec une stabilité orbitale établie, les propriétés qualitatives de sa généralisation en 2D restent largement inexplorées. Le système 2D n'est pas intégrable, et il n'existe pas de résultats rigoureux concernant le comportement à long terme, l'existence globale ou l'explosion en temps fini. L'objectif principal est de déterminer numériquement l'existence de structures localisées stables (telles que des solitons 2D ou des solitons de ligne) dans le système KBK 2D et de caractériser l'évolution de données initiales localisées générales.

Méthodologie
Les auteurs emploient une approche numérique de haute précision combinant des méthodes spectrales dans l'espace avec des schémas de différences exponentielles (ETD) dans le temps.

• Formulation : Le système est reformulé en utilisant un potentiel de vitesse $V$ (où $v = \nabla V$) et transformé en un système d'équations de Schrödinger non linéaires (NLS) fractionnaires dans l'espace de Fourier. Cette formulation facilite l'utilisation d'intégrateurs exponentiels, spécifiquement le schéma ETD d'ordre quatre de Cox et Matthews, nécessaire pour gérer la rigidité du système résultant de la dispersion d'ordre élevé.
• Solutions statiques : Pour construire des solutions statiques localisées, les auteurs résolvent d'abord l'équation différentielle ordinaire (ODE) radiale associée en coordonnées polaires en utilisant une discrétisation spectrale de Chebyshev et une méthode de Newton-Krylov avec relaxation. Cette solution radiale 1D est ensuite interpolée sur une grille de Fourier 2D pour servir d'initialisation pour la résolution du système non linéaire complet en 2D.
• Évolution temporelle : L'évolution temporelle est simulée sur un domaine doublement périodique (2-tore) pour approximer des fonctions à décroissance rapide. Les auteurs surveillent la conservation relative de l'énergie ($\Delta E$) pour garantir la précision numérique, visant des erreurs nettement inférieures à $10^{-3}$.
• Analyse d'explosion : Pour investiguer les singularités potentielles, les auteurs utilisent une technique de redimensionnement dynamique. Cela implique de redimensionner les variables spatiales et le temps pour suivre le comportement des solutions près d'un temps d'explosion potentiel $t^*$. Le facteur d'échelle $L(t)$ et les exposants des normes $L^\infty$ et $L^2$ sont ajustés à des lois de puissance de la forme $(t^* - t)^\alpha$ pour déterminer la nature de l'explosion (algébrique vs exponentielle).

Contributions et résultats clés

1. Existence d'une solution statique : Les auteurs construisent avec succès une solution statique localisée à symétrie radiale pour le système KBK 2D. Cette solution correspond à un état fondamental d'une équation NLS vectorielle.
2. Instabilité des solutions statiques : Les perturbations numériques de cette solution statique révèlent une dichotomie dans le comportement :
   • Dispersion : Les perturbations avec une masse (ou norme $L^2$) légèrement inférieure à celle de la solution statique conduisent à la dispersion de l'onde à l'infini.
   • Explosion en temps fini : Les perturbations avec une masse légèrement supérieure entraînent un processus rapide d'auto-focalisation menant à une explosion en temps fini. Les données numériques suggèrent que l'explosion est auto-similaire, avec le facteur d'échelle $L$ se comportant comme $L \propto (t^* - t)^{1/2}$, et les normes croissant comme $\|\eta\| \sim (t^* - t)^\alpha$ avec $\alpha \approx -1$. Ce comportement est analogue au régime de focalisation dans les équations NLS surcritiques en $L^2$.
3. Instabilité transversale des solitons de ligne : L'étude examine les solitons de ligne 1D (solutions indépendantes de la coordonnée transverse $y$) étendus trivialement en 2D. Ces structures se révèlent fortement instables face aux perturbations transverses. Les perturbations positives et négatives conduisent toutes deux à une explosion en temps fini, caractérisée par la formation de pics aigus et un scaling auto-similaire similaire au cas de la solution statique 2D. Fait notable, l'explosion se produit même lorsque la masse perturbée est inférieure à celle de la solution statique construite, un comportement distinct de l'équation NLS critique en $L^2$.
4. Données localisées générales : Les simulations avec des données initiales gaussiennes générales (où aucune solution statique spécifique n'est utilisée comme base) montrent que le système ne supporte pas de structures localisées stables. Selon les paramètres, les données se dispersent ou subissent une explosion.

Signification et affirmations
L'article conclut que le système KBK 2D ne supporte pas de structures localisées stables. Les auteurs affirment que :

• Les solutions pour des données initiales lisses et localisées peuvent exploser en temps fini via un mécanisme auto-similaire, reflétant le comportement des équations NLS surcritiques.
• Les ondes solitaires de ligne sont fortement instables et sujettes à l'explosion sous l'effet de perturbations transverses.
• Il n'existe aucune structure stable localisée en deux dimensions pour ce système.

Les auteurs déclarent explicitement que les critères précis pour qu'une donnée explose et les profils exacts des solutions d'explosion ($U^\infty$ et $V^\infty$) restent inconnus. Ce travail sert de caractérisation numérique de la dynamique du système, mettant en évidence l'absence de stabilité en 2D malgré l'intégrabilité et les propriétés de stabilité de son homologue 1D.