With the aim of solving in a four dimensional phase space a multi-scale Vlasov-Poisson system, we propose in a Particle-In-Cell framework a robust time-stepping method that works uniformly when the small parameter vanishes. As an exponential integrator, the scheme is able to use large time steps with respect to the typical size of the solution’s fast oscillations. In addition, we show numerically that the method has accurate long time behaviour and that it is asymptotic preserving with respect to the limiting Guiding Center system.

In plasma physics domain, the electron transport is described with the Fokker-Planck-Landau equation. The direct numerical solution of the kinetic equation is usually intractable due to the large number of independent variables. That is why we propose in this paper a new model whose derivation is based on an angular closure in the phase space and retains only the energy of particles as kinetic dimension. To find a solution compatible with physics conditions, the closure of the moment system is obtained under a minimum entropy principle. This model is proved to satisfy the fundamental properties like a H theorem. Moreover an entropic discretization in the velocity variable is proposed on the semi-discrete model. Finally, we validate on numerical test cases the fundamental properties of the full discrete model.

We construct an efficient numerical scheme for the quantum Fokker-Planck- Landau (FPL) equation that works uniformly from kinetic to fluid regimes. Such a scheme inevitably needs an implicit discretization of the nonlinear collision operator, which is difficult to invert. Inspired by work [9] we seek a linear operator to penalize the quantum FPL collision term QqFPL in order to remove the stiffness induced by the small Knudsen number. However, there is no suitable simple quantum operator serving the purpose and for this kind of operators one has to solve the complicated quantum Maxwellians (Bose-Einstein or Fermi-Dirac distribution). In this paper, we propose to penalize QqFPL by the "classical" linear Fokker-Planck operator. It is based on the observation that the classical Maxwellian, with the temperature replaced by the internal energy, has the same first five moments as the quantum Maxwellian. Numerical results for Bose and Fermi gases are presented to illustrate the efficiency of the scheme in both fluid and kinetic regimes.