Classification of singular limits for free boundary and singularly perturbed elliptic problems: the Dancer-Yan spikes revisited

Énoncé du Problème
Cet article traite de la classification des limites singulières pour un problème à frontière libre issu de la physique des plasmas, modélisant spécifiquement les équilibres de plasma dans une section transversale de Tokamak. Le problème est formulé comme une équation elliptique singulièrement perturbée en dimension $d=2$ :
$$\begin{cases} -\Delta v = [v]_+^p & \text{dans } \Omega, \\ v = \gamma & \text{sur } \partial\Omega, \\ \int_\Omega [v]_+^p = I, \end{cases}$$
où $\Omega \subset \mathbb{R}^2$ est un domaine borné de classe $C^{2,\beta}$, $p \in (1, +\infty)$, et $I > 0$ représente un courant total fixé. Les auteurs étudient le comportement asymptotique des solutions lorsque le courant total $I \to +\infty$. Dans cette limite, la « région de plasma » $\Omega_+ = \{v > 0\}$ est censée développer une structure non triviale, pouvant se concentrer en « pics » (régions où $v$ est significativement plus grand que la valeur à la frontière).

Alors que le comportement de tels problèmes en dimensions supérieures ($d \ge 3$) a été établi comme une superposition de plusieurs pics finis (spécifiquement des pics de Dancer-Yan), le cas bidimensionnel présente des difficultés uniques. En $d=2$, le problème global limite standard $-\Delta w = [w-1]_+^p$ dans $\mathbb{R}^2$ avec une masse finie n'admet aucune solution, empêchant l'application directe des arguments de blow-up des dimensions supérieures.

Méthodologie
Les auteurs emploient une analyse de blow-up raffinée adaptée au cadre bidimensionnel, s'appuyant sur le cadre de Dancer et Yan mais introduisant de nouvelles techniques de rescalage pour gérer l'absence de solutions globales standards.

1. Transformation de Variables : Le problème est reformulé à l'aide des variables $(\alpha, \psi)$ où $v = \lambda^{1/(p-1)}(\alpha + \lambda \psi)$ et $\lambda = I^{1/q}$ ($q$ étant le conjugué de Hölder de $p$). Cela transforme le problème en une forme singulièrement perturbée impliquant un petit paramètre $\varepsilon_n \to 0$.
2. Bornes Intégrales : L'analyse repose sur deux bornes intégrales naturelles, notées (H1) et (H2), qui contrôlent respectivement la masse $(p-1)$ et la masse $p$ des solutions. Ces bornes sont cruciales pour exclure les solutions à « masse infinie » et assurer la compacité.
3. Rescalage Raffiné : Contrairement au cas $d \ge 3$ où les maxima locaux restent éloignés de 1, en $d=2$, les maxima $v_n(x_n)$ convergent vers 1. Les auteurs introduisent un rescalage délicat à plusieurs échelles impliquant les paramètres $s_n$ et $\theta_n$ (dérivés de la solution de l'équation d'Emden $\phi$) pour capturer la structure « microscopique » des pics.
   • Ils définissent une fonction rescalée $\tilde{v}_n$ et une fonction normalisée supplémentaire $u_n$ pour analyser le profil près des maximisateurs locaux.
4. Classification des Solutions Entières : Une étape technique clé consiste en la classification des solutions à masse finie de l'équation planaire $-\Delta w = [w]_+^p$ dans $\mathbb{R}^2$. Les auteurs prouvent que de telles solutions sont soit triviales, soit radiales et prennent une forme spécifique impliquant la solution d'Emden $\phi$ à l'intérieur d'un disque et une décroissance logarithmique à l'extérieur.
5. Extraction de Pics Multiples : Pour gérer le regroupement de plusieurs pics au même point, les auteurs utilisent un argument inductif. Ils prouvent que deux pics distincts ne peuvent pas être arbitrairement proches par rapport à leurs paramètres de rescalage et que la présence d'un pic de Type I ou de Type II exclut l'existence d'un pic « s'estompant » (Fading) au même emplacement.
6. Identité de Pohozaev : Pour la structure globale, les auteurs utilisent l'identité de Pohozaev pour dériver une contrainte sur les emplacements des pics, montée que ceux-ci doivent être des points critiques d'un Hamilton de Kirchhoff-Routh.

Contributions Principales et Résultats

1. Classification des Limites Singulières (Théorème 1.1) :
   Les auteurs établissent que sous les bornes intégrales (H1) et (H2), la limite singulière des solutions tombe dans l'une des trois catégories suivantes autour de toute séquence de maximisateurs locaux :

   • Vanissement (Vanishing) : La solution reste en dessous du seuil 1 dans les sous-ensembles compacts.
   • Pics de Type I : La solution converge vers un profil de « pic de Dancer-Yan ». Ici, la fonction rescalée converge vers la solution entière standard $w^*$, et la masse est quantifiée.
   • Pics de Type II : Un phénomène nouveau spécifique à $d=2$. La fonction rescalée converge également vers $w^*$, mais les paramètres de rescalage divergent de telle sorte que la masse $p$ tend vers zéro dans la limite, tandis que la masse $(p-1)$ reste quantifiée.
   • Pics S'estompant (Fading Spikes) : L'amplitude du pic décroît si rapidement qu'il disparaît effectivement dans la limite, la solution se comportant comme une fonction harmonique perturbée par un terme de faible amplitude.

   Crucialement, l'article démontre qu'contrairement aux dimensions supérieures, toute solution n'est pas une simple superposition de pics de Dancer-Yan. La présence de pics de Type II et de pics s'estompant indique une structure singulière plus riche en 2D.

2. Structure Globale avec Conditions de Dirichlet (Théorème 1.2) :
   En imposant des conditions aux limites de Dirichlet ($v=0$ sur $\partial\Omega$) et une condition de non-annulation (NVp) sur la masse $p$, les auteurs fournissent une classification détaillée pour la solution globale, tout en notant qu'une description complète de l'interaction entre différents types de pics se regroupant au même point reste un problème ouvert.

   • L'ensemble singulier $\Sigma$ consiste en un nombre fini de points intérieurs.
   • La solution est une superposition de pics de Type I et de Type II centrés en ces points.
   • La masse $(p-1)$ est quantifiée : elle est égale à $(N_I + N_{II})I_{p-1}$, où $N_I$ et $N_{II}$ sont les comptages des pics de Type I et de Type II.
   • La masse $p$ n'est pas quantifiée en général ; elle dépend des limites de rescalage spécifiques des pics de Type I.
   • La région de plasma $\Omega_{n,+}$ consiste en des composantes asymptotiquement rondes autour des centres des pics.
   • Les emplacements des pics de Type I sont des points critiques d'un Hamilton de Kirchhoff-Routh impliquant la fonction de Green du domaine.
   • Note : Bien que les auteurs parviennent à exclure les pics « s'estompant » en présence de pics de Type I ou de Type II sous ces conditions, ils déclarent explicitement qu'une description complète de l'interaction locale lorsque plusieurs types de pics (Type I et Type II) se regroupent au même point est actuellement manquante.

3. Application au Problème de Plasma Original (Théorème 1.3) :
   Les auteurs traduisent ces résultats vers les variables originales $(v, \gamma)$ pour le problème du plasma. Ils montrent que sous des contraintes asymptotiques spécifiques sur les paramètres $\lambda_n$ et $\alpha_n$, les solutions présentent les structures de pics classifiées. Ils fournissent des formules asymptotiques précises pour la valeur à la frontière $\gamma$ et le courant total en fonction du nombre et du type de pics.

Signification et Revendications
L'article prétend fournir la première description détaillée de la limite singulière pour ce problème à frontière libre en dimension $d=2$. Sa signification primaire réside dans la révélation que le cas bidimensionnel est fondamentalement différent du cas de dimension supérieure :

• Richesse de la Structure : L'existence de pics de Type II et de pics s'estompant montre que l'espace des solutions est plus complexe qu'un simple collage de profils standards.
• Phénomène de Quantification : L'article met en évidence un comportement de quantification distinct où la masse $(p-1)$ est quantifiée (similaire aux équations de type Liouville), mais la masse $p$ ne l'est pas, défiant l'intuition héritée des dimensions $d \ge 3$.
• Avancée Méthodologique : Le travail introduit un raffinement nécessaire des techniques de rescalage pour gérer la disparition de la valeur maximale en 2D, surmontant l'obstacle du fait que le problème global limite standard ne possède pas de solutions.

Les auteurs déclarent explicitement qu'une description complète de l'interaction entre plusieurs pics de types différents se regroupant au même point reste un problème ouvert, bien qu'ils parviennent à exclure certaines configurations (comme la coexistence de pics s'estompant avec des pics de Type I/II) sous des conditions de Dirichlet. Ils conjecturent également que les pics de Type II pourraient ne pas exister sur des domaines convexes, suggérant une dépendance vis-à-vis de la géométrie du domaine.