Using a local gyrokinetic code to study global ion temperature gradient modes in tokamaks

作者：Abdoul P A^*; Dickinson D; Roach C M; Wilson H R

来源：Plasma Physics and Controlled Fusion, 2015, 57(6): 065004.

DOI：10.1088/0741-3335/57/6/065004

• local gyrokinetics
• profile shearing
• flow shear
• global gyrokinetics

In this paper the global eigenmode structures of linear ion temperature gradient (ITG) modes in tokamak plasmas are obtained using a novel technique which combines results from the local gyrokinetic code GS2 with analytical theory to reconstruct global properties. Local gyrokinetic calculations are performed for a range of radial flux surfaces, x, and ballooning phase angles, p, to map out the local complex mode frequency, Omega(0)(x, p) = omega(0)(x, p) + i gamma(0)(x, p) for a single toroidal mode number, n. Taylor expanding Omega(0) about a reference surface at x = 0, and employing the Fourier-ballooning representation leads to a second order ODE for the amplitude envelope, A(p), which describes how the local results are combined to form the global mode. The equilibrium profiles impact on the variation of Omega(0)(x, p) and hence influence the global mode structure. The simulations presented here are based upon a global extension to the CYCLONE base case and employ the circular Miller equilibrium model. In an equilibrium with radially varying profiles of a/L-T and a/L-n, peaked at x = 0, and with all other equilibrium profiles held constant, including eta(i) = L-n/L-T, Omega(0)(x, p) is found to have a stationary point. The reconstructed global mode sits at the outboard mid-plane of the tokamak, with global growth rate, gamma similar to Max[gamma(0)]. Including the radial variation of other equilibrium profiles like safety factor and magnetic shear, leads to a mode that peaks away from the outboard mid-plane, with a reduced global growth rate. Finally, the influence of toroidal flow shear has also been investigated through the introduction of a Doppler shift, omega(0) --> omega(0) - n Omega(phi)' x, where Omega(phi) is the equilibrium toroidal flow, and a prime denotes the radial derivative. The equilibrium profile variations introduce an asymmetry into the global growth rate spectrum with respect to the sign of Omega(phi)', such that the maximum growth rate is achieved with non-zero shearing, consistent with recent global gyrokinetic calculations.

• 出版日期2015-6