Conservative DG Method for the Micro-Macro Decomposition of the Vlasov-Poisson-Lenard-Bernstein Model

Abstract

The micro-macro (mM) decomposition approach is considered for the numerical solution of the Vlasov--Poisson--Lenard--Bernstein (VPLB) system, which is relevant for plasma physics applications. In the mM approach, the kinetic distribution function is decomposed as f=E[f]+g, where E is a local equilibrium distribution, depending on the macroscopic moments f=∫Re fdv=e fR, where e=(1,v,12v2)T, and g, the microscopic distribution, is defined such that e gR=0. We aim to design numerical methods for the mM decomposition of the VPLB system, which consists of coupled equations for f and g. To this end, we use the discontinuous Galerkin (DG) method for phase-space discretization, and implicit-explicit (IMEX) time integration, where the phase-space advection terms are integrated explicitly and the collision operator is integrated implicitly. We give special consideration to ensure that the resulting mM method maintains the e gR=0 constraint, which may be necessary for obtaining (i) satisfactory results in the collision dominated regime with coarse velocity resolution, and (ii) unambiguous conservation properties. The constraint-preserving property is achieved through a consistent discretization of the equations governing the micro and macro components. We present numerical results that demonstrate the performance of the mM method. The mM method is also compared against a corresponding DG-IMEX method solving directly for f.

0

Turn this paper into a lesson

ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…