Gyrokinetic simulations of the kinetic electron effects on the electrostatic instabilities on the ITER baseline scenario

【摘要】 The linear and nonlinear simulations are carried out using the gyrokinetic code NLT for the electrostatic instabilities in the core region of a deuterium plasma based on the International Thermonuclear Experimental Reactor（ITER） baseline scenario.The kinetic electron effects on the linear frequency and nonlinear transport are studied by adopting the adiabatic electron model and the fully drift-kinetic electron model in the NLT code,respectively.The linear simulations focus on the dependence of linear frequency on the plasma parameters,such as the ion and electron temperature gradients κ_Li,e≡R/L_Ti,e,the density gradient K_n≡R/L_n and the ion-electron temperature ratio τ=T_e/T_i.Here,R is the major radius,and T_e and T_i denote the electron and ion temperatures,respectively.L_A=-（?_rlnA）^-1 is the gradient scale length,with A denoting the density,the ion and electron temperatures,respectively.In the kinetic electron model,the ion temperature gradient（ITG） instability and the trapped electron mode（TEM） dominate in the small and large k_θ region,respectively,where k_θ is the poloidal wavenumber.The TEMdominant region becomes wider by increasing（decreasing） κ_Te(κ_Ti) or by decreasing K_n.For the nominal parameters of the ITER baseline scenario,the maximum growth rate of dominant ITG instability in the kinetic electron model is about three times larger than that in the adiabatic electron model.The normalized linear frequency depends on the value of τ,rather than the value of T_e or T_i,in both the adiabatic and kinetic electron models.The nonlinear simulation results show that the ion heat diffusivity in the kinetic electron model is quite a lot larger than that in the adiabatic electron model,the radial structure is finer and the time oscillation is more rapid.In addition,the magnitude of the fluctuated potential at the saturated stage peaks in the ITGdominated region,and contributions from the TEM（dominating in the higher k_θ region） to the nonlinear transport can be neglected.In the adiabatic electron model,the zonal radial electric field is found to be mainly driven by the turbulent energy flux,and the contribution of turbulent poloidal Reynolds stress is quite small due to the toroidal shielding effect.However,in the kinetic electron model,the turbulent energy flux is not strong enough to drive the zonal radial electric field in the nonlinear saturated stage.The kinetic electron effects on the mechanism of the turbulence-driven zonal radial electric field should be further investigated.更多还原

【Abstract】 The linear and nonlinear simulations are carried out using the gyrokinetic code NLT for the electrostatic instabilities in the core region of a deuterium plasma based on the International Thermonuclear Experimental Reactor（ITER） baseline scenario.The kinetic electron effects on the linear frequency and nonlinear transport are studied by adopting the adiabatic electron model and the fully drift-kinetic electron model in the NLT code,respectively.The linear simulations focus on the dependence of linear frequency on the plasma parameters,such as the ion and electron temperature gradients κ_Li,e≡R/L_Ti,e,the density gradient K_n≡R/L_n and the ion-electron temperature ratio τ=T_e/T_i.Here,R is the major radius,and T_e and T_i denote the electron and ion temperatures,respectively.L_A=-（?_rlnA）^-1 is the gradient scale length,with A denoting the density,the ion and electron temperatures,respectively.In the kinetic electron model,the ion temperature gradient（ITG） instability and the trapped electron mode（TEM） dominate in the small and large k_θ region,respectively,where k_θ is the poloidal wavenumber.The TEMdominant region becomes wider by increasing（decreasing） κ_Te(κ_Ti) or by decreasing K_n.For the nominal parameters of the ITER baseline scenario,the maximum growth rate of dominant ITG instability in the kinetic electron model is about three times larger than that in the adiabatic electron model.The normalized linear frequency depends on the value of τ,rather than the value of T_e or T_i,in both the adiabatic and kinetic electron models.The nonlinear simulation results show that the ion heat diffusivity in the kinetic electron model is quite a lot larger than that in the adiabatic electron model,the radial structure is finer and the time oscillation is more rapid.In addition,the magnitude of the fluctuated potential at the saturated stage peaks in the ITGdominated region,and contributions from the TEM（dominating in the higher k_θ region） to the nonlinear transport can be neglected.In the adiabatic electron model,the zonal radial electric field is found to be mainly driven by the turbulent energy flux,and the contribution of turbulent poloidal Reynolds stress is quite small due to the toroidal shielding effect.However,in the kinetic electron model,the turbulent energy flux is not strong enough to drive the zonal radial electric field in the nonlinear saturated stage.The kinetic electron effects on the mechanism of the turbulence-driven zonal radial electric field should be further investigated.更多还原