Plasma and Fusion Research

Volume 15, 1403086 (2020)

Regular Articles

Study of Plasma Shaping Effects on ITG Instability Using Global Gyrokinetic Code GKNET with Analytical Magnetic Equilibrium
Kenji IMADERA, Jhih-Yi LIN, Daichi NAKAJIMA and Yasuaki KISHIMOTO
Department of Fundamental Energy Science, Graduate School of Energy Science, Kyoto University, Gokasho, Uji, Kyoto 611-0011, Japan
(Received 12 February 2020 / Accepted 13 October 2020 / Published 11 December 2020)


We extended the global gyrokinetic code GKNET (GyroKinetic Numerical Experiment of Tokamak) to non-circular shaped plasmas with analytical magnetic equilibria, which satisfy the Grad-Shafranov equation up to the second order with respect to aspect ratio. The extended version allows us to set the equilibria with non-unity elongation and non-zero triangularity, where the finite Shafranov shift is consistently determined. The allocated mesh follows the magnetic field line with periodic boundary conditions along the poloidal and toroidal directions so that the calculation cost for solving the gyrokinetic quasi-neutrality condition can be reduced by utilizing 1D FFT and MPI_ALLtoALL transpose technique. Based on the developed code, we studied the effect of elongation and triangularity on linear Ion Temperature Gradient (ITG) instability with adiabatic electrons in a non-circular shaped Tokamak by paying attention to the symmetry breaking of mode structure due to global profile effects. It is found that elongation reduces ITG instability owing to the effective reduction of flux-surface averaged ion temperature gradient and increases the asymmetry characterized by the Bloch angle θb. On the other hand, when elongation is approximately unity, triangularity weakly affects the growth rate, while negative triangularity stabilizes ITG modes and increases the asymmetry in a large elongation regime.


global gyrokinetic simulation, elongation, negative triangularity, symmetry breaking, Bloch angle

DOI: 10.1585/pfr.15.1403086


  • [1] J.G. Cordey et al., Nucl. Fusion 45, 1078 (2005).
  • [2] F. Troyon et al., Plasma Phys. Control. Fusion 26, 209 (1984).
  • [3] J.E. Kinsey et al., Phys. Plasmas 14, 102306 (2007).
  • [4] E.A. Belli et al., Phys. Plasmas 15, 092303 (2008).
  • [5] P. Angelino et al., Phys. Rev. Lett. 102, 195002 (2009).
  • [6] Y. Xiao and P.J. Catto, Phys. Plasmas 13, 082307 (2006).
  • [7] K. Imadera et al., Proc. of 26th IAEA Fusion Energy Conf., TH/P3-3 (2016).
  • [8] K. Obrejan et al., Plasma Fusion Res. 10, 3403042 (2015).
  • [9] K. Obrejan et al., Comput. Phys. Commum. 216, 8 (2017).
  • [10] R.L. Miller et al., Phys. Plasmas 5, 973 (1998).
  • [11] J.P. Graves, Plasma Phys. Control. Fusion 55, 074009 (2013).
  • [12] A. Fukuyama et al., Proc. of 20th IAEA Fusion Energy Conf., IAEA-CSP-25/TH/P2-3 (2004).
  • [13] V. Grandgirard et al., Comput. Phys. Commun. 207, 35 (2016).
  • [14] Y. Morinishi et al., J. Comput. Phys. 143, 90 (1998).
  • [15] Y. Idomura et al., J. Comput. Phys. 226, 224 (2007).
  • [16] Y. Kishimoto et al., Plasma Phys. Control. Fusion 41, A663 (1999).
  • [17] Y. Camenen et al., Nucl. Fusion 51, 073039 (2011).
  • [18] A. Ishizawa et al., Phys. Rev. Lett. 123, 025003 (2019).
  • [19] G.M. Staebler et al., Nucl. Fusion 58, 115001 (2018).