Verification and validation of linear gyrokinetic and kinetic-MHD simulations for internal kink instability in DIII-D tokamak

Verification and linear validation of the internal kink instability in tokamak have been performed for both gyrokinetic (GTC) and kinetic-MHD codes (GAM-solver, M3D-C1-K, NOVA, XTOR-K). Using realistic magnetic geometry and plasma profiles from the same equilibrium reconstruction of the DIII-D shot #141216, these codes exhibit excellent agreement for the growth rate and mode structure of the internal kink mode when all kinetic effects are suppressed. The simulated radial mode structures, obtained from linear simulations, are in reasonable agreement with the normalised electron cyclotron emission measurement after adjusting, within the experimental uncertainty, the safety factor q=1 flux-surface location in the equilibrium reconstruction. Compressible magnetic perturbations strongly destabilize the kink, while poloidal variations of the equilibrium current density reduce the growth rate of the kink. Furthermore, kinetic effects of thermal ions are found to decrease the kink growth rate in kinetic-MHD simulations, but increase the kink growth rate in gyrokinetic simulations, due to the additional drive of the ion temperature gradient and parallel electric field. Kinetic thermal electrons are found to have negligible effects on the internal kink instability.

New Identification Approach and Methods for Plasma Equilibrium Reconstruction in D-Shaped Tokamaks

The paper deals with the identification of plasma equilibrium reconstruction in D-shaped tokamaks on the base of plasma external magnetic measurements. The methods of such identification are directed to increase their speed of response when plasma discharges are relatively short, like in the spherical Globus-M2 tokamak (Ioffe Inst., St. Petersburg, Russia). The new approach is first to apply to the plasma discharges data the off-line equilibrium reconstruction algorithm based on the Picard iterations, and obtain the gaps between the plasma boundary and the first wall, and the second is to apply new identification methods to the gap values, producing plasma shape models operating in real time. The inputs for on-line robust identification algorithms are the measurements of magnetic fluxes on magnetic loops, plasma current, and currents in the poloidal field coils measured by the Rogowski loops. The novel on-line high-performance identification algorithms are designed on the base of (i) full-order observer synthesized by linear matrix inequality (LMI) methodology, (ii) static matrix obtained by the least square technique, and (iii) deep neural network. The robust observer is constructed on the base of the LPV plant models which have the novelty that the state vector contains the gaps which are estimated by the observer, using input and output signals. The results of the simulation of the identification systems on the base of experimental data of the Globus-M2 tokamak are presented.