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

Resumen Técnico: Estudio Numérico del Sistema de Boussinesq Kaup-Broer-Kuperschmidt Bidimensional

Planteamiento del Problema
Este trabajo investiga la versión bien planteada del sistema bidimensional (2D) Kaup-Broer-Kuperschmidt (KBK), un caso específico de los sistemas de Boussinesq 'abcd'. Mientras que el sistema KBK unidimensional (1D) es completamente integrable y posee soluciones de ondas solitarias conocidas con estabilidad orbital establecida, las propiedades cualitativas de su generalización 2D permanecen en gran medida inexploradas. El sistema 2D no es integrable, y carecen resultados rigurosos sobre el comportamiento a largo plazo, la existencia global o la explosión en tiempo finito. El objetivo principal es determinar numéricamente si existen estructuras localizadas estables (tales como solitones 2D o solitones de línea) en el sistema KBK 2D y caracterizar la evolución de datos iniciales localizados generales.

Metodología
Los autores emplean un enfoque numérico de alta precisión que combina métodos espectrales en el espacio con esquemas de diferenciación exponencial en el tiempo (ETD).

• Formulación: El sistema se reformula utilizando un potencial de velocidad $V$ (donde $v = \nabla V$) y se transforma en un sistema de ecuaciones de Schrödinger no lineales (NLS) fraccionarias en el espacio de Fourier. Esta formulación facilita el uso de integradores exponenciales, específicamente el esquema ETD de cuarto orden de Cox y Matthews, necesario para manejar la rigidez del sistema derivada de la dispersión de alto orden.
• Soluciones Estáticas: Para construir soluciones estáticas localizadas, los autores resuelven primero la ecuación diferencial ordinaria (ODE) radial asociada en coordenadas polares utilizando una discretización espectral de Chebyshev y un método Newton-Krylov con relajación. Esta solución radial 1D se interpola luego en una cuadrícula de Fourier 2D para servir como una estimación inicial para resolver el sistema no lineal completo 2D.
• Evolución Temporal: La evolución temporal se simula en un dominio doblemente periódico (2-toro) para aproximar funciones de decaimiento rápido. Los autores monitorean la conservación relativa de la energía ($\Delta E$) para garantizar la precisión numérica, apuntando a errores significativamente por debajo de $10^{-3}$.
• Análisis de Explosión: Para investigar posibles singularidades, los autores utilizan una técnica de reescalado dinámico. Esto implica reescalar las variables espaciales y el tiempo para rastrear el comportamiento de las soluciones cerca de un tiempo potencial de explosión $t^*$. El factor de escala $L(t)$ y los exponentes de las normas $L^\infty$ y $L^2$ se ajustan a leyes de potencia de la forma $(t^* - t)^\alpha$ para determinar la naturaleza de la explosión (algebraica vs. exponencial).

Contribuciones y Resultados Clave

1. Existencia de una Solución Estática: Los autores construyen exitosamente una solución estática localizada y simétrica radialmente para el sistema KBK 2D. Esta solución corresponde a un estado fundamental de una ecuación NLS vectorial.
2. Inestabilidad de las Soluciones Estáticas: Las perturbaciones numéricas de esta solución estática revelan una dicotomía en el comportamiento:
   • Dispersión: Las perturbaciones con una masa (o norma $L^2$) ligeramente menor que la de la solución estática conducen a la dispersión de la onda hacia el infinito.
   • Explosión en Tiempo Finito: Las perturbaciones con una masa ligeramente mayor resultan en un proceso rápido de autoenfoque que conduce a una explosión en tiempo finito. Los datos numéricos sugieren que la explosión es autosimilar, con el factor de escala $L$ comportándose como $L \propto (t^* - t)^{1/2}$, y las normas creciendo como $\|\eta\| \sim (t^* - t)^\alpha$ con $\alpha \approx -1$. Este comportamiento es análogo al régimen de enfoque en ecuaciones NLS supercríticas en $L^2$.
3. Inestabilidad Transversal de Solitones de Línea: El estudio examina solitones de línea 1D (soluciones independientes de la coordenada transversal $y$) extendidos trivialmente a 2D. Se encuentra que estas estructuras son fuertemente inestables frente a perturbaciones transversales. Tanto las perturbaciones positivas como las negativas conducen a una explosión en tiempo finito, caracterizada por la formación de picos agudos y un escalado autosimilar similar al caso de la solución estática 2D. Cabe destacar que la explosión ocurre incluso cuando la masa perturbada es menor que la de la solución estática construida, un comportamiento distinto al de la ecuación NLS crítica en $L^2$.
4. Datos Localizados Generales: Las simulaciones con datos iniciales gaussianos generales (donde no se utiliza ninguna solución estática específica como base) muestran que el sistema no soporta estructuras localizadas estables. Dependiendo de los parámetros, los datos se dispersan o sufren una explosión.

Significado y Afirmaciones
El artículo concluye que el sistema KBK 2D no soporta estructuras localizadas estables. Los autores afirman que:

• Las soluciones para datos iniciales suaves localizados pueden explotar en tiempo finito mediante un mecanismo autosimilar, reflejando el comportamiento de las ecuaciones NLS supercríticas.
• Las ondas solitarias de línea son fuertemente inestables y propensas a la explosión bajo perturbaciones transversales.
• No existen estructuras estables localizadas en dos dimensiones para este sistema.

Los autores declaran explícitamente que los criterios precisos para que los datos exploten y los perfiles exactos de las soluciones de explosión ($U^\infty$ y $V^\infty$) permanecen desconocidos. El trabajo sirve como una caracterización numérica de la dinámica del sistema, destacando la ausencia de estabilidad en 2D a pesar de las propiedades de integrabilidad y estabilidad de su contraparte 1D.