Gyrokinetic Studies of Turbulence in Steep Gradient Region: Role of Turbulence Spreading and E x B Shear
Authors: T.S. Hahm, Z. Lin, P.H. Diamond, G. Rewoldt, W.X. Wang, S. Ethier, O. Gurcan, W.W. Lee, and W.M. Tang
Date of PPPL Report: December 2004
Presented at: the 20th IAEA Fusion Energy Conference, 1-6 November 2004, Vilamoura, Portugal. The papers will be published by the IAEA as unedited proceedings in electronic format on CD-ROM and on the IAEA Physics Section web site as soon as possible after the conference.
An integrated program of gyrokinetic particle simulation and theory has been developed to investigate several outstanding issues in both turbulence and neoclassical physics. Gyrokinetic particle simulations of toroidal ion temperature gradient (ITG) turbulence spreading using the GTC code and its related dynamical model have been extended to the case with radially increasing ion temperature gradient, to study the inward spreading of edge turbulence toward the core. Due to turbulence spreading from the edge, the turbulence intensity in the core region is significantly enhanced over the value obtained from simulations of the core region only. Even when the core gradient is within the Dimits shift regime (i.e., self-generated zonal flows reduce the transport to a negligible value), a significant level of turbulence and transport is observed in the core due to spreading from the edge. The scaling of the turbulent front propagation speed is closer to the prediction from our nonlinear diffusion model than one based on linear toroidal coupling. A calculation of ion poloidal rotation in the presence of sharp density and toroidal angular rotation frequency gradients from the GTC-Neo particle simulation code shows that the results are significantly different from the conventional neoclassical theory predictions. An energy conserving set of a fully electromagnetic nonlinear gyrokinetic Vlasov equation and Maxwell's equations, which is applicable to edge turbulence, is being derived via the phase-space action variational Lie perturbation method. Our generalized ordering takes the ion poloidal gyroradius to be on the order of the radial electric field gradient length.