Effet de pincement dans le plasma : comment un tokamak attire lui-même les particules vers le centre
Contexte
Ce travail se situe en marge des recherches sur la turbulence et le transport dans le plasma de bord des tokamaks. La plupart des travaux voisins étudient le transport turbulent de chaleur et de particules vers l'extérieur, en fixant la densité à la frontière interne ou en utilisant des modèles simplifiés. Ici, pour la première fois dans une modélisation fluide globale avec un profil de densité auto-cohérent, on montre que la pince se produit naturellement, sans hypothèses artificielles. Cela relie ce travail à des questions anciennes de la physique du confinement, comme la formation du piedestal de densité dans le mode H.
Méthodes
On utilise le code GDB, un modèle global tridimensionnel électromagnétique basé sur les équations de Braginskii réduites par dérive. On simule un plasma avec des paramètres proches de ceux d'Alcator C-Mod : grand rayon de 68,5 cm, petit rayon de 22 cm, champ de 5,3 T. La zone de calcul couvre toute la couche toroïdale annulaire, de 5,3 cm à l'intérieur de la dernière surface fermée jusqu'à 3,3 cm dans la couche de bord (scrape-off layer). Le profil de densité initial est plat, la source de particules est localisée près de la séparatrice, et les sources de chaleur se trouvent à la frontière interne. La simulation dure 8,4 ms de temps de transport, ce qui permet au profil de densité d'atteindre un quasi-équilibre. Pour analyser les mécanismes, on effectue des calculs supplémentaires en désactivant la conductivité thermique électronique et les flux thermiques transversaux.
Résultats
La simulation montre qu'à partir d'un profil initial plat, un profil avec un pic de densité au centre de la zone de calcul se forme en quelques millisecondes. Le flux de particules vers l'intérieur existe pendant toute l'évolution, mais sa nature change. Au début, lorsque le gradient de densité est positif (la densité augmente vers l'extérieur), le flux turbulent domine, lié à la thermodiffusion électronique. Une analyse linéaire d'un modèle simplifié montre que lorsque le rapport des gradients de densité et de température est négatif (η<0), les ondes de dérive créent un déphasage entre les fluctuations de densité et de potentiel, ce qui donne un flux négatif (vers l'intérieur). À un stade ultérieur, lorsque le profil de densité devient plat, la contribution turbulente s'affaiblit, mais un flux E×B d'équilibre apparaît, dû à l'asymétrie poloïdale de la densité et du potentiel. Cette asymétrie est engendrée par le flux de chaleur de Pfirsch-Schlüter, mais la théorie néoclassique classique prédit un flux total nul en raison de la réponse boltzmannienne des électrons. Dans le modèle GDB, on tient compte des effets non idéaux (résistance, inertie des électrons, induction électromagnétique) qui créent un petit déphasage poloïdal δs ≈ −0.08π entre la densité et le potentiel. Ce déphasage conduit à un flux fini vers l'intérieur, ce qui explique quantitativement l'accumulation de densité. Les calculs de contrôle avec désactivation des flux thermiques transversaux montrent que le pic de densité se dissipe, confirmant le rôle de l'asymétrie poloïdale.
Implications
Ce travail se distingue des recherches voisines en ce qu'il n'impose pas le profil de densité de l'extérieur, mais lui permet d'évoluer de manière auto-cohérente. C'est la première démonstration dans un modèle fluide global du bord que la pince de particules se produit naturellement et conduit à un pic de densité. On montre que le tableau néoclassique classique est incomplet : la prise en compte des non-idéalités dans la loi d'Ohm ouvre un canal de transport supplémentaire. Cela change la compréhension du bilan des particules dans le plasma de bord et fournit une base pour des prédictions plus précises des profils de densité.
Développements futurs
Le sujet pourrait se développer vers la prise en compte d'une géométrie plus réaliste (par exemple, l'élongation, le divertor) et des effets cinétiques pour les plasmas chauds. Il est important de vérifier si la pince E×B non idéale découverte persiste dans les modèles gyrocinétiques et à plus faible collisionalité. Il serait également intéressant d'étudier le rôle de l'ionisation des neutres asymétrique poloïdalement, qui dans les tokamaks réels pourrait renforcer ou affaiblir la pince.
Impact
Les résultats influenceront la modélisation du confinement du plasma dans les tokamaks et les stellarators, en particulier la prédiction des profils de densité et la transition vers les régimes de confinement amélioré. Ils sont également utiles pour l'interprétation des expériences d'injection de gaz et d'injection de pastilles.
Prochaines étapes
Les prochaines étapes comprennent la comparaison avec les profils de densité expérimentaux sur Alcator C-Mod et d'autres installations, ainsi que l'extension du modèle pour inclure les corrections cinétiques et la géométrie réelle.
Problèmes ouverts clés
Ce travail est directement lié au problème non résolu de la formation du piedestal de densité dans le mode H et à la question de savoir quels mécanismes déterminent le pic de densité dans les différents régimes de confinement. Il éclaire également une énigme ancienne : pourquoi dans le mode I il n'y a pas de barrière pour les particules alors qu'il y en a une pour la chaleur.
Dans la simulation, les particules « nagent à contre-courant » : elles se déplacent du bord vers le centre, alors que le gradient de densité est initialement dirigé vers l'extérieur. Ce transport contre-intuitif est possible grâce à la corrélation entre les fluctuations de densité et le champ électrique, et non à la diffusion ordinaire.
L'idée d'une compression spontanée du plasma vers le centre fait écho au concept de « pince » dans la science-fiction, par exemple dans le roman de Larry Niven « L'Anneau-Monde », où des champs magnétiques confinent le plasma dans un réacteur.
Chiffres clés
- durée de la simulation: 8.4 ms
- déphasage δs: ≈ −0.08π
- pic de densité: 1.6×1020 m^-3
- vitesse de croissance de la densité: en 4 ms le profil devient plat, en 8 ms il est piqué
- rapport des gradients η: de −2 à +2