W7-X and the sawtooth instability: towards realistic simulations of current-driven magnetic reconnection. (24th June 2021)

• Record Type: Journal Article
• Authors: Zocco, Alessandro; Mishchenko, Alexey; Nührenberg, Carolin; Könies, Axel; Kleiber, Ralf; Borchardt, Matthias; Slaby, Christoph; Zanini, Marco; Stange, Torsten; Laqua, Heinrich Peter; Rahbarnia, Kian; Thomsen, Henning; Wolf, R.C.; Helander, Per; Hatzky, Roman; Cole, Michael D. J.
• Is Part Of: Nuclear fusion. Volume 61:Number 8(2021)
• Publication Date: 2021-06-24
• Subjects: magnetic reconnection -- stellarator -- sawtooth instability; Nuclear fusion -- Periodicals; 621.48405
• DOI: 10.1088/1741-4326/ac077d
• Languages: English
• ISSNs: 0029-5515

Abstract

Abstract: Magnetic reconnection in W7-X is studied by means of global numerical simulations in a series of models of increasing complexity. The magnetic geometry ranges from that of a cylinder to the full three-dimensional field of W7-X, and the equations solved range from ideal magnetohydrodynamics (MHD) to gyrokinetics. We simulate plasmas from the first operation phase with electron cyclotron current drive (ECCD). These are characterized by an equilibrium magnetic field featuring an ECCD-distorted 'humped' profile of the rotational transform ι, with ι = 1 in two radial locations. Such plasmas generally show sawtooth activity, hence motivating the present study. We pay particular attention to the role of equilibrium current density gradients in the destabilization of reconnecting modes. When the equilibrium temperature and density gradients are artificially suppressed (to eliminate the pressure gradient drive), the perturbed electrostatic potential is radially localized between the locations at which ι = 1. This is shown with a purely collisionless gyrokinetic model, in cylindrical geometry. In the real toroidal geometry of W7-X, for a non-ideal MHD model including a uniform resistivity, electron inertia and (numerical) viscosity, the same qualitative behaviour is observed. In particular, even if a resonant ( m, n ) = (1, −1) perturbation is initialized, the most unstable mode is the ( m, n ) = (−4, 4), where m and n are the poloidal and toroidal mode numbers, respectively. Other modes are destabilized due to geometric coupling. The growth rate of this instability scales as η 1/3, where η is the plasma resistivity, thus suggesting that ECCD drives ideal MHD stable W7-X plasmas towards non-ideal marginality. An ideal magnetohydrodynamic analysis confirms the result. A fluid-kinetic hybrid version of the EUTERPE code shows that gyrokinetic ions have a stabilizing effect on these modes. For W7-X relevant collisionalities, the growth rate scales linearly with the electron skin depth, d e . Implications of our results for sawtoothing W7-X operation are discussed.