After twenty years of development, first-principle gyrokinetic simulation has become a major tool to study turbulent transport in tokamaks. In this work, we carried out gyrokinetic simulations using GTC code for several recent shots of HL-2A tokamak in H-mode. The major electrostatic instability with experimental parameters is identified by the GTC simulation to be trapped electron mode (TEM), with good agreement in both nonlinear frequency and poloidal mode number. At the linear stage, the electrostatic modes driven by steep grdient with the H-mode edge parameters show un-conventional mode structures, which are not localized at the outboard mid-plane of the poloidal cross section, while the conventional ballooning mode structures are linked with L-mode parameters. This un-conventional structures can peak at any poloidal angle between 0 and 2pi and can also have multi-peaks. The linear mode frequency is observed to jump to another branch when the gradient exceeds a critical value. The poloidal mode spectrum cascades inversely into longer wavelength region during the nonlinear saturation of turbulence, which are found more likely to come from the nonlinear evolution of each single mode instead of the mode-mode coupling process. The instability problem is mimicked by a model equation where a series of unstable eigen solutions are found. Under weak gradient, the most unstable modes are found to be in the ground state, which shows the conventional ballooning structure. However, under strong gradient, the most unstable solution jumps into the non-ground state, which shows aforementioned un-conventional mode structure. Thus, the L-H transition could be analogous to the transition between eigenstates. By nonlinear gyrokinetic simulations, we found that zonal flow would be less important to suppress turbulent transport in the strong gradient regime. More interestingly, we discovered that the turbulent transport coefficient would decrease with the gradient increasing when the gradient exceed a critical value. This provide a new route for the L to H transition without invoking shear flow or zonal flow. Physically, we would argue that the change of linear eigenmode structure leads to change of transport characteristics for the H-mode edge plasma, e.g., the reduction of radial correlation length, or change from diffusion and convection mixed transport to mainly diffusion transport.
[mehr]
