Extension of Gkeyll Discontinuous Galerkin Kinetic Code to 2D

Gkeyll is a discontinuous Galerkin (DG) code under development for modeling the edge plasma in fusion devices and basic plasma experiments. High-order accurate, energy-conserving numerical algorithms for general Hamiltonian systems are implemented in Gkeyll. Details of the recent extension of the code dimensionality to 2D2V will be presented. Since DG schemes allow for flexibility in the choice of basis functions, we will discuss how various types of basis functions affect code accuracy and efficiency. Test problems in 2D, such as toroidal ITG instabilities and turbulence in a local limit, will be presented. We will also show initial results from 2D kinetic simulations of transport in a scrape-off layer plasma, using a specified diffusion coefficient to model radial transport.