Search

Semi-Lagrangian solvers for the Vlasov system offer noiseless solutions compared to Lagrangian particle methods and can handle larger time steps compared to Eulerian methods. In order to reduce the computational complexity of the interpolation steps, it is common to use a directional splitting. However, this typically yields the wrong angular velocity. In this paper, we analyze a semi-Lagrangian method that treats the $v \times B$ term with a rotational grid and combines this with a directional splitting for the remaining terms. We analyze the convergence properties of the scheme both analytically and numerically. The favorable numerical properties of the rotating grid solution are demonstrated for the case of ion Bernstein waves.

With the increase in computational capabilities over the last years it becomes possible to simulate more and more complex and accurate physical models. Gyrokinetic theory has been introduced in the 1960s and 1970s in the need of describing a plasma with more accurate models than fluid equations, but eliminating the complexity of the fast gyration about the magnetic field lines. Although results from current gyrokinetic computer simulations are in fair agreement with experimental results in core physics, crucial assumptions made in the derivation make it unreliable in regimes of higher fluctuations and stronger gradient, such as the tokamak edge. With our novel optimized and scalable semi-Lagrangian solver we are able to simulate ion-temperature gradient modes with the 6D kinetic model including the turbulent saturation. After thoroughly testing our simulation code against analytical computations and gyrokinetic simulations (with the gyrokinetic code GYRO), it has been possible to show first plasma properties that go beyond standard gyrokinetic simulations. This includes the explicit description of the complete perpendicular energy fluxes and the excitation of high frequency waves (around the Larmor frequency) in the nonlinear saturation phase.