Gyrokinetic Semi-Lagrangian Parallel Simulation using a Hybrid OpenMP/MPI programming

Guillaume Latu (1, 2) , Nicolas Crouseilles (1, 3) , Virginie Grandgirard (4) , Eric Sonnendrücker (1, 3)

1 CALVI - Scientific computation and visualization
2 LaBRI - Laboratoire Bordelais de Recherche en Informatique
3 IRMA - Institut de Recherche Mathématique Avancée
4 IRFM - Institut de Recherche sur la Fusion par confinement Magnétique

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.

Domaines: Systèmes dynamiques [math.DS]

Guillaume Latu, Nicolas Crouseilles, Virginie Grandgirard, Eric Sonnendrücker. Gyrokinetic Semi-Lagrangian Parallel Simulation using a Hybrid OpenMP/MPI programming. 14th European PVM/MPI User's Group Meeting, Sep 2007, Paris, France. pp.356-364, ⟨10.1007/978-3-540-75416-9⟩. ⟨hal-00605748⟩