A new model for runaway electron transport based on chaotic Hamiltonian systems

Technische Samenvatting: Een nieuw model voor transport van runaway elektronen gebaseerd op chaotische Hamiltoniaanse systemen

Probleemstelling
Het transport van runaway elektronen (RE's) in ergodische magnetische geometrieën, met name tijdens tokamak-disrupties, vormt een kritieke uitdaging voor gecontroleerde kernfusie. Hoewel de generatie van RE-bundels via avalanche-vermenigvuldiging goed begrepen is, blijft het overleven van de initiële zaadpopulatie tijdens de quench-fase onzeker. Huidige modellering leunt vaak op de Rechester-Rosenbluth (RR) diffusiebenadering om transport in stochastische velden te beschrijven. Echter, het RR-model gaat uit van een volledig afgesloten ergodisch volume en puur diffusief (exponentieel) verval. Simulatiestudies geven aan dat in open ergodische volumes of complexe magnetische topologieën de RR-benadering ontoereikend is. Specifiek houdt het er geen rekening mee dat er "sticky regio's" (plakkerige regio's) zijn—trappinglagen rond magnetische eilanden waar de ontsnapping van deeltjes een power-law verval volgt in plaats van een exponentieel verval. Deze discrepantie beperkt de nauwkeurigheid van gereduceerde simulaties die worden gebruikt voor het interpreteren van experimentele data en het valideren van RE-modellen.

Methodologie
De auteurs stellen een semi-analytisch transportmodel voor, afgeleid van de theorie van chaotische Hamiltoniaanse systemen. De methodologie berust op de correspondentie tussen de topologie van magnetische veldlijnen en de faseruimte van dynamische systemen.

1. Theoretisch Kader: Het model behandelt gemengde faseruimten die zowel chaotische zeeën als magnetische eilanden bevatten. Het stelt dat de overlevingskans van deeltjes, $P(t)$, een superpositie is van twee verschillende vervalmechanismen die gelijktijdig optreden:
   • Exponentieel Verval: Gedomineerd door deeltjes in de chaotische zee, beschreven door $e^{-\kappa t}$, waarbij $\kappa$ de ontsnappingssnelheid is.
   • Power-Law Verval: Gedomineerd door deeltjes die gevangen zitten in "sticky regio's" (restanten van gebroken eilanden), beschreven door een term proportioneel aan $(1 + (t/\tau)^\alpha)^{-1}$.
2. Modelformulering: De totale overlevingskans wordt geformuleerd als:
   $$P(t) = A e^{-\kappa t} + \frac{1-A}{1 + (t/\tau)^\alpha}$$
   waarbij $A$ het initiële fractie deeltjes in het chaotische regime vertegenwoordigt, $\tau$ de onderscheidende tijdschaal voor het power-law verval is, en $\alpha$ de power-law exponent is. In tegen tegenstelling tot eerdere modellen die specifieke tijdsintervallen aan deze regimes toewijzen, gaat dit model ervan uit dat beide termen op elk moment aanwezig zijn.
3. Validatiegevallen: Het model wordt getest tegen twee verschillende scenario's:
   • Ullmann-Caldas Map: Een vereenvoudigde 2D-mapping die magnetische veldlijnen in een TBR-1 tokamak-achtige geometrie representeert.
   • JET Disruption Simulatie: Een complexe 3D niet-lineaire magnetohydrodynamische (MHD) simulatie van een JET-disruptie (ontlading #95135) met behulp van de JOREK-code, gevolgd door relativistische deeltjestracing.

Belangrijkste Resultaten

• Ullmann-Caldas Map: Het model heeft de numerieke overlevingskans van 300.000 deeltjes succesvol gefit. De fit onthulde een initiële exponentiële verliesfase gedomineerd door de chaotische zee, gevolgd door een langdurige staart die wordt beheerst door het power-law verval. De crossover-tijd ($t_x$) tussen deze regimes werd berekend op ongeveer 158 $\mu$s. De onderscheidende tijdschalen werden bevestigd: de exponentiële vervaltijd ($1/\kappa \approx 45,45$ $\mu$s) was aanzienlijk verschillend van de power-law tijdschaal ($\tau \approx 22,38$ $\mu$s).
• JET Disruption Scenario: Toegepast op een JOREK-simulatie van een JET-disruptie met 1,8 miljoen tracer-deeltjes, bereikte het model een "opmerkelijk goede fit" bij de transportcurve. De simulatie toonde een duidelijke chaotische zee met prominente sticky regio's rond eilandketens. Het model legde de overgang van snel initiële verlies naar de tragere power-law staart accuraat vast, met een crossover-tijd van ongeveer 125 $\mu$s.
• Parameterconsistentie: In beide gevallen waren de gefitte parameters ($A$, $\kappa$, $\tau$, $\alpha$) fysiek consistent, met $\alpha$-waarden nabij 1, wat de theoretische afleiding ondersteunt.

Betekenis en Claims
Het artikel claimt dat het voorgestelde model een significante verbetering biedt ten opzichte van de Rechester-Rosenbluth-benadering door expliciet rekening te houden met de effecten van sticky regio's in gemengde faseruimten. De primaire bijdrage is de functionele vorm die simultaan exponentieel en power-law verval meeneemt met onderscheidende tijdschalen, in plaats van deze als opeenvolgende of wederzijds uitsluitende regimes te behandelen.

De auteurs demonstreren dat dit eenvoudige transportmodel accuraat de temporele evolutie van de deeltjesconfinement-kansen kan beschrijven in zowel geïdealiseerde mapping-systemen als complexe, realistische 3D MHD-simulaties. Door een betere fit te bieden bij directe deeltjestraject-simulaties, biedt het model een robuuster instrument voor het integreren van transporteffecten in gereduceerde simulaties van runaway elektron-dynamica. Deze verbeterde nauwkeurigheid wordt gepresenteerd als cruciaal voor de interpretatie van huidige RE-experimenten en voor de validatie van modellen bedoeld voor toekomstige apparaten zoals ITER, waar het begrijpen van het transport van de zaadpopulatie essentieel is voor het voorspellen van RE-bundelvorming en -mitigatie.