Resistive reduced MHD modeling of multi-edge-localized-mode

The full dynamics of a multi-edge-localized-mode (ELM) cycle is modeled for the first time in realistic tokamak X-point geometry with the nonlinear reduced MHD code JOREK. The diamagnetic rotation is found to be instrumental to stabilize the plasma after an ELM crash and to model the cyclic reconstruction and collapse of the plasma pressure profile. ELM relaxations are cyclically initiated each time the pedestal gradient crosses a triggering threshold. Diamagnetic drifts are also found to yield a near-symmetric ELM power deposition on the inner and outer divertor target plates, consistent with experimental measurements.