VEQ: a fast parametric Grad--Shafranov solver for fixed-boundary tokamak equilibria with flexible source profiles

L'article présente VEQ, un solveur paramétrique rapide pour les équilibres de tokamaks à frontières fixes qui utilise des profils de sources flexibles et des développements en harmoniques pour obtenir une reconstruction d'équilibre à faible latence et haute précision, adaptée aux requêtes répétées dans les flux de travail de couplage de géométrie de transport.

L'essentiel

Imaginez un tokamak (un réacteur de fusion nucléaire en forme de doughnut) comme un immense ballon invisible rempli de plasma surchauffé. Pour empêcher ce ballon d'éclater ou de s'effondrer, les scientifiques doivent savoir exactement comment la « peau » magnétique qui le maintient ensemble est façonnée et comment les forces à l'intérieur sont équilibrées. C'est ce qu'on appelle trouver l'« équilibre ».

Habituellement, calculer cet équilibre revient à essayer de résoudre un puzzle géant en 3D où chaque pièce change de forme dès qu'on la touche. C'est précis, mais cela prend beaucoup de temps — trop longtemps si vous devez vérifier la forme des centaines de fois par seconde pour contrôler le réacteur.

Ce document présente VEQ (Veloce EQuilibrium), un nouvel outil conçu pour être le bouton « avance rapide » de ces calculs. Voici comment il fonctionne, en utilisant des analogies simples :

1. Le « Métamorphe » vs le « Peintre de Pixels »

Les méthodes traditionnelles sont comme un Peintre de Pixels. Elles tentent de cartographier l'ensemble du champ magnétique en calculant la valeur de millions de minuscules points de grille (pixels). C'est détaillé, mais lourd et lent.

VEQ ressemble davantage à un Métamorphe. Au lieu de peindre chaque pixel, il décrit la forme du plasma à l'aide d'un « squelette » mathématique flexible composé de quelques boutons et curseurs clés (appelés paramètres).

• Imaginez la forme du plasma comme un morceau d'argile.
• Les solveurs traditionnels essaient de sculpter chaque minuscule bosse sur l'argile.
• VEQ utilise un ensemble d'outils de « déformation » et de « courbure » prédéfinis (harmoniques et polynômes mathématiques). Il vous suffit de tourner quelques boutons pour étirer l'argile afin de lui donner la bonne forme. Cela rend le calcul incroyablement rapide car il y a beaucoup moins de boutons à tourner que de pixels à peindre.

2. Le « Traducteur Universel »

L'un des plus grands casse-têtes de la recherche sur la fusion est que différents programmes informatiques parlent des langues différentes. Un programme peut vous donner la « pression », un autre le « courant », et un autre le « facteur de sécurité » (une mesure de stabilité).

VEQ agit comme un Traducteur Universel. Le document montre que VEQ possède six « routes d'entrée » différentes (comme des ports USB différents). Vous pouvez brancher des données provenant de n'importe laquelle de ces sources différentes, et VEQ les traduit toutes dans son propre langage interne pour résoudre le problème.

• L'affirmation : Les auteurs ont testé cela en injectant les mêmes données parfaites dans les six ports. Ils ont constaté que, quel que soit le port utilisé, VEQ produisait exactement le même résultat. Cela prouve que le traducteur fonctionne parfaitement et n'introduit pas d'erreurs simplement parce que vous avez utilisé un câble d'entrée différent.

3. Le compromis « Vitesse vs Précision »

Le document ne se contente pas de dire qu'il est rapide ; il explore le compromis entre vitesse et précision, comme le choix entre un croquis et une photographie.

• Le Croquis (Peu de paramètres) : Vous utilisez très peu de boutons. C'est instantané (millisecondes) et suffisant pour un coup d'œil rapide, mais cela peut manquer de petits détails.
• La Photo (Nombreux paramètres) : Vous utilisez beaucoup de boutons. Cela prend un peu plus de temps (toujours très rapide, environ 19 millisecondes pour des formes complexes), mais cela capture la forme avec une haute précision.
• Le Résultat : Les auteurs ont testé cela sur trois types différents de formes de plasma (une forme en « D » standard, une forme « H-mode » de haute performance, et une forme complexe avec un « point X »). Ils ont découvert que même avec un petit nombre de boutons, VEQ pouvait recréer les formes complexes avec une erreur si petite qu'elle est presque invisible (moins de 0,2 % de la taille du réacteur).

4. Le « Test de Stress » (Où les fissures apparaissent)

Les auteurs ont été honnêtes sur les limites. Ils ont vérifié là où l'« équilibre des forces » (la tension qui maintient le plasma) était parfait et là où il ne l'était pas.

• L'Intérieur : Au milieu du plasma, VEQ est excellent. Les forces sont équilibrées presque parfaitement.
• Le Bord : Près de la peau extérieure très fine (la limite), les erreurs sont légèrement plus élevées. C'est parce que VEQ utilise un squelette lisse et flexible, alors que les limites réelles du plasma peuvent être dentelées ou présenter des coins tranchants (comme un point X).
• La Conclusion : VEQ est excellent pour le « milieu » du réacteur. Si vous avez besoin de savoir exactement ce qui se passe à la limite extrême, vous devrez peut-être effectuer une double vérification avec un outil plus lent et plus détaillé. Mais pour la plupart des tâches de contrôle, VEQ est assez rapide et précis.

5. Le Test de « Transport »

Enfin, ils ont testé si les petites erreurs de forme de VEQ causeraient de gros problèmes si vous essayiez de l'utiliser pour prédire comment la chaleur se déplace à travers le plasma (un test de « transport »).

• Le Résultat : Les erreurs étaient infimes (moins de 1 %). C'est comme si vous mesuriez une pièce avec une règle légèrement tordue ; l'erreur sur la taille de la pièce ne changerait pas votre décision sur la taille du tapis à acheter.

Résumé

VEQ est un nouveau calculateur ultra-rapide pour les réacteurs de fusion. Au lieu de cartographier chaque point individuel, il utilise un squelette mathématique flexible pour décrire la forme du plasma.

• Il est rapide : Il peut résoudre des formes complexes en quelques millisecondes.
• Il est flexible : Il accepte des données provenant de nombreuses sources différentes.
• Il est fiable : Il fonctionne bien pour l'intérieur du réacteur et est suffisamment précis pour la plupart des systèmes de contrôle, à condition de surveiller de près les bords extérieurs.

Les auteurs concluent que VEQ est parfait pour les systèmes qui ont besoin de demander « À quoi ressemble le plasma en ce moment ? » de manière répétée, comme dans le contrôle en temps réel d'un réacteur ou lors de l'exécution de milliers de simulations pour trouver les meilleures conditions de fonctionnement.

Résumé technique

Résumé Technique : VEQ – Un solveur de Grad–Shafranov paramétrique rapide pour les équilibres de tokamaks à frontières fixes

Énoncé du Problème
Les flux de travail de modélisation intégrée pour les tokamaks (par exemple, IM mass, RAPTOR, SuperCode) nécessitent des requêtes répétées de la géométrie de l'équilibre, des coefficients métriques et des profils de diagnostic. Ces flux de travail diffèrent des reconstructions de haute fidélité ponctuelles ; ils exigent une représentation suffisamment compacte pour les balayages de paramètres et les itérations orientées vers le contrôle, tout en conservant une information de mise en forme 2D et d'équilibre des forces suffisante pour être utile au-delà des simples fermetures analytiques. Les outils existants se situent souvent aux deux extrêmes : les fermetures d'ingénierie à bas ordre (par exemple, les modèles Miller) sont rapides mais manquent de diagnostics résiduels, tandis que les solveurs d'équilibre complets (par exemple, CHEASE, EFIT, DESC) fournissent des sorties natives de haute fidélité mais sont souvent trop coûteux en termes de calcul ou lourds en termes d'infrastructure pour un usage dans les boucles internes. Il existe un fossé pour une représentation continue à frontières fixes qui préserve une mise en forme 2D utile et des diagnostics résiduels tout en restant peu coûteuse à régénérer.

Méthodologie
Le document présente VEQ (Veloce EQuilibrium), un cadre paramétrique compact conçu pour combler ce régime intermédiaire. L'implémentation spécifique testée, VEQPy, est un solveur de Grad–Shafranov axisymétrique à frontières fixes.

• Représentation Paramétrique : L'équilibre est défini par une application de surface de flux paramétrique $(\rho, \theta) \to [R(\rho, \theta), Z(\rho, \theta)]$. La géométrie est paramétrée à l'aide d'harmoniques de surface de flux de type MXH (angles poloïdaux déformés) et de polynômes de Chebyshev décalés pour les profils radiaux. Les inconnues actives sont les coefficients de ces harmoniques et de ces profils, plutôt que les valeurs nodales d'une carte de flux poloidal sur une grille.
• Formulation Variationnelle : Le solveur impose un résiduel de Grad–Shafranov projeté par voie variationnelle. Au lieu de résoudre l'équation sous forme forte point par point, la méthode restreint les variations admissibles à celles générées par les coefficients actifs de la représentation paramétrique. Cela aboutit à un système de dimension finie où le bloc de forme est une projection de structure variationnelle (condition de Petrov–Galerkin) et où les blocs de profils sont fermés par des projections de moments.
• Voies d'Entrée Flexibles : Pour assurer la compatibilité avec diverses entrées provenant d'études analytiques, de fichiers de reconstruction et de modèles de contrôle, VEQ prend en charge six voies d'entrée (PF, PP, PI, PJ1, PJ2, PQ). Ces voies acceptent des gradients de pression, des fonctions de champ toroïdal, des gradients de flux poloidal, du courant toroïdal enfermé, de la densité de courant ou des profils de facteur de sécurité. Une voie de fermeture spécifique cartographie ces entrées hétérogènes vers une représentation de source normalisée commune, permettant à toutes les voies d'alimenter le même opérateur résiduel de dimension finie.
• Implémentation : Le solveur utilise un solveur non linéaire (la méthode hybride de Powell de SciPy) avec un backend accéléré (Numba). Il sépare la configuration (précalcul dépendant de la grille, données de l'opérateur) de la phase de résolution pour permettre des requêtes répétées à faible latence.

Contributions Clés

1. Représentation Paramétrique de Dimension Finie : Définition d'une représentation de surface de flux compacte de type MXH–Chebyshev et dérivation du résiduel de Grad–Shafranov projeté par voie variationnelle, incluant les variations convectives associées au mouvement des surfaces de flux.
2. Fermetures de Source Spécifiques aux Voies : Introduction d'un cadre qui cartographie des entrées 1D hétérogènes (pression, courant, facteur de sécurité, etc.) vers un système résiduel commun, distinguant la cohérence de la voie d'entrée de l'erreur d'approximation géométrique.
3. Approche de Diagnostic Duale : Rapport de la norme du résiduel projeté (minimisée par le solveur) et des résiduels de forme forte échantillonnés par points. Cela distingue la convergence du système réduit du criblage de l'équilibre des forces ponctuel.
4. Évaluation Comparative (Benchmarking) : Tests complets sur des cas Solov'ev, CHEASE et dérivés de G-EQDSK d'EFIT, quantifiant la cohérence des voies d'entrée, l'erreur de reconstruction de la surface de flux, les compromis précision–coût et la latence de résolution répétée.

Résultats

• Cohérence des Voies : Des tests contrôlés utilisant des entrées lisses et mutuellement compatibles générées à partir d'un équilibre de référence commun ont démontré que toutes les six voies d'entrée convergent vers le même équilibre paramétrique à frontières fixes, les différences étant attribuables uniquement aux chemins de récupération numérique (par exemple, étapes d'intégration/différentiation).
• Précision et Latence : Le solveur a été testé sur trois cas représentatifs : un cas Solov'ev en forme de D, un cas CHEASE en mode H, et un cas EFIT avec point X.
  • Forme de D : A atteint une erreur de forme normalisée par le rayon mineur de $1,4 \times 10^{-3}$ avec 9 paramètres actifs et un temps de résolution seule médian de 1,6 ms.
  • Mode H : A atteint une erreur de $1,1 \times 10^{-3}$ avec 65 paramètres et un temps médian de 19 ms.
  • Point X : A atteint une erreur de $1,9 \times 10^{-3}$ avec 94 paramètres et un temps médian de 15 ms.
• Diagnostics de Résiduel : Les diagnostics de forme forte ponctuels ont révélé qu'enrichir la représentation active améliore principalement l'équilibre des forces intérieur. Dans les cas de mode H et de point X, les valeurs RMS globales et maximales des résiduels restaient dominées par les contributions proches de la frontière, indiquant que l'ansatz à frontières fixes ne peut pas compenser entièrement les désaccords d'ajustement de la frontière ou les structures singulières (comme les points X) par les seules variations de forme intérieure.
• Couplage de Transport : Dans un test de couplage géométrie-transport 1D isolé, la réponse du profil de température aux erreurs de géométrie de VEQ est restée inférieure à environ 1 %.

Signification et Revendications
L'article affirme que VEQ est adapté pour des requêtes répétées de géométrie d'équilibre dans les flux de travail de modélisation intégrée et de contrôle, à condition de conserver des diagnostics ponctuels pour filtrer les cas nécessitant un raffinement de la frontière, une correction locale ou des solveurs de plus haute fidélité.

Les auteurs positionnent VEQ comme complémentaire aux solveurs de reconstruction de haute fidélité, et non comme un remplacement universel. Sa valeur primaire réside dans la fourniture d'une représentation compacte et continue avec des résiduels projetés explicites et des diagnostics échantillonnés qui sont peu coûteux à régénérer. L'article stipule explicitement que les résultats de temps de calcul reflètent la latence de « résolution seule » après la configuration et ne doivent pas être interprétés comme un remplacement des temps de pipeline complets incluant le prétraitement, les E/S et la discrétisation. La méthode est destinée à combler le fossé entre les fermetures analytiques et les pipelines natifs de solveurs complets, offrant un équilibre entre richesse géométrique et efficacité computationnelle.