Non-linear dynamics of the double tearing mode

3D non-linear magnetohydrodynamics simulations of a double tearing mode with the JOREK code are presented in the context of trying to better understand the benign termination of runaway electron beams observed in some experiments. It is shown that the non-linear behavior qualitatively depends on the resistivity η via its effect on how fast secondary, non-linearly destabilized, tearing modes grow relative to the primary mode. Within a certain range of η, a violent and global relaxation is observed, consistent with the “Kadomtsev-predicted” reconnection region extending from almost the very center up to the edge of the plasma.