Global existence and optimal decay for a three-dimensional penalized Navier--Stokes system with biharmonic damping

Resumo Técnico: Existência Global e Decaimento Ótimo para um Sistema de Navier–Stokes Penalizado Tridimensional com Amortecimento Biharmônico

Enunciado do Problema
Este artigo investiga um sistema parabólico tridimensional que serve como uma aproximação hiperviscosa e penalizada das equações de Navier–Stokes incompressíveis. O modelo aborda os desafios matemáticos das equações de Navier–Stokes clássicas (1.1), particularmente o problema em aberto da regularidade global e unicidade para dados iniciais arbitrários em três dimensões. O sistema proposto (1.3) combina três mecanismos:

1. Difusão viscosa clássica ($-\nu\Delta u$).
2. Amortecimento biharmônico (hiperviscoso) ($\beta\Delta^2 u$), fornecendo dissipação de ordem superior.
3. Penalização da divergência ($-\frac{1}{\varepsilon}\nabla \text{div} u$), que relaxa a restrição estrita de incompressibilidade ($\text{div} u = 0$) ao impor quase-incompressibilidade conforme o parâmetro $\varepsilon \to 0$.

Crucialmente, o termo de convecção não linear é modificado usando uma correção de tipo Temam de simetria anti-simétrica ($N(u) = (u \cdot \nabla)u + \frac{1}{2}u \text{div} u$). Esta correção compensa a compressibilidade fraca introduzida pela penalização, restaurando a propriedade de cancelamento $L^2$ do termo de transporte, que é tipicamente perdida em formulações compressíveis.

Metodologia
A análise baseia-se em uma combinação de métodos de energia, teoria de semigrupos e técnicas de divisão de Fourier, com foco central na derivação de estimativas que são uniformes em relação ao parâmetro de penalização $\varepsilon > 0$.

• Soluções Fracas (Teorema 1.3): A existência de soluções fracas globais para quaisquer dados iniciais $L^2$ é estabelecida via uma aproximação de tipo Galerkin de Friedrichs. Os autores projetam a equação sobre bolas de Fourier usando operadores de corte $J_n$ em vez do projetor de Leray. Esta distinção é crítica; usar o projetor de Leray anularia o termo de penalização, alterando o problema limite. Estimativas a priori uniformes são derivadas da identidade de energia, e a compacidade é obtida via o lema de Aubin–Lions–Simon.
• Soluções Fortes (Teorema 1.4): Para dados iniciais pequenos em $H^2(\mathbb{R}^3)$, soluções fortes globais são provadas usando um argumento de ponto fixo na formulação tênue associada ao semigrupo linear $S_\varepsilon(t) = e^{tL_\varepsilon}$. O operador linear $L_\varepsilon = \nu\Delta - \beta\Delta^2 + \frac{1}{\varepsilon}\nabla\text{div}$ gera um semigrupo com propriedades de suavização. A prova baseia-se em uma estimativa de energia $H^2$ uniforme, onde os termos não lineares são absorvidos pela dissipação viscosa e biharmônica.
• Estimativas de Decaimento (Teorema 1.6): Taxas de decaimento de longo tempo ótimas são derivadas usando o método de divisão de Fourier (Schonbek). A análise separa as frequências baixas e altas. O comportamento de baixa frequência é controlado via a fórmula de Duhamel e a decomposição específica da não linearidade, enquanto as frequências altas são geridas através de desigualdades de energia diferenciadas e interpolação de Gagliardo–Nirenberg.

Principais Contribuições e Resultados

1. Existência Global Fraca: O artigo prova a existência global de soluções fracas para quaisquer dados iniciais $u_0 \in L^2(\mathbb{R}^3)$. As soluções satisfazem classes de regularidade específicas, incluindo $L^\infty(0, T; L^2)$, $L^2(0, T; H^2)$ e $L^4(0, T; L^4)$.
2. Soluções Fortes Globais com Pequenez Uniforme: Para dados iniciais suficientemente pequenos em $H^2(\mathbb{R}^3)$, o sistema admite uma solução forte global única. Uma contribuição significativa é que o limiar de pequena magnitude $\delta_0$ e as estimativas de energia associadas são independentes de $\varepsilon$. Esta uniformidade garante que os resultados se apliquem a toda a família de sistemas penalizados, não apenas a uma instância fixa.
3. Taxas de Decaimento Ótimas: Para dados iniciais em $L^1(\mathbb{R}^3) \cap H^2(\mathbb{R}^3)$, os autores derivam estimativas de decaimento ótimas:
   $$\|u(t)\|_{L^2} \leq C(1+t)^{-3/4}$$
   $$\|\nabla^k u(t)\|_{L^2} \leq C(1+t)^{-3/4 - k/2}, \quad k=1, 2$$
   Estas taxas coincidem com o decaimento assintótico da equação do calor clássica em três dimensões. A análise confirma que, embora o termo biharmônico domine a dissipação de alta frequência, o decaimento de longo tempo é governado pelo comportamento de fluxo de calor das baixas frequências.
4. Uniformidade em $\varepsilon$: Uma característica definidora da análise é que todas as constantes e estimativas permanecem uniformes conforme $\varepsilon \to 0$. O termo de penalização $\frac{1}{\varepsilon}\|\text{div} u\|^2$ é tratado como um termo de dissipação não negativo que é ou retido ou descartado, nunca limitado superiormente por uma quantidade envolvendo $1/\varepsilon$.

Significância e Alegações
O artigo afirma fornecer uma "fundação analítica estável e rigorosa" para a aproximação penalizada de fluxos incompressíveis. Ao incorporar o amortecimento biharmônico e a correção de Temam, o modelo alcança resultados de bem-postura melhorados em comparação com as equações de Navier–Stokes clássicas, especificamente no que diz respeito à disponibilidade de estimativas de energia $L^2$ incondicionais e suavização de ordem superior.

Os autores observam que suas estimativas uniformes sugerem que, conforme $\varepsilon \to 0$, as soluções fortes convergem para um limite que resolve um sistema de Navier–Stokes hiperviscoso incompressível. No entanto, eles declaram explicitamente que identificar o sistema limite e provar as taxas de convergência em $\varepsilon$ são direções deixadas para trabalhos futuros. Além disso, embora a regularidade global para dados arbitrários seja conhecida para o sistema de Navier–Stokes hiperviscoso incompressível (devido a Lions), os autores não alegam tal resultado para o sistema penalizado com dados grandes, observando que a estrutura compressível anti-simétrica modifica as estimativas de energia de ordem superior de uma forma que impede uma extensão direta do resultado de dados grandes.

O trabalho distingue-se de estudos anteriores sobre equações de Navier–Stokes amortecidas (por exemplo, com termos de absorção) ao focar em um mecanismo de amortecimento de ordem superior e um método de penalização que elimina o termo de pressão, simplificando assim a formulação fraca ao remover a restrição de divergência nula sobre as funções de teste.