A Discontinuous Galerkin Method for the Micro-Macro Formulation of the Vlasov-Poisson-Lenard-Bernstein System
The Vlasov-Poisson-Lenard-Bernstein (VPLB) system can serve as an approximate model to study transport processes in fusion and space plasmas. These systems are often characterized by multiple spatial and temporal scales, depending on the degree of collisionality. On the one hand, when inter-particle collisions are frequent, the plasma can be described quite well by fluid models. On the other hand, when collisions are infrequent, a fully kinetic description is needed to capture non-equilibrium effects. The micro-Macro (mM) formulation decomposes the particle distribution function into fluid and kinetic components, which are evolved separately with a coupled system of equations. One motivation for adopting the mM model is to gain computational efficiency in collision-dominated regimes. However, maintaining structural properties (e.g., particle, momentum, and energy conservation) is challenging. In this talk, we present a numerical method for the mM formulation of the VPLB system. This method is based on the discontinuous Galerkin method for phase-space discretization and implicit-explicit time stepping. We focus on the design of consistent discretization of the micro and macro components, which helps recover conservation properties of the method. Numerical results, demonstrating conservation properties and accuracy with coarse phase-space resolution, are also presented.