Modelling of NSTX hot vertical displacement events using M3D-C1

The main results of an intense vertical displacement event (VDE) modelling activity using the implicit 3D extended MHD code M3D-C1 are presented. A pair of nonlinear 3D simulations are performed using realistic transport coefficients based on the reconstruction of a so-called NSTX frozen VDE where the feedback control was purposely switched off to trigger a vertical instability. The vertical drift phase is solved assuming axisymmetry until the plasma contacts the first wall, at which point the intricate evolution of the plasma, decaying to large extent in force-balance with induced halo/wall currents, is carefully resolved via 3D nonlinear simulations. The faster 2D nonlinear runs allow to assess the sensitivity of the simulations to parameter changes. In the limit of perfectly conducting wall, the expected linear relation between vertical growth rate and wall resistivity is recovered. For intermediate wall resistivities, the halo region contributes to slowing the plasma down, and the characteristic VDE time depends on the choice of halo temperature. The evolution of the current quench and the onset of 3D halo/eddy currents are diagnosed in detail. The 3D simulations highlight a rich structure of toroidal modes, penetrating inwards from edge to core and cascading from high-n to low-n mode numbers. The break-up of flux-surfaces results in a progressive stochastisation of field-lines precipitating the thermalisation of the plasma with the wall. The plasma current then decays rapidly, inducing large currents in the halo region and the wall. Analysis of normal currents flowing in and out of the divertor plate reveals rich time-varying patterns. © 2018 Author(s).