The global well-posedness of the multi-dimensional compressible Euler system with damping in the $L^p$ critical Besov spaces for $p<2$

Résumé Technique : Globalité bien posée du système d'Euler compressible avec amortissement en espaces de Besov critiques $L^p$ pour $p < 2$

1. Problématique et Contexte

L'article étudie le problème de Cauchy pour le système d'Euler compressible avec amortissement linéaire, un modèle décrivant l'écoulement d'un gaz à travers un milieu poreux où la matrice solide exerce une force de friction proportionnelle et opposée à la quantité de mouvement du fluide. Le système est régi par les équations :
$$\begin{cases} \partial_t\rho + \text{div}(\rho u) = 0, \\ \partial_t(\rho u) + \nabla\cdot (\rho u \otimes u) + \nabla P(\rho) + \alpha \rho u = 0, \end{cases}$$
où $\rho$ est la densité, $u$ la vitesse, et $\alpha > 0$ le coefficient d'amortissement.

Contexte mathématique :
Des résultats antérieurs ont établi l'existence et le comportement asymptotique de solutions dans des espaces de Sobolev et des espaces de Besov critiques de type $L^2$ (notamment les travaux de Fang, Xu, Crin-Barat, Danchin et al.). Récemment, des extensions ont été faites vers des espaces hybrides $L^2-L^p$ pour $p \ge 2$.
La question centrale : Peut-on étendre ces résultats de bien-posé global au cas $1 \le p < 2$ ? Ce cas est particulièrement difficile car les systèmes hyperboliques sans termes d'ordre zéro sont généralement mal posés dans les espaces $L^p$ pour $p \neq 2$. L'amortissement joue ici un rôle crucial pour stabiliser le système.

2. Méthodologie et Outils Techniques

Pour aborder ce problème, les auteurs adoptent une approche basée sur l'analyse de Fourier et les espaces de Besov critiques, en divisant le spectre de fréquences en basses et hautes fréquences.

A. Reformulation du système
Les auteurs introduisent une nouvelle variable inconnue, l'enthalpie $n(\rho)$, définie par $n(\rho) \triangleq \int_1^\rho \frac{P'(s)}{s} ds$. Le système est réécrit en termes de $(n, v)$ (où $v$ est la vitesse) :
$$\begin{cases} \partial_t n + v \cdot \nabla n + P'(\bar{\rho})\text{div}v + G(n)\text{div}v = 0, \\ \partial_t v + v \cdot \nabla v + \nabla n = -\frac{1}{\varepsilon}v, \end{cases}$$
avec $\varepsilon$ lié au coefficient d'amortissement.

B. Cadre fonctionnel : Espaces de Besov Hybrides
L'analyse se déroule dans un espace fonctionnel $E^J_\varepsilon$ défini par une séparation des fréquences à un seuil $J_\varepsilon \approx -\log_2 \varepsilon$.

• Basses fréquences ($j \le J_\varepsilon$) : Analysées dans l'espace $L^p$ critique $\dot{B}^{\frac{d}{p}}_{p,1}$.
• Hautes fréquences ($j > J_\varepsilon$) : Analysées dans l'espace $L^2$ critique $\dot{B}^{\frac{d}{2}+1}_{2,1}$.
  Cette structure hybride permet de traiter la singularité des données initiales (plus concentrées pour $p<2$) tout en conservant la régularité nécessaire pour les termes non linéaires.

C. Estimations a priori et Lemme de Produit Clé
La difficulté majeure réside dans le contrôle des termes non linéaires (convection, pression) dans le cadre $L^p$ pour $p<2$.

• Nouvelle estimation de produit (Lemme 2.1) : Les auteurs établissent une nouvelle inégalité de produit dans les espaces hybrides $L^p-L^2$. Pour $1 \le p < 2$ et $0 < s_1 < d/2$, ils prouvent :
  $$\|ab\|_{\ell, J_\varepsilon}^{\dot{B}^{s_1}_{p,1}} \lesssim \left( \|b\|_{\ell, J_\varepsilon}^{\dot{B}^{s_1}_{p,1}} + 2^{(s_1 - d/p)J_\varepsilon} \|b\|_{h, J_\varepsilon}^{\dot{B}^{d/2}_{2,1}} \right) \|a\|_{\dot{B}^{d/2}_{2,1}}$$
  Cette estimation permet de gérer l'interaction entre les basses fréquences (où la norme $L^p$ est faible) et les hautes fréquences (où la norme $L^2$ domine).

D. Stratégie de preuve
La preuve de l'estimation a priori globale se décompose en deux étapes :

1. Estimation des basses fréquences ($L^p$) : En introduisant une "vitesse effective" $z = \varepsilon^{-1}v + \nabla n$, l'équation pour $n$ devient une équation de chaleur et celle pour $z$ une équation amortie. Les auteurs utilisent l'estimation de produit ci-dessus pour contrôler les termes non linéaires, obtenant une estimation de type $\|X(T)\| \lesssim X(0) + X(T)^2$.
2. Estimation des hautes fréquences ($L^2$) : Une méthode d'énergie pondérée (Lyapunov) est utilisée, similaire aux travaux sur les systèmes dissipatifs. Cela permet de contrôler la dissipation à haute fréquence et de fermer l'estimation globale.

3. Résultats Principaux

Théorème 1.1 (Existence et Unicité Globale) :
Sous l'hypothèse $1 \le p < 2$ et $d \ge 1$, il existe une constante $\delta_1 > 0$ telle que si les données initiales $(n_0, v_0)$ satisfont :
$$\|(n_0, v_0)\|_{\ell, J_\varepsilon}^{\dot{B}^{d/p}_{p,1}} + \varepsilon \|(n_0, v_0)\|_{h, J_\varepsilon}^{\dot{B}^{d/2+1}_{2,1}} \le \delta_1,$$
alors le problème de Cauchy admet une solution unique globale dans l'espace $E^\infty_{J_\varepsilon}$ pour tout $\varepsilon > 0$.

Points clés du résultat :

• Extension du cadre $p$ : Le résultat valide la bien-posé global pour $p < 2$, comblant une lacune par rapport aux travaux antérieurs limités à $p \ge 2$.
• Espace plus petit : L'espace $\dot{B}^{d/p}_{p,1}$ pour $p < 2$ est strictement inclus dans $\dot{B}^{d/2}_{2,1}$. Cela signifie que l'existence est prouvée pour des données initiales plus singulières et plus concentrées (pertinente pour des sources ponctuelles ou des filaments de vortex).
• Limite de relaxation : Le résultat implique que la limite de relaxation vers l'équation du milieu poreux peut être justifiée pour $p < 2$ (Remark 1.3).

4. Signification et Impact

• Avancée théorique : Ce travail surmonte l'obstacle technique majeur lié à l'ill-posé des systèmes hyperboliques purs dans $L^p$ ($p \neq 2$) en exploitant l'amortissement et une décomposition hybride des fréquences.
• Nouveaux outils : Le lemme de produit établi (Lemme 2.1) est un outil technique puissant qui peut être appliqué à d'autres systèmes, tels que les systèmes de chimiotaxie hyperbolique-parabolique ou les équations de Navier-Stokes dans des cadres $L^p-L^2$ pour $p < 2$.
• Pertinence physique : En permettant l'étude de données initiales plus singulières (concentrées), le modèle devient plus applicable à des phénomènes physiques réels où les conditions initiales ne sont pas nécessairement lisses ou distribuées uniformément.

En conclusion, cet article étend significativement la théorie de la bien-posé globale pour les fluides compressibles amortis, ouvrant la voie à l'analyse de régimes physiques plus extrêmes et fournissant de nouveaux outils d'analyse harmonique pour les EDP non linéaires.