Gradient blowup of smooth vacuum solutions to 1D compressible Euler equations

1. Énoncé du problème

L'article étudie les équations d'Euler compressibles isentropiques en une dimension pour un gaz polytropique ($p = \rho^\gamma$) dans le domaine de la demi-droite $\Omega = \{x > 0\}$. Le fluide est en contact avec une frontière de vide stationnaire en $x=0$.

• Équations gouvernantes : Le système est écrit en termes de densité $\rho$, de vitesse $u$ et de vitesse du son $c$ :
  $$\rho_t + (\rho u)_x = 0, \quad (\rho u)_t + (\rho u^2 + p)_x = 0$$
  Dans la région non vide, cela est équivalent à un système pour la vitesse $u$ et la vitesse du son $c$.
• Conditions aux limites : La frontière de vide est stationnaire ($x=0$) avec $c(t,0)=0$ et $u(t,0)=0$.
• Le phénomène : Les auteurs étudient l'explosion du gradient (perte de régularité) se produisant à la frontière de vide en temps fini. Plus précisément, ils construisent des solutions initialement lisses (spécifiquement $C^{\lfloor \beta \rfloor, \beta-\lfloor \beta \rfloor}$) et de carré intégrable, mais qui évoluent pour devenir seulement Hölder-continues de classe $C^{1-\mu}$ près de la frontière au moment de l'explosion $t=0$, où les premières dérivées de la vitesse et de la vitesse du son divergent.
• Paramètre de régularité : La régularité est caractérisée par $\mu \in [1/2, 1)$ et $\beta = \frac{1}{1-\mu} \geq 2$. Le cas $\mu = 1/2$ correspond au profil le plus stable (lié à la singularité de rebord $C^{1/3}$ dans l'équation de Burgers).

2. Méthodologie

La preuve repose sur une analyse de stabilité des solutions auto-similaires construites dans un travail antérieur [27]. La méthodologie implique plusieurs étapes sophistiquées :

A. Reformulation auto-similaire

Les auteurs utilisent la symétrie d'échelle des équations d'Euler pour introduire des variables auto-similaires $(\tau, y)$ :
$$\tau = -\log(\mu(-t)), \quad y = \frac{x}{(-t)^\delta}$$
où $\delta = 1/\mu$. Le système est transformé en un ensemble d'équations pour les invariants de Riemann $(\hat{z}, \hat{w})$. Les solutions auto-similaires construites dans [27] apparaissent comme des états stationnaires $(\bar{z}, 0)$ dans ces variables.

B. Localisation et construction du profil

Puisque les solutions auto-similaires sont non bornées à l'infini, les auteurs construisent un profil localisé $\mathring{z} = \chi \bar{z}$ en utilisant une fonction de coupure dépendante du temps $\chi$. Cette coupure est transportée par la vitesse caractéristique de l'écoulement de fond pour garantir que les termes d'erreur décroissent de manière appropriée.

C. Décomposition des perturbations

La solution est décomposée en trois parties :

1. Le profil : La solution auto-similaire localisée $\mathring{z}$.
2. Termes de modulation : Une somme finie de polynômes en $y$ avec des coefficients dépendant du temps $(z_i(\tau), w_i(\tau))$. Ces coefficients correspondent au développement de Taylor de la solution à l'origine.
3. Reste : Composantes haute fréquence $(Z, W)$ qui s'annulent à la frontière jusqu'à un ordre élevé ($n_\mu - 1$).

La décomposition est :
$$\hat{z} = \mathring{z} + \sum_{i=2}^{n_\mu-1} \phi(y) \frac{z_i(\tau)}{i!} y^i + Z, \quad \hat{w} = \sum_{i=2}^{n_\mu-1} \phi(y) \frac{w_i(\tau)}{i!} y^i + W$$

D. Analyse de modulation (Système d'EDO)

Les coefficients $(z_i, w_i)$ satisfont un système d'équations différentielles ordinaires (EDO) dérivé en imposant des conditions d'annulation sur les dérivées d'ordre supérieur à la frontière.

• Stabilité : Les auteurs identifient une variété stable de dimension finie pour ces EDO. Pour empêcher une croissance exponentielle (instabilité) entraînée par des termes anti-amortis, les données initiales doivent satisfaire des contraintes algébriques spécifiques (relations polynomiales entre $z_i$ et $w_i$).
• Symétries : Les contraintes $z_1 = w_1 = 0$ éliminent les symétries de translation temporelle et spatiale, fixant le temps d'explosion à $t=0$.

E. Estimations d'énergie pour le reste

Pour le reste $(Z, W)$, les auteurs effectuent une analyse de stabilité non linéaire rigoureuse en utilisant des arguments de bootstrap et des estimations d'énergie :

• Poids singuliers : Ils emploient des espaces $L^2$ pondérés avec des poids singuliers $y^{-n_\mu}$ pour gérer la dégénérescence à la frontière de vide.
• Estimations d'ordre élevé : Ils dérivent des estimations dans $\dot{H}^{n_\mu}$ (espaces de Sobolev) pour contrôler la régularité.
• Dissipation : L'opérateur linéarisé autour du profil présente une dissipativité stricte pour les hautes fréquences en raison de la vitesse de transport et du choix spécifique de l'indice de régularité $n_\mu$.
• Gestion de $\beta$ non entier : Une attention particulière est portée lorsque $\beta$ n'est pas un entier (c'est-à-dire $\mu \notin S$), où le profil a une régularité limitée. Dans les cas où $\beta$ est un demi-entier, une contrainte supplémentaire $w_2 = 0$ est imposée pour éviter des singularités non intégrables dans les termes d'interaction.

3. Contributions et résultats clés

Théorème principal (Théorèmes 1.1 & 2.7)

Les auteurs prouvent l'existence d'une large classe de données initiales lisses, de carré intégrable et non auto-similaires qui évoluent en des solutions présentant une explosion de gradient à la frontière de vide stationnaire.

• Transition de régularité : Les solutions commencent comme $C^{\lfloor \beta \rfloor, \beta-\lfloor \beta \rfloor}$ et passent à $C^{1-\mu}$ à $t=0$.
• Taux d'explosion : Près du point d'explosion $(0,0)$, la vitesse et la vitesse du son se comportent asymptotiquement comme :
  $$u(t, x) \sim -\tilde{a} x^{1-\mu}, \quad c(t, x) \sim \tilde{b} x^{1-\mu}$$
  conduisant à des gradients non bornés $\partial_x u, \partial_x c \to \infty$.
• Cofinitude : L'ensemble des données initiales menant à ce scénario d'explosion spécifique forme un sous-ensemble de codimension finie de l'espace de Sobolev $H^{n_\mu}$.
  • La codimension est $\lfloor \beta + 1 \rfloor$ pour $\mu$ générique.
  • Pour $\mu = 1/2$, la codimension est 3, correspondant aux trois symétries (translation temporelle, boost galiléen, mise à l'échelle).
  • Pour $\mu \in H$ ($\beta$ demi-entier), une contrainte supplémentaire $w_2=0$ est requise, augmentant la codimension.

Innovations techniques

1. Stabilité des solutions auto-similaires : Il s'agit de la première preuve rigoureuse de la stabilité des solutions auto-similaires de "temps d'attente" pour les équations d'Euler compressibles avec une frontière de vide stationnaire, étendant les résultats antérieurs sur l'équation de Burgers.
2. Gestion des frontières stationnaires : Contrairement aux problèmes de frontière libre où le vide se déplace, ce travail traite le cas où la frontière est fixe ($u=0$ en $x=0$) mais où le gradient explose, nécessitant une analyse minutieuse de l'interaction entre la condition aux limites et le profil auto-similaire.
3. Modulation avec poids singuliers : Le développement d'estimations d'énergie avec des poids singuliers $y^{-n_\mu}$ combiné à la théorie de la modulation pour contrôler les modes basse fréquence près de la singularité.

4. Importance

• Physique mathématique : L'article résout une question de longue date concernant la généralité de l'explosion de gradient dans les écoulements d'Euler avec vide. Il démontre que l'explosion de gradient n'est pas seulement une caractéristique de solutions auto-similaires spéciales, mais qu'elle est stable sous les perturbations d'une large classe de données initiales lisses.
• Lien avec l'équation de Burgers : Les résultats établissent un lien rigoureux entre la formation de singularités de rebord dans l'équation de Burgers et le système plus complexe d'Euler compressible, validant l'utilisation des profils de Burgers comme modèles de base pour la formation de singularités en dynamique des gaz.
• Dynamique de la frontière de vide : Il clarifie le comportement du phénomène de "temps d'attente", montrant que même avec une condition aux limites stationnaire, le système peut évoluer naturellement vers un état de vide physique (où la frontière commencerait à bouger) si la régularité est perdue, ou maintenir l'état stationnaire avec des gradients infinis.
• Perspectives futures : Ce travail ouvre la voie à l'étude de la continuation de ces solutions au-delà du temps d'explosion, menant potentiellement à la formation de frontières libres ou de discontinuités faibles, un problème laissé pour des études ultérieures.

En résumé, cet article fournit un cadre mathématique rigoureux pour comprendre comment des données initiales lisses et physiquement pertinentes en dynamique des gaz peuvent évoluer vers des états singuliers avec des gradients infinis à une frontière de vide, établissant la stabilité de ces singularités face aux perturbations.