Mercier--Cotsaftis and Grad--Shafranov equations for anisotropic plasma

Cette brève revue examine les aspects historiques de la généralisation de l'équation de Grad-Shafranov au cas d'un plasma anisotrope.

L'essentiel

🌌 L'histoire d'une équation qui a perdu son nom (et qu'il faut retrouver)

Imaginez que vous êtes un architecte qui construit des maisons. Pour construire une maison solide, vous avez besoin d'un plan de base, une règle d'or qui dit comment les murs et le toit doivent s'assembler pour ne pas s'effondrer.

En physique des plasmas (ce gaz super chaud et magnétique qu'on utilise pour essayer de créer l'énergie des étoiles, comme dans les réacteurs à fusion), il existe un "plan d'architecte" très célèbre appelé l'équation de Grad-Shafranov. C'est la règle qui permet de calculer comment le plasma se comporte dans un réacteur en forme de beignet (un tore).

Mais il y a un problème : ce plan fonctionne parfaitement pour un gaz "normal" (isotrope), mais il devient flou quand le gaz est "tordu" ou "anisotrope" (c'est-à-dire quand la pression du gaz n'est pas la même dans toutes les directions, un peu comme si vous poussiez sur un ballon qui résiste plus fort d'un côté que de l'autre).

C'est ici que l'auteur, Igor Kotelnikov, nous raconte une petite histoire de science un peu confuse, un peu comme un jeu de "téléphone arabe".

1. Le mystère du "diamagnétisme"

Récemment, des chercheurs ont essayé de simuler un piège à plasma très spécial (un "piège diamagnétique") où le plasma chasse presque totalement le champ magnétique de son intérieur. Pour le faire, ils ont utilisé une équation qu'ils ont appelée "l'équation de Grad-Shafranov".

Mais l'auteur dit : "Attendez une minute !".
Il compare cela à appeler un avion un "bateau" juste parce qu'ils ont tous les deux des ailes. L'équation utilisée par ces chercheurs est mathématiquement différente de la vraie équation de Grad-Shafranov. C'est comme si le plan de construction avait changé de nature : il est devenu une équation "intégro-différentielle" (un mot compliqué pour dire qu'elle dépend du passé et du futur en même temps, pas juste de l'instant présent).

2. Le casse-tête des noms

L'auteur a alors décidé de faire un peu de détective. Il a regardé dans les livres et les articles scientifiques pour voir comment on nomme cette équation pour les plasmas "tordus" (anisotropes).

Le résultat est fouillis ! On trouve quatre ou cinq noms différents pour la même chose :

• L'équation de Grad-Shafranov généralisée.
• L'équation de Grad-Shafranov modifiée.
• L'équation de Grad-Shafranov anisotrope.
• L'équation de Mercier-Cotsaftis.

C'est comme si on appelait la même voiture "La Ford", "L'Automobile", "Le Véhicule" ou "La Taxis" selon le pays où vous vous trouvez. Cela crée de la confusion et fait perdre du temps aux scientifiques.

3. L'histoire oubliée de Mercier et Cotsaftis

Voici le cœur de l'histoire. En 1961, deux physiciens français, Claude Mercier et Michel Cotsaftis, ont écrit la première version correcte de cette équation pour les plasmas tordus. C'était leur "invention".

Mais quelques années plus tard, un grand physicien américain, Harold Grad, a redécouvert cette équation en 1967. Il l'a expliquée plus clairement, mais il n'a pas mentionné le nom de Mercier et Cotsaftis. Comme Grad était très célèbre et très respecté, tout le monde a commencé à dire : "Ah, c'est l'équation de Grad !" (ou plutôt "l'équation de Grad-Shafranov généralisée").

L'auteur compare cela à l'histoire de l'astronomie :

• Pendant longtemps, on a appelé la loi de l'expansion de l'univers la "Loi de Hubble".
• Mais en réalité, un prêtre belge, Georges Lemaître, l'avait découverte et prouvée mathématiquement deux ans avant Hubble !
• Ce n'est qu'en 2018 que l'Union Astronomique Internationale a décidé de corriger l'histoire et d'appeler cela la "Loi de Hubble-Lemaître".

4. La conclusion de l'auteur

Igor Kotelnikov nous dit : "Il est temps de corriger l'histoire !"

Il suggère que pour éviter la confusion, nous devrions arrêter d'appeler cette équation "Généralisée" ou "Modifiée" (ce qui est ennuyeux et vague) et lui donner le nom de ses véritables inventeurs français : l'équation de Mercier-Cotsaftis.

C'est une façon de rendre justice à ceux qui ont posé les fondations, tout en simplifiant la vie de tous ceux qui utilisent cette équation aujourd'hui.

En résumé, avec une analogie culinaire 🍲

Imaginez que vous avez une recette de gâteau (l'équation de base).

1. Grad et Shafranov ont écrit la recette originale pour un gâteau simple.
2. Mercier et Cotsaftis ont inventé la recette pour un gâteau avec des couches de fruits exotiques (le plasma anisotrope).
3. Plus tard, un grand chef célèbre (Grad à nouveau) a repris la recette des fruits, l'a expliquée très bien, mais a oublié de citer les inventeurs des fruits.
4. Aujourd'hui, tout le monde appelle ça "La recette de Grad aux fruits", alors que c'est en fait "La recette de Mercier-Cotsaftis".

L'auteur dit : "Arrêtons de mélanger les noms. Appelons-le par son vrai nom, c'est plus juste et plus clair pour tout le monde !"

Résumé technique

1. Problématique

L'article aborde une confusion terminologique et conceptuelle persistante dans la littérature sur la physique des plasmas, concernant l'équation d'équilibre pour les plasmas à pression anisotrope.

• Le contexte : Les équations de magnétohydrodynamique (MHD) classiques, notamment l'équation de Grad–Shafranov (GSE), supposent une pression scalaire isotrope. Cependant, dans les plasmas modernes soumis à un chauffage intense (injection de neutres, ondes radiofréquences), la pression devient anisotrope ($p_\parallel \neq p_\perp$), nécessitant l'utilisation d'un tenseur de pression.
• La confusion : Des travaux récents (notamment sur les pièges diamagnétiques) utilisent des équations d'équilibre pour des plasmas anisotropes en les appelant simplement « équation de Grad–Shafranov ». L'auteur souligne que cette désignation est inexacte car l'équation utilisée dans ces cas spécifiques devient une équation intégro-différentielle (non locale) plutôt qu'une équation aux dérivées partielles (EDP) non linéaire classique.
• L'objectif : Clarifier l'historique, la terminologie et la structure mathématique de ces équations, en mettant en lumière le rôle souvent sous-estimé de Claude Mercier et Michel Cotsaftis par rapport à Harold Grad et V. Shafranov.

2. Méthodologie

L'auteur adopte une approche historique et analytique :

• Revue de la littérature : Analyse statistique du nombre de publications mentionnant différentes variantes de l'équation (GSE généralisée, modifiée, anisotrope, ou équation de Mercier–Cotsaftis) dans des bases de données et des journaux majeurs (PoP, JPP, NF, PPCF).
• Analyse mathématique comparative : Comparaison rigoureuse des formes canoniques et généralisées des équations d'équilibre. L'auteur examine la dépendance des termes sources (courant, pression) par rapport au flux magnétique $\psi$ et au champ magnétique $B$.
• Examen historique : Investigation des publications originales de Mercier et Cotsaftis (1961) et de Grad (1967) pour retracer la genèse des équations et les hypothèses sous-jacentes (notamment la nature des fonctions de surface).

3. Contributions Clés

A. Clarification Terminologique et Historique

• L'article identifie que l'équation d'équilibre pour un plasma anisotrope porte plusieurs noms : Équation de Grad–Shafranov Généralisée (GGSE), Modifiée (MGSE), Anisotrope (AGSE), ou Équation de Mercier–Cotsaftis (MCE).
• Il démontre que l'équation a été formulée par Mercier et Cotsaftis en 1961, puis réintroduite par Harold Grad en 1967 sans citer les précurseurs.
• Les statistiques montrent que le terme « Grad–Shafranov » (généralisé) est beaucoup plus utilisé (plus de 10 000 références estimées) que « Mercier–Cotsaftis » (quelques centaines), bien que ce dernier soit historiquement prioritaire pour le cas anisotrope.

B. Distinction Mathématique Fondamentale

L'auteur distingue deux types d'équations pour les plasmas anisotropes :

1. L'Équation de Grad–Shafranov Généralisée (GGSE / MCE) :

   • Forme : $\nabla \cdot \left( \frac{\sigma}{R^2} \nabla \psi \right) = -\mu_0 \frac{\partial p_\parallel}{\partial \psi} - \mu_0^2 \frac{F}{\sigma R^2} \frac{dF}{d\psi}$ (avec $\sigma = 1 - \mu_0(p_\parallel - p_\perp)/B^2$).
   • C'est une EDP non linéaire elliptique.
   • Elle repose sur l'hypothèse que le flux de courant modifié $F = \sigma i_p$ est une fonction de surface (dépendant uniquement de $\psi$).
   • Elle permet d'analyser les instabilités « miroir » et « firehose » via la condition $\sigma > 0$.

2. L'Équation des Pièges Diamagnétiques (Cas non local) :

   • Dans les pièges diamagnétiques où le champ est expulsé, le courant diamagnétique $J_\theta$ dépend de manière non locale du flux $\psi$ (liée à la taille des orbites des ions).
   • L'équation devient intégro-différentielle (le terme de droite dépend des dérivées ou de l'histoire du flux, pas seulement de $\psi$).
   • L'auteur critique l'appellation « Grad–Shafranov » pour ce cas, suggérant qu'il s'agit d'une classe d'équations différente, potentiellement nommée « équation de Beklemishev–Khristo ».

C. Application aux Pièges à Miroir

• L'article montre que pour les pièges à miroir axiaux (GDT, GDMT), l'équation se simplifie considérablement (pas de courant poloidal $i_p=0$), réduisant la MCE à une forme plus simple.
• Il note que l'approximation paraxiale a longtemps suffi pour ces systèmes, retardant l'application des équations d'équilibre complètes, sauf pour les configurations de pièges diamagnétiques où l'approximation échoue.

4. Résultats

• Validation de la priorité historique : Confirmation que Mercier et Cotsaftis ont établi les bases mathématiques de l'équilibre anisotrope avant Grad.
• Propriétés mathématiques : Démonstration que la condition de stabilité (absence d'instabilité firehose) est intrinsèquement liée à la nature elliptique de l'équation de Mercier–Cotsaftis ($\sigma > 0$).
• Statistiques d'usage : Mise en évidence du décalage entre l'usage courant (prédominance du nom Grad–Shafranov) et la réalité historique et technique (l'équation de Mercier–Cotsaftis est le terme le plus précis pour le cas anisotrope historique).
• Critique des modèles récents : Identification d'une erreur de catégorisation dans les simulations récentes de pièges diamagnétiques qui traitent des équations intégro-différentielles comme des EDP classiques.

5. Signification et Conclusion

• Justification historique : L'auteur plaide pour une reconnaissance accrue du nom « Équation de Mercier–Cotsaftis », comparant la situation à celle de la loi de Hubble-Lemaître, où la justice historique a été rétablie tardivement.
• Utilité pratique : L'utilisation du terme correct est cruciale pour éviter la confusion entre les équations différentielles locales (valables pour la plupart des tokamaks et pièges à miroir) et les équations intégro-différentielles (nécessaires pour les pièges diamagnétiques à fort courant).
• Impact scientifique : En clarifiant ces termes, l'article vise à améliorer la compréhension des problèmes de stabilité locale (critère de Mercier) et à guider les développeurs de codes numériques vers les formulations mathématiques appropriées selon le régime de plasma étudié.

En résumé, cet article est une mise au point nécessaire qui réhabilite le travail de Mercier et Cotsaftis, tout en fournissant une distinction mathématique rigoureuse entre les équations d'équilibre classiques pour plasmas anisotropes et les nouvelles équations nécessaires pour décrire les états diamagnétiques extrêmes.