Self-consistent modelling and qualitative comparison of mildly relativistic runaway electron dynamics with a closed flux surface formation model during tokamak startup

Cet article présente un modèle cinétique réduit et auto-cohérent pour les électrons de fuite légèrement relativistes, intégré dans le code prédictif DYON (DYON-RE), qui reproduit avec succès les observations du démarrage ohmique de KSTAR et offre un cadre pour concevoir des scénarios de démarrage sans électrons de fuite pour les futurs dispositifs tels que ITER.

L'essentiel

Imaginez un tokamak (un réacteur à fusion nucléaire en forme de beignet) tentant de démarrer comme un moteur de voiture. Il doit transformer un vide froid et vide en une boule chaude et tourbillonnante de plasma. Mais il y a un effet secondaire dangereux : parfois, quelques électrons sont projetés si violemment qu'ils deviennent des particules « fuyardes », se déplaçant à une vitesse proche de celle de la lumière. Si trop de ces électrons fuyards se forment, ils peuvent agir comme un faisceau laser haute puissance, faisant fondre les parois du réacteur et mettant fin à l'expérience.

Ce document traite de la construction d'une meilleure carte pour prédire quand et comment ces électrons fuyards apparaissent lors de cette phase délicate de « démarrage ». Les auteurs, travaillant avec le réacteur à fusion KSTAR en Corée du Sud, ont développé un nouveau modèle appelé DYON-RE.

Voici une décomposition de leur travail utilisant des analogies simples :

1. Le Problème : L'Erreur de la « Vitesse de la Lumière »

Par le passé, les scientifiques tentaient de prédire ces électrons fuyards en supposant qu'ils voyageaient déjà à la vitesse de la lumière dès le moment où ils commençaient à fuir.

• L'Analogie : Imaginez un pilote de voiture de course. Les anciens modèles supposaient que dès que le pilote appuyait sur l'accélérateur, la voiture atteignait instantanément 320 km/h.
• La Réalité : Dans les premières étapes du démarrage du réacteur, les électrons sont « modérément relativistes ». Ils sont rapides, mais ils n'ont pas encore atteint leur vitesse maximale. Ils ressemblent davantage à une voiture accélérant de 0 à 100 km/h.
• La Solution : Les auteurs ont créé un nouveau modèle qui prend en compte cette phase d'accélération. En réalisant que les électrons ne sont pas instantanément à vitesse maximale, leur modèle cesse de surestimer l'ampleur du courant dangereux généré par ces électrons. C'est comme réaliser que la voiture ne va qu'à 65 km/h, et non 320 km/h, ce qui modifie l'ampleur des dégâts qu'elle pourrait causer.

2. Le Défi : Le Piège « Ouvert vs Fermé »

Pendant le démarrage, les champs magnétiques qui maintiennent le plasma en place changent de forme.

• L'Analogie : Pensez au champ magnétique comme à une clôture.
  • Champ Ouvert : Au tout début, la clôture présente des brèches. Si un électron fuyard tente de s'enfuir, il heurte une brèche et s'échappe (comme un chien sortant par une porte ouverte).
  • Champ Fermé : À mesure que le réacteur chauffe, la clôture se referme pour former un cercle parfait (une surface de flux fermée). Désormais, l'électron fuyard est piégé à l'intérieur d'une cage et ne peut pas s'échapper.
• L'Ancienne Méthode : Les modèles précédents traitaient la clôture comme étant soit toujours ouverte, soit toujours fermée, ou utilisaient une moyenne floue des deux.
• La Nouvelle Méthode : Le modèle DYON-RE agit comme un système de sécurité intelligent qui sait exactement quand la clôture se referme. Il suit les électrons séparément : ceux qui courent dans le « champ ouvert » (où ils se perdent rapidement) et ceux piégés dans le « champ fermé » (où ils s'accumulent). Ceci est crucial car le moment où la clôture se referme est celui où le danger commence vraiment à s'accumuler.

3. L'Expérience : Observer le « Thermomètre à Rayonnement »

L'équipe a testé son nouveau modèle contre des données réelles provenant du réacteur KSTAR. Ils ne pouvaient pas voir directement les électrons fuyards, alors ils ont cherché des indices.

• L'Analogie : Imaginez essayer de déterminer si une pièce est remplie de gens en écoutant le niveau de bruit.
• L'Indice : Ils ont utilisé un outil appelé Émission Cyclotronique Électronique (ECE), qui agit comme un « thermomètre à rayonnement ». Lorsque les électrons fuyards sont excités, ils émettent un rayonnement qui fait afficher une température très élevée sur ce thermomètre.
• Le Résultat : Ils ont examiné deux tentatives de démarrage différentes :
  1. Le Tir « Riche en Fuyards » : Le réacteur contenait beaucoup d'électrons fuyards. Le modèle l'avait prédit, et le « thermomètre » a montré une énorme pointe de température, exactement comme le modèle l'indiquait.
  2. Le Tir « Pauvre en Fuyards » : Le réacteur contenait très peu d'électrons fuyards. Le modèle l'avait également prédit, et le thermomètre est resté relativement calme, avec seulement de petits « pics » rythmiques (comme un battement de cœur) au lieu d'une énorme pointe.

4. L'Ingrédient Secret : Les Murs

L'une des découvertes clés de l'article est que les parois du réacteur jouent un rôle plus important que précédemment supposé.

• L'Analogie : Imaginez essayer de remplir un seau avec un tuyau (injection de gaz). Si le seau présente une fuite cachée (les parois absorbant le gaz), vous devez ouvrir le tuyau plus fort pour obtenir la même quantité d'eau.
• La Découverte : Les chercheurs ont constaté que même lorsqu'ils utilisaient exactement les mêmes paramètres de gaz, le réacteur se comportait différemment parce que les « parois » agissaient différemment (absorbant ou libérant le gaz à des rythmes différents). Pour que leur modèle fonctionne, ils ont dû ajuster ces conditions de paroi. Sans prendre en compte les parois, le modèle ne pouvait pas prédire correctement la densité électronique.

Résumé

L'article ne prétend pas avoir résolu le problème des électrons fuyards pour toujours, mais il a construit un simulateur meilleur et plus réaliste.

• Il cesse de supposer que les électrons sont instantanément à vitesse maximale.
• Il suit exactement quand la « clôture » magnétique se referme pour les piéger.
• Il prédit avec succès les « pointes de température » observées dans les expériences réelles.

Cela offre aux scientifiques un outil plus fiable pour concevoir les futurs réacteurs (comme ITER) afin qu'ils puissent démarrer en toute sécurité sans créer accidentellement un faisceau d'électrons susceptible d'endommager la machine.

Résumé technique

1. Énoncé du problème

Les électrons de fuite (RE) constituent une préoccupation critique dans les tokamaks, traditionnellement associée aux disruptions. Cependant, les RE peuvent également se former durant la phase de démarrage (spécifiquement lors de la pénétration thermique et de la montée précoce du courant), entravant potentiellement l'allumage du plasma ou endommageant l'appareil.

Les modèles existants pour les RE de démarrage souffrent de plusieurs limitations :

• Hypothèse de vitesse : De nombreux modèles supposent que les électrons de fuite voyagent à la vitesse de la lumière ($v=c$). Cela surestime considérablement la densité de courant de fuite lorsque les électrons se trouvent dans le régime « légèrement relativiste » (accélérés mais pas encore ultra-relativistes).
• Configuration magnétique : Les modèles échouent souvent à distinguer entre les lignes de champ ouvertes (début du démarrage) et les surfaces de flux fermées (formées ultérieurement). La transition entre ces configurations modifie drastiquement le confinement des RE, ce que les modèles précédents à population unique ne pouvaient pas capturer avec précision.
• Validation : Il manque une validation auto-cohérente par rapport aux données expérimentales, en particulier concernant le comportement qualitatif des RE (par exemple, les instabilités cinétiques) plutôt que simplement le courant total.

2. Méthodologie

Les auteurs ont développé un nouveau cadre prédictif intégré au code d'allumage du plasma DYON, désigné sous le nom de DYON-RE.

A. Le modèle d'électrons de fuite légèrement relativistes

• Approche cinétique réduite : Au lieu de supposer $v=c$, le modèle calcule la vitesse moyenne ($\beta = v/c$) de la population de RE.
• Réduction de multi-fluide à mono-fluide : Les auteurs ont d'abord formulé un modèle multi-fluide suivant les RE générés à différents pas de temps. Ils ont ensuite réduit cela à une forme mono-fluide en supposant que la vitesse moyenne caractérise l'évolution globale du courant. Cela préserve la précision tout en maintenant l'efficacité numérique.
• Définition de la frontière : Le modèle définit la frontière de fuite non seulement à l'impulsion critique ($p_c$) mais à une frontière basée sur une région ($p_V \approx 1.3 p_c$), ce qui s'aligne mieux avec les simulations cinétiques du taux de génération Dreicer.
• Champs effectifs : Le modèle intègre un champ électrique critique effectif et un champ Dreicer effectif pour tenir compte des collisions inélastiques (ionisation) avec les particules neutres et les impuretés, qui sont abondantes durant le démarrage.

B. Configuration magnétique double (modèle CFSF)

Le modèle sépare explicitement les populations de RE en deux groupes basés sur la topologie magnétique :

1. Lignes de champ ouvertes ($n_{op}^{RE}$) : Les RE sont perdus principalement par écoulement parallèle le long des lignes de champ.
2. Surfaces de flux fermées ($n_{cl}^{RE}$) : Les RE sont confinés radialement, les pertes étant gouvernées par la diffusion Rechester-Rosenbluth.

• Transition : Le modèle suit dynamiquement la formation de Surfaces de Flux Fermées (CFSF). À mesure que le courant plasma augmente et qu'une région locale fermée se forme (même avant l'existence d'une surface fermée globale), les RE générés dans cette région bénéficient d'un confinement nettement amélioré.

C. Couplage auto-cohérent

Le modèle de RE est couplé de manière auto-cohérente avec :

• Circuit électromagnétique : Le courant de fuite est soustrait du courant plasma total dans l'équation du circuit pour calculer avec précision le champ électrique résistif.
• Allumage du plasma : Le code simule le processus complet de pénétration thermique, incluant l'ionisation du gaz, le recyclage et l'afflux d'impuretés.

D. Outils de validation

Pour comparer avec les expériences KSTAR, les auteurs ont utilisé :

• KIAT (Kinetic Instability Analysis Tool) : Pour prédire le début des instabilités cinétiques (ondes siffleuses) entraînées par les RE.
• SYNO (Synthetic Non-thermal ECE Reconstruction) : Pour reconstruire la température radiative ($T_{rad}$) attendue des RE et la comparer avec les diagnostics d'émission cyclotronique électronique (ECE).

3. Contributions clés

1. Correction légèrement relativiste : Démonstration que l'hypothèse $v=c$ surestime les courants de RE de démarrage. Le nouveau modèle prend en compte la distribution de vitesse réelle, fournissant une prédiction de densité de courant plus précise.
2. Transport à double topologie : Introduction d'un modèle distinguant les configurations magnétiques ouvertes et fermées. Cela capture l'effet de confinement des « graines » où les RE générés dans des surfaces de flux fermées locales survivent plus longtemps, impactant significativement la population totale de RE durant la phase de transition.
3. Cadre de validation qualitative : Établissement d'une méthode pour valider les modèles de RE non seulement par l'amplitude du courant (difficile à mesurer directement dans KSTAR) mais en comparant les caractéristiques du rayonnement ECE (motifs de dents de scie, pics de température) et les seuils d'instabilité cinétique.
4. Sensibilité aux conditions de paroi : Mise en évidence du fait que la prédiction précise des RE de démarrage nécessite une modélisation précise des conditions de paroi (coefficients de recyclage) et de la pression de gaz de remplissage préalable, car celles-ci dictent l'évolution de la densité électronique qui entraîne la génération de RE.

4. Résultats

Le modèle a été testé contre deux décharges ohmiques KSTAR :

• Tir #26031 (Riche en RE) :
  • Observation : Signaux HXR élevés et pic rapide de température radiative ECE ($>12$ keV), indiquant une forte instabilité cinétique.
  • Prédiction DYON-RE : A reproduit avec succès le courant et la densité du plasma. Il a prédit un courant de RE initial élevé (~24 kA) dû au confinement des graines dans des surfaces de flux fermées locales. Le modèle a montré que la densité de RE était suffisante pour entraîner de fortes instabilités cinétiques ($\gamma_{kin} \gg \gamma_{damp}$), correspondant au comportement ECE observé.
• Tir #27340 (Pauvre en RE) :
  • Observation : Signaux HXR plus faibles et éclatements « de type dents de scie » dans la température ECE, suggérant une instabilité marginale.
  • Prédiction DYON-RE : A prédit un courant de RE plus faible (~16 kA). Le modèle a montré que le taux de croissance cinétique était légèrement supérieur au taux d'amortissement, cohérent avec le comportement d'éclatement observé.
• Comparaison : Sans le modèle à double configuration (en utilisant une simple interpolation), le courant de RE initial prédit était d'environ 1 kA seulement, échouant à expliquer les observations expérimentales. Le nouveau modèle a expliqué avec succès les différences entre les deux tirs comme résultant de conditions de paroi variables (recyclage) affectant l'évolution de la densité, plutôt que simplement de la pression de remplissage préalable.

5. Importance

• Capacité prédictive pour les futurs dispositifs : Le code DYON-RE fournit un outil robuste pour concevoir des scénarios de démarrage sans fuite pour les futurs réacteurs comme ITER et CPD (Compact Fusion Reactor).
• Sécurité et exploitation : En comprenant que les RE peuvent se former lors du démarrage et persister en raison des surfaces de flux fermées locales, les opérateurs peuvent mieux atténuer les risques de tremblement des aimants ou d'endommagement des composants.
• Insight physique : Ce travail clarifie que la nature « légèrement relativiste » des RE de démarrage et la topologie magnétique spécifique durant la transition CFSF sont des facteurs critiques souvent négligés dans les modèles simplifiés précédents.
• Orientation expérimentale : L'étude souligne que le contrôle des conditions de paroi (recyclage) est aussi crucial que l'injection de gaz pour la gestion des RE de démarrage, suggérant que les simulations prédictives doivent prendre en compte des coefficients de recyclage variables dans le temps pour être fiables.

En conclusion, cet article présente une avancée significative en physique du démarrage des tokamaks en intégrant un modèle de RE physiquement rigoureux et légèrement relativiste avec des changements de topologie magnétique dynamique, validé qualitativement par des données expérimentales KSTAR.