The aim of the talk is to introduce a nonlinear Discrete Duality Finite Volume scheme to approximate the solutions of drift-diffusion equations. The scheme is built to preserve at the discrete level even on severely distorted meshes the energy / energy dissipation relation. This relation is of paramount importance to capture the long-time behavior of the problem in an accurate way. To enforce it, the linear convection diffusion equation is rewritten in a nonlinear form before being discretized. This is a joint work with Clément Cancès (Lille) and Stella Krell (Nice).[-]

The dynamical low-rank approximation is a low-rank factorization updating technique. It leads to differential equations for factors in a decomposition of the solution, which need to be solved numerically. The dynamical low-rank method seems particularly suitable for solving kinetic equations, because in many relevant cases the effective dynamics takes place on a lower-dimensional manifold and thus the solution has low rank. In this way, the 5-dimensional (3 space, 2 angle) radiation transport problem is reduced, both in computational cost as well as in memory footprint. We show several numerical examples.[-]

We present an efficient algorithm for the long time behavior of plasma simulations. We will focus on 4D drift-kinetic model, where the plasma's motion occurs in the plane perpendicular to the magnetic field and can be governed by the 2D guiding-center model. Hermite WENO reconstructions, already proposed in [1], are applied for solving the Vlasov equation. Here we consider an arbitrary computational domain with an appropriate numerical method for the treatment of boundary conditions. Then we apply this algorithm for plasma turbulence simulations. We first solve the 2D guiding-center model in a D-shape domain and investigate the numerical stability of the steady state. Then, the 4D drift-kinetic model is studied with a mixed method, i.e. the semi-Lagrangian method in linear phase and finite difference method during the nonlinear phase. Numerical results show that the mixed method is efficient and accurate in linear phase and it is much stable during the nonlinear phase. Moreover, in practice it has better conservation properties.

Keywords: Cartesian mesh - semi-Lagrangian method - Hermite WENO reconstruction - guiding-center - drift-kinetic model[-]