Stability and instability of a one-dimensional MHD model

1. Problématique et Contexte

L'article étudie la stabilité et l'instabilité d'un modèle de magnétohydrodynamique (MHD) en une dimension, introduit par Dai et al. (2023). Ce modèle, posé sur le tore $\mathbb{T} = [-\pi, \pi]$, généralise le modèle de De Gregorio pour les équations d'Euler incompressibles.

Le système considéré (équations 1.1-1.2) décrit l'évolution de deux champs de vorticité $\omega_+$ et $\omega_-$ (liés aux variables d'Elsässer $u_\pm = v \pm b$) sous l'effet d'un transport et d'un étirement de tourbillon, avec des paramètres $a, p, q$.

• Cas limite : En l'absence de champ magnétique ($b=0$, donc $\omega_+ = \omega_-$), le système se réduit au modèle de De Gregorio ($a=1$).
• Objectif : L'analyse se concentre sur le comportement des solutions près du premier état excité (ou premier état stationnaire non trivial) sur le tore, défini par :
  $$\Omega_2(\theta) = -\sin(2\theta), \quad U_2(\theta) = \frac{1}{2}\sin(2\theta).$$
  Ce travail généralise les résultats récents de Guo et Jiu (2025) obtenus pour le modèle de De Gregorio au cas MHD couplé.

2. Méthodologie

Les auteurs adoptent une approche analytique rigoureuse combinant analyse spectrale, méthodes de Galerkin et estimées d'énergie dans des espaces de Sobolev pondérés.

A. Cadre Fonctionnel

L'analyse est menée dans l'espace de Hilbert $H_2$, inspiré des travaux de Lai et al. (2020) et Guo-Jiu (2025) :
$$H_2 := \left\{ \eta \in H^1(\mathbb{T}) : \eta \text{ est impair et } \int_{\mathbb{T}} \frac{|\partial_\theta \eta|^2}{\sin^2(\theta)} d\theta < \infty \right\}.$$
La norme est définie par une pondération singulière $\rho_2 = (4\pi \sin^2 \theta)^{-1}$. Une base orthonormée $\{e_{2,k}\}_{k \ge 1}$ est utilisée pour décomposer les perturbations, où :
$$e_{2,k} = \frac{\sin((k+2)\theta)}{k+2} - \frac{\sin(k\theta)}{k}.$$

B. Linéarisation et Opérateurs

Les auteurs étudient les perturbations $\eta_\pm$ autour de l'état excité $\Omega_2$. Le système linéarisé (2.7) s'écrit :
$$\partial_t \eta_+ = L_+ \eta_+, \quad \partial_t \eta_- = L_- \eta_-,$$
où $L_+$ et $L_-$ sont des opérateurs linéaires non locaux.

• Observation clé : La structure de l'opérateur $L_+$ est identique à celle du modèle de De Gregorio. C'est cet opérateur qui porte la source d'instabilité.
• L'opérateur $L_-$, en revanche, possède des propriétés de dissipation (décroissance exponentielle).

C. Stratégie de Preuve

1. Bien-posé global (Linéaire) : Utilisation de la méthode de Galerkin avec la base $\{e_{2,k}\}$ pour établir l'existence et l'unicité de solutions classiques pour le système linéarisé.
2. Analyse d'instabilité :
   • Développement des coefficients de Fourier $\eta_{2,k}(t)$.
   • Dérivation d'une équation différentielle ordinaire (EDO) du second ordre pour la norme $L^2$ pondérée des coefficients.
   • Analyse d'une forme quadratique positive définie issue de cette EDO pour prouver une croissance exponentielle pour certaines classes de données initiales.
3. Instabilité Non Linéaire : Utilisation d'un argument de contradiction (inspiré de [JJN13]). On construit une famille de solutions non linéaires à partir de données initiales petites, on les rescale et on montre que, par compacité, elles convergent vers la solution linéaire instable, contredisant ainsi la stabilité supposée.
4. Stabilité Non Linéaire : Pour une sous-classe spécifique de données initiales (coefficients pairs uniquement dans le développement de Fourier), les termes instables disparaissent, permettant de prouver la stabilité globale via un argument de continuité et des estimées d'énergie fines.

3. Résultats Principaux

Les résultats sont synthétisés dans les théorèmes 3.1 à 3.10 :

1. Bien-posé global du système linéarisé (Théorèmes 3.1 & 3.2) :
   Pour des données initiales dans $H_2$ (avec régularité suffisante $m > 3$), il existe une unique solution classique globale pour le système linéarisé.

2. Instabilité Linéaire (Théorème 3.3) :
   Pour une large classe de données initiales (notamment celles ayant une composante non nulle sur le premier mode impair $e_{2,1}$), la solution du système linéarisé croît exponentiellement :
   $$\|\eta_+(t)\|_{H_2} + \|\eta_-(t)\|_{H_2} \ge C e^{\lambda t}.$$
   Le mécanisme d'instabilité est identifié comme étant l'opérateur $L_+$, dont la structure spectrale permet une croissance similaire à celle du modèle de De Gregorio.

3. Instabilité Non Linéaire (Théorème 3.7) :
   Le système non linéaire est instable au sens de Lipschitz. Il existe des données initiales arbitrairement petites dans $H_2$ pour lesquelles la solution s'éloigne exponentiellement de l'état d'équilibre, confirmant que l'instabilité linéaire persiste dans le régime non linéaire.

4. Stabilité Globale et Décroissance (Théorème 3.10) :
   Pour une sous-classe spécifique de données initiales (où $\eta_+$ ne contient que des modes pairs dans la base $\{e_{2,k}\}$), le système est globalement bien posé et stable. La solution converge exponentiellement vers l'état d'équilibre :
   $$I_0(t) \le C e^{-t/2} I_0(0).$$
   Cela démontre que l'instabilité n'est pas universelle mais dépend de la structure spectrale de la perturbation.

4. Contributions Clés

• Généralisation au MHD 1D : Extension des résultats de stabilité/instabilité du modèle de De Gregorio (Euler) au modèle MHD couplé, en traitant la complexité supplémentaire introduite par le couplage des champs magnétiques et fluides.
• Identification du mécanisme d'instabilité : Démonstration que l'instabilité provient spécifiquement de la composante $\eta_+$ et de son opérateur linéarisé, qui conserve la structure critique du modèle de De Gregorio, tandis que la composante $\eta_-$ est dissipative.
• Analyse fine des états excités : Contrairement aux études sur l'état fondamental (ground state), l'analyse des états excités ($\kappa \ge 2$) révèle une sensibilité critique aux modes de Fourier spécifiques (pairs vs impairs).
• Traitement des paramètres $q \neq 0$ : L'article aborde la difficulté technique liée aux termes non locaux supplémentaires lorsque $q \neq 0$ (via l'opérateur $Q$), en prouvant sa bornitude dans l'espace pondéré $H_2$.

5. Signification et Impact

Ce travail est significatif pour plusieurs raisons :

• Compréhension de la singularité : Les modèles 1D MHD servent de banc d'essai pour comprendre la formation de singularités dans les équations MHD 3D. La démonstration d'instabilité pour les états excités suggère que la turbulence ou la formation de singularités pourrait être déclenchée par des perturbations spécifiques de ces états.
• Nuance entre stabilité et instabilité : L'article met en lumière la dualité du système : il est instable pour une large classe de perturbations, mais stable pour une autre classe (modes pairs). Cela souligne l'importance de la structure spectrale initiale dans la dynamique des fluides magnétisés.
• Méthodologie : La combinaison de l'analyse spectrale fine (formes quadratiques, EDO du second ordre) avec des estimées non linéaires dans des espaces pondérés fournit un cadre robuste pour l'étude de modèles MHD non locaux.

En conclusion, l'article établit que le modèle MHD 1D est fondamentalement instable autour du premier état excité pour des perturbations génériques, généralisant ainsi les phénomènes d'instabilité connus pour les équations d'Euler, tout en identifiant des sous-espaces de stabilité pour des données structurées.