Overview
Research introduces a transport model for runaway electrons (RE) described as applicable to chaotic Hamiltonian systems. This model extends beyond conventional diffusion approximations, such as the Rechester-Rosenbluth model, by incorporating specific effects observed in particle transport within ergodic magnetic geometries. A key distinction of this new model is its consideration of 'sticky regions' that form around magnetic islands.
Research Context
The transport of runaway electrons in ergodic magnetic geometries is an active area of study. Direct simulation of particle trajectories to compute this transport is computationally intensive. To address this, diffusion models are frequently employed to integrate transport effects into reduced simulations. One such model is the Rechester-Rosenbluth approximation.
However, comparisons between diffusion-based simulations and direct simulations indicate that the observed transport is not consistently purely diffusive. This discrepancy suggests a need for transport models that can account for non-diffusive behaviors, particularly in the context of chaotic systems.
Approach
The proposed model is rooted in chaos theory. It aims to improve upon the Rechester-Rosenbluth approximation by including the influence of 'sticky regions'. These sticky regions are defined as trapping layers that surround magnetic islands. Within these regions, the mechanism of particle escape deviates from the exponential decay characteristic of purely diffusive processes, exhibiting instead a power-law decay.
To assess its applicability, the model was tested in two distinct scenarios:
- The Ullmann-Caldas map, using parameters relevant to the TBR-1 tokamak.
- A JOREK simulation, specifically modeling a JET disruption scenario.
Findings
The new transport model demonstrates the ability to account for both chaotic diffusion and the effects of sticky regions. A primary finding is the observation that particle escape within these sticky regions follows a power-law decay, in contrast to the exponential decay implied by purely diffusive models. This characteristic directly impacts the representation of slower particle escape from these trapped layers.
The applicability of the model was tested and confirmed in both the Ullmann-Caldas map, under conditions corresponding to the TBR-1 tokamak, and within a JOREK simulation simulating a JET disruption. In both evaluated scenarios, the model achieved good fits.
Why This Matters
The development of this model addresses the limitations of purely diffusion-based approximations for runaway electron transport, which often fail to fully capture the complexities of transport in ergodic magnetic geometries. By incorporating phenomena like sticky regions and their associated power-law decay of particle escape, the model contributes to a more accurate representation of runaway electron behavior in fusion devices. This enhanced understanding is relevant for the study of plasma physics and fusion energy research.