Global existence of classical solutions for the multi-dimensional compressible Navier-Stokes-Poisson equations on solid balls for arbitrary spherically symmetric large initial data

Título: Existência Global de Soluções Clássicas para as Equações Compressíveis de Navier-Stokes-Poisson Multi-dimensionais em Esferas Sólidas com Dados Iniciais Grandes Arbitrários

1. O Problema Investigado

O artigo aborda um problema aberto de longa data na análise matemática de fluidos: a existência global de soluções clássicas para as equações de Navier-Stokes-Poisson (NSP) compressíveis em três dimensões (e também em duas), considerando dados iniciais grandes (arbitrariamente grandes) e não-vácuo.

O sistema modela o movimento de estrelas gasosas ou plasmas em um campo gravitacional ou elétrico auto-consistente. As equações governam a densidade ($\rho$), a velocidade ($u$) e o potencial ($\Phi$). O desafio principal reside na combinação de:

• Não-linearidade forte das equações de Navier-Stokes.
• O acoplamento com a equação de Poisson (termo de fonte $\rho \nabla \Phi$).
• A presença de singularidades no centro da esfera (origem $r=0$) devido à simetria esférica e aos termos de viscosidade dependentes da densidade.
• A necessidade de lidar com dados iniciais grandes, onde métodos de perturbação de pequenos dados não são aplicáveis.

O foco é em domínios limitados (esferas sólidas) com condições de contorno de deslizamento e viscosidade dependente da densidade do tipo BD (Bresch-Desjardins), onde os coeficientes de viscosidade cisalhante ($\mu$) e volumétrica ($\lambda$) satisfazem $\mu(\rho) = \rho^\alpha$ e $\lambda(\rho) = (\alpha-1)\rho^\alpha$.

2. Metodologia e Estratégias Chave

Os autores desenvolvem uma análise de estimativas a priori robustas para estender soluções locais para soluções globais. A metodologia baseia-se em:

• Coordenadas Lagrangianas: O sistema é transformado para coordenadas lagrangianas para lidar com o domínio variável e simplificar a estrutura das equações, fixando o domínio espacial em $[0, 1]$.
• Velocidade Efetiva ($w$): Introdução da variável $w = u + r^{N-1}(\rho^\alpha)_x$ (ou forma equivalente). Esta variável é crucial para explorar a estrutura dissipativa do sistema e obter estimativas de energia adicionais.
• Estimativas de Entropia BD: Utilização da desigualdade de entropia de Bresch-Desjardins para obter controle sobre o gradiente da densidade, especificamente termos da forma $\nabla \rho^{\alpha-1/2}$.
• Tratamento da Singularidade Central:
  • O maior obstáculo é o termo de Poisson $\frac{x}{r^k}$ (onde $k=1$ em 2D e $k=2$ em 3D), que gera singularidades fortes na origem.
  • Os autores estabelecem estimativas ponderadas para a densidade perto da origem (ex: $\|\rho^{\alpha-1/2} r^\xi\|_{L^\infty}$).
  • Eles derivam relações explícitas entre a coordenada lagrangiana $x$ e a coordenada radial $r$ (ex: $x \leq C r^\beta$), permitindo controlar integrais singulares como $\int \frac{x^2}{r^2} dx$ através de normas de energia da velocidade.
• Estimativas de Alta Integrabilidade:
  • Para fechar as estimativas, os autores demonstram que a velocidade $u$ e a velocidade efetiva $w$ possuem alta integrabilidade ($L^{2k}$ ou $L^{2m}$).
  • Isso é alcançado escolhendo parâmetros específicos de $\alpha$ e $\gamma$ (expoente adiabático) e utilizando o conjunto denso $A_{set}$ para garantir a existência de expoentes adequados.
  • Em 3D, o tratamento é mais delicado devido à singularidade mais forte ($\frac{x}{r^2}$), exigindo condições mais restritivas em $\gamma$ (ex: $\gamma > 4/3$ para estrelas gasosas).
• Limites Superiores e Inferiores para a Densidade:
  • Limite Superior: Obtido via desigualdade de Sobolev unidimensional aplicada a potências da densidade, utilizando as estimativas de alta integrabilidade de $u$ e $w$.
  • Limite Inferior (Ausência de Vácuo): Obtido analisando a evolução do volume específico $v = 1/\rho$ e utilizando uma desigualdade de Gronwall acoplada para as normas $L^\infty$ de $u$ e $w$.

3. Principais Resultados

O artigo estabelece a existência global de soluções clássicas sfericamente simétricas para dados iniciais grandes, sem vácuo, sob as seguintes condições:

• Caso Bidimensional ($N=2$):

  • Para $\alpha \in (1/2, 1]$, existe solução global para qualquer $\gamma > 2 - \alpha$ (se $\alpha < 1$) ou $\gamma > 1$ (se $\alpha = 1$).
  • A densidade permanece estritamente positiva e limitada superiormente: $C(T)^{-1} \leq \rho(x,t) \leq C(T)$.

• Caso Tridimensional ($N=3$):

  • Para $\alpha \in (5/6, 1]$, existe solução global se:
    • $\gamma > 4/3$ (necessário para estrelas gasosas, $\kappa=1$, devido à singularidade do Poisson).
    • Para plasmas ($\kappa=-1$), a restrição é menos severa, mas ainda requer $\gamma > 1$.
  • As condições em $\alpha$ e $\gamma$ são otimizadas através de funções auxiliares ($n_2(\alpha)$ e $n_3(\alpha)$) definidas no artigo.

• Regularidade: As soluções obtidas pertencem a espaços de Sobolev de alta regularidade ($H^3$ para densidade, $H^4$ para velocidade em tempos positivos), garantindo que são soluções clássicas.

4. Contribuições e Significância

• Primeiro Resultado em Domínio Multi-dimensional Real: Este é considerado o primeiro resultado que prova a existência global de soluções clássicas para dados iniciais arbitrariamente grandes em um domínio verdadeiramente multi-dimensional (2D e 3D) com características de alta dimensão, superando limitações de modelos 1D ou simetrias planas.
• Resolução da Singularidade Central: O trabalho fornece uma técnica inovadora para controlar a singularidade no centro da esfera ($r=0$) através de estimativas ponderadas de densidade derivadas da entropia BD, algo que métodos anteriores não conseguiam fazer de forma completa para dados grandes.
• Generalidade dos Dados: Ao lidar com dados iniciais grandes (não perturbativos), o resultado elimina a necessidade de assumir que o sistema está próximo de um estado estacionário, tornando a teoria aplicável a cenários físicos mais realistas e extremos.
• Distinção entre Estrelas e Plasmas: O artigo trata separadamente os casos de estrelas gasosas ($\kappa > 0$) e plasmas ($\kappa < 0$), identificando como a natureza do acoplamento afeta as condições necessárias para a existência global (especialmente a restrição $\gamma > 4/3$ em 3D para estrelas).

Em resumo, o artigo representa um avanço significativo na teoria de equações diferenciais parciais não lineares, provando que, sob simetria esférica e com viscosidade do tipo BD, o sistema de Navier-Stokes-Poisson é bem-comportado globalmente mesmo para configurações de energia e densidade iniciais muito elevadas.