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

Resumo Técnico: Estudo Numérico do Sistema Boussinesq Kaup-Broer-Kuperschmidt Bidimensional

Formulação do Problema
Este trabalho investiga a versão bem-posta do sistema bidimensional (2D) Kaup-Broer-Kuperschmidt (KBK), um caso específico dos sistemas Boussinesq 'abcd'. Enquanto o sistema KBK unidimensional (1D) é completamente integrável e possui soluções de onda solitária conhecidas com estabilidade orbital estabelecida, as propriedades qualitativas de sua generalização 2D permanecem amplamente inexploradas. O sistema 2D não é integrável, e resultados rigorosos sobre comportamento de longo prazo, existência global ou explosão em tempo finito estão ausentes. O objetivo principal é determinar numericamente se estruturas localizadas estáveis (como sólitons 2D ou sólitons de linha) existem no sistema KBK 2D e caracterizar a evolução de dados iniciais localizados gerais.

Metodologia
Os autores empregam uma abordagem numérica de alta precisão que combina métodos espectrais no espaço com esquemas de diferenciação exponencial no tempo (ETD).

• Formulação: O sistema é reformulado usando um potencial de velocidade $V$ (onde $v = \nabla V$) e transformado em um sistema de equações de Schrödinger não lineares (NLS) fracionárias no espaço de Fourier. Esta formulação facilita o uso de integradores exponenciais, especificamente o esquema ETD de quarta ordem de Cox e Matthews, necessário para lidar com a rigidez do sistema decorrente da dispersão de alta ordem.
• Soluções Estáticas: Para construir soluções estáticas localizadas, os autores primeiro resolvem a equação diferencial ordinária (EDO) radial associada em coordenadas polares usando uma discretização espectral de Chebyshev e um método Newton-Krylov com relaxação. Esta solução radial 1D é então interpolada em uma grade de Fourier 2D para servir como um chute inicial para resolver o sistema não linear 2D completo.
• Evolução Temporal: A evolução temporal é simulada em um domínio duplamente periódico (2-toro) para aproximar funções de decaimento rápido. Os autores monitoram a conservação relativa de energia ($\Delta E$) para garantir precisão numérica, visando erros significativamente abaixo de $10^{-3}$.
• Análise de Explosão: Para investigar potenciais singularidades, os autores utilizam uma técnica de redimensionamento dinâmico. Isso envolve redimensionar variáveis espaciais e o tempo para rastrear o comportamento das soluções perto de um tempo potencial de explosão $t^*$. O fator de escala $L(t)$ e os expoentes das normas $L^\infty$ e $L^2$ são ajustados a leis de potência da forma $(t^* - t)^\alpha$ para determinar a natureza da explosão (algébrica vs. exponencial).

Principais Contribuições e Resultados

1. Existência de uma Solução Estática: Os autores constroem com sucesso uma solução estática localizada e radialmente simétrica para o sistema KBK 2D. Esta solução corresponde a um estado fundamental de uma equação NLS vetorial.
2. Instabilidade das Soluções Estáticas: Perturbações numéricas desta solução estática revelam uma dicotomia no comportamento:
   • Dispersão: Perturbações com uma massa (ou norma $L^2$) ligeiramente menor que a da solução estática levam à dispersão da onda para o infinito.
   • Explosão em Tempo Finito: Perturbações com uma massa ligeiramente maior resultam em um processo rápido de auto-focalização levando à explosão em tempo finito. Os dados numéricos sugerem que a explosão é auto-similar, com o fator de escala $L$ comportando-se como $L \propto (t^* - t)^{1/2}$, e as normas crescendo como $\|\eta\| \sim (t^* - t)^\alpha$ com $\alpha \approx -1$. Este comportamento é análogo ao regime de focalização em equações NLS supercríticas em $L^2$.
3. Instabilidade Transversal de Sólitons de Linha: O estudo examina sólitons de linha 1D (soluções independentes da coordenada transversal $y$) estendidos trivialmente para 2D. Estas estruturas são encontradas como fortemente instáveis contra perturbações transversais. Tanto perturbações positivas quanto negativas levam à explosão em tempo finito, caracterizada pela formação de picos agudos e escala auto-similar semelhante ao caso da solução estática 2D. Notavelmente, a explosão ocorre mesmo quando a massa perturbada é menor que a da solução estática construída, um comportamento distinto da equação NLS crítica em $L^2$.
4. Dados Localizados Gerais: Simulações com dados iniciais gaussianos gerais (onde nenhuma solução estática específica é usada como base) mostram que o sistema não suporta estruturas localizadas estáveis. Dependendo dos parâmetros, os dados dispersam ou sofrem explosão.

Significado e Afirmações
O artigo conclui que o sistema KBK 2D não suporta estruturas localizadas estáveis. Os autores afirmam que:

• Soluções para dados iniciais suaves e localizados podem explodir em tempo finito através de um mecanismo auto-similar, espelhando o comportamento de equações NLS supercríticas.
• Ondas solitárias de linha são fortemente instáveis e propensas à explosão sob perturbações transversais.
• Não há estruturas estáveis localizadas em duas dimensões para este sistema.

Os autores afirmam explicitamente que os critérios precisos para que os dados explodam e os perfis exatos das soluções de explosão ($U^\infty$ e $V^\infty$) permanecem desconhecidos. O trabalho serve como uma caracterização numérica da dinâmica do sistema, destacando a ausência de estabilidade em 2D apesar das propriedades de integrabilidade e estabilidade de sua contraparte 1D.