3D modelling of thermal loads during unmitigated vertical displacement events in ITER and JET

Deze paper presenteert een nieuwe physics-gebaseerde workflow voor het voorspellen van driedimensionale thermische belastingen tijdens verticale verplaatsingsgebeurtenissen, die is gevalideerd met JET-gegevens en wordt toegepast om de veerkracht van de ITER-wolfraamwand tegen disrupties te beoordelen.

De kern

Stel je voor dat je een supermoderne, gigantische kookpan hebt (de ITER-tokamak) die bedoeld is om de zon op aarde na te bootsen. Deze pan moet extreem heet worden om energie op te wekken. Maar er is een probleem: soms gaat het "koken" helemaal mis. De vlammen slaan plotseling alle kanten op en de pan raakt uit balans. Dit noemen wetenschappers een disruptie.

Dit wetenschappelijke artikel gaat over hoe we kunnen voorspellen wat er gebeurt als die "pan" uit balans raakt en de enorme hitte tegen de wanden van de machine slaat.

Hier is de uitleg in begrijpelijke taal:

1. Het probleem: De "Ongeleide Raket" (Vertical Displacement Event)

Normaal gesproken zweeft het gloeiend hete plasma (de brandstof) netjes in het midden van de machine, zwevend op magnetische velden. Maar bij een bepaald type fout — een Vertical Displacement Event (VDE) — verliest het plasma zijn evenwicht.

Zie het als een ongeleide raket die plotseling uit koers raakt en recht op de wand van de machine af stort. In plaats van dat de hitte netjes wordt verdeeld, slaat de energie als een hittegolf op één specifieke plek tegen de wand aan. Als dat gebeurt, kan het materiaal van de wand (het pantser) letterlijk smelten.

2. Wat hebben de onderzoekers gedaan? (De Digitale Tweeling)

Omdat je niet zomaar een miljardenmachine als ITER kapot wilt laten smelten om te kijken wat er gebeurt, hebben de onderzoekers een supercomputer-simulatie gemaakt.

Ze hebben een soort "digitale tweeling" van de machine gebouwd. Ze gebruikten een programma genaamd JOREK. Dit programma is als een extreem geavanceerde weersvoorspelling, maar dan voor magnetische velden en hitte. Ze hebben eerst getest of hun digitale model klopte door het te vergelijken met echte experimenten in een kleinere machine (JET). Het werkte: de computer voorspelde precies waar het smelten gebeurde, net als in het echt.

3. De belangrijkste ontdekkingen

A. De 3D-factor: Het is geen platte cirkel
Oude modellen gingen ervan uit dat de hitte zich als een perfecte, gelijkmatige cirkel over de wand zou verspreiden (2D). Maar deze onderzoekers lieten zien dat de wereld 3D is. Door de vorm van de machine en de wiebelende magnetische velden, slaat de hitte als een laserstraal op bepaalde plekken, terwijl andere plekken relatief koel blijven. Dit noemen ze "lokalisatie".

B. Het "Twee-Stappen-Plan" van de hitte
De hitte komt in twee golven:

1. De Thermische Quench (De flits): Een plotselinge, korte flits van hitte.
2. De Stroom Quench (De langdurige hitte): Een langere periode waarin de magnetische energie wordt omgezet in hitte.

De onderzoekers ontdekten dat de eerste flits de wand "voorverwarmt" (zoals een pan die je eerst heet maakt), waardoor de tweede golf de wand veel makkelijker kan laten smelten. Het is de combinatie van de twee die de genadeslag geeft.

C. Goed nieuws voor ITER: Het nieuwe pantser
In het verleden gebruikten machines materiaal zoals beryllium. Maar voor de nieuwe ITER-plannen gebruiken ze wolfraam (tungsten). Wolfraam is een soort "super-metaal" dat een extreem hoog smeltpunt heeft. De simulaties laten zien dat dit nieuwe pantser veel beter bestand is tegen deze "raketinslagen" dan het oude materiaal.

Samenvatting in één beeld

Stel je voor dat je met een zaklamp in het donker schijnt. Een oude simulatie dacht dat de zaklamp een zachte, gelijkmatige gloed over de hele muur zou werpen. Deze nieuwe studie laat zien dat de zaklamp eigenlijk een fel, bewegend punt is dat op specifieke plekken de verf van de muur kan wegbranden.

Waarom is dit belangrijk?
Nu we weten waar en hoe hard die hitte gaat slaan, kunnen ingenieurs de wanden van de ITER-machine precies op de juiste plekken versterken. Zo kunnen we de zon op aarde temmen zonder dat de pan smelt!

Technische samenvatting

Technische Samenvatting: 3D-modellering van thermische belastingen tijdens ongedempte verticale verplaatsingsgebeurtenissen (VDE) in ITER en JET

1. Probleemstelling

Tijdens een disruptie in een tokamak (zoals ITER) treden twee kritieke fasen op: de Thermal Quench (TQ) en de Current Quench (CQ). Bij ongedempte verticale verplaatsingsgebeurtenissen (UVDE) gaat een groot deel van de plasma-energie verloren via parallel transport langs open veldlijnen, wat leidt tot extreem hoge lokale energiedichtheden op de plasma-facing components (PFC's).

Bestaande voorspellende modellen zijn vaak gebaseerd op axiale symmetrie (2D). Dit is problematisch omdat UVDE's gepaard gaan met grote externe kink-modi die de toroidale symmetrie doorbreken. Het ontbreken van 3D-modellen maakt het onmogelijk om de werkelijke lokalisatie van de hitteflux en het risico op smelten van de wand (bijv. wolfraam in de nieuwe ITER-baseline) nauwkeurig te voorspellen.

2. Methodologie

De auteurs presenteren een physics-gebaseerde workflow die uit drie gekoppelde stadia bestaat:

1. MHD-simulaties: Gebruik van de JOREK-code voor niet-lineaire resistieve magnetohydrodynamische (MHD) simulaties. Dit levert tijd-opgeloste 3D-magnetische velden en brontermen op (zoals parallelle hitteflux en halo-stromen).
2. Veldlijn-tracing: Een post-processor brengt de plasma-belastingen in kaart op een hoogwaardig 3D-model van de eerste wand (First Wall, FW) met behulp van veldlijn-tracing om te bepalen welke wandelementen daadwerkelijk "nat" (verbonden met het plasma) zijn.
3. Thermische respons: Voor elk wandelement wordt een eendimensionaal transient warmtegeleidingsprobleem opgelost om de piek-oppervlaktetemperatuur te berekenen en te beoordelen of smelten optreedt.

De methode is gevalideerd tegen experimentele data van JET (met beryllium-armouring) voordat deze is toegepast op ITER-scenario's (met wolfraam-armouring).

3. Belangrijkste Resultaten

Validatie in JET:

• De simulaties reproduceerden de globale dynamiek (plasma-stroom $I_p$, verticale verplaatsing van het stroomcentrum) en de toroidale asymmetrie van de halo-stromen uit de experimenten zeer nauwkeurig.
• Smeltgedrag: Voor de hoge-stroom JET-puls (#84832) voorspelde het model correct het smelten van beryllium. Een cruciale bevinding is dat smelten pas optreedt door de gecombineerde actie van de TQ (die de wand voorverwarmt) en de CQ (die de temperatuur verder verhoogt).

Voorspellingen voor ITER (15 MA scenario):

• Resilience (Veerkracht): De ITER wolfraam (W) wand is aanzienlijk veerkrachtiger dan de eerdere beryllium wand. Hoewel er lokaal smelten wordt voorspeld, is dit beperkt tot de randen van de panelen (bij de bovenste poorten) waar de veldlijnen bijna loodrecht op de wand vallen.
• 3D vs. 2D: Axiale symmetrische modellen onderschatten de lokale hittebelasting aanzienlijk. De 3D-asymmetrie zorgt voor een "peaking factor" van respectievelijk ~3,5 (tijdens TQ) en ~2 (tijdens CQ) in de temperatuurverhoging.
• Verbreding van de footprint: In tegenstelling tot 2D-modellen (zoals TOKES), die de hitte zeer lokaal voorspellen, laten de 3D-simulaties zien dat de hitteflux-footprint veel breder is. Dit komt doordat de 3D MHD-perturbaties de magnetische topologie vervormen, waardoor veldlijnen diep in de kern van het plasma doordringen.
• Proxy voor hitteflux: Er is een sterke correlatie tussen de verdeling van de halo-stroom en de verdeling van de afgevoerde energie, wat betekent dat halo-stroommetingen gebruikt kunnen worden als een proxy voor de thermische belasting.

4. Betekenis en Conclusie

Dit onderzoek levert een essentieel instrument voor het beheer van de "disruption budget" van toekomstige reactoren. De belangrijkste conclusies zijn:

• De overstap van beryllium naar wolfraam in de ITER-baseline verhoogt de veiligheidsmarge tegen smelten aanzienlijk.
• 3D-effecten zijn niet optioneel; ze zijn noodzakelijk om zowel de extreme lokalisatie (piektemperaturen) als de bredere spreiding (totale energieverdeling) van de thermische belasting correct te begrijpen.
• De breedte van de thermische belasting wordt primair bepaald door de magnetische verbinding met de plasma-kern, een inzicht dat cruciaal is voor het ontwerpen van beschermingssystemen.