Development of a Gyrokinetic Particle-in-Cell Code for Whole-Volume Modeling of Stellarators

Abstract

We present initial results in the development of a gyrokinetic particle-in-cell code for the whole-volume modeling of stellarators. This is achieved through two modifications to the X-point Gyrokinetic Code (XGC), originally developed for tokamaks. One is an extension to three-dimensional geometries with an interface to Variational Moments Equilibrium Code (VMEC) data. The other is a connection between core and edge regions that have quite different field-line structures. The VMEC equilibrium is smoothly extended to the edge region by using a virtual casing method. Non-axisymmetric triangular meshes in which triangle nodes follow magnetic field lines in the toroidal direction are generated for field calculation using a finite-element method in the entire region of the extended VMEC equilibrium. These schemes are validated by basic benchmark tests relevant to each part of the calculation cycle, that is, particle push, particle-mesh interpolation, and field solver in a magnetic field equilibrium of Large Helical Device including the edge region. The developed code also demonstrates collisionless damping of geodesic acoustic modes and steady states with residual zonal flow in the core region.