Researchers have developed a new hybrid kinetic-magnetohydrodynamic (MHD) simulation to model the dynamics of runaway electrons (REs) during tokamak disruptions. The model integrates a full-f particle-in-cell (PIC) treatment of the RE population with the 3D nonlinear extended MHD framework of the NIMROD code. This approach aims to provide a more physically complete picture of the formation of the RE current plateau, a critical phase in a disruption where the plasma current is converted into a beam of highly energetic electrons. The work, detailed in a preprint on arXiv, focuses on capturing kinetic effects that are often simplified or omitted in purely fluid-based models, which is essential for predicting RE behavior in high-current machines like ITER. Source: arXiv

The simulation advances the RE population by tracking their motion along guiding-center (GC) orbits, while RE generation is handled through analytical source terms. This hybrid methodology allows the model to capture the kinetic nature of the relativistic electron beam while still leveraging the robust and widely-used MHD capabilities of NIMROD for the bulk plasma. The coupling provides a self-consistent framework where the evolution of the background plasma and the RE beam are interdependent. The model's performance was benchmarked against established fluid-based RE codes, demonstrating strong agreement when the RE guiding-center drifts were intentionally excluded, thereby validating its foundational physics implementation. Source: arXiv

The simulation advances the RE population by tracking their motion along guiding-center (GC) orbits, while RE generation is handled through analytical source terms.

A key finding from the simulation is the significant role of guiding-center drifts for highly relativistic electrons. The model shows that for these electrons, the grad-B and curvature drifts, which are often neglected, become important factors in their generation and subsequent motion. This effect is amplified by the evolution of the safety factor profile during the current quench phase of a disruption. As the plasma current collapses and is replaced by the RE beam, the magnetic field structure changes, making these drift effects more pronounced. This insight highlights a limitation of purely fluid descriptions and underscores the necessity of a kinetic treatment for accurately predicting the spatial distribution and confinement of the RE beam. Source: arXiv

The implications of this work are directly relevant to the operational safety of future high-current tokamaks. Unmitigated runaway electron beams can deposit immense localized energy onto plasma-facing components, causing significant damage. Developing effective RE mitigation strategies, such as massive material injection, depends on predictive models that can accurately capture the formation, energy spectrum, and spatial evolution of the RE plateau. By incorporating previously overlooked kinetic drift effects, this new NIMROD-based model offers a higher-fidelity tool for designing and validating these critical mitigation systems. Further development could see the model used to simulate specific disruption scenarios and test the efficacy of various intervention techniques. Source: arXiv