Improved n=1 Empirical Error Field Penetration Threshold Scaling with Ohmic and L-Mode Conventional Tokamak Plasma Discharges

Ce papier présente une mise à l'échelle améliorée du seuil de pénétration du champ d'erreur empirique n=1, dérivée d'une base de données élargie de décharges tokamak ohmiques et en mode L conventionnelles, offrant une qualité d'ajustement supérieure et une incertitude réduite pour mieux éclairer les tolérances d'ingénierie et la conception des tokamaks conventionnels futurs.

L'essentiel

Imaginez un tokamak (une machine conçue pour produire de l'énergie de fusion) comme un tourbillon géant et invisible de gaz surchauffé, maintenu en place par de puissants champs magnétiques. Idéalement, ce tourbillon est parfaitement symétrique, tel une toupie lisse en rotation. Cependant, dans le monde réel, les aimants qui le retiennent ne sont pas parfaits. Ils présentent de minuscules inclinaisons, des déplacements et des imperfections. Ces imperfections créent des « champs d'erreur » — de minuscules ondulations magnétiques indésirables qui peuvent perturber la rotation fluide du gaz.

Si ces ondulations deviennent trop fortes, elles peuvent amener le tourbillon à développer un « nœud » (une île magnétique) qui reste bloqué sur place. Une fois bloqué, l'ensemble du système peut s'effondrer, entraînant un arrêt soudain appelé « disruption ». C'est une mauvaise nouvelle pour la machine et pour les personnes qui la construisent.

Le Problème : Quelle est la limite ?
Les ingénieurs doivent connaître la limite exacte : Jusqu'à quelle taille ces imperfections magnétiques peuvent-elles grandir avant que la machine ne cède ? S'ils fixent la limite trop bas, ils doivent construire la machine avec une précision impossible, la rendant astronomiquement chère et lente à construire. S'ils la fixent trop haut, la machine pourrait s'effondrer.

Pendant des années, les scientifiques ont tenté de créer une « règle empirique » (une loi d'échelle) pour prédire cette limite en fonction de la taille de la machine et du comportement du gaz. Mais les anciennes règles étaient un peu fragiles, comme une carte aux bords flous.

La Solution : Une carte plus précise
Cet article présente une nouvelle « carte » améliorée (une loi d'échelle empirique) beaucoup plus claire et plus fiable. Voici comment ils l'ont faite, en utilisant des analogies simples :

• Nettoyage des données : Les chercheurs ont passé en revue une vaste base de données d'expériences passées provenant de tokamaks du monde entier (comme DIII-D, JET et KSTAR). Ils ont décidé de se concentrer uniquement sur des types spécifiques de « météo » dans la machine : les conditions « Ohmiques » et « L-mode ». Ils ont exclu le « H-mode » car cet état est comme une forteresse solide — il est très difficile à briser, il n'aide donc pas à comprendre le point le plus faible de la machine. En se concentrant sur les états vulnérables, ils ont trouvé la véritable zone de danger.
• Ajout de nouveaux ingrédients : Ils ont ajouté de nouvelles données provenant de deux machines spécifiques : J-TEXT (qui est plus petite et fonctionne à des courants plus faibles) et davantage de données de JET (qui est énorme, similaire à la future machine ITER). Pensez-y comme à l'ajout de nouveaux essais routiers à une base de données de sécurité automobile. Vous avez besoin de petites voitures et de camions géants dans les données pour savoir comment les règles de sécurité s'appliquent à n'importe quel véhicule que vous pourriez construire à l'avenir.
• Meilleure mathématique : Ils ont utilisé une méthode mathématique plus sophistiquée pour trouver la relation entre la taille de la machine, la force du champ magnétique, la densité du gaz et le courant électrique. Ils ont découvert que le courant du plasma (la quantité d'électricité qui circule dans le gaz) est un facteur critique qu'ils n'avaient pas pleinement pris en compte auparavant.

Les nouvelles découvertes
La nouvelle « règle empirique » nous indique que :

1. Une densité plus élevée est votre alliée : Empaqueter plus de gaz dans la machine rend plus difficile pour les champs d'erreur de provoquer un effondrement.
2. Les machines plus grandes sont étonnamment résilientes : Les machines plus grandes (comme le futur ITER) peuvent supporter de plus grandes imperfections magnétiques que nous ne le pensions auparavant.
3. Le courant compte : La quantité de courant circulant dans le plasma modifie la façon dont la machine réagit à ces erreurs.

Pourquoi cela compte pour l'avenir
L'article examine spécifiquement le projet ITER, une expérience de fusion internationale massive actuellement en construction. En utilisant leur nouvelle carte plus précise, les chercheurs ont effectué des millions de simulations (comme exécuter une prévision météorologique un million de fois avec des conditions de départ légèrement différentes).

Le Résultat : Ils ont découvert que ITER est dans bien meilleure posture que nous ne le pensions. La « zone de danger » pour les erreurs magnétiques est beaucoup plus éloignée que les imperfections réelles que ITER devrait avoir.

• L'ancienne carte : Suggérait qu'il y avait une chance décente que ITER trébuche sur ses propres lacets (modes bloqués).
• La nouvelle carte : Montre que la probabilité que cela se produise est incroyablement faible (moins de 1 sur un million pour le scénario le plus probable).

Le fond du problème
Cet article ne dit pas simplement « la fusion est difficile ». Il donne aux ingénieurs une règle de mesure beaucoup plus confiante et précise pour évaluer les tolérances de leurs machines. Parce que les nouvelles règles montrent que les machines sont plus robustes face aux erreurs magnétiques, les ingénieurs pourraient ne pas avoir besoin de construire les aimants avec une précision extrême et coûteuse. Cela pourrait économiser du temps et de l'argent tout en maintenant la machine en sécurité.

En bref : Ils ont pris une carte floue et confuse des limites de sécurité magnétique, l'ont nettoyée avec de meilleures données et des mathématiques plus intelligentes, et ont découvert que l'avenir des centrales à fusion est plus sûr et plus réalisable que nous ne le croyions auparavant.

Résumé technique

1. Énoncé du problème

Les tokamaks nécessitent une axisymétrie précise pour fonctionner de manière stable. Cependant, les imperfections de fabrication (inclinaisons, décalages et désalignements des bobines) créent des champs d'erreur (EF) non axisymétriques. Si ces champs dépassent un seuil de pénétration critique, ils peuvent déclencher une reconnexion magnétique, former des îlots magnétiques et provoquer le « verrouillage » de ces îlots sur la paroi du réacteur. Cela conduit à des modes verrouillés et à des disruptions ultérieures du plasma.

La prédiction de ce seuil est cruciale pour :

• Les tolérances d'ingénierie : Déterminer les désalignements de bobines admissibles lors de la construction.
• Les coûts et les délais : Des tolérances trop strictes augmentent les coûts et les délais de construction ; des tolérances trop laxistes risquent la défaillance du dispositif.
• Les futurs dispositifs : Des projets comme ITER et SPARC possèdent une énergie stockée élevée, rendant les disruptions particulièrement dangereuses.

Les lois d'échelle empiriques précédentes souffraient de :

• Une faible qualité d'ajustement : Des coefficients de détermination ($R^2$) aussi bas que 0,38.
• Une incertitude élevée : Jusqu'à 48 % d'incertitude dans les projections pour ITER.
• Des données incohérentes : Inclusion mixte de plasmas en mode H, mode L et Ohmiques, de tokamaks sphériques (NSTX) et de niveaux variables de modélisation des champs d'erreur intrinsèques à travers différents dispositifs.
• Un espace de paramètres limité : Manque de données pour les dispositifs à grand rayon majeur (comme ITER) et pour des régimes de densité spécifiques.

2. Méthodologie

Les auteurs ont développé une loi d'échelle empirique mise à jour en affinant la base de données et la méthodologie d'ajustement.

A. Affinement de la base de données

• Sélection du régime : L'étude a restreint l'ensemble de données aux tokamaks conventionnels (non sphériques) opérant dans les régimes Ohmique et L-mode. Le mode H a été exclu car il est plus résilient à la pénétration des champs d'erreur, et les auteurs cherchaient à modéliser le régime « le plus dangereux ». Les tirs à confinement Ohmique saturé (SOC) ont été exclus pour éviter des dépendances de densité confondantes.
• Inclusion des dispositifs :
  • Ajout : De nouvelles données provenant de J-TEXT (fournissant des points de faible courant/densité et une nouvelle taille de rayon majeur) et un ensemble de données élargi de JET (21 tirs contre 4 précédemment, couvrant un grand rayon majeur et un courant de plasma élevé).
  • Exclusion : NSTX (tokamak sphérique) et les tirs en mode H.
• Normalisation :
  • DIII-D : Passage d'une approximation d'équilibre groupée à des reconstructions d'équilibre individuelles pour chaque tir, utilisant le code SURFMN pour modéliser plus précisément les champs d'erreur intrinsèques.
  • J-TEXT : Utilisation du code TokaMaker pour reconstruire les équilibres (les reconstructions magnétiques de routine n'étant pas disponibles), assurant la compatibilité avec le GPEC (General Perturbed Equilibrium Code) utilisé pour tous les autres dispositifs.
  • Métrique : Utilisation du Recouvrement du mode dominant ($\delta$), une métrique indépendante du système de coordonnées qui prend en compte la réponse du plasma, calculée via GPEC.

B. Méthodologie de mise à l'échelle

• Technique de régression : Utilisation d'une régression par moindres carrés avec sélection progressive pour identifier la loi de puissance log-linéaire optimale.
• Sélection des paramètres : Expansion de l'espace de paramètres au-delà des études précédentes pour inclure des dépendances explicites sur le Courant de plasma ($I_p$), le Champ toroïdal ($B_T$), le Rayon majeur ($R_0$), la Densité ($n_e$) et la Pression normalisée ($\beta_n/\ell_i$).
• Validation : Utilisation des Nombres de condition pour détecter la multicolinéarité (arrêt à 5 variables pour maintenir le nombre de condition < 20) et comparaison des Moindres Carrés Ordinaires (OLS) avec les Moindres Carrés Pondérés (WLS) utilisant l'estimation de densité de noyau (KDE) pour éviter la suréchantillonnage des régions riches en données.
• Quantification de l'incertitude : Hypothèse d'une incertitude de base de 10 % sur le temps d'identification du mode verrouillé, propagée via des simulations de Monte Carlo pour générer des fonctions de densité de probabilité (PDF) et des fonctions de distribution cumulative (CDF) pour les prédictions.

3. Contributions clés

1. Nouvelles lois d'échelle : Dérivation de deux nouvelles équations d'échelle (Éq. 9 pour OLS et Éq. 10 pour WLS) avec des métriques statistiques considérablement améliorées.
   • Équation 9 (OLS) : $\delta \propto (\beta_n/\ell_i)^{0,25} I_p^{-0,97} R_0^{1,88} n_e^{0,77} B_T^{0,20}$
   • Équation 10 (WLS) : $\delta \propto (\beta_n/\ell_i)^{0,13} I_p^{-1,01} R_0^{1,57} n_e^{0,56} B_T^{0,30}$
2. Inclusion d'une dépendance explicite en $I_p$ : Contrairement aux échelles précédentes qui reposaient sur des dépendances implicites, les nouveaux modèles incluent explicitement le courant de plasma, révélant qu'un courant de plasma plus faible augmente la résilience aux champs d'erreur.
3. Changement de dépendance du champ toroïdal : L'échelle pour $B_T$ est passée d'une forte dépendance négative (dans les modèles précédents) à une faible dépendance positive, entraînée par l'inclusion de $I_p$ et la suppression des données de tokamaks sphériques.
4. Tests de robustesse : Démonstration que l'échelle est robuste face à la suppression de dispositifs individuels (par exemple, DIII-D, J-TEXT, C-MOD), tandis que des sous-ensembles « minimaux » de dispositifs (par exemple, JET+KSTAR+C-MOD) ont conduit à un surajustement sévère.

4. Résultats

• Qualité d'ajustement améliorée : La nouvelle échelle a atteint un $R^2$ de 0,63 (OLS) et 0,66 (WLS), une amélioration significative par rapport aux valeurs précédentes (0,38–0,58).
• Réduction de l'incertitude :
  • Chi-carré réduit : Amélioration de 17,0 (modèle 2020) à 12,99.
  • Erreur logarithmique absolue moyenne : Amélioration de 1,336 à 1,306 (équivalent à une réduction de l'incertitude de prédiction de $\pm33,6\%$ à $\pm30,6\%$).
• Dépendances des paramètres :
  • Densité ($n_e$) : Confirmation d'une dépendance positive ($n_e^{0,56-0,77}$), cohérente avec de nombreux dispositifs individuels mais distincte de certains modèles théoriques.
  • Rayon majeur ($R_0$) : Une forte dépendance positive ($R_0^{1,57-1,88}$), indiquant que les dispositifs plus grands sont plus résilients.
  • Courant de plasma ($I_p$) : Une forte dépendance négative ($I_p^{-0,97}$ à $-1,01$), indiquant qu'un courant plus élevé rend le plasma plus susceptible au verrouillage.
• Valeurs aberrantes spécifiques aux dispositifs :
  • NSTX : La nouvelle échelle sous-estime considérablement la tolérance aux champs d'erreur de NSTX (prévoyant $\delta \sim 10^{-5}$ contre $\sim 10^{-4}$ réel), probablement en raison de sa nature sphérique et de son profil de $q_{95}$ élevé, confirmant qu'il devrait être exclu des échelles de tokamaks conventionnels.
  • C-MOD : Les données à haute densité de C-MOD ont été cruciales pour capturer les comportements d'échelle de densité non monotones.

5. Signification et implications pour ITER

• Projection pour ITER :
  • La nouvelle échelle prédit un seuil de pénétration des champs d'erreur beaucoup plus élevé pour ITER par rapport aux modèles précédents.
  • $\delta$ intrinsèque le plus probable pour ITER : $0,38 \times 10^{-4}$.
  • Seuil prédit (Éq. 10) : $4,01 \times 10^{-4}$ (avec une plage d'incertitude à 1 sigma).
  • Probabilité de verrouillage : La probabilité que ITER se verrouille dans son scénario de base est calculée comme extrêmement faible (0,0001 % pour le champ intrinsèque le plus probable, et 22,3 % même pour le champ intrinsèque maximal possible).
• Impact sur l'ingénierie :
  • Les résultats suggèrent que les tolérances d'ingénierie actuelles d'ITER pour l'alignement des bobines sont conservatrices et suffisantes.
  • L'échelle fournit un outil plus fiable pour concevoir de futurs dispositifs (comme SPARC) et optimiser les opérations sur les machines existantes en prédisant la probabilité de verrouillage de mode basée sur des simulations de Monte Carlo des champs d'erreur intrinsèques.
• Travaux futurs : Les auteurs suggèrent que l'ajout de données provenant de JT-60SA (grand rayon) et d'ASDEX-U (haut $\beta_n$) affinerait davantage l'échelle. De plus, développer une échelle basée sur des tirs avec des mesures de rotation explicites pourrait améliorer encore le pouvoir prédictif, car la rotation est un facteur dominant dans la physique des champs d'erreur.

En conclusion, cet article établit une échelle empirique plus robuste, cohérente avec la physique et statistiquement fiable pour la pénétration des champs d'erreur n=1, augmentant considérablement la confiance dans la conception et l'exploitation des tokamaks conventionnels de nouvelle génération.