Axisymmetric reduction of the adjoint operator in ideal MHD equilibrium

Résumé technique : Réduction axisymétrique de l'opérateur adjoint dans l'équilibre MHD idéal

Énoncé du problème
Les calculs d'équilibre de la magnétohydrodynamique (MHD) idéale sont fondamentaux pour la conception de dispositifs de confinement magnétique et la reconstruction des états de plasma à partir de données de diagnostic. Si les équilibres axisymétriques (par exemple, les tokamaks) sont régis par l'équation de Grad–Shafranov, les équilibres tridimensionnels (3D) (par exemple, les stellarators) nécessitent des solveurs numériques. Les deux régimes utilisent de plus en plus des approches basées sur l'optimisation, où une analyse de sensibilité par rapport à de nombreux paramètres est requise. La méthode de l'adjoint offre un cadre efficace pour une telle analyse de sensibilité. Cependant, une relation mathématique précise entre la formulation de l'adjoint du problème complet d'équilibre 3D et l'adjoint de l'équation de Grad–Shafranov axisymétrique n'a pas été explicitement démontrée. Plus précisément, il reste à déterminer comment l'opérateur adjoint 3D se réduit sous l'hypothèse de l'axisymétrie et si la structure résultante s'aligne avec les formulations de l'adjoint axisymétrique établies, particulièrement concernant les facteurs de pondération géométrique inhérents à l'énergie magnétique axisymétrique.

Méthodologie
Les auteurs emploient une approche théorique des opérateurs pour dériver l'adjoint des équations d'équilibre de la MHD idéale sans supposer initialement la symétrie.

1. Formulation générale de l'adjoint : L'équilibre MHD est traité comme un problème d'optimisation sous contraintes où les variables d'état sont la densité de courant ($J$), le champ magnétique ($B$) et la pression ($p$). Un lagrangien est construit en utilisant les contraintes d'équilibre et une fonctionnelle objectif.
2. Dérivation de l'adjoint 3D : En prenant les variations du lagrangien par rapport aux variables d'état et en appliquant le théorème de la divergence sous des conditions de bord fixes ($B \cdot n = 0$), les auteurs dérivent les équations de l'adjoint 3D. Les variables adjointes ($\alpha, \beta, \gamma$) sont introduites comme des multiplicateurs de Lagrange duaux aux contraintes d'équilibre statique, distinctes des déplacements lagrangiens physiques.
3. Réduction axisymétrique : Les auteurs imposent l'axisymétrie sur le système adjoint 3D dérivé. Ils utilisent des représentations standards pour les champs vectoriels sans divergence en géométrie toroïdale, exprimant le champ magnétique et le champ adjoint $\alpha$ en termes de fonctions de flux poloidal ($\psi$) et d'angles toroïdaux ($\phi$).
4. Projection et simplification : Le système réduit est projeté selon et transversalement à la direction toroïdale pour isoler les composantes scalaires. Les auteurs analysent spécifiquement la quantité scalaire $\eta = -\alpha \cdot \nabla \psi$, l'interprétant comme une perturbation de type flux.
5. Comparaison des opérateurs : L'équation scalaire résultante est comparée à l'opérateur de Grad–Shafranov perturbé. Les auteurs examinent également la structure du produit scalaire, spécifiquement le rôle du poids géométrique $1/R^2$ provenant de la variation de l'énergie magnétique poloidale.

Contributions clés et résultats

• Réduction explicite : L'article démontre que la réduction axisymétrique du système adjoint homogène 3D complet produit un opérateur scalaire qui coïncide exactement avec l'opérateur de Grad–Shafranov perturbé (Zakharov & Pletzer 1999). Cet opérateur régit les variations infinitésimales de la fonction de flux d'équilibre.
• Correspondance des variables : L'étude établit une correspondance directe entre les variables adjointes 3D et la variable scalaire réduite. La variable scalaire de l'adjoint réduite est identifiée comme $\eta = -\alpha \cdot \nabla \psi$, où $\alpha$ est le champ adjoint incompressible associé à l'équation de force équilibrée.
• Structure du produit scalaire pondéré : L'analyse clarifie que l'opérateur scalaire réduit est formellement auto-adjoint uniquement par rapport à un produit scalaire pondéré $(u, v)_w = \int u v R^{-2} dV$. Ce facteur de poids $1/R^2$ est montré comme émanant naturellement de l'expression axisymétrique de l'énergie magnétique poloidale, plutôt que d'une imposition arbitraire.
• Cohérence structurelle : Ce travail résout les divergences avec les formulations précédentes (par exemple, Sanpei et al. 2021, 2023) qui omettaient le poids géométrique dans le produit scalaire. Les auteurs montrent que, bien que les formulations non pondérées soient pratiquement viables pour l'optimisation numérique, la formulation pondérée est nécessaire pour maintenir la compatibilité structurelle et l'auto-adjonction entre l'adjoint 3D et sa réduction axisymétrique.

Signification
Cet article fournit un énoncé de cohérence rigoureux au niveau des opérateurs, reliant la formulation de l'adjoint 3D complète au problème de l'adjoint de Grad–Shafranov axisymétrique. En clarifiant comment la structure de l'adjoint axisymétrique est intégrée dans le cadre 3D, ce travail offre une base théorique unifiée pour l'analyse de sensibilité. Cela est particulièrement pertinent pour :

• Configurations faiblement non-axisymétriques : La formulation de l'adjoint 3D sert de cadre « parent » pour analyser les équilibres qui sont de petites perturbations d'un état axisymétrique, permettant l'organisation systématique des réponses axisymétriques et des corrections 3D.
• Reconstruction et optimisation basées sur l'adjoint : La correspondance explicite entre les variables 3D et axisymétriques aide à interpréter les résultats de sensibilité et peut aider à l'analyse de bifurcation ou à la caractérisation des réponses singulières dans les systèmes d'équilibre linéarisés.

Les auteurs concluent que cette réduction clarifie la relation mathématique entre les problèmes d'adjoint MHD idéale 3D et axisymétriques, garantissant que les propriétés structurelles (telles que l'auto-adjonction) sont préservées à travers les dimensions lorsque les poids géométriques appropriés sont appliqués.