Global existence of classical solutions for the multi-dimensional compressible Navier-Stokes-Poisson equations on solid balls for arbitrary spherically symmetric large initial data

Résumé Technique : Existence Globale de Solutions Classiques pour les Équations de Navier-Stokes-Poisson Compressibles

Titre de l'article : Global existence of classical solutions for the multi-dimensional compressible Navier-Stokes-Poisson equations on solid balls for arbitrary spherically symmetric large initial data
Auteurs : Jie Fan, Xiangdi Huang, Muxi Lei.

1. Le Problème Étudié

L'article s'attaque à un problème mathématique ouvert de longue date : l'existence de solutions classiques globales pour les équations de Navier-Stokes-Poisson (NSP) compressibles en dimension 3 (et 2) pour des données initiales arbitrairement grandes, sans hypothèse de petite perturbation.

Le système modélise le mouvement d'étoiles gazeuses ou de plasmas dans un champ gravitationnel ou électrique auto-cohérent. Les équations sont données par :
$$\begin{cases} \rho_t + \text{div}(\rho u) = 0, \\ \rho(u_t + u \cdot \nabla u) + \nabla P(\rho) = -\rho\nabla\Phi + \text{div}(\mu(\rho)Du) + \nabla(\lambda(\rho)\text{div}u), \\ \Delta\Phi = \kappa(\rho - \bar{\rho}), \end{cases}$$
où $\rho$ est la densité, $u$ la vitesse, $\Phi$ le potentiel, et $\mu, \lambda$ les coefficients de viscosité dépendant de la densité.

Contraintes spécifiques :

• Domaine : Boule solide (domaine borné contenant l'origine).
• Symétrie : Solutions sphériquement symétriques.
• Viscosité : Coefficients de type Bresch-Desjardins (BD) : $\mu(\rho) = \rho^\alpha$ et $\lambda(\rho) = (\alpha-1)\rho^\alpha$.
• Données initiales : Grandes données, loin du vide ($0 < \rho_* \le \rho_0 \le \rho^*$).
• Difficulté majeure : La singularité à l'origine ($r=0$) due au terme de Poisson ($\frac{x}{r}$ en 2D, $\frac{x}{r^2}$ en 3D) et la dégénérescence de la viscosité lorsque $\rho \to 0$ (bien que le vide soit exclu ici, la structure de la viscosité pose des défis d'estimation).

2. Méthodologie

Les auteurs développent une série d'estimations a priori sophistiquées pour étendre la solution locale (dont l'existence est connue) à une solution globale. La stratégie repose sur trois piliers principaux :

A. Utilisation de l'Entropie BD et de la Vitesse Effective

• Introduction de la vitesse effective $w = u + r^{N-1}(\rho^\alpha)_x$ (en coordonnées lagrangiennes).
• Exploitation de l'inégalité d'entropie de Bresch-Desjardins pour obtenir des estimations de dissipation supplémentaires sur le gradient de densité, cruciales pour contrôler les termes non-linéaires.

B. Gestion de la Singularité Centrale

• Le cœur de la difficulté réside dans le terme de Poisson qui devient singulier à l'origine.
• Les auteurs établissent des estimations pondérées de la densité près de l'origine (de la forme $\|\rho^{\alpha-1/2}r^\xi\|_{L^\infty}$).
• Ces estimations permettent de convertir les intégrales singulières (comme $\int \frac{x^2}{r^2} dx$) en quantités contrôlées par les normes $L^2$ de la vitesse et de la vitesse effective.
• Une relation clé est établie entre la coordonnée lagrangienne $x$ et la coordonnée eulérienne $r$, montrant que $x \le C r^\beta$, ce qui permet d'éliminer la singularité dans les estimations de haute intégrabilité.

C. Estimations de Haute Intégrabilité et Inégalités de Sobolev

• Pour obtenir des bornes supérieures et inférieures uniformes sur la densité (nécessaires pour éviter le vide et l'effondrement), les auteurs dérivent des estimations de haute intégrabilité pour $u$ et $w$ (normes $L^{2k}$).
• Ils utilisent des inégalités de type Gronwall couplées et des embeddings de Sobolev unidimensionnels adaptés au contexte radial.
• Le choix des paramètres $\alpha$ et $\gamma$ (exposant adiabatique) est crucial pour fermer les estimations, notamment via l'ensemble dense $A_{set}$ défini dans l'article.

3. Résultats Principaux

L'article établit l'existence et l'unicité de solutions classiques globales sphériquement symétriques pour des données initiales grandes, sous des conditions spécifiques sur les paramètres $\alpha$ et $\gamma$.

Cas 1 : $N=2$ (Dimension 2)

• Condition : $\frac{1}{2} < \alpha < 1$ et $\gamma > 2 - \alpha$.
• Résultat : Existence globale de solutions classiques avec des bornes uniformes sur la densité : $(C(T))^{-1} \le \rho(x,t) \le C(T)$.

Cas 2 : $N=3$ (Dimension 3)

• Condition : $\frac{5}{6} < \alpha < 1$ et $\frac{4}{3} < \gamma < 4\alpha - 1 - \alpha^2$.
• Note : Pour le cas des étoiles gazeuses ($\kappa=1$), la condition $\gamma > 4/3$ est nécessaire pour contrôler le terme de Poisson dans l'estimation d'énergie de base. Pour les plasmas ($\kappa=-1$), la condition est moins restrictive.
• Résultat : Existence globale de solutions classiques avec des bornes uniformes sur la densité.

Cas 3 : Le cas limite $\alpha = 1$

• Condition : $\gamma > 1$ en 2D, et $\frac{4}{3} < \gamma < 3$ en 3D.
• Résultat : L'existence globale est également prouvée pour ce cas de viscosité constante (ou linéaire), étendant les résultats précédents.

4. Contributions Clés et Innovations

1. Résolution d'un problème ouvert en dimension réelle : C'est le premier résultat établissant l'existence de solutions classiques globales pour des données initiales arbitrairement grandes dans un domaine multidimensionnel véritable (boule solide) pour les équations NSP compressibles avec viscosité dépendante de la densité.
2. Traitement de la singularité à l'origine : La méthode proposée pour gérer la singularité du terme de Poisson à $r=0$ via des estimations pondérées spécifiques et la relation $x-r$ est une avancée technique majeure, permettant de contourner les obstacles qui avaient jusqu'alors limité les résultats aux données petites ou aux domaines sans cœur.
3. Optimisation des conditions sur $\gamma$ : Les auteurs précisent les plages admissibles pour l'exposant adiabatique $\gamma$ en fonction de $\alpha$, montrant que des solutions globales existent même pour des gaz avec des équations d'état complexes, tant que la viscosité suit la relation BD.
4. Unification des cas $\kappa=1$ et $\kappa=-1$ : Bien que traités séparément en raison de la nature du terme de Poisson (gravité vs électricité), la méthodologie couvre les deux cas physiques (étoiles et plasmas) dans le même cadre théorique.

5. Signification Scientifique

Ce travail représente une percée significative dans la théorie mathématique des fluides compressibles.

• Physique : Il valide rigoureusement l'existence de solutions pour des modèles d'étoiles et de plasmas soumis à de fortes perturbations initiales, sans risque de formation de singularités en temps fini (dans le cadre de la symétrie sphérique).
• Mathématiques : Il démontre la puissance de la structure d'entropie de Bresch-Desjardins combinée à l'analyse de symétrie pour surmonter les difficultés liées aux singularités géométriques et aux viscosités dégénérées.
• Perspectives : Ces résultats ouvrent la voie à l'étude de problèmes similaires dans des géométries plus complexes ou pour des systèmes couplés plus généraux, en fournissant un cadre robuste pour le contrôle des singularités centrales.

En conclusion, l'article répond affirmativement à la question de l'existence globale pour des données grandes dans un cadre multidimensionnel réaliste, en surmontant les obstacles analytiques liés à la singularité centrale et à la non-linéarité forte du système.