Unstable Interface Dynamics for Gravity Stokes Flow

Résumé Technique : Dynamique d'Interface Instable pour l'Écoulement de Stokes sous Gravité

Énoncé du Problème
Cet article étudie la dynamique à long terme et les mécanismes d'instabilité du problème d'interface nette pour un écoulement de Stokes en deux dimensions sous l'effet de la gravité. Le système modélise deux fluides incompressibles et homogènes possédant des densités constantes distinctes ($\rho_+$ et $\rho_-$) évoluant selon la loi de Stokes sous une force gravitationnelle. L'interface $\Gamma(t)$ sépare la région supérieure $\Omega_+(t)$ de la région inférieure $\Omega_-(t)$ dans un domaine périodique horizontal ($\mathbb{T} \times \mathbb{R}$).

L'étude se concentre sur deux régimes primaires :

1. Régime Instable (Rayleigh-Taylor) : Le fluide le plus dense se trouve au-dessus du fluide le plus léger ($\rho_+ < \rho_-$). Bien que la bien-posée globale pour des données initiales lisses soit connue (ex. [12]), le comportement à long terme des quantités géométriques (périmètre, courbure) reste largement inexploré.
2. Régime Stable : Le fluide le plus dense se trouve en dessous du fluide le plus léger ($\rho_+ > \rho_-$). L'article examine si des configurations initiales stables peuvent évoluer vers des états instables, en étudiant spécifiquement les « instabilités de retournement » (turning instabilities) où l'interface cesse d'être un graphe en un temps fini.

Le système est analysé via l'Équation de Dynamique de Contour (CDE) dérivée de [11] :
$$\partial_t z(\alpha, t) = (\rho_- - \rho_+) \int_{\mathbb{T}} S(z(\alpha, t) - z(\beta, t)) \dot{z}^\perp(\beta, t) z^2(\beta, t) d\beta$$
où $S$ est le stokeslet $\mathbb{T}$-périodique en $x_1$.

Méthodologie
Les auteurs emploient une combinaison de techniques analytiques et de simulations numériques :

• Analyse de l'Énergie Potentielle : Un outil central est un fonctionnel d'énergie potentielle modulée $E(t) = \int_{\mathbb{T}\times\mathbb{R}} x_2(\rho_s(x_2) - \rho(x,t)) dx$, où $\rho_s$ représente l'état stratifié instable. Les auteurs prouvent que $E(t)$ est monotone croissante et que sa dérivée temporelle est liée à la norme $L^2$ d'un opérateur non local appliqué à la densité.
• Exploitation des Symétries : L'analyse repose fortement sur des symétries spécifiques des données initiales : la symétrie centrale ($z(\alpha) = -z(-\alpha)$) et la symétrie paire par rapport à $x_1 = \pm \pi/2$. Ces symétries sont prouvées comme étant conservées par le flux et permettent la réduction du domaine à une bande restreinte pour les estimations géométriques.
• Dichotomie Géométrique : Pour prouver les résultats de croissance, les auteurs établissent une dichotomie basée sur la quantité $\delta(t) = \|( -\Delta)^{-1} \partial_1 \rho(t) \|_{L^2}^2$. Ils démontrent que si le périmètre $L(t)$ ne croît pas suffisamment vite, la quantité $\delta(t)$ doit devenir petite, ce qui force en retour la courbure maximale $K(t)$ à croître rapidement.
• Analyticité et Fingering (Doigtement) : Pour l'analyse du doigtement, les auteurs utilisent l'analyticité des solutions dans une bande complexe. Ils appliquent un lemme concernant la distribution des points critiques de fonctions analytiques pour estimer le nombre de « doigts » (extrema locaux) pouvant se développer.
• Construction de l'Instabilité de Retournement : Pour prouver la rupture en temps fini dans le régime stable, les auteurs construisent des données initiales spécifiques satisfaisant des conditions de symétrie paire/impaire. Ils analysent la dérivée temporelle de la pente à l'origine pour montrer qu'une interface stable peut développer une tangente verticale et se retourner.
• Simulations Numériques : Les auteurs implémentent un schéma numérique de type spectral utilisant des différences centrales et la règle de Simpson pour les intégrales singulières, régularisées par une viscosité artificielle, afin de visualiser la croissance de la longueur de l'interface et de la courbure.

Contributions Clés et Résultats

1. Croissance à Temps Infini dans le Régime Instable (Théorème 1.1) :
   L'article prouve que pour des données initiales instables satisfaisant des conditions de symétrie spécifiques, l'interface présente une croissance non bornée d'au moins une quantité géométrique (périmètre $L(t)$ ou courbure maximale $K(t)$) quand $t \to \infty$. Plus précisément, pour tout $\epsilon > 0$ :
   $$\limsup_{t \to \infty} t^{-\frac{1}{17} + \epsilon} (L(t) + K(t)) = \infty$$
   Ce résultat établit que les solutions ne se stabilisent pas vers un état stationnaire mais présentent plutôt une croissance générique à long terme, contrastant avec le résultat d'existence globale qui est indépendant de la stratification de densité.

2. Estimation sur le Phénomène de Doigtement (Théorème 1.3) :
   Les auteurs fournissent une borne inférieure sur la croissance du nombre de doigts (points critiques de la pente de l'interface) pour une classe de solutions analytiques. Ils montrent que le nombre de doigts $F(t)$ croît au moins comme $C(1 + \sqrt{t})$. Cela relie la croissance analytique de la solution au phénomène physique du doigtement visqueux.

3. Instabilité de Retournement en Temps Fini dans le Régime Stable (Théorème 1.4) :
   Contrairement à l'intuition selon laquelle la stratification stable empêche l'instabilité, les auteurs prouvent l'existence de données initiales lisses dans le régime stable ($\rho_+ > \rho_-$) qui évoluent vers une instabilité de retournement en temps fini. L'interface cesse d'être un graphe, et la condition de Rayleigh-Taylor s'effondre. L'article construit une famille de telles courbes initiales et prouve que le retournement se produit en des points symétriques spécifiques.

4. Validation Numérique :
   Les simulations numériques confirment les prédictions théoriques, montrant une croissance significative de la longueur de l'interface et une courbure bornée à court terme pour les configurations instables, ainsi que le développement de doigts dans l'évolution de l'interface.

Signification et Revendications
L'article affirme traiter une lacune dans la compréhension de la dynamique à long terme du système de transport de Stokes. Bien que la bien-posée globale ait été établie dans des travaux antérieurs [12], le comportement de ces solutions globales était inconnu. Les auteurs soutiennent que leurs résultats démontrent que :

• L'existence globale n'implique pas la stabilité ; les configurations instables mènent à une croissance géométrique non bornée.
• Les configurations stables ne sont pas à l'abri de l'instabilité ; des données initiales lisses spécifiques peuvent évoluer vers des singularités de retournement, brisant la propriété de graphe.
• La croissance des quantités géométriques (périmètre et courbure) est une caractéristique générique du régime instable, soutenant le lien avec les phénomènes de doigtement visqueux observés dans les écoulements en milieux poreux.

Les auteurs notent explicitement que le taux de croissance dérivé dans le Théorème 1.1 ($t^{1/17}$) n'est probablement pas optimal, les preuves numériques suggérant un taux de croissance linéaire, mais la preuve analytique fournit une borne inférieure rigoureuse pour l'instabilité. Le travail met en évidence la distinction entre la dynamique de l'écoulement de Stokes et celle, plus singulière, du problème de Muskat, notamment concernant la nature de l'amortissement linéaire et les mécanismes spécifiques d'instabilité.