Single exponential $H^1$-upper bounds for the primitive equations

Résumé technique : Bornes supérieures exponentielles simples $H^1$ pour les équations primitives

Énoncé du problème
L'article étudie les équations primitives tridimensionnelles régissant les écoulements géophysiques à grande échelle (atmosphère et océan) dans un domaine de boîte périodique horizontale $\Omega = G \times (-h, 0)$. Le système inclut une viscosité complète et est soumis soit à des conditions aux limites de Neumann homogènes, soit de Dirichlet sur les limites supérieures et inférieures, avec des conditions périodiques sur les limites latérales. L'objectif principal est d'établir des bornes a priori pour les solutions fortes $v$ dans l'espace $L^\infty(0, \infty; H^1(\Omega)) \cap L^2(0, \infty; \dot{H}^2(\Omega))$.

Plus précisément, ce travail porte sur le taux de croissance de ces bornes par rapport aux données initiales $a$ et au temps $t$. La littérature existante, notamment les résultats de bien la posée globale par Cao et Titi [3], fournit des bornes supérieures qui croissent de manière double exponentielle par rapport aux données initiales et au temps. Des travaux ultérieurs [9, 10, 7] ont amélioré la dépendance temporelle ou traité des conditions aux limites spécifiques, mais ont généralement conservé une dépendance au moins double exponentielle vis-à-vis des données initiales ou du temps, s'appuyant souvent sur des inégalités de Sobolev logarithmiques ou des hypothèses de petite taille spécifiques.

Méthodologie
L'auteur emploie une série d'estimations d'énergie et d'inégalités d'interpolation pour dériver les bornes, en distinguant l'approche selon les conditions aux limites :

1. Conditions aux limites de Neumann (Proposition 1) :

   • La preuve utilise des estimations $L^q$ pour la vitesse horizontale $v$ (spécifiquement $q=4, 8$) et des estimations $L^r$ pour la dérivée verticale $v_z$ (spécifiquement $r=2, 4$).
   • La méthode repose sur les propriétés de l'annulation du terme de convection et sur l'inégalité de Calderón–Zygmund pour traiter le terme de pression.
   • Une difficulté clé dans ce cadre est l'absence de l'inégalité de Poincaré pour $v$ (puisque la moyenne de $v$ n'est pas nécessairement nulle). Par conséquent, la dérivation initiale produit une borne dépendante du temps ($e^{C\sqrt{t}}$).

2. Conditions aux limites de Dirichlet (Théorème 1) :

   • Pour obtenir une borne indépendante du temps, l'auteur déplace la stratégie d'estimation vers des estimations $L^6$ pour $v$ et des estimations $L^3$ pour $v_z$.
   • Un défi technique critique dans le cas de Dirichlet est le couplage entre le terme de pression $\pi$ et les valeurs aux limites. Le gradient de pression $\nabla_H \pi$ est décomposé en $\nabla_H \pi_1$ (dépendant de la non-linéarité de volume) et $\nabla_H \pi_2$ (dépendant de la trace de bord de $v_z$).
   • La preuve résout l'interdépendance entre l'estimation $L^q$ pour $v$ et l'estimation $L^r$ pour $v_z$ en choisissant soigneusement la paire $(q, r) = (6, 3)$. Cela permet de contrôler les termes de bord issus d'une estimation par les termes de dissipation de l'autre.

3. Affinement pour le cas de Neumann (Théorème 1) :

   • Pour éliminer la dépendance temporelle trouvée dans la Proposition 1, l'auteur introduit la variable de fluctuation $\tilde{v} = v - \bar{v}$ (où $\bar{v}$ est la moyenne verticale).
   • Puisque $\tilde{v}$ satisfait une inégalité de type Poincaré ($\|\tilde{v}\|_{L^q} \leq C \|v_z\|_{L^q}$), l'analyse se concentre sur les estimations $L^6$ pour $\tilde{v}$ et les estimations $L^3$ pour $v_z$.
   • Cette approche traite efficacement les « termes d'ordre inférieur » qui causaient auparavant la croissance temporelle, en exploitant le fait que la moyenne de la dérivée verticale est nulle.

Résultats clés
L'article établit le théorème principal suivant pour les conditions aux limites de Neumann et de Dirichlet :

Pour des données initiales $a \in H^1(\Omega)$ avec $a_z \in L^3(\Omega)$ (et satisfaisant des conditions de compatibilité), la solution forte $v$ satisfait :
$$\|v(t)\|_{H^1(\Omega)}^2 + \int_0^t \|\nabla^2 v(s)\|_{L^2(\Omega)}^2 \, ds \leq B(a) e^{C\|a\|_{L^2(\Omega)}^2}$$
pour tout $t \in [0, T]$, où $B(a)$ est une fonction polynomiale des normes des données initiales, et $C$ est une constante dépendant uniquement du domaine $\Omega$.

Crucialement, le facteur exponentiel dépend uniquement du carré de la norme $L^2$ des données initiales, $\|a\|_{L^2(\Omega)}^2$, et est indépendant du temps $t$.

Signification et revendications
L'article revendique deux contributions primaires :

1. Réduction de la croissance exponentielle : La borne supérieure dérivée est exponentielle simple par rapport aux données initiales ($e^{C\|a\|_{L^2}^2}$), contrastant nettement avec les bornes double exponentielles ($e^{e^{C\|a\|^2}}$) trouvées dans la littérature précédente (ex. [3, 9, 10, 7]). Cela représente une amélioration significative de la compréhension quantitative de la stabilité des équations primitives.
2. Estimations uniformes dans le temps pour le cas de Neumann : L'article fournit ce qu'il prétend être la première borne supérieure $H^1$ uniforme dans le temps pour le cas de la condition de Neumann. Dans ce cadre, l'inégalité de Poincaré standard n'est pas disponible pour le champ de vitesse $v$, rendant la dérivation de bornes indépendantes du temps non triviale et auparavant non réalisée dans la littérature.

L'auteur note que bien que le cas de Dirichlet permette une décroissance exponentielle lorsque $t \to \infty$ si l'inégalité de Poincaré est utilisée, ce comportement de décroissance spécifique n'est pas l'objet principal du présent travail. Les résultats sont présentés comme des estimations a priori pour des solutions fortes, supposant leur existence, laquelle est soutenue par des résultats de existence locale antérieurs.