Implementation of multiple species in the GT5D code

Jolliet, S.; Idomura, Yasuhiro   

In this work, the gyrokinetic model of the GT5D code is extended to include an arbitrary number of ion species. First, the implementation of a collisionless, multiple-species gyrokinetic model is presented and compared to the single species model. Linear simulations are then presented and compared to analytical and numerical results. Then, the implementation of a recently proposed multiple-species collision operator is presented. This approximated linear operator has all the physical (conservation of density, momentum and energy, H-theorem) and mathematical (self-adjointness) properties desired to accurately describe neoclassical processes, while the usual linearized Landau operator does not satisfy H-theorem and is not self-adjoint when species have different temperatures. The implementation is presented and compared to the single-species case. Finally, several numerical tests are presented.