Free boundary regularity and well-posedness of physical solutions to the supercooled Stefan problem

1. Problématique et Contexte

L'article étudie le problème de Stefan surfondu (supercooled Stefan problem) en une dimension, qui modélise le gel de l'eau surfondue en glace. Le système est décrit par les équations suivantes :
$$\begin{cases} u_t - \frac{1}{2}u_{xx} = 0 & t > 0, x > \Lambda_t, \\ \dot{\Lambda}_t = \frac{\alpha}{2} u_x(\Lambda_t, t) & t > 0, \\ u(x, 0) = u_0(x), \quad u(\Lambda_t, t) = 0 & x \ge 0, t \ge 0, \end{cases}$$
où $u$ est la température (négative), $\Lambda_t$ est la frontière libre (front de gel), et $\alpha > 0$.

Défis principaux :

• Contrairement au problème de Stefan classique, ce système est mal posé : les solutions classiques peuvent exploser en temps fini, entraînant des discontinuités de saut dans la frontière $\Lambda$ et des pics de température instantanés.
• La définition de la solution au-delà de l'explosion nécessite une reformulation probabiliste. On considère une population de particules browniennes interagissant via la frontière mobile.
• Une notion de solution physique a été introduite pour sélectionner une solution unique parmi les solutions probabilistes possibles. Elle impose que la taille du saut de la frontière soit minimale (condition de saut physique).
• Des questions ouvertes subsistaient concernant l'accumulation des temps de saut, la régularité globale de la frontière, et l'unicité globale pour des données initiales très générales (notamment oscillantes).

2. Méthodologie

L'auteur développe une approche analytique combinant l'analyse stochastique, la théorie des équations aux dérivées partielles (EDP) et l'analyse des problèmes d'obstacle.

A. Reformulation en Problème d'Obstacle Pondéré
L'analyse repose sur la transformation de la densité de probabilité $u$ en un potentiel $w$ défini par :
$$w(x, t) = \int_t^\infty u(x, s) ds.$$
Ce potentiel satisfait un problème d'obstacle parabolique pondéré :
$$w_t - \frac{1}{2}w_{xx} = -\nu(x)\chi_{\{w>0\}},$$
où $\nu$ est la densité de probabilité des particules arrêtées. Cette transformation permet d'utiliser la théorie classique de la régularité des problèmes d'obstacle.

B. Estimations Clés
L'article établit des estimations cruciales pour les solutions probabilistes et physiques :

1. Estimations d'ordre deux : Bornes supérieures sur la dérivée temporelle $u_t$ et la dérivée spatiale $u_{xx}$ (dans le sens des distributions).
2. Estimation de non-dégénérescence linéaire : Pour les solutions physiques, la densité $u$ croît linéairement par rapport à la distance à la frontière : $u(x, t) \ge c(x - \Lambda_t)_+$. Cela implique que le temps de gel $s(x) = \inf\{t : \Lambda_t > x\}$ est localement lipschitzien.

C. Propagation Rétrograde des Oscillations
C'est un outil technique majeur. L'auteur démontre que tout changement de signe de la dérivée spatiale $u_x$ à un temps futur $t^*$ peut être "transporté" vers l'arrière le long de courbes continues jusqu'à un temps initial $t_0$, sans toucher la frontière libre.

• Lien avec les sauts : Chaque saut de la frontière $\Lambda$ à un temps $t^*$ force une oscillation supplémentaire de $u$ près de la frontière à ce même instant.
• En propageant ces oscillations vers $t_0$, on peut compter le nombre de sauts futurs en fonction du nombre d'oscillations passées.

3. Résultats Principaux

A. Absence d'accumulation des temps de saut (Théorème 1.3)

L'auteur répond à une question ouverte de [8] : Les temps de discontinuité de $\Lambda$ ne peuvent pas avoir de points d'accumulation dans $(0, \infty)$.

• Grâce à la propriété de propagation rétrograde, si les sauts s'accumulaient, cela impliquerait un nombre infini d'oscillations de $u$ à un temps initial fini, ce qui contredit les estimations de régularité locale.
• Conséquence : Les oscillations initiales sont résolues instantanément, et les sauts ne peuvent s'accumuler qu'à $t=0$.

B. Régularité de la Frontière Libre (Théorèmes 1.4 et 1.6)

L'article fournit une classification complète et des résultats de régularité sans hypothèse de monotonie ou de bornitude sur la donnée initiale $u_0$ :

1. Régularité $C^1$ globale : La fonction temps de gel $s(x)$ est de classe $C^1$ sur $(\Lambda_0, \alpha)$.
2. Régularité $C^\infty$ : Il existe un ouvert $R \subset (\Lambda_0, \alpha)$ dont le complémentaire est dénombrable, tel que $s \in C^\infty(R)$.
3. Classification des points :
   • Points réguliers : Là où $u$ s'annule continûment ou à l'intérieur strict d'un saut. La frontière est lisse ($C^\infty$) ici.
   • Points singuliers : Les extrémités des sauts ou les points où $u$ chute de la valeur critique $\alpha^{-1}$ à 0. L'ensemble des points singuliers est dénombrable.
   • À l'inverse des résultats précédents ([8]), le cas intermédiaire (régularité $C^{1/2}$) n'apparaît jamais pour $t > 0$.

C. Unicité Globale (Théorème 1.7 et Corollaires)

L'auteur prouve que l'unicité locale implique l'unicité globale pour les solutions physiques.

• Mécanisme : Toute solution physique devient localement monotone (en nombre de changements de monotonie finis sur les compacts) à tout temps $t > 0$, même si la donnée initiale $u_0$ est très oscillante.
• Applications :
  • Corollaire 1.8 : Unicité globale pour des données initiales satisfaisant une condition locale près de $x=0$ (pas d'oscillation excessive autour de $\alpha^{-1}$), répondant à une question ouverte de [21].
  • Corollaire 1.9 : Unicité pour une famille spécifique de données initiales oscillantes (cas critique de [29]), résolvant la question de l'unicité globale pour ces cas.

4. Signification et Impact

Ce travail représente une avancée majeure dans la théorie du problème de Stefan surfondu :

1. Généralité des hypothèses : Les résultats de régularité et d'unicité s'appliquent à des données initiales très générales (intégrables, non bornées, fortement oscillantes), dépassant les régimes de monotonie locale ou de données bornées traités précédemment.
2. Structure fine de la singularité : La démonstration que l'ensemble des points singuliers est dénombrable et que la régularité est $C^\infty$ presque partout (sauf sur un ensemble dénombrable) offre une compréhension beaucoup plus précise de la géométrie de la frontière libre que les résultats antérieurs.
3. Résolution de conjectures : L'article résout définitivement la question de l'accumulation des sauts et de l'unicité globale pour des classes larges de données, unifiant ainsi les régimes de solutions précédemment étudiés séparément.
4. Méthodologie innovante : L'utilisation de la propagation rétrograde des oscillations et de l'analyse de blow-up sur le problème d'obstacle pondéré fournit de nouveaux outils puissants pour l'étude des EDP non linéaires avec conditions aux limites libres et discontinuités.

En résumé, Sebastian Munoz établit un cadre robuste pour la régularité et l'unicité des solutions physiques du problème de Stefan surfondu, démontrant que malgré la nature potentiellement explosive du système, la structure de la solution et de sa frontière libre reste hautement régulière et bien comportée pour tout temps strictement positif.