Towards iterative 2D and 3D kinetic modelling of RMP interaction with tokamak plasmas

Suppressing edge-localized modes (ELMs) with help of external resonant magnetic field perturbations (RMPs) is an experimentally proven technique for the reduction of tokamak first wall damage by large, type-I ELMs. Linear and nonlinear modelling of the interaction of RMPs with tokamak plasma is presently developing on the basis of MHD and kinetic theory [1]. In this report, we extend an iterative approach proposed in Ref. [2] and realized in the code MEPHIT for the ideal MHD plasma response in Ref. [3] to the approach combining the ideal MHD and kinetic plasma response. In this approach, plasma response current density computed for the given perturbation magnetic field and computations of this field by a finite element Maxwell solver in realistic tokamak geometry with a given response current coupled together within an
iterative procedure.