Kinetic Equilibrium Prediction at TCV using RAPTOR and FBT

Dit artikel presenteert een nieuwe workflow voor TCV die RAPTOR-transportsimulaties koppelt aan FBT-evenwichtsberekeningen om nauwkeurige voorspellingen te doen voor plasma-eigenschappen en spoelstromen, waardoor tokamak-exploitanten hun schotplannen realistischer kunnen aanpassen.

De kern

Stel je voor dat je een heel complexe, futuristische auto gaat bouwen: een tokamak. Dit is geen gewone auto, maar een gigantische oven die probeert dezelfde energie te maken als de zon (kernfusie). Het probleem? De "brandstof" in deze oven is een plasma, een wolk van geladen deeltjes die zich gedraagt als een levend, onvoorspelbaar dier. Als je het niet perfect in de gaten houdt, ontsnapt het of stort het in.

Dit artikel beschrijft een nieuwe manier om deze "auto" voor te bereiden voordat hij zelfs maar start. Het team van het TCV-experiment (in Zwitserland) heeft twee slimme computersystemen gekoppeld om te voorspellen hoe het plasma zich gaat gedragen.

Hier is de uitleg in simpele taal, met wat creatieve vergelijkingen:

1. Het Probleem: Gokken met een onzichtbare bal

Normaal gesproken moeten operators (de bestuurders van de tokamak) raden hoe ze de elektromagneten moeten instellen om het plasma in vorm te houden. Ze gebruiken vaak een "gokje" gebaseerd op ervaring.

• De analogie: Stel je voor dat je een ballon moet vasthouden met een touw, maar je kunt de ballon niet zien. Je moet raden hoe hard je moet trekken. Als je te zacht trekt, vliegt de ballon weg. Als je te hard trekt, knapt hij.

2. De Oplossing: Twee slimme vrienden die samenwerken

De auteurs hebben twee softwareprogramma's laten praten met elkaar: RAPTOR en FBT.

• RAPTOR (De Voorspeller):
  Dit programma is als een snel voorspellingsmodel voor het weer, maar dan voor de binnenkant van de oven. Het kijkt naar de plannen (hoeveel stroom, hoeveel verwarming) en zegt: "Als we dit doen, zal de temperatuur hier stijgen en de druk daar dalen." Het simuleert hoe het plasma zich gedraagt, van het moment dat het start tot het moment dat het stopt.

  • Vergelijking: Het is alsof je een video maakt van hoe de ballon zal bewegen voordat je hem überhaupt loslaat.

• FBT (De Ontwerper):
  Dit programma is de architect die de magneten instelt. Het berekent precies welke stroom door welke magneet moet lopen om de gewenste vorm te houden.

  • Vergelijking: De architect tekent de blauwdruk voor de touwen die de ballon vasthouden.

De Magie (KEP):
Vroeger gaf de architect (FBT) een standaardontwerp aan de bestuurder, gebaseerd op een simpel idee van hoe de ballon eruit zou zien. Nu geven ze eerst de voorspelling van RAPTOR aan FBT.

• Het resultaat: De architect krijgt nu te zien: "O, de ballon wordt zwaarder aan de linkerkant en warmer aan de rechterkant." Hierdoor past hij de touwen (de magneten) direct aan voor die specifieke situatie. Dit heet Kinetic Equilibrium Prediction (KEP).

3. Wat levert dit op?

Door deze twee systemen te koppelen, krijgen de operators een veel realistischer plan.

• Betere controle: Ze weten precies hoeveel stroom ze nodig hebben om de vorm van het plasma (bijvoorbeeld een sneeuwvlok-vorm of een gewone bol) vast te houden.
• Veiligheid: Het voorkomt dat het plasma "ontsnapt" of instabiel wordt, wat zou kunnen leiden tot een onderbreking van het experiment (een "crash").
• Snelheid: Het systeem is zo snel dat het binnen enkele minuten een voorspelling kan doen voor een experiment dat pas over een uur plaatsvindt.

4. De Test: Het bewijs in de praktijk

De auteurs hebben dit getest met 211 verschillende experimenten op de TCV-tokamak.

• Ze zagen dat hun voorspellingen over de temperatuur en dichtheid van het plasma zeer nauwkeurig waren (binnen een marge van 20%).
• Ze konden zelfs voorspellen wanneer het plasma zou switchen van een "luie" modus (L-mode) naar een "energieke" modus (H-mode), net zoals een auto die van een stadsweg op de snelweg gaat.
• Ze testten het zelfs met een lastige vorm: een sneeuwvlok-divertor (een complexe vorm met meerdere punten). Zonder de nieuwe voorspelling was dit moeilijk te controleren; met de nieuwe methosis bleef de vorm stabiel en precies op de juiste plek.

Samenvatting

In essentie hebben de wetenschappers een simulatie-systeem gebouwd dat de tokamak "in de toekomst" laat kijken. In plaats van te gokken met de instellingen, gebruiken ze een digitale tweeling van het plasma om te zien wat er gaat gebeuren, en passen ze de echte instellingen daarop aan.

Het is alsof je voor het eerst niet meer blindelings een auto bestuurt in de mist, maar een navigatiesysteem hebt dat je precies vertelt waar de bochten komen en hoe je het stuur moet draaien om veilig aan te komen. Dit maakt het onderzoek naar schone energie (fusie) veiliger, sneller en succesvoller.

Technische samenvatting

Titel: Kinetisch Evenwichtvoorspelling bij TCV met RAPTOR en FBT

Auteurs: C. Contré et al. (EPFL/Swiss Plasma Center en partners)
Doel: Presentatie van een nieuwe workflow voor het voorspellen van het plasma-evenwicht en de benodigde spoelstromen vóór een experiment op de Tokamak Configuration Variable (TCV).

---

1. Het Probleem

Het optimaliseren van tokamak-prestaties vereist het vermijden van operationele limieten die kunnen leiden tot storingen (disruptions) of vroegtijdige beëindiging van een puls. Traditionele methoden voor het voorbereiden van experimenten (shot preparation) vertrouwen vaak op vereenvoudigde aannames over de interne plasma-profielen (zoals druk en stroomdichtheid) bij het berekenen van de poloidale veldspoelstromen (PF-coils).

• Beperkingen: Bestaande methoden gebruiken vaak polynomen voor de druk- ($p'$) en stroomfuncties ($TT'$) die gebaseerd zijn op operatorervaring in plaats van fysica-gebaseerde voorspellingen.
• Gevolg: Dit kan leiden tot onnauwkeurige schattingen van kritieke grootheden zoals de interne inductie ($l_i$) en de poloidale beta ($\beta_{pol}$). Onjuiste voorspellingen kunnen resulteren in een slechte vormcontrole, misalignement van het plasma of verticale instabiliteiten.
• Uitdaging: Er is behoefte aan een "predict-first" workflow die snel en automatisch de kinetische evolutie (temperatuur, dichtheid, stroom) voorspelt op basis van de ingevoerde engineeringparameters, zodat het evenwicht nauwkeuriger kan worden berekend voordat het experiment start.

2. Methodologie

De auteurs hebben een gekoppelde workflow ontwikkeld, genaamd Kinetic Equilibrium Prediction (KEP), die twee codes combineert:

1. RAPTOR (Transport Code):

   • Een snelle 1.5D transportcode die de evolutie van stroomdichtheid, elektronen- en ionentemperatuur en deeltjesdichtheid simuleert.
   • Voorspellingsmodel: Gebruikt een op gradiënten gebaseerd transportmodel (gradient-based model) dat het gedrag van het plasma in zowel L-mode als H-mode simuleert.
   • Nieuwe ontwikkelingen:
     • Uitbreiding van de Martin-schaling voor de L-H overgang, aangepast voor gunstige en ongunstige magnetische configuraties (zoals de richting van de ionen-drift).
     • Aanpassing voor Negatieve Triangulariteit (NT) plasmas, die bekend staan om hun hoge prestaties in L-mode zonder over te gaan naar H-mode.
     • Voorspelling van de lijngemiddelde dichtheid ($n_{e,l}$) op basis van gasinvoerprogramma's (hoewel dit nog een beperking is).
   • Snelheid: Geoptimaliseerd voor snelle simulaties (compilatie naar C-code), waardoor volledige puls-simulaties in minuten kunnen worden uitgevoerd.

2. FBT (Inverse Evenwicht Code):

   • Een statische free-boundary solver die gebruikt wordt om de stromen in de poloidale veldspoelen te berekenen die nodig zijn om een gewenste plasma-vorm en positie te behouden.
   • FBT lost de Grad-Shafranov vergelijking op, waarbij de interne profielen ($p'$ en $TT'$) als randvoorwaarden dienen.

3. De Koppeling (KEP Workflow):

   • Stap 1: RAPTOR voert een "predict-first" simulatie uit op basis van het pulse schedule (stroom, verwarming, vorm) om de tijdsevolutie van $p'$ en $TT'$ te voorspellen.
   • Stap 2: Deze voorspelde profielen worden doorgegeven aan FBT.
   • Stap 3: FBT berekent de benodigde spoelstromen ($I_a$) om de gewenste vorm te behouden met deze realistische interne profielen.
   • Iteratie: Het proces wordt een paar keer herhaald (meestal 2-3 iteraties) totdat de oplossingen van RAPTOR en FBT convergeren (verschil in $\beta_{pol}$ en $l_i$ < 5%).

3. Belangrijkste Bijdragen

• Integratie van RAPTOR en FBT: Eerste toepassing van een gekoppelde transport-evenwicht workflow voor volledige TCV-discharges in de shot preparation fase.
• Uitbreiding van H-mode voorspelling: Implementatie van een aangepaste Martin-schaling die rekening houdt met de divertor-geometrie (single-null, double-null) en de richting van de ionen-drift, inclusief een specifieke behandeling voor Negatieve Triangulariteit (NT) plasmas.
• Validatie op grote schaal: Benchmarking van het model op 211 experimenten (shot range 82000-83000), inclusief diverse vormen (LSN, USN, Snowflake) en modi (L-mode, H-mode, NT).
• Experimentele validatie: Toepassing van de nieuwe voorspellingen op twee specifieke experimenten (#81882 en #83575) om de verbetering in vormcontrole te demonstreren.

4. Resultaten

• Voorspellingsnauwkeurigheid:
  • De totale energie (elektronen en ionen) wordt binnen een marge van 20% voorspeld voor de meeste shots.
  • De timing van de L-H overgang wordt correct voorspeld binnen een venster van 100ms voor 48 van de 51 geteste H-mode shots.
  • Voor NT-plasmas wordt de hoge L-mode confinement goed gereproduceerd, terwijl de overgang naar H-mode correct wordt onderdrukt (zoals verwacht voor NT).
• Invloed op Evenwichtsberekening:
  • De koppeling leidt tot significante correcties in de voorspelde poloidale veldspoelstromen. Correcties van enkele tientallen tot honderden ampères zijn mogelijk, afhankelijk van de fase van de puls en de nauwkeurigheid van de initiële aannames.
  • De voorspelde interne inductie ($l_i$) en poloidale beta ($\beta_{pol}$) komen beter overeen met de fysica dan bij traditionele methoden.
• Experimentele Uitkomsten:
  • Shot #81882 (Standaard H-mode): De met KEP voorbereide shot (#83740) toonde een betere uitlijning van het X-punt met het doelwit vergeleken met een shot voorbereid met standaard FBT-instellingen.
  • Shot #83575 (Negatieve Triangularity met Snowflake): NT-plasmas zijn bekend om hun moeilijkheid in verticale stabiliteit. De KEP-bereiding resulteerde in een aanzienlijk stabielere vorm, waarbij de primaire en secundaire X-punten langer dicht bij hun doelwaarden bleven. Dit bevestigt dat nauwkeurige voorspelling van $l_i$ en $\beta_{pol}$ cruciaal is voor complexe vormen.

5. Betekenis en Toekomstperspectief

Deze studie toont aan dat het koppelen van snelle transportvoorspellingen aan evenwichtsberekeningen de kwaliteit van de shot preparation op TCV aanzienlijk verbetert.

• Betere Stabiliteit: Door realistischere schattingen van de interne profielen te gebruiken, kunnen operators de spoelstromen nauwkeuriger instellen, wat het risico op verticale instabiliteiten en vormmisalignement verlaagt.
• Efficiëntie: De workflow is snel genoeg om gebruikt te worden in de voorbereidingstijd tussen twee experimenten (typisch 5-15 minuten).
• Toekomst: De auteurs zien mogelijkheden voor verdere verbetering, zoals het integreren van geavanceerdere dichtheidsvoorspellingsmodellen, het meenemen van stralingsverliezen (impurity transport), en het ontwikkelen van een "tight coupling" waarbij verwarmings- en stroomdrijfsystemen (zoals NBI en ECRH) direct gekoppeld worden aan de transport- en evenwichtsberekening in één iteratieve loop.

Kortom, de KEP-workflow biedt een robuust, fysica-gebaseerd instrument dat de brug slaat tussen theoretische simulatie en praktische tokamak-operatie, wat essentieel is voor de voorbereiding van toekomstige experimenten op grotere machines zoals ITER.