Gyroaveraging operations using adaptive matrix operators

A new adaptive scheme to be used in Particle-In-Cell codes for carrying out gyroaveraging operations with matrices is presented. This new scheme uses an intermediate velocity grid whose resolution is adapted to the local thermal Larmor radius. The charge density is computed by projecting marker weights in a field-line following manner while preserving the adiabatic magnetic moment $\mu$. These choices permit to improve the accuracy of the gyroaveraging operations performed with matrices even when strong spatial variation of temperature and magnetic field is present. Accuracy of the scheme in different geometries from simple 2d slab geometry to realistic 3d toroidal equilibrium has been studied. A successful implementation in the grokinetic code XGC is presented in the delta-f limit.

[1]  Patrick H. Worley,et al.  A fast low-to-high confinement mode bifurcation dynamics in the boundary-plasma gyrokinetic code XGC1 , 2018 .

[2]  Hans-Joachim Bungartz,et al.  Block-structured grids in full velocity space for Eulerian gyrokinetic simulations , 2017, Comput. Phys. Commun..

[3]  Laurent Villard,et al.  An arbitrary wavelength solver for global gyrokinetic simulations. Application to the study of fine radial structures on microturbulence due to non-adiabatic passing electron dynamics , 2017 .

[4]  Robert Hager,et al.  A new hybrid-Lagrangian numerical scheme for gyrokinetic simulation of tokamak edge plasma , 2016, J. Comput. Phys..

[5]  M. Nakata,et al.  Cross-Scale Interactions between Electron and Ion Scale Turbulence in a Tokamak Plasma. , 2015, Physical review letters.

[6]  F. Jenko,et al.  How non-adiabatic passing electron layers of linear microinstabilities affect turbulent transport , 2015 .

[7]  V. Grandgirard,et al.  Gyroaverage operator for a polar mesh , 2015 .

[8]  Martin Greenwald,et al.  Synergistic cross-scale coupling of turbulence in a tokamak plasma , 2014 .

[9]  Frank Jenko,et al.  The global version of the gyrokinetic turbulence code GENE , 2011, J. Comput. Phys..

[10]  L. Gavrilov,et al.  Journal of Physics: Conference Series 279 (2011) 012019 , 2011 .

[11]  O. Sauter,et al.  On the effects of the equilibrium model in gyrokinetic simulations: from s-α to diverted MHD equilibrium , 2010 .

[12]  Laurent Villard,et al.  Clarifications to the limitations of the s-α equilibrium model for gyrokinetic computations of turbulence , 2009 .

[13]  Shinji Tokuda,et al.  Conservative global gyrokinetic toroidal full-f five-dimensional Vlasov simulation , 2008, Comput. Phys. Commun..

[14]  X. Garbet,et al.  Global full-f gyrokinetic simulations of plasma turbulence , 2007 .

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

[16]  S. Lanthaler,et al.  Theory of Fusion Plasmas , 2007 .

[17]  Roman Hatzky,et al.  Particle simulations with a generalized gyrokinetic solver , 2005 .

[18]  Scott E. Parker,et al.  A δ f particle method for gyrokinetic simulations with kinetic electrons and electromagnetic perturbations , 2003 .

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

[20]  Laurent Villard,et al.  Finite element approach to global gyrokinetic Particle-In-Cell simulations using magnetic coordinates , 1998 .

[21]  Lin,et al.  Method for solving the gyrokinetic Poisson equation in general geometry. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

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