Effect of convective transport in edge/SOL plasmas of ADITYA-U tokamak

Met behulp van de UEDGE-code is aangetoond dat er in de randregio van de ADITYA-U tokamak naast diffusie ook een constante inwaartse convectiesnelheid nodig is om de gemeten elektronen dichtheidsprofielen nauwkeurig te modelleren.

De kern

De "Heetste Rand van de Zon" in een Magnetische Kooi: Een Uitleg

Stel je voor dat je probeert een gloeiend hete, vloeibare zon in een kamer te houden. Het probleem? De zon is zo heet dat hij alles wat hij aanraakt direct laat smelten. Je kunt de zon dus niet in een gewone bak of pot doen. Wat doe je dan? Je bouwt een "magnetische kooi" (een tokamak) om de hete deeltjes in de lucht te laten zweven, zodat ze de wanden van de machine niet raken.

In dit onderzoek kijken wetenschappers naar de "rand van de kooi" (de edge/SOL region). Dit is de gevaarlijkste zone: het is de plek waar de magnetische kooi de wand van de machine (de limiter) bijna raakt.

Het Probleem: De Onzichtbare Wind

De onderzoekers gebruikten een supercomputerprogramma genaamd UEDGE om te simuleren wat er precies gebeurt in die randzone. Ze wilden weten hoe de deeltjes (elektronen) zich verplaatsen.

In de wetenschap denken we vaak dat deeltjes zich gedragen als mensen in een drukke winkelstraat: ze verspreiden zich gewoon een beetje alle kanten op (dat noemen we diffusie). Als je een druppel inkt in water gooit, verspreidt die zich langzaam. Dat is wat de onderzoekers eerst dachten dat er gebeurde.

Maar toen kwamen ze erachter: de simulatie klopte niet met de werkelijkheid. De deeltjes bewogen niet alleen "een beetje alle kanten op", er was iets anders aan de hand.

De Ontdekking: De "Magnetische Stofzuiger"

De onderzoekers ontdekten dat er naast die normale verspreiding ook een convectie aanwezig is.

De Metafoor:
Stel je voor dat je in een drukke winkelstraat loopt (diffusie). Je loopt een beetje willekeurig rond. Maar stel je nu voor dat er plotseling een enorme, krachtige stofzuiger aan de rand van de straat wordt aangezet die iedereen met een constante kracht naar binnen zuigt. Dat is wat er in de plasma-rand gebeurt!

Deze "magnetische stofzuiger" (de convectieve snelheid) trekt de deeltjes met een snelheid van 1,5 meter per seconde naar binnen, richting het midden van de plasma-wolk. Zonder deze "zuigkracht" in hun computerprogramma, konden ze de metingen van de echte machine niet verklaren.

Waarom is dit belangrijk?

Waarom maken we ons druk om een beetje "wind" of "zuigkracht" aan de rand?

1. Bescherming van de machine: Als we precies weten hoe de hitte en de deeltjes naar de rand bewegen, kunnen we voorkomen dat de machine zichzelf kapot smelt. Het is alsowen we de thermostaat van een oven heel precies moeten afstellen om te voorkomen dat de binnenkant verbrandt.
2. De weg naar schone energie: Tokamaks zijn de sleutel tot fusie-energie (de energie van de zon). Als we de rand van de plasma-wolk perfect kunnen beheersen, kunnen we in de toekomst veilige, onbeperkte schone energie opwekken.

Samenvatting in drie zinnen:

Wetenschappers hebben een digitale kopie gemaakt van een magnetische plasma-machine om te zien hoe deeltjes zich aan de rand gedragen. Ze ontdekten dat de deeltjes niet alleen willekeurig rondzweven, maar ook door een soort constante "magnetische wind" naar binnen worden getrokken. Deze ontdekking helpt ons om toekomstige energiecentrales beter te bouwen en te beschermen.

Technische samenvatting

Technische Samenvatting: Effect van convectief transport in de edge/SOL-plasma's van de ADITYA-U tokamak

1. Probleemstelling

Het begrijpen van de plasma-edge en de Scrape-Off Layer (SOL) is cruciaal voor de werking van tokamaks, omdat dit gebied de hete kern scheidt van de fysieke wanden van de machine. In deze regio bepalen transportprocessen (het verplaatsen van deeltjes en energie) hoe de wanden worden belast en hoe efficiënt de kern wordt geïsoleerd.

Bij de ADITYA-U tokamak wordt gebruikgemaakt van een limiter-configuratie (in plaats van een divertor). De uitdaging is dat bestaande simulatiecodes, zoals UEDGE, primair zijn ontwikkeld voor divertor-geometrieën met magnetische null-punten (X-punten). Er was behoefte aan een nauwkeurige modellering van de deeltjes- en energietransporten in de specifieke limiter-geometrie van ADITYA-U om de experimentele metingen (via Langmuir-probes) te kunnen verklaren.

2. Methodologie

De onderzoekers hebben de volgende stappen ondernomen:

• Code-modificatie: De 2D-vloeistoftransportcode UEDGE is aangepast om limiter-geometrieën te ondersteunen.
• Mesh-generatie: Er is een nieuwe routine ontwikkeld in MATLAB/Python die gebruikmaakt van magnetische fluxoppervlakken (verkregen via de IPREQ-code) om een computationeel rooster (mesh) te genereren dat specifiek is voor de ADITYA-U limiter-configuratie. Dit rooster is niet-uniform, met een zeer fijne resolutie nabij de limiter-tip om de scherpe gradiënten in dichtheid en temperatuur op te vangen.
• Simulatieparameters: De simulaties werden uitgevoerd met behulp van een Krylov-Newton solver om een steady-state oplossing te bereiken. De randvoorwaarden werden ingesteld op basis van de bekende plasma-eigenschappen van de ADITYA-U (zoals plasma-stroom en dichtheid).
• Vergelijking: De gesimuleerde profielen van de elektronendichtheid ($n_e$) en temperatuur ($T_e$) werden vergeleken met experimentele data van Langmuir-probes.

3. Belangrijkste Bijdragen

• Geometrische integratie: Het succesvol integreren van een limiter-specifieke mesh in de UEDGE-code, wat een belangrijke stap is voor het modelleren van machines die niet in divertor-modus werken.
• Identificatie van transportmechanismen: Het aantonen dat louter diffusief transport (het verspreiden van deeltjes) onvoldoende is om de experimentele data te verklaren; er is expliciet een convectief component nodig.
• Validatie van transportcoëfficiënten: Het kwantificeren van de effectieve diffusiviteit en convectiesnelheid in de edge-regio van een limiter-tokamak.

4. Resultaten

• Convectie en Diffusie: Om de gemeten elektronendichtheidsprofielen nauwkeurig te reproduceren, is een combinatie nodig van:
  • Een constante radiale inwaartse convectiesnelheid ($v_{conv}$) van 1,5 m/s.
  • Een constante perpendiculaire diffusiecoëfficiënt ($D_{\perp}$) van ongeveer 0,2 m²/s.
• Vergelijking met andere modellen: De gevonden $D_{\perp}$ (~0,2 m²/s) is groter dan de neoclassicistische diffusiviteit, maar aanzienlijk kleiner dan de Bohm-diffusiviteit en de diffusiviteit veroorzaakt door fluctuaties (zoals ITG-modi of resistieve ballooning-modi). Dit suggereert dat fluctuatie-geïnduceerd transport in dit specifieke gebied van ADITYA-U minder dominant is.
• Warmteflux: De onderzoekers stelden vast dat de radiale elektronwarmteflux maximaal is nabij de limiter-tip. De geleidingsflux (conductive heat flux) domineert in de regio dieper in de edge, maar de convectieve warmteflux draagt maximaal bij aan de totale warmteoverdracht nabij de limiter-tip.
• Neutrale deeltjes: De gesimuleerde neutrale waterstofdichtheid vertoonde een goede overeenstemming met de verwachte fysica, waarbij de dichtheid snel afneemt naarmate men de limiter-tip nadert.

5. Betekenis en Conclusie

Dit onderzoek bewijst dat voor een correcte modellering van de plasma-rand in limiter-tokamaks, zowel diffusie als convectie (zoals de Ware-pinch) meegenomen moeten worden. De succesvolle aanpassing van de UEDGE-code biedt een krachtig instrument voor toekomstig onderzoek naar de interactie tussen plasma en wanden. Bovendien vormt deze methodologie een belangrijke basis voor het koppelen van UEDGE aan de DEGAS2-code, wat essentieel is voor het simuleren van complexe neutrale deeltjesdynamica in toekomstige fusie-experimenten.