Development of semi-Lagrangian gyrokinetic code for full-f turbulence simulation in general tokamak geometry

In this work, we report the development of a new gyrokinetic code for full-f simulation of electrostatic turbulence in general tokamak geometry. Backward semi-Lagrangian scheme is employed for noise-free simulation of the gyrokinetic Vlasov equation with finite Larmor radius effects. Grid systems and numerical interpolations are implemented to deal with arbitrary equilibrium information given in a GEQDSK format file. In particular, we introduce an adaptive interpolation technique for fluctuating quantities, which have elongated structures along equilibrium magnetic fields. This field-aligned interpolation can reduce the required number of grid points to represent the fluctuating quantities. Also, it is shown that the new interpolation allows us to choose bigger time sizes with better simulation accuracy. Several benchmark simulation results are presented for comparison with previously known cases. It is demonstrated that the new code can reproduce the well known results of zonal flow and linear ITG instabilities in concentric circular equilibrium. It is also shown that the code can capture the effects of plasma shaping on the zonal flow and ITG instabilities, and the stabilization effects of the shaping result in significant up-shifts of the threshold of ion temperature gradient for ITG in nonlinear simulation.

[1]  Laurent Villard,et al.  The effect of plasma triangularity on turbulent transport: modeling TCV experiments by linear and non-linear gyrokinetic simulations , 2009 .

[2]  Özgür D. Gürcan,et al.  Analysis of symmetry breaking mechanisms and the role of turbulence self-regulation in intrinsic rotation , 2011 .

[3]  Peter J. Catto,et al.  Plasma shaping effects on the collisionless residual zonal flow level , 2006 .

[4]  Srinath Vadlamani,et al.  The particle-continuum method: an algorithmic unification of particle-in-cell and continuum methods , 2004, Comput. Phys. Commun..

[5]  Marshall N. Rosenbluth,et al.  POLOIDAL FLOW DRIVEN BY ION-TEMPERATURE-GRADIENT TURBULENCE IN TOKAMAKS , 1998 .

[6]  Charlson C. Kim,et al.  Large-scale gyrokinetic turbulence simulations: Effects of profile variation* , 1999 .

[7]  Shinji Tokuda,et al.  Study of ion turbulent transport and profile formations using global gyrokinetic full-f Vlasov simulation , 2009 .

[8]  Laurent Villard,et al.  Global gyrokinetic ion temperature gradient turbulence simulations of ITER , 2013 .

[9]  Choong-Seock Chang,et al.  Full-f gyrokinetic particle simulation of centrally heated global ITG turbulence from magnetic axis to edge pedestal top in a realistic tokamak geometry , 2009 .

[10]  E. Sonnendrücker,et al.  Comparison of Eulerian Vlasov solvers , 2003 .

[11]  Frank Jenko,et al.  On the role of numerical dissipation in gyrokinetic Vlasov simulations of plasma microturbulence , 2010, Comput. Phys. Commun..

[12]  R. E. Waltz,et al.  An Eulerian gyrokinetic-Maxwell solver , 2003 .

[13]  Timothy A. Davis,et al.  A column pre-ordering strategy for the unsymmetric-pattern multifrontal method , 2004, TOMS.

[14]  Shinji Tokuda,et al.  Global gyrokinetic simulation of ion temperature gradient driven turbulence in plasmas using a canonical Maxwellian distribution , 2003 .

[15]  M. Mehrenberger,et al.  Numerical Solution of the Gyroaverage Operator for the Finite Gyroradius Guiding-center Model , 2010 .

[16]  Jean-Michel Dischler,et al.  Beyond scale separation in gyrokinetic turbulence , 2007 .

[17]  T. Arber,et al.  A Critical Comparison of Eulerian-Grid-Based Vlasov Solvers , 2002 .

[18]  Volker Naulin,et al.  Transport of parallel momentum by collisionless drift wave turbulence , 2008 .

[19]  W. Lee,et al.  Gyrokinetic Particle Simulation Model , 1987 .

[20]  E. Sonnendrücker,et al.  The Semi-Lagrangian Method for the Numerical Resolution of the Vlasov Equation , 1999 .

[21]  Y. Camenen,et al.  Transport of parallel momentum induced by current-symmetry breaking in toroidal plasmas. , 2009, Physical review letters.

[22]  T. H. Stix,et al.  Decay of poloidal rotation in a tokamak plasma , 1973 .

[23]  Charlson C. Kim,et al.  Comparisons and physics basis of tokamak transport models and turbulence simulations , 2000 .

[24]  X. Garbet,et al.  Computing ITG turbulence with a full-f semi-Lagrangian code , 2008 .

[25]  Alain J. Brizard,et al.  Foundations of Nonlinear Gyrokinetic Theory , 2007 .

[26]  Guillaume Latu,et al.  Conservation equations and calculation of mean flows in gyrokinetics , 2011 .

[27]  Eric Sonnendrücker,et al.  A forward semi-Lagrangian method for the numerical solution of the Vlasov equation , 2008, Comput. Phys. Commun..

[28]  Robert D. Falgout,et al.  hypre: A Library of High Performance Preconditioners , 2002, International Conference on Computational Science.

[29]  Virginie Grandgirard,et al.  Defining an equilibrium state in global full-f gyrokinetic models , 2008 .

[30]  T. S. Hahm,et al.  Nonlinear gyrokinetic equations for tokamak microturbulence , 1988 .

[31]  Frank Jenko,et al.  Electron temperature gradient driven turbulence , 1999 .

[32]  A. Staniforth,et al.  Semi-Lagrangian integration schemes for atmospheric models - A review , 1991 .