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

Sintesi Tecnica: Stabilità di Stati Stazionari Puramente Convettivi delle Equazioni di Boussinesq Frazionarie in un Dominio Esterno

Formulazione del Problema
Il lavoro indaga la stabilità globale e il comportamento asintotico dei flussi di convezione termica in un dominio illimitato tridimensionale $\Omega = \{x \in \mathbb{R}^3; |x| > 1\}$, esterno a una sfera unitaria. Il fluido è riscaldato dal confine della sfera, generando una forza di galleggiamento guidata dalla gravità. Lo studio modifica il sistema classico di Boussinesq sostituendo la forza viscosa standard (determinata dall'operatore di Stokes $-\Delta$) con una potenza frazionaria dell'operatore di Stokes, indicata come $A^\kappa$ (dove $0 < \kappa \leq 1$).

Il sistema governante per le variabili di perturbazione $(u, \theta)$, che rappresentano le deviazioni di velocità e temperatura rispetto a uno stato stazionario puramente conduttivo $(u_0, T_0)$, è dato da:
$$\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}$$
dove $P$ è la proiezione di Helmholtz, $g$ è il vettore dell'accelerazione di gravità e $\alpha$ è una costante legata alla temperatura al confine. L'obiettivo primario è stabilire la stabilità globale $L^2$ dello stato stazionario e derivare tassi di decadimento algebrico per le perturbazioni al tendere di $t \to \infty$.

Metodologia
L'analisi procede attraverso quattro fasi principali:

1. Teoria degli Operatori e Generazione di Semigruppi: Gli autori definiscono un operatore lineare $L$ associato al sistema linearizzato. Utilizzando il teorema di Kato–Rellich, dimostrano che $L$ è autoaggiunto e positivo sotto la condizione $0 < \alpha < 2^{-\kappa}$. Ciò permette la costruzione di un semigruppo analitico fortemente continuo $e^{-tL}$, che governa l'evoluzione lineare del sistema.
2. Esistenza di Soluzioni Deboli: Per gestire i termini non lineari, gli autori impiegano una tecnica di mollificazione utilizzando gli operatori $J_n$ e $I_n$. Costruiscono soluzioni approssimate alla forma integrale dell'equazione del sistema. Stabilendo stime energetiche uniformi e utilizzando argomenti di compattezza (in particolare risultati di tipo lemma di Aubin-Lions), dimostrano l'esistenza di una soluzione debole globale nello spazio $L^\infty(0, \tau; L^2) \cap L^2(0, \tau; D(A^{\kappa/2}) \times W^{1,2}_0)$.
3. Stime di Decadimento $L^2$: Il decadimento della soluzione è analizzato suddividendo il termine integrale non lineare nella formulazione della soluzione debole. Gli autori utilizzano la rappresentazione spettrale dell'operatore $L$ e derivano specifici tassi di decadimento per il semigruppo $e^{-tL}$ che agisce sui dati iniziali. Crucialmente, impiegano disuguaglianze di interpolazione (Lemma 2.2 e Lemma 2.3) e immersioni di Sobolev per limitare i termini convettivi non lineari rispetto ai tassi di decadimento lineari.
4. Iterazione Bootstrap: Per stabilire il tasso di decadimento algebrico ottimale, il lavoro utilizza uno schema di iterazione bootstrap. Partendo da un tasso di decadimento di base, gli autori migliorano iterativamente la stima alimentando il precedente limite di decadimento indietro nelle stime integrali non lineari, convergendo infine al tasso di decadimento target determinato dal problema linearizzato.

Contributi e Risultati Chiave
Il lavoro stabilisce i seguenti risultati principali sotto le condizioni $3/4 < \kappa \leq 1$, $0 < \beta < 1/2$ e $0 < \alpha < 2^{-\kappa}$:

• Esistenza Globale: Per qualsiasi dato iniziale $(a, b) \in L^2_\sigma(\Omega)^3 \times L^2(\Omega)$, il sistema di Boussinesq frazionario ammette una soluzione debole globale.
• Stabilità Globale: La soluzione $(u(t), \theta(t))$ converge a zero nella norma $L^2$ al tendere di $t \to \infty$, confermando la stabilità globale dello stato stazionario puramente conduttivo.
• Decadimento Algebrico: Se il sistema linearizzato esibisce un tasso di decadimento algebrico di $t^{-\beta}$, il sistema non lineare eredita tale tasso. Nello specifico, la soluzione soddisfa:
  $$\|(u(t), \theta(t))\|_2 \leq C t^{-\beta}, \quad t > 0.$$
  Il tasso di decadimento è vincolato dalla potenza frazionaria $\kappa$ e dalle proprietà dei dati iniziali.

Significato e Ambito
Il lavoro si posiziona come il primo tentativo di comprendere la stabilità dei flussi di Boussinesq esterni frazionari. Estende lavori precedenti sulle equazioni di Navier-Stokes classiche (dove $\kappa=1$) e sulle equazioni di Boussinesq standard al caso frazionario in un dominio esterno.

Gli autori notano esplicitamente che il risultato di stabilità si basa sull'assunzione che il parametro della forza di galleggiamento $\alpha$ sia sufficientemente piccolo. Riconoscono che se la differenza di temperatura è grande, possono sorgere instabilità (come quelle osservate nella convezione di Rayleigh-Bénard), rendendo la stabilità non fattibile. Lo studio non propone nuove applicazioni o setup sperimentali, ma fornisce piuttosto una fondazione matematica rigorosa per il comportamento asintotico di questi specifici modelli fluidi frazionari, colmando una lacuna nella letteratura riguardante i flussi esterni frazionari.