Abstract:Recently, the collisionless pitch-angle scattering for relativistic runaway electrons (REs) in toroidal geometries such as tokamaks was discovered through a full orbit simulation approach (Liu et al., Nucl. Fusion, vol. 56, 2016, p. 064002), and it was then theoretically investigated that a new expression for the magnetic moment, including the second-order corrections, could essentially reproduce the so-called collisionless pitch-angle scattering process (Liu et al., Nucl. Fusion, vol. 58, 2018, p. 106018). In this paper, with synchrotron radiation, extensive numerical verification of the validity of the high-order guiding-centre theory is given for simulations involving REs by incorporating such an expression for the magnetic moment into our particle tracing code. A high-order guiding-centre simulation approach with synchrotron radiation (HGSA) is applied. Synchrotron radiation plays an essential role in the life cycle of REs. The energy of REs first increases and then becomes saturated until the electric field acceleration is balanced by the radiation dissipation. Unfortunately, the process cannot be simulated accurately with the standard guiding-centre model, i.e. the first-order guiding-centre model. Remarkably, it is found that the HGSA can effectively produce the fundamental process of REs. Since the time scale of the energy saturation of REs is close to seconds, the computational cost becomes significant. In order to save costs, it is necessary to estimate the time of energy saturation. An analytical estimate is derived for the time it takes for synchrotron drag to balance an accelerating electric field and the provided formula has been numerically verified. Test calculations reveal that HGSA is favourable for exploiting the dynamics of REs in tokamak plasmas.

Higher order Larmor radius corrections to guiding-centre equations and application to fast ion equilibrium distributions

VENUS-LEVIS and its spline-Fourier interpolation of 3D toroidal magnetic field representation for guiding-centre and full-orbit simulations of charged energetic particles

Monte Carlo Orbit-Following Simulations Including the Finite Larmor Radius Effect Based on a Phase-Space Coordinate Transform Method

Impact of RMP magnetic field simulation models on fast ion losses

Numerical integration of particle orbits in discontinuous fields using VENUS-LEVIS and SPEC

Energetic particle tracing in optimized quasi-symmetric stellarator equilibria

Energetic Particle Tracing in Optimized Quasisymmetric Stellarator Equilibria

Equilibrium ion distribution in the presence of clearing electrodes and its influence on electron dynamics

Accurate Transfer Maps for Realistic Beamline Elements: Part I, Straight Elements

Problem and solution with the longitudinal tracking of the ORBIT code

Guiding center and gyrokinetic orbit theory for large electric field gradients and strong shear flows

Variational Formulation of Higher-order Guiding-center Vlasov-Maxwell Theory

Calculation of Collisionless Pitch-Angle Scattering of Runaway Electrons with Synchrotron Radiation Via High-Order Guiding-Centre Equation

From the orbit theory to a guiding center parametric distribution function

Conservative magnetic moment of runaway electrons and collisionless pitch-angle scattering

Improved boundary conditions for coupled geospace models: An application in modeling spacecraft surface charging environment

Electron acceleration driven by the lower-hybrid-drift instability: an extended quasilinear model

Guiding-center transformation of the radiation-reaction force in a nonuniform magnetic field

Higher-order moment convergent method in weakly anisotropic plasma and the NLVFP code for solution of the 0D-2V Vlasov-Fokker-Planck equation

Optimizing Beam-Ion Confinement in ITER by Adjusting the Toroidal Phase of the 3D Magnetic Fields Applied for ELM Control

Kinetic Corrections to Heat-flow and Nernst advection for Laser Heated Plasmas