Yanjun Liu

Fishbone modes and Alfvén eigenmodes (AEs) are two of the most significant core magnetohydrodynamic (MHD) events observed in advanced scenarios on the Experimental Advanced Superconducting Tokamak. Microwave reflectometry, with its excellent spatial and temporal resolution, has been utilized to study their behaviors. During the fishbone event, the density profile and its gradient ( ∇ n e ), measured by density profile reflectometry, in the vicinity of the fishbone mode are modulated by the rotation of internal kink modes. Furthermore, the position of the maximum density gradient ( ∇ n e,max ) aligns with the fishbone location, and its value exhibits a positive correlation with the amplitude of the internal kink modes. Additionally, the fishbone can be detected by monitoring the center of gravity oscillations in a conventional reflectometry called multi-channel poloidal correlation reflectometry (PCR). Furthermore, the spectrogram of density fluctuation amplitude, measured by the PCR, reveals AEs more distinctly than commonly used diagnostics such as electron cyclotron emission and high frequency Mirnov coils. These findings not only help to better understand the characteristics of these MHD behaviors, but also contribute to the q-profile reconstruction and investigating the related physical processes.

A module with self-consistent evolution of driven current is developed and coupled with the resistive-MHD equations in the three-dimensional, toroidal, and nonlinear simulation code (CLT). The driven current equation is solved with a second-order accuracy symmetric scheme, which exhibits good conservation properties. With the new module, we find that the driven current can self-consistently concentrate inside the magnetic island when the parallel diffusion of the driven current is sufficiently large. The efficiency of the driven current on tearing mode suppression will then be much higher than those with stationary distributions. With the new module, the influence of electron cyclotron current drive (ECCD) on the nonlinear evolution of the 2/1 double tearing modes (DTMs) is investigated. When co-ECCD deposits on the outer resonant surface, the local magnetic shear is reduced, and the growth rates of the DTMs decrease; if ctr-ECCD deposits on the outer resonant surface, the local magnetic shear increases, and the DTMs become more unstable. However, things will be different if ECCD deposits on the inner resonant surface since the local magnetic shear is negative. The co-ECCD deposited on the inner resonant surface increases the negative shear and then promotes the growth of the DTMs; while the ctr-ECCD suppresses the DTMs. It is also found that the off-axis and central pressure crashes associated with the 2/1 DTMs can be converted to each other by properly depositing the driven current. To convert a central crash to an off-axis crash, the co-ECCD should be deposited on the outer resonant surface, or the ctr-ECCD deposited on the inner resonant surface. While, the co-ECCD should be deposited on the inner rational surface, or the ctr-ECCD deposited on the outer rational surface to convert an off-axis crash to a central crash. The co- or ctr-ECCD should be larger than a threshold for such transitions, and the threshold value is mainly determined by the location of the inner resonant surface.

This paper reports the recent observation of a weakly coherent mode (WCM) within a conventional reflectometer on EAST and successfully determines its poloidal wavenumber range. During the transition from the L-mode to I-mode, the line-averaged density remains nearly unchanged while a significant change is observed in the electron cyclotron emission (ECE) signals at the boundary. The difference between the signals for the two channels at the edge increased, coinciding with the appearance of the WCM and a simultaneous rise in the boundary electron temperature. Further investigation unveiled the modulating role of edge temperature ring oscillation (ETRO) (Liu et al 2020 Nucl. Fusion 60 126016) on high-frequency density fluctuations. Statistical results unveil an inverse relationship between the centeral frequency of the WCM and q ₉₅. Simulation results provide additional insights, demonstrating that the simulated ‘WCM’ in the density fluctuations aligns with experimental data in terms of center frequency. Additionally, the radial distribution of the simulated ‘WCM’ closely corresponds to regions with the strongest electron temperature gradients. Finally, through a cross-correlation analysis of the simulated fluctuations, the following phase relationship for the wavenumber range of ‘WCM’ was observed: α T ~ e > α n ~ i ∼ α ϕ ~ > α T ~ i .

The improved energy confinement mode (I-mode) is widely considered as an important operation regime for ITER. I-mode implementation depends on the specified basic plasma parameters and certain operation conditions, which are discovered by statistical plasma characteristics from a large number of I-mode discharges on a tokamak. The extraction process of I-mode plasma characteristics is complicated, time-consuming, and limited to the sampling rate of the measured signals. Experimental observation of the I-mode is accompanied by the appearance of a weakly coherent mode (WCM). However, it takes much time to accurately scan and quantify WCM characteristics when analyzing many I-mode discharges. Recently, a neural network identification method was developed as an I-mode detector to traverse a whole database as a replacement for manual identification. Two fully connected neural network models were trained with the spectrum of propagation velocity of density perturbation from Doppler backward scattering and the electron density measured by a polarimeter-interferometer system with the experimental advanced superconducting tokamak I-mode database. An accuracy of 98.30% in identifying WCMs in I-mode discharges is achieved with the WCM classification model. In addition, the regime classification model was also utilized to successfully distinguish between the low confinement mode (L-mode), I-mode, and high confinement mode (H-mode) with 96.03% accuracy. Finally, ablation experiments were performed on the regime classifiers, showing that there is potential for further performance improvement with future use of RNN model.

The pressure crashes observed in shot No. 71326 in the Experimental Advanced Superconducting Tokamak are numerically investigated with the three-dimensional, toroidal, and full resistive-magnetohydrodynamics code (CLT). Based on the experimental observations, the pressure crash is caused by the nonlinear evolution of the m/n = 2/1 double tearing mode (DTM), where n and m are the toroidal and poloidal mode numbers, respectively. However, we find that the m/n = 2/1 DTM is stable based on the safety factor (q) profile from the equilibrium code EFIT, which indicates that the original q profile is somewhat inconsistent with the actual profile due to q measurement uncertainty. Since there is no motional Stark effect diagnostic for this shot, the local information of the magnetic field is missing, which leads to the largest contribution to the discrepancy. If other information is perfectly known and the q profile is the only uncertainty, then we could provide some information for the EFIT reconstruction by comparing our simulation results with electron cyclotron emission signals to constrain the uncertainty of the q profile to a much smaller region. The influence of plasma rotation and the two-fluids effect is also discussed.

The I-mode is a promising operation mode for fusion in the future, featuring high-temperature and low-density confinement, but the reason why the temperature and density are decoupled remains an important aspect to be explored. The experimental results from the experimental advanced superconducting tokamak (EAST) showed that the weakly coherent mode (WCM) is directly related to sustaining the I-mode and that the peak amplitude of the WCM is proportional to the temperature in the pedestal. Simulating the experimental data from EAST with the six-field model of BOUT++, we find a density perturbation close to the frequency of the WCM observed in the experiment. By testing all the physical terms in this model, we find that the density perturbation and particle transport are directly related to the drift Alfvén wave (DAW) mode. We also use the SymPIC program (Xiao et al Plasma Sci. Technol. 20 110501; Phys. Plasmas 22 112504; Plasma Sci. Technol. 23 055102) to simulate the same experimental data and find that the frequency range of the WCM is close to both experimental and BOUT++ results. Therefore, the WCM of the I-mode can be considered to be driven by the DAW, which helps improve the transport of the I-mode.