Stability of purely convective steady-states of fractional Boussinesq equations in an exterior domain

Resumo Técnico: Estabilidade de Estados Estacionários Puramente Convectivos de Equações de Boussinesq Fracionárias em um Domínio Externo

Formulação do Problema
O artigo investiga a estabilidade global e o comportamento assintótico de fluxos de convecção térmica em um domínio ilimitado tridimensional $\Omega = \{x \in \mathbb{R}^3; |x| > 1\}$, exterior a uma esfera unitária. O fluido é aquecido a partir da fronteira da esfera, criando uma força de empuxo impulsionada pela gravidade. O estudo modifica o sistema clássico de Boussinesq substituindo a força viscosa padrão (determinada pelo operador de Stokes $-\Delta$) por uma potência fracionária do operador de Stokes, denotada como $A^\kappa$ (onde $0 < \kappa \leq 1$).

O sistema governante para as variáveis de perturbação $(u, \theta)$, representando desvios de velocidade e temperatura de um estado estacionário puramente condutivo $(u_0, T_0)$, é dado por:
$$\begin{cases} \partial_t u + A^\kappa u + P(\alpha g \theta) = -P(u \cdot \nabla u), \\ \partial_t \theta - \Delta \theta + u \cdot \nabla T_0 = -u \cdot \nabla \theta, \\ u|_{\partial \Omega} = 0, \quad \theta|_{\partial \Omega} = 0, \end{cases}$$
onde $P$ é a projeção de Helmholtz, $g$ é o vetor de aceleração gravitacional e $\alpha$ é uma constante relacionada à temperatura de fronteira. O objetivo principal é estabelecer a estabilidade global $L^2$ do estado estacionário e derivar taxas de decaimento algébrico para as perturbações quando $t \to \infty$.

Metodologia
A análise procede através de quatro etapas principais:

1. Teoria de Operadores e Geração de Semigrupos: Os autores definem um operador linear $L$ associado ao sistema linearizado. Usando o teorema de Kato–Rellich, provam que $L$ é autoadjunto e positivo sob a condição $0 < \alpha < 2^{-\kappa}$. Isso permite a construção de um semigrupo analítico fortemente contínuo $e^{-tL}$, que governa a evolução linear do sistema.
2. Existência de Soluções Fracas: Para lidar com os termos não lineares, os autores empregam uma técnica de molificação usando operadores $J_n$ e $I_n$. Eles constroem soluções aproximadas para a forma integral da equação do sistema. Ao estabelecer estimativas de energia uniformes e utilizar argumentos de compacidade (especificamente resultados do tipo lema de Aubin-Lions), provam a existência de uma solução fraca global no espaço $L^\infty(0, \tau; L^2) \cap L^2(0, \tau; D(A^{\kappa/2}) \times W^{1,2}_0)$.
3. Estimativas de Decaimento $L^2$: O decaimento da solução é analisado dividindo o termo integral não linear na formulação da solução suave. Os autores utilizam a representação espectral do operador $L$ e derivam taxas de decaimento específicas para o semigrupo $e^{-tL}$ atuando sobre os dados iniciais. Crucialmente, eles empregam desigualdades de interpolação (Lema 2.2 e Lema 2.3) e imersões de Sobolev para limitar os termos convectivos não lineares contra as taxas de decaimento lineares.
4. Iteração Bootstrap: Para estabelecer a taxa de decaimento algébrico ótima, o artigo utiliza um esquema de iteração bootstrap. Começando com uma taxa de decaimento base, os autores melhoram iterativamente a estimativa alimentando o limite de decaimento anterior de volta nas estimativas integrais não lineares, convergindo eventualmente para a taxa de decaimento alvo determinada pelo problema linearizado.

Principais Contribuições e Resultados
O artigo estabelece os seguintes resultados principais sob as condições $3/4 < \kappa \leq 1$, $0 < \beta < 1/2$ e $0 < \alpha < 2^{-\kappa}$:

• Existência Global: Para quaisquer dados iniciais $(a, b) \in L^2_\sigma(\Omega)^3 \times L^2(\Omega)$, o sistema de Boussinesq fracionário admite uma solução fraca global.
• Estabilidade Global: A solução $(u(t), \theta(t))$ converge para zero na norma $L^2$ quando $t \to \infty$, confirmando a estabilidade global do estado estacionário puramente condutivo.
• Decaimento Algébrico: Se o sistema linearizado exibe uma taxa de decaimento algébrico de $t^{-\beta}$, o sistema não linear herda essa taxa. Especificamente, a solução satisfaz:
  $$\|(u(t), \theta(t))\|_2 \leq C t^{-\beta}, \quad t > 0.$$
  A taxa de decaimento é restringida pela potência fracionária $\kappa$ e pelas propriedades dos dados iniciais.

Significado e Escopo
O artigo posiciona-se como a primeira tentativa de compreender a estabilidade de fluxos de Boussinesq exteriores fracionários. Ele estende trabalhos anteriores sobre as equações de Navier-Stokes clássicas (onde $\kappa=1$) e as equações de Boussinesq padrão para o caso fracionário em um domínio externo.

Os autores observam explicitamente que o resultado de estabilidade depende da suposição de que o parâmetro da força de empuxo $\alpha$ seja suficientemente pequeno. Eles reconhecem que, se a diferença de temperatura for grande, instabilidades (como as observadas na convecção de Rayleigh-Bénard) podem surgir, tornando a estabilidade inviável. O estudo não propõe novas aplicações ou configurações experimentais, mas fornece uma base matemática rigorosa para o comportamento assintótico desses modelos específicos de fluidos fracionários, preenchendo uma lacuna na literatura relativa a fluxos exteriores fracionários.