The linear and nonlinear effects of hydrogen isotope species on ion temperature gradient (ITG) instability in the Large Helical Device (LHD) stellarator are studied using radially global gyrokinetic simulation. We found that the coupling range of linear toroidal harmonics depends on the ion mass of the hydrogen isotope. The growth rates of ITG mode are almost the same for H, D, and T plasmas, indicating a gyro-Bohm scaling of ion-mass dependence. The nonlinear electrostatic simulations show that the zonal flow breaks the radially elongated eigenmode structures and reduces the size of the turbulence eddies, which suppresses the turbulence and the ion heat transport in the LHD. The turbulence amplitude without the zonal flow is almost the same for H, D, and T plasmas, while it decreases with increasing ion mass of the hydrogen isotope when the zonal flow is present. The reduction of the turbulent transport with larger ion mass is mostly due to the enhancement of zonal flows by larger ion mass. The ion heat conductivity deviates from the gyro-Bohm scaling for both cases with and without the zonal flow.
On the role of ion temperature gradient turbulence in driving ion thermal transport in neutral beam injection-heated L-mode plasmas in a superconducting tokamak