This paper describes a parallel implementation of a numerical solver for the Vlasov equation. The solver is based on a kinetic model describing the motion of charged particles in a plasma. The evolution of the distribution of particles in phase space is computed with an explicit method, and we take into account the self-consistent electric field through the coupling with a Poisson type equation. In this paper, we focus on a recently developed 5D parallel numerical application dedicated to gyrokinetic simulation of tokamak systems and ITG turbulence simulation. A semi-Lagrangian Vlasov solver is used. A specific cubic spline interpolation allows us to formulate a domain decomposition method. A hybrid MPI/OpenMP paradigm was used to benefit from a large number of processors while reducing communication costs.
[1]
C. Birdsall,et al.
Plasma Physics via Computer Simulation
,
2018
.
[2]
P. Bertrand,et al.
Conservative numerical schemes for the Vlasov equation
,
2001
.
[3]
Laurent Villard,et al.
A drift-kinetic Semi-Lagrangian 4D code for ion turbulence simulation
,
2006,
J. Comput. Phys..
[4]
Guillaume Latu,et al.
Hermite Spline Interpolation on Patches for Parallelly Solving the Vlasov-Poisson Equation
,
2007,
Int. J. Appl. Math. Comput. Sci..
[5]
G. Knorr,et al.
The integration of the vlasov equation in configuration space
,
1976
.