Section: New Results

DG scheme for Drift-Kinetic equation

Participants : Laurent Navoret [correspondent] , Philippe Helluy, Nhung Pham.

Using the discontinuous Galerkin solver of Schnaps, we have implemented a numerical scheme for the drift-kinetic model (in a cylinder geometry). The equation is written as an hyperbolic system after reduction in velocity (using spectral finite element). The code is parallelized on a multi-CPU or GPU architecture using OpenCL instructions. To solve the quasineutral equation (for the electric potential), the elliptic solver (already present in Schnaps) has been extended to be used slice by slice (of the cylinder). We have started by validating the code on the 2D guiding-center model and the diocotron instability test-case: we observe that the geometry approximation of the computational domain has a major impact on the precision of the numerical simulations.