Modeling of plasma transport during edge-localized mode in tokamak using a kinetic Vlasov-Poisson code

Binnen de enorme, donutvormige machines die tokamak's worden genoemd, proberen wetenschappers de kracht van de sterren te recreëren om schone energie op te wekken. Deze apparaten houden superverhit gas, plasma genoemd, op zijn plaats met behulp van krachtige magnetische velden. Dit plasma is echter niet perfect kalm; het ondergaat af en toe plotselinge, gewelddadige uitbarstingen die bekend staan als edge-localized modes. Stel je het plasma voor als een kokende pan waarvan het oppervlak plotseling barst, waardoor een golf van intense hitte en deeltjes tegen de binnenwanden van de machine slaat. Deze wanden zijn bekleed met speciale materialen, vaak wolfraam, die ontworpen zijn om de hitte te weerstaan, maar herhaalde klappen kunnen ze doen smelten of stukjes ervan laten afbreken. Als deze geërodeerde stukjes terug naar de kern van het plasma drijven, kunnen ze de reactie vergiftigen en stilleggen. Begrijpen hoe precies deze uitbarstingen van het centrum van de machine naar de wanden reizen, is cruciaal voor het bouwen van een veilige, langdurige fusiereactor, maar de betrokken fysica is ongelooflijk complex, waarbij deeltjes met verschillende snelheden bewegen en op manieren interageren die eenvoudige vloeistofmodellen niet volledig kunnen vatten.

Om dit probleem aan te pakken, hebben onderzoekers Ce Wang, Sven Van Loo en Geert Verdoolaege aan de Universiteit Gent een nieuwe computersimulatietool genaamd KOBRA ontwikkeld. In plaats van het plasma te behandelen als een gladde vloeistof, wat een gebruikelijke maar soms onnauwkeurige benadering is voor deze hogesnelheidsgebeurtenissen, volgt hun code het gedrag van individuele deeltjes terwijl ze door de ruimte en tijd bewegen. Deze methode, bekend als kinetische modellering, stelt hen in staat om fijne details te zien, zoals hoe snel bewegende elektronen zich anders gedragen dan zwaardere ionen. Het specifieke probleem waarmee zij werden geconfronteerd, was een kwestie van schaal. De afstand die het plasma aflegt van het centrum van de machine naar de wanden is enorm, ongeveer tien tot twintig meter, terwijl de minuscule afstand waarover elektrische krachten werken slechts ongeveer de breedte van een menselijke haar is. Het simuleren van elk klein detail over deze enorme afstand met een standaard, uniform rooster zou een computer zoveel berekeningen laten uitvoeren dat het onbetaalbaar en traag zou zijn.

Het team loste dit op door een techniek te gebruiken die adaptive mesh refinement wordt genoemd. Denk hierbij aan een digitale camera die automatisch inzoomt om een foto met hoge resolutie te maken op de plek waar de actie plaatsvindt, terwijl de achtergrond met een lagere resolutie behouden blijft om ruimte te besparen. In hun simulatie wordt het computernetwerk zeer fijn en gedetailleerd in de regio's waar de deeltjesverdeling snel verandert, zoals de voorrand van de plasmaburst, en blijft het grover in de gladdere, kalmere regio's. Deze aanpak stelde hen in staat om een volledige simulatie van het plasmavervoer van de mid-plane van de machine naar de divertor targets uit te voeren, wat de specifieke platen zijn die ontworpen zijn om de hitte te absorberen. Door hun adaptieve methode te vergelijken met een traditionele simulatie die overal een hoogresolutie rooster gebruikte, ontdekten ze dat de nieuwe methode bijna identieke resultaten produceerde met betrekking tot de stroming van deeltjes en hitte, met verschillen van minder dan één procent.

De simulaties onthulden de specifieke mechanica van hoe deze uitbarstingen zich ontvouwen. In de allereerste momenten raast een kleine groep hoogenergetische elektronen voor de rest van het plasma uit, en bereikt de wand voordat de zwaardere ionen hen kunnen inhalen. Dit vroege aankomst creëert een tijdelijke scheiding van lading, wat een sterk elektrisch veld genereert. Dit veld werkt vervolgens als een rem op de snelle elektronen en een duw op de ionen, waardoor de hele plasmawolk snel van vorm verandert. De onderzoekers observeerden dat de deeltjesflux bij de wand scherp stijgt en daarna langzaam afneemt, een patroon dat overeenkomt met eerdere experimentele waarnemingen. Cruciaal was dat de adaptieve mesh-methode deze vroege piek en de daaropvolgende evolutie net zo nauwkeurig vastlegde als de veel zwaardere, uniforme rooster-methode, wat bewees dat de "inzoom"-strategie geen fysieke waarheid opoffert voor snelheid.

Naast nauwkeurigheid benadrukte de studie de dramatische efficiëntiewinsten van deze nieuwe aanpak. Hoewel de adaptieve simulatie aanvankelijk trager was omdat deze constant zijn rooster moest reorganiseren, werd de methode veel sneller naarmate het plasma tot rust kwam en de behoefte aan extreem detail afnam. Aan het einde van de simulatie voltooide de adaptieve methode de taak in ongeveer 68 uur, terwijl de uniforme rooster-methode 113 uur in beslag nam, waardoor de rekentijd bijna werd gehalveerd. Bovendien had de adaptieve methode aanzienlijk minder computergeheugen nodig, waarbij aan het einde van de run slechts ongeveer 18 procent van het geheugen van het uniforme model werd gebruikt. Deze efficiëntie is van vitaal belang omdat het de deur opent naar complexere, hogere-dimensie simulaties die voorheen onmogelijk uit te voeren waren. Het werk demonstreert dat door de resolutie van de simulatie intelligent aan te passen aan het fysieke gedrag van het plasma, wetenschappers deze gevaarlijke uitbarstingen met grotere snelheid en minder hulpbronnen kunnen modelleren, wat een duidelijk pad biedt naar het begrijpen en beperken van de schade die deze gebeurtenissen veroorzaken aan toekomstige fusiereactoren.