Formation of Implosion Singularities in 3D Compressible Navier-Stokes-Korteweg Equation

Resumo Técnico: Formação de Singularidades de Implosão em Equações de Navier–Stokes–Korteweg de Compressíveis 3D

Enunciado do Problema
Este artigo investiga a formação de singularidades em tempo finito para o sistema de Navier–Stokes–Korteweg (NSK) compressível tridimensional com coeficientes de viscosidade e capilaridade dependentes da densidade. O sistema é definido por:
$$\begin{cases} \rho_t + \text{div}(\rho u) = 0, \\ (\rho u)_t + \text{div}(\rho u \otimes u) + \nabla P = \text{div}(2\mu(\rho)D(u)) + \nabla(\lambda(\rho)\text{div} u) + \text{div} K, \end{cases}$$
onde os coeficientes satisfazem as relações específicas de lei de potência:
$$\mu(\rho) = \nu\rho^\alpha, \quad \lambda(\rho) = 2\nu(\alpha - 1)\rho^\alpha, \quad \kappa(\rho) = \varepsilon^2\alpha^2\rho^{2\alpha-3}.$$
Trabalhos anteriores de Gu–Huang–Meng–Zhou [19] e Huang–Lei–Zhou [21] estabeleceram soluções fortes globais para dados iniciais arbitrariamente grandes (longe do vácuo) quando o expoente $\alpha$ reside em um intervalo adequado específico (por exemplo, $\alpha \in (\frac{9\sqrt{3}-4\sqrt{2}}{9\sqrt{3}-2\sqrt{2}}, 1)$ para 3D). A questão central abordada aqui é se as soluções existem globalmente ou sofrem explosão (blow-up) em tempo finito quando $\alpha$ cai fora deste intervalo, especificamente para expoentes positivos pequenos ($\alpha < 1/2$). Neste regime, a viscosidade volumétrica efetiva torna-se negativa ($\lambda(\rho) + \frac{2}{3}\mu(\rho) < 0$), uma característica associada a meios acústicos fora do equilíbrio e matéria ativa, o que desafia as teorias padrão de existência global.

Metodologia
Os autores empregam uma abordagem perturbativa baseada na construção de perfis de implosão auto-similares suaves derivados das equações de Euler compressíveis. A metodologia procede através de várias etapas principais:

1. Reformulação da Velocidade Efetiva: Para lidar com o tensor de terceira ordem de Korteweg, o sistema é reformulado usando uma velocidade efetiva $v = u + d\rho^{\alpha-2}\nabla\rho$ (onde $d = \nu - \sqrt{\nu^2-\varepsilon^2}$). Esta transformação recascula o sistema em uma estrutura parabólica-parabólica para $(\rho, v)$, revelando um termo de difusão degenerado oculto $d\Delta(\rho^\alpha)$ na equação de massa.
2. Escalonamento Auto-similar: O problema é transformado em coordenadas auto-similares $(\tau, y)$ onde $\tau = -\log(T-t)/\Lambda$ e $y = x/(T-t)^{1/\Lambda}$. Nestas coordenadas, os efeitos viscosos e capilares aparecem como perturbações exponencialmente decrescentes ($e^{-f_1\tau}$) desde que o parâmetro de escalonamento $\Lambda$ satisfaça uma condição específica ($f_1 > 0$).
3. Análise de Estabilidade de Perturbações: O núcleo da prova envolve a análise da estabilidade não linear do sistema em torno de um perfil de Euler auto-similar, esféricamente simétrico e suave $(\bar{S}, \bar{V})$. As equações de perturbação são derivadas, e um argumento de bootstrap é construído para controlar a solução.
4. Estimativas de Energia e Decomposição Espectral:
   • Estimativas de Energia Ponderadas: Estimativas de energia ponderadas de alta ordem são estabelecidas para controlar o decaimento das perturbações. Os autores utilizam uma função de peso específica $\phi$ e uma hierarquia de parâmetros para gerenciar a interação entre transporte e dissipação.
   • Argumento de Bootstrap: Um esquema de bootstrap rigoroso é desenhado para fechar as estimativas para derivadas de ordem inferior (especificamente o limite inferior da densidade) e derivadas de ordem superior. Uma inovação crítica é o uso de estruturas de transporte e de dissipação para estabelecer um limite inferior para a densidade, superando a dificuldade de que a equação da densidade possui uma estrutura de dissipação além da de transporte.
   • Controle de Modos Instáveis: O operador linearizado é decomposto em subespaços estáveis e instáveis. Os componentes instáveis, que correspondem a um espaço de dimensão finita, são controlados pela seleção dos coeficientes dos dados iniciais. Esta seleção garante que a solução permaneça em uma variedade de codimensão finita que leva à implosão.

Principais Contribuições e Resultados
O artigo prova a existência de dados iniciais suaves com densidade uniformemente separada do vácuo que evoluem para singularidades de implosão em tempo finito para o sistema NSK 3D quando $\alpha$ é suficientemente pequeno ($\alpha < 1/2$).

• Teorema Principal (Teorema 1.1): Para parâmetros satisfazendo ou:
  • (P1) $1 < \gamma < 1 + 2/\sqrt{3}$, $0 < \alpha < 1/2$, e $f_1 > 0$; ou
  • (P2) $\gamma \in (1, 1 + 2/\sqrt{3}) \setminus \mathcal{J}$ (onde $\mathcal{J}$ é um conjunto excepcional enumerável) e $0 < \alpha < \alpha^*(\gamma)$,
    existem dados iniciais $C^\infty$ $(\rho_0, u_0)$ tais que a solução correspondente $(\rho, u)$ existe em $[0, T)$ e desenvolve uma singularidade de implosão em $t=T$.
• Características da Singularidade: No tempo de explosão $T$:
  • A densidade torna-se infinita na origem: $\lim_{t \to T^-} \rho(t, 0) = +\infty$.
  • A velocidade efetiva $v$ é ilimitada em toda vizinhança da origem.
  • A solução converge, após o redimensionamento, para o perfil de implosão auto-similar prescrito.
• Corolário 1.1: O sistema original desenvolve uma singularidade de tempo finito, especificamente caracterizada pela explosão da densidade ($\lim_{t \to T^-} \rho(t, 0) = +\infty$). No entanto, a explosão da velocidade física $u$ não pode ser diretamente concluída a partir da explosão da velocidade efetiva, pois o termo de gradiente $\nabla\rho$ também pode tornar-se singular próximo à origem, e um possível cancelamento na relação $v = u + d\rho^{\alpha-2}\nabla\rho$ não pode ser descartado.

Significância
O artigo afirma complementar a teoria de existência global existente para o sistema NSK. Enquanto resultados anteriores demonstraram o bem-posto global para um intervalo específico de $\alpha$, este trabalho identifica um mecanismo distinto de explosão em tempo finito para o regime de "$\alpha$ pequeno". Este regime é caracterizado por uma estrutura de viscosidade volumétrica efetiva negativa, um fenômeno relevante para a física fora do equilíbrio. Os resultados demonstram que, mesmo com dados iniciais suaves longe do vácuo, a interação entre viscosidade dependente da densidade, capilaridade e viscosidade volumétrica negativa pode levar a singularidades de implosão onde os próprios campos de densidade e velocidade efetiva explodem, em vez de apenas suas derivadas (formação de choque). A construção baseia-se na estabilidade de perfis de Euler auto-similares, estendendo o mecanismo de implosão conhecido para as equações de Euler para o cenário viscoso e capilar de NSK.