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

Résumé technique : Existence globale et décroissance optimale pour un système de Navier–Stokes pénalisé tridimensionnel avec amortissement biharmonique

Énoncé du problème
Cet article étudie un système parabolique tridimensionnel servant d'approximation hypervisqueuse et pénalisée des équations de Navier–Stokes incompressibles. Le modèle aborde les défis mathématiques des équations de Navier–Stokes classiques (1.1), en particulier le problème ouvert de la régularité globale et de l'unicité pour des données initiales arbitraires en trois dimensions. Le système proposé (1.3) combine trois mécanismes :

1. Diffusion visqueuse classique ($-\nu\Delta u$).
2. Amortissement biharmonique (hypervisqueux) ($\beta\Delta^2 u$), fournissant une dissipation d'ordre supérieur.
3. Pénalisation de la divergence ($-\frac{1}{\varepsilon}\nabla \text{div} u$), qui relâche la contrainte stricte d'incompressibilité ($\text{div} u = 0$) en imposant une quasi-incompressibilité lorsque le paramètre $\varepsilon \to 0$.

Crucialement, le terme de convection non linéaire est modifié par une correction antisymétrique de type Temam ($N(u) = (u \cdot \nabla)u + \frac{1}{2}u \text{div} u$). Cette correction compense la faible compressibilité introduite par la pénalisation, restaurant la propriété de l'annulation en $L^2$ du terme de transport, qui est typiquement perdue dans les formulations compressibles.

Méthodologie
L'analyse repose sur une combinaison de méthodes d'énergie, de la théorie des semi-groupes et des techniques de division de Fourier (Fourier splitting), avec un accent central mis sur la dérivation d'estimations qui sont uniformes par rapport au paramètre de pénalisation $\varepsilon > 0$.

• Solutions faibles (Théorème 1.3) : L'existence de solutions faibles globales pour des données initiales $L^2$ arbitraires est établie via une approximation de type Galerkin de Friedrichs. Les auteurs projettent l'équation sur des boules de Fourier en utilisant des opérateurs de coupure $J_n$ plutôt que le projecteur de Leray. Cette distinction est critique ; l'utilisation du projecteur de Leray annulerait le terme de pénalisation, modifiant le problème limite. Des estimations a priori uniformes sont dérivées de l'identité énergétique, et la compacité est obtenue via le lemme d'Aubin–Lions–Simon.
• Solutions fortes (Théorème 1.4) : Pour de petites données initiales dans $H^2(\mathbb{R}^3)$, des solutions fortes globales sont prouvées en utilisant un argument de point fixe sur la formulation douce associée au semi-groupe linéaire $S_\varepsilon(t) = e^{tL_\varepsilon}$. L'opérateur linéaire $L_\varepsilon = \nu\Delta - \beta\Delta^2 + \frac{1}{\varepsilon}\nabla\text{div}$ génère un semi-groupe possédant des propriétés de lissage. La preuve repose sur une estimation d'énergie $H^2$ uniforme où les termes non linéaires sont absorbés par la dissipation visqueuse et biharmonique.
• Estimations de décroissance (Théorème 1.6) : Des taux de décroissance optimale à grand temps sont dérivés en utilisant la méthode de division de Fourier (Schonbek). L'analyse sépare les basses et hautes fréquences. Le comportement aux basses fréquences est contrôlé via la formule de Duhamel et la décomposition spécifique de la non-linéarité, tandis que les hautes fréquences sont gérées par des inégalités d'énergie différenciées et l'interpolation de Gagliardo–Nirenberg.

Contributions clés et résultats

1. Existence globale faible : L'article prouve l'existence globale de solutions faibles pour toute donnée initiale $u_0 \in L^2(\mathbb{R}^3)$. Les solutions satisfont des classes de régularité spécifiques, incluant $L^\infty(0, T; L^2)$, $L^2(0, T; H^2)$, et $L^4(0, T; L^4)$.
2. Solutions fortes globales avec petiteur de données uniforme : Pour des données initiales suffisamment petites dans $H^2(\mathbb{R}^3)$, le système admet une solution forte globale unique. Une contribution significative est que le seuil de petiteur $\delta_0$ et les bornes d'énergie associées sont indépendants de $\varepsilon$. Cette uniformité garantit que les résultats s'appliquent à toute la famille de systèmes pénalisés, et pas seulement à une instance fixée.
3. Taux de décroissance optimaux : Pour des données initiales dans $L^1(\mathbb{R}^3) \cap H^2(\mathbb{R}^3)$, les auteurs dérivent des estimations de décroissance optimales :
   $$\|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$$
   Ces taux correspondent à la décroissance asymptotique de l'équation de la chaleur classique en trois dimensions. L'analyse confirme que si le terme biharmonique domine la dissipation des hautes fréquences, la décroissance à grand temps est gouvernée par le comportement de l'écoulement de la chaleur des basses fréquences.
4. Uniformité en $\varepsilon$ : Une caractéristique définissante de l'analyse est que toutes les constantes et estimations restent uniformes lorsque $\varepsilon \to 0$. Le terme de pénalisation $\frac{1}{\varepsilon}\|\text{div} u\|^2$ est traité comme un terme de dissipation non négatif qui est soit conservé, soit écarté, mais jamais borné de manière supérieure par une quantité impliquant $1/\varepsilon$.

Signification et affirmations
L'article affirme fournir un « fondement analytique stable et rigoureux » pour l'approximation pénalisée des écoulements incompressibles. En incorporant l'amortissement biharmonique et la correction de Temam, le modèle atteint des résultats de bien-poséité améliorés par rapport aux équations de Navier–Stokes classiques, spécifiquement concernant la disponibilité d'estimations d'énergie $L^2$ inconditionnelles et le lissage d'ordre supérieur.

Les auteurs notent que leurs estimations uniformes suggèrent qu'à mesure que $\varepsilon \to 0$, les solutions fortes convergent vers une solution d'un système de Navier–Stokes hypervisqueux incompressible. Cependant, ils précisent explicitement que l'identification du système limite et la preuve des taux de convergence en $\varepsilon$ sont des directions laissées pour des travaux futurs. De plus, bien que la régularité globale pour des données arbitraires soit connue pour le système de Navier–Stokes hypervisqueux incompressible (grâce à Lions), les auteurs ne prétendent pas un tel résultat pour le système pénalisé avec de grandes données, notant que la structure compressible asymétrique de Temam modifie les estimations d'énergie d'ordre supérieur d'une manière qui empêche une extension directe du résultat pour grandes données.

Le travail se distingue des études précédentes sur les équations de Navier–Stokes amorties (par exemple, avec des termes d'absorption) en se concentrant sur un mécanisme d'amortissement d'ordre supérieur et une méthode de pénalisation qui élimine le terme de pression, simplifiant ainsi la formulation faible en supprimant la contrainte de divergence sur les fonctions tests.