An asymptotic preserving scheme based on a micro-macro decomposition for collisional Vlasov equations: diffusion and high-field scaling limits (Q549878)

In the present work, the micro-macro decomposition based numerical schemes developed in [\textit{R. Belaouar, N. Crouseilles, P. Degond} and \textit{E. Sonnendrücker}, J. Sci. Comput. 41, No. 3, 341--365 (2009; Zbl 1203.65217)] are extended to the solution of the kinetic equation of a plasma with a linear collision operator, including the Poisson equation. This system of equations is called ``collisional Vlasov-Poisson model''. The authors first derive the Vlasov-Poisson model as a system that couples the macroscopic (equilibrium) part with the fluctuations of the system. A suitable discretization of this micro-macro model enables to derive an asymptotic preserving scheme in the diffusion and high-field asymptotics. In case of the diffusion limit, the collision operator is the Bhatnagar-Gross-Krook one, for the high-field limit the Fokker-Planck collision operator is taken into account. In addition, two main improvements are presented: On the one hand, a self-consistent electric field is introduced, which induces a specific discretization in the velocity direction and represents a wide range of applications in plasma physics. On the other hand, as suggested in [\textit{M. Lemou}, C. R., Math., Acad. Sci. Paris 348, No. 7--8, 455--460 (2010; Zbl 1188.65120)], a suitable reformulation of the micro-macro scheme is introduced, which leads to the following asymptotic preserving property: it degenerates into an implicit scheme for the diffusion limit model when \(\varepsilon\rightarrow 0\) (\(\varepsilon\) is a small parameter), which makes it free from the usual diffusion constraint \(\Delta t= O(\Delta x^2)\) (\(\Delta t\) -- time step, \(\Delta x\) -- spatial step) in all regimes. Numerical examples are used to demonstrate the efficiency and applicability of the schemes for both regimes. There, the periodic boundary conditions and Landau damping, the asymptotic equivalence of the proposed scheme to the limit scheme up to the order \(O(\varepsilon^2)\), and the Dirichlet boundary condition and boundary layer (imposed by a non-equilibrium boundary condition) are considered.