Section: New Software and Platforms
FBGKI (Full Braginskii)
Functional Description The Full Braginskii solver considers the equations proposed by Braginskii (1965), in order to describe the plasma turbulent transport in the edge part of tokamaks. These equations rely on a two fluid (ion - electron) description of the plasma and on the electroneutrality and electrostatic assumptions. One has then a set of 10 coupled non-linear and strongly anisotropic PDEs. FBGKI makes use in space of high order methods: Fourier in the toroidal periodic direction and spectral elements in the poloidal plane. The spectral vanishing viscosity (SVV) technique is implemented for stabilization. Static condensation is used to reduce the computational cost. In its sequential version, a matrix free solver is used to compute the potential. The parallel version of FBGKI presents two layers of parallelization: The first one corresponds to the poloidal plane and the second one to the toroidal direction. In the poloidal plane, the domain decomposition is achieved using the software METIS. For the parallel linear algebra, one uses the software PETSC (Portable Extensible Toolkit for Scientific Computation). The time discretization makes use of a Strang splitting, that decouples the explicit treatment of the advection and Braginskii terms, from the implicit treatment of the Lorentz forces and the computation of the electric potential. Whereas the explicit part is easily parallelized, the implicit one requires solving a strongly anisotropic elliptic problem for the potential. In the parallel version of FBGKI the system matrix is assembled in sparse manner, in order to allow using the multigrid HYPRE preconditionner implemented in PETSC. Till now results have only been obtained for computations done on a few tens of processors. Both the weak and strong scalings look satisfactory. Numerical experiments are still required to go up to hundreds or thousands of processors.