On decay of solutions to the anisotropic Boussinesq equations near the hydrostatic balance in half space $\mathbb{R}_+^3$

Sintesi Tecnica: Decadimento delle Soluzioni delle Equazioni di Boussinesq Anisotrope vicino all'Equilibrio Idrostatico in Semispazio

Enunciato del Problema
Questo articolo investiga il problema iniziale-valore di bordo per il sistema di Boussinesq anisotropo tridimensionale incomprimibile nel semispazio $\mathbb{R}^3_+$. Il sistema modella i fluidi geofisici e incorpora la dissipazione orizzontale (viscosità e diffusività termica che agiscono solo nelle direzioni orizzontali, $\Delta_h = \partial_{x_1}^2 + \partial_{x_2}^2$) pur mancando di dissipazione verticale. Lo studio si concentra sulla stabilità e sul comportamento a lungo termine delle perturbazioni vicino a uno stato di equilibrio idrostatico specifico definito da $u^{(0)} = (0,0,0)$, $\Theta^{(0)} = x_3$, e $P^{(0)} = \frac{1}{2}x_3^2$. Il sistema è soggetto a condizioni al contorno di scivolamento di Navier per la velocità ($u_3=0, \partial_3 u_1 = \partial_3 u_2 = 0$) e condizioni di Dirichlet per la temperatura ($\Theta=0$) al bordo $x_3=0$.

Metodologia
Gli autori utilizzano una combinazione di metodi energetici, argomenti di bootstrapping e rappresentazioni integrali tramite analisi di Fourier per affrontare la stabilità globale e i tassi di decadimento.

1. Stabilità Nonlineare (Metodi Energetici): Per dimostrare l'esistenza globale e la stabilità nello spazio Sobolev $H^3(\mathbb{R}^3_+)$, gli autori utilizzano un argomento di bootstrapping. Definiscono funzionali di energia che coinvolgono le norme $H^3$ e i termini di dissipazione orizzontale. Un componente critico di questa dimostrazione è la gestione dei termini al contorno derivanti dall'integrazione per parti. Nello specifico, gli autori stabiliscono che, sotto le date di condizioni al contorno fornite (in particolare $\partial_3^2 \theta_0 = 0$ in $x_3=0$), la derivata verticale di terzo ordine della pressione si annulla al contorno ($\partial_3^3 p = 0$), il che permette di chiudere le stime di energia senza perdita di regolarità. Si affidano a disuguaglianze del calcolo negli spazi Sobolev e a disuguaglianze di interpolazione anisotropa (simili a quelle nel Lemma 2.2) per limitare i termini non lineari.

2. Tassi di Decadimento (Analisi di Fourier e Rappresentazione Integrale): Per derivare i tassi di decadimento per il sistema non lineare, gli autori vanno oltre le dirette stime energetiche, che notano essere insufficienti per catturare un decadimento ottimale. Invece, essi:

   • Linearizzano il sistema attorno all'equilibrio idrostatico.
   • Applicano una trasformata di Fourier nelle variabili orizzontali ($x_h = (x_1, x_2)$) e trasformate di Fourier coseno/seno nella variabile verticale ($x_3$) alle equazioni linearizzate.
   • Derivano una formula esplicita per la soluzione del problema lineare che coinvolge sette funzioni kernel nello spazio delle frequenze. Questi kernel mostrano degenerazione e anisotropia.
   • Utilizzano il principio di Duhamel per esprimere la soluzione del sistema non lineare come una somma dell'evoluzione lineare e di termini integrali che coinvolgono le non linearità.
   • Stimano questi termini integrali combinando le proprietà di decadimento dei kernel lineari con i tassi di decadimento assunti della soluzione (bootstrapping del decadimento).

Contributi Chiave e Risultati

• Stabilità Globale (Teorema 1.1): L'articolo stabilisce l'esistenza globale e l'unicità delle soluzioni in $H^3(\mathbb{R}^3_+)$ per dati iniziali sufficientemente piccoli. Le soluzioni soddisfano un limite uniforme in $H^3$ ed esibiscono dissipazione nei gradienti orizzontali ($\nabla_h u, \nabla_h \theta$). Questo risultato conferma la stabilità dell'equilibrio idrostatico sotto la specificata dissipazione anisotropa e le condizioni al contorno.

• Tassi di Decadimento (Teorema 1.2): Gli autori derivano espliciti tassi di decadimento algebrico per la soluzione globale e le sue derivate. Assumendo che i dati iniziali possiedano una certa regolarità a bassa frequenza (misurata dall'operatore $\Lambda_h^{-\sigma}$ con $9/10 < \sigma < 1$), la soluzione decade come segue per tempo grande $t$:

  • Le norme $L^2$ di velocità e temperatura decadono come $(1+t)^{-\sigma/2 + \delta}$.
  • Le derivate verticali ($\partial_3 u, \partial_3 \theta$) decadono come $(1+t)^{-\sigma/2 + 3\delta}$.
  • I gradienti orizzontali ($\nabla_h u, \nabla_h \theta$) e la derivata verticale della componente della velocità verticale ($\partial_3 u_3$) decadono come $(1+t)^{-\sigma + 1/2 + \delta}$.
    Qui, $\delta$ è un piccolo parametro positivo vincolato dalla regolarità dei dati iniziali.

Significato e Rivendicazioni
L'articolo afferma che, mentre i tassi di decadimento per le equazioni linearizzate corrispondono a quelli trovati negli studi precedenti sul problema di Cauchy (specificamente riferendosi a Ji, Yan e Wu [10]), il comportamento a lungo termine per il problema iniziale-valore non lineare al contorno è distinto. Gli autori sottolineano che, in presenza di confini e dissipazione anisotropa, i tassi di decadimento dipendono fortemente dall'alta regolarità delle soluzioni ($H^3$) e dal particolare accoppiamento all'interno del sistema.

Il lavoro evidenzia le sfide poste dalla dissipazione anisotropa in domini limitati, dove i metodi energetici standard sono insufficienti per catturare i precisi tassi di decadimento. Costruendo una rappresentazione integrale con kernel degeneri e anisotropi, gli autori dimostrano che l'accoppiamento tra velocità e temperatura fornisce un ulteriore smoothing, influenzando la struttura del decadimento. I risultati estendono la comprensione dei sistemi di Boussinesq dai problemi di Cauchy in spazio intero ai problemi al valore iniziale-bordo in semispazio con dissipazione parziale, un contesto rilevante per la dinamica dei fluidi geofisici dove il mixing verticale è spesso trascurabile rispetto al mixing orizzontale.