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

Het artikel introduceert VEQ, een snelle parametrische solver voor tokamak-evenwichten met vaste grenzen die flexibele bronprofielen en harmonische expansies gebruikt om lage latentie en hoge nauwkeurigheid van evenwichtsreconstructie te bereiken, geschikt voor herhaalde queries in transport-geometrie koppelingsworkflows.

De kern

Stel je een tokamak (een donutvormige kernfusiereactor) voor als een gigantische, onzichtbare ballon gevuld met superheet plasma. Om deze ballon van uiteen te laten knappen of in te laten storten, moeten wetenschappers precies weten hoe de magnetische "huid" die het bij elkaar houdt gevormd is en hoe de krachten binnenin in evenwicht zijn. Dit wordt het "evenwicht" genoemd.

Normaal gesproken is het berekenen van dit evenwicht alsof je een enorme, 3D-puzzel probeert op te lossen waarbij elk stukje van vorm verandert zodra je het aanraakt. Het is accuraat, maar het duurt lang—te lang als je de vorm honderden keren per seconde moet controleren om de reactor aan te sturen.

Dit artikel introduceert VEQ (Veloce EQuilibrium), een nieuwe tool die is ontworpen om de "fast-forward"-knop te zijn voor deze berekeningen. Zo werkt het, met behulp van eenvoudige analogieën:

1. De "Vormveranderaar" versus de "Pixel-schilder"

Traditionele methoden zijn als een Pixel-schilder. Ze proberen het hele magnetische veld in kaart te brengen door de waarde op miljoenen piepkleine rasterpunten (pixels) te berekenen. Het is gedetailleerd, maar zwaar en traag.

VEQ is meer als een Vormveranderaar. In plaats van elke pixel te schilderen, beschrijft VEQ de plasma vorm met een flexibel, wiskundig "skelet" gemaakt van een paar belangrijke knoppen en schuiven (parameters).

• Denk aan de plasma vorm als een stuk klei.
• Traditionele oplossers proberen elk klein bultje op de klei te boetseren.
• VEQ gebruikt een set vooraf gedefinieerde "rekkende" en "buigende" instrumenten (wiskundige harmonieken en polynomen). Je draait gewoon een paar knoppen om de klei in de juiste vorm te trekken. Dit maakt de berekening ongelooflijk snel omdat er veel minder knoppen gedraaid hoeven te worden dan pixels geschilderd moeten worden.

2. De "Universele Vertaler"

Een van de grootste hoofdpijndossiers in fusieonderzoek is dat verschillende computerprogramma's verschillende talen spreken. Het ene programma geeft je misschien de "druk", een ander de "stroom", en een derde de "veiligheidsfactor" (een maatstaf voor stabiliteit).

VEQ fungeert als een Universele Vertaler. Het artikel laat zien dat VEQ zes verschillende "inputroutes" heeft (zoals verschillende USB-poorten). Je kunt gegevens van al deze verschillende bronnen inpluggen, en VEQ vertaalt ze allemaal naar zijn eigen interne taal om het probleem op te lossen.

• De claim: De auteurs testten dit door dezelfde perfecte data in alle zes de poorten te voeren. Ze ontdekten dat ongeacht welke poort werd gebruikt, VEQ exact hetzelfde resultaat produceerde. Dit bewijst dat de vertaler perfect werkt en geen fouten introduceert enkel omdat je een andere inputkabel gebruikte.

3. De "Snelheid versus Nauwkeurigheid" afweging

Het artikel zegt niet alleen "het is snel"; het onderzoekt de afweging tussen snelheid en precisie, zoals het kiezen tussen een schets en een foto.

• De Schets (Lage parameters): Je gebruikt heel weinig knoppen. Het is direct beschikbaar (milliseconden) en goed genoeg voor een snelle blik, maar het kan kleine details missen.
• De Foto (Hoge parameters): Je gebruikt veel knoppen. Het duurt iets langer (nog steeds zeer snel, rond de 19 milliseconden voor complexe vormen), maar legt de vorm met hoge precisie vast.
• Het resultaat: De auteurs testten dit op drie verschillende soorten plasma vormen (een standaard "D"-vorm, een hoogwaardige "H-mode", en een complexe vorm met een "X-punt"). Ze ontdekten dat zelfs met een klein aantal knoppen, VEQ de complexe vormen kon reproduceren met een fout die zo klein is dat hij bijna onzichtbaar is (minder dan 0,2% van de grootte van de reactor).

4. De "Stress-test" (Waar de barsten zichtbaar worden)

De auteurs waren eerlijk over de beperkingen. Ze controleerden waar de "krachtbalans" (de spanning die het plasma bij elkaar houdt) perfect was en waar niet.

• Het Interieur: In het midden van het plasma is VEQ uitstekend. De krachten zijn bijna perfect in evenwicht.
• De Rand: Nabij de uiterste buitenste huid (de grens), zijn de fouten iets hoger. Dit komt omdat VEQ een glad, flexibel skelet gebruikt, terwijl echte plasma grenzen grillig kunnen zijn of scherpe hoeken kunnen hebben (zoals een X-punt).
• De Conclusie: VEQ is geweldig voor het "midden" van de reactor. Als je precies wilt weten wat er aan de uiterste rand gebeurt, heb je misschien een tragere, meer gedetailleerde tool nodig om het te controleren. Maar voor de meeste besturingstaken is VEQ snel en nauwkeurig genoeg.

5. De "Transport"-test

Ten slotte hebben ze getest of de kleine fouten in de vorm van VEQ grote problemen zouden veroorzaken als je het zou gebruiken om te voorspellen hoe warmte door het plasma beweegt (een "transport"-test).

• Het resultaat: De fouten waren minuscuul (minder dan 1%). Het is alsof je een kamer meet met een licht verbogen liniaal; de fout in de afmeting van de kamer zou je beslissing over de grootte van het tapijt niet veranderen.

Samenvatting

VEQ is een nieuwe, supersnelle rekenmachine voor fusiereactoren. In plaats van elk enkel punt in kaart te brengen, gebruikt het een flexibel wiskundig skelet om de plasma vorm te beschrijven.

• Het is snel: Het kan complexe vormen oplossen in milliseconden.
• Het is flexibel: Het accepteert gegevens van veel verschillende bronnen.
• Het is betrouwbaar: Het werkt goed voor het binnenste van de reactor en is nauwkeurig genoeg voor de meeste besturingssystemen, mits men de uiterste randen in de gaten houdt.

De auteurs concluderen dat VEQ perfect is voor systemen die herhaaldelijk de vraag moeten stellen: "Hoe ziet het plasma er op dit moment uit?", zoals bij real-time reactorbesturing of het draaien van duizenden simulaties om de beste operationele condities te vinden.

Technische samenvatting

Technische Samenvatting: VEQ – Een Snelle Parametrische Grad–Shafranov Solver voor Fixed-Boundary Tokamak-evenwichten

Probleemstelling
Geïntegreerde modelleringsworkflows voor tokamaks (bijv. IMAS, RAPTOR, SuperCode) vereisen herhaalde queries van evenwichtsgeometrie, metriek-coëfficiënten en diagnostische profielen. Deze workflows verschillen van eenmalige high-fidelity reconstructies; ze vereisen een representatie die compact genoeg is voor parametrische scans en controle-georiënteerde iteraties, maar tegelijkertijd voldoende twee-dimensionale vormgeving en krachtbalans-informatie behoudt om bruikbaar te zijn buiten eenvoudige analytische sluitingen. Bestaande tools vallen vaak uiteen in twee extremen: laag-orde engineering-sluitingen (bijv. Miller-modellen) zijn snel maar missen residuele diagnostiek, terwijl volledige evenwichtssolvers (bijv. CHEASE, EFIT, DESC) high-fidelity native outputs bieden maar vaak te rekenintensief of infrastructuur-zwaar zijn voor gebruik in de binnenste lus. Er is een gat voor een continue, fixed-boundary representatie die nuttige 2D-vormgeving en residuele diagnostiek behoudt terwijl deze goedkoop te regenereren blijft.

Methodologie
Het artikel introduceert VEQ (Veloce EQuilibrium), een compact parametrisch framework dat ontworpen is om dit tussenliggende regime op te vullen. De specifieke implementatie die wordt getest, VEQPy, is een axiaal symmetrische, fixed-boundary Grad–Shafranov solver.

• Parametrische Representatie: Het evenwicht wordt gedefinieerd door een parametrische flux-oppervlakte kaart $(\rho, \theta) \to [R(\rho, \theta), Z(\rho, \theta)]$. De geometrie wordt geparametriseerd met behulp van MXH-type flux-oppervlakte harmonieken (vervormde polaire hoeken) en verschoven-Chebyshev polynomen voor radiale profielen. De actieve onbekenden zijn de coëfficiënten van deze harmonieken en profielen, in plaats van de nodalwaarden van een grid-gebaseerde polaire fluxkaart.
• Variationele Formulering: De solver dwingt een variationeel geïnduceerde geprojecteerde Grad–Shafranov residu af. In plaats van de sterk-vorm vergelijking puntgewijs op te lossen, beperkt de methode de toegestane variaties tot die gegenereerd door de actieve coëfficiënten van de parametrische representatie. Dit resulteert in een eind-dimensionaal systeem waarbij de vorm-block een variationeel gestructureerde projectie is (Petrov–Galerkin conditie) en de profiel-blocks worden gesloten door moment-projecties.
• Flexibele Inputroutes: Om compatibiliteit te garanderen met diverse inputs uit analytische studies, reconstructiebestanden en controlemodellen, ondersteunt VEQ zes inputroutes (PF, PP, PI, PJ1, PJ2, PQ). Deze routes accepteren drukgradiënten, toroidale-veld functies, polaire-flux gradiënten, ingesloten toroidale stroom, stroomdichtheid of veiligheidsfactor-profielen. Een route-specifieke sluiting brengt deze heterogene inputs naar een gemeenschappelijke genormaliseerde bronrepresentatie, waardoor alle routes in dezelfde eind-dimensionale residu-operator kunnen voeden.
• Implementatie: De solver gebruikt een niet-lineaire solver (SciPy's Powell hybrid methode) met een versnelde backend (Numba). Het scheidt de setup (grid-afhankelijke precomputatie, operator data) van de solve-fase om lage-latentie herhaalde queries mogelijk te maken.

Belangrijkste Bijdragen

1. Eind-dimensionale Parametrische Representatie: Definitie van een compacte MXH–Chebyshev flux-oppervlakte representatie en de afleiding van de variationeel geïnduceerde geprojecteerde Grad–Shafranov residu, inclusclusief convectieve variaties geassocieerd met flux-oppervlakte beweging.
2. Route-Specifieke Bron-Sluitingen: Introductie van een framework dat heterogene 1D inputs (druk, stroom, veiligheidsfactor, etc.) mapt naar een gemeenschappelijke residu-systeem, waarbij onderscheid wordt gemaakt tussen input-route consistentie en geometrische benaderingsfout.
3. Duale Diagnostische Aanpak: Rapportage van zowel de geprojecteerde residu-norm (geminimaliseerd door de solver) als gesamplede puntgewijze strong-form Grad–Shafranov residuen. Dit onderscheidt de convergentie van het gereduceerde systeem van puntgewijze krachtbalans-screening.
4. Benchmarking: Uitgebreide testen op Solov'ev, CHEASE, en EFIT-afgeleide G-EQDSK gevallen, waarbij input-route consistentie, flux-oppervlakte reconstructiefout, accuratesse–kosten trade-offs en herhaalde solve-latentie worden gekwantificeerd.

Resultaten

• Route Consistentie: Gecontroleerde tests met gladde, onderling compatibele inputs gegenereerd uit een gemeenschappelijk referentie-evenwicht toonden aan dat alle zes de inputroutes naar hetzelfde fixed-boundary parametrische evenwicht convergeren, waarbij verschillen uitsluitend toe te schrijven zijn aan numerieke recovery paden (bijv. integratie/differentiatie stappen).
• Accuratesse en Latentie: De solver werd getest op drie representatieve gevallen: een D-vormige Solov'ev case, een H-mode CHEASE case, en een X-point EFIT case.
  • D-vormig: Bereikte een minor-radius-genormaliseerde vormfout van $1.4 \times 10^{-3}$ met 9 actieve parameters en een solve-only median tijd van 1.6 ms.
  • H-mode: Bereikte een fout van $1.1 \times 10^{-3}$ met 65 parameters en een median tijd van 19 ms.
  • X-point: Bereikte een fout van $1.9 \times 10^{-3}$ met 94 parameters en een median tijd van 15 ms.
• Residu Diagnostiek: Puntgewijze strong-form diagnostiek onthulde dat het verrijken van de actieve representatie primair de interne krachtbalans verbetert. In H-mode en X-point gevallen bleven globale RMS en maximale residu waarden gedomineerd door nabij-de-rand bijdragen, wat aangeeft dat de fixed-boundary ansatz de boundary-fit mismatches of singulariteiten (zoals X-punten) niet volledig kan compenseren via interne vormvariaties alleen.
• Transport Koppeling: In een geïsoleerde 1D transport-geometrie koppeling test, bleef de temperatuur-profiel respons op VEQ geometrie fouten onder ongeveer 1%.

Significantie en Claims
Het artikel claimt dat VEQ geschikt is voor herhaalde evenwichts-geometrie queries in geïntegreerde modellering en controle workflows, mits puntgewijze diagnostiek behouden blijft om gevallen te screenen die verdere boundary verfijning, lokale correctie of hogere-fidelity solvers vereisen.

De auteurs positioneren VE_Q als complementair aan high-fidelity reconstructie solvers, niet als een universele vervanging. De primaire waarde ligt in het bieden van een continue, compacte representatie met expliciete geprojecteerde residuen en gesamplede diagnostiek die goedkoop te regenereren is. Het artikel stelt expliciet dat de timing resultaten "solve-only" latentie reflecteren na de setup en niet geïnterpreteerd moeten worden als een vervanging voor volledige pipeline timings die preprocessing, I/O en discretisatie bevatten. De methode is bedoeld om de kloof te overbruggen tussen analytische sluitingen en volledige solver-native pipelines, waarbij een balans wordt geboden tussen geometrische rijkdom en computationele efficiëntie.