Axisymmetric reduction of the adjoint operator in ideal MHD equilibrium

Technische Samenvatting: Axiale Reductie van de Adjoint-operator in Ideale MHD-evenwichten

Probleemstelling
Berekeningen van ideale magnetohydrodynamische (MHD) evenwichten zijn fundamenteel voor het ontwerp van magnetische opsluitingsapparaten en de reconstructie van plasmatoestanden uit diagnostische gegevens. Terwijl axiale evenwichten (bijv. tokamaks) worden beheerst door de Grad–Shafranov-vergelijking, vereisen driedimensionale (3D) evenwichten (bijv. stellaratoren) numerieke oplosers. Beide regimes maken steeds vaker gebruik van optimalisatiegestuurde benaderingen waarbij gevoeligheidsanalyse met betrekking tot talrijke parameters vereist is. De adjoint-methode biedt een efficiënt kader voor een dergelijke gevoeligheidsanalyse. Er bestaat echter geen expliciet aangetoonde wiskundige relatie tussen de adjoint-formulering van het volledige 3D-evenwichtsprobleem en de adjoint van de axiale Grad–Shafranov-vergelijking. Specifiek is het onduidelijk hoe de 3D adjoint-operator reduceert onder de aanname van axialiteit en of de resulterende structuur overeenkomt met gevestigde axiale adjoint-formuleringen, met name met betrekking tot de geometrische weegfactoren die inherent zijn aan axiale magnetische energie.

Methodologie
De auteurs hanteren een operator-theoretische benadering om de adjoint van de ideale MHD-evenwichtsvergelijkingen af te leiden zonder aanvankelijk symmetrie aan te nemen.

1. Algemene Adjoint-formulering: Het MHD-evenwicht wordt behandeld als een geconstrueerd optimalisatieprobleem waarbij de toestandsvariabelen de stroomdichtheid ($J$), het magnetisch veld ($B$) en de druk ($p$) zijn. Een Lagrangiaan wordt geconstrueerd met behulp van de evenwichtsvoorwaarden en een objectieve functionaal.
2. 3D Adjoint-afleiding: Door variaties van de Lagrangiaan te nemen met betrekking tot de toestandsvariabelen en de divergentiestelling toe te passen onder vaste randvoorwaarden ($B \cdot n = 0$), leiden de auteurs de 3D adjoint-vergelijkingen af. De adjoint-variabelen ($\alpha, \beta, \gamma$) worden geïntroduceerd als Lagrange-multiplicatoren dual aan de statische evenwichtsbeperkingen, die verschillen van fysieke Lagrangiaanse verplaatsingen.
3. Axiale Reductie: De auteurs leggen axialiteit op aan de afgeleide 3D adjoint-systeem. Zij maken gebruik van standaardrepresentaties voor divergentievrije vectoren in toroidale geometrie, waarbij het magnetische veld en de adjoint-veld $\alpha$ worden uitgedrukt in termen van poloidale fluxfuncties ($\psi$) en toroidale hoeken ($\phi$).
4. Projectie en Vereenvoudiging: Het gereduceerde systeem wordt langs en dwars op de toroidale richting geprojecteerd om scalaire componenten te isoleren. De auteurs analyseren specifiek de scalaire grootheid $\eta = -\alpha \cdot \nabla \psi$, die wordt geïnterpreteerd als een flux-achtige perturbatie.
5. Operatorvergelijking: De resulterende scalaire vergelijking wordt vergeleken met de geperturbeerde Grad–Shafranov-operator. De auteurs onderzoeken ook de inproductstructuur, specifiek de rol van de geometrische weegfactor $1/R^2$ die voortvloeit uit de variatie van de poloidale magnetische energie.

Belangrijkste Bijdragen en Resultaten

• Expliciete Reductie: Het artikel toont aan dat de axiale reductie van de volledige 3D homogene adjoint-systeem een scalaire operator oplevert die exact overeenkomt met de geperturbeerde Grad–Shafranov-operator (Zakharov & Pletzer 1999). Deze operator beheerst de infinitesimale variaties van de evenwichtsfluxfunctie.
• Variabele Correspondentie: De studie vestigt een directe correspondentie tussen de 3D adjoint-variabelen en de gereduceerde scalaire variabele. De gereduceerde scalaire adjoint-variabele wordt geïdentificeerd als $\eta = -\alpha \cdot \nabla \psi$, waarbij $\alpha$ het incompressibele adjoint-veld is dat geassocieerd wordt met de krachtbalansvergelijking.
• Gewogen Inproductstructuur: De analyse verheldert dat de gereduceerde scalaire operator formeel zelf-adjoint is met betrekking tot een gewogen inproduct $(u, v)_w = \int u v R^{-2} dV$. Deze weegfactor $1/R^2$ blijkt natuurlijk voort te vloeien uit de axiale expressie van de poloidale magnetische energie, in plaats van een willekeurige oplegging te zijn.
• Structurele Consistentie: Het werk lost discrepanties op met eerdere formuleringen (bijv. Sanpei et al. 2021, 2023) die de geometrische weegfactor in het inproduct weglieten. De auteurs tonen aan dat, hoewel ongegewogen formuleringen praktisch levensvatbaar zijn voor numerieke optimalisatie, de gewogen formulering noodzakelijk is om de structurele compatibiliteit en zelf-adjointheid tussen de 3D adjoint en de axiale reductie ervan te handhaven.

Betekenis
Dit artikel biedt een rigoureuze operator-niveau consistentieverklaring die de volledige 3D adjoint-formulering verbindt met het axiale Grad–Shafranov adjoint-probleem. Door te verhelderen hoe de axiale adjoint-structuur binnen het 3D-raamwerk is ingebed, biedt het werk een verenigde theoretische basis voor gevoeligheidsanalyse. Dit is bijzonder relevant voor:

• Zwak Axiaal Niet-Symmetrische Configuraties: De 3D adjoint-formulering dient als een "parent" raamwerk voor het analyseren van evenwichten die kleine perturbaties zijn van een axiale toestand, waardoor de systematische organisatie van axiale responsen en 3D-correcties mogelijk wordt.
• Adjoint-gebaseerde Reconstructie en Optimalisatie: De expliciete correspondentie tussen 3D en axiale variabelen helpt bij het interpreteren van gevoeligheidsresultaten en kan assisteren bij bifurcatie-analyse of het karakteriseren van singuliere responsen in geliniariseerde evenwichtssystemen.

De auteurs concluderen dat deze reductie de wiskundige relatie tussen 3D en axiale ideale MHD adjoint-problemen verheldert, waardoor de structurele eigenschappen (zoals zelf-adjointheid) behouden blijven over dimensies wanneer de juiste geometrische gewichten worden toegepast.