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=\mathcal{E}[{\mathbf{\rho }}_f]+g$, where $\mathcal{E}$ is a local equilibrium distribution, depending on the macroscopic moments ${\mathbf{\rho }}_f={\int }_{\mathbb{R}}\mathbf{e}fdv=⟨\mathbf{e}f{⟩}_{\mathbb{R}}$, where $\mathbf{e}=(1,v,\frac{1}{2}{v}^2{)}^{\mathrm{T}}$, and $g$, the microscopic distribution, is defined such that $⟨\mathbf{e}g{⟩}_{\mathbb{R}}=0$. We aim to design numerical methods for the mM decomposition of the VPLB system, which consists of coupled equations for ${\mathbf{\rho }}_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 $⟨\mathbf{e}g{⟩}_{\mathbb{R}}=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$.

arXiv commentSubmitted to Journal of Computational Physics

References in corpus (15)

• A class of asymptotic preserving schemes for kinetic equations and related problems with stiff sources
• A positivity-preserving high-order semi-Lagrangian discontinuous Galerkin scheme for the Vlasov-Poisson equations
• Discontinuous Galerkin algorithms for fully kinetic plasmas
• Energy-conserving discontinuous Galerkin methods for the Vlasov-Ampère system
• Gyrokinetic continuum simulation of turbulence in a straight open-field-line plasma
• An asymptotic-preserving dynamical low-rank method for the multi-scale multi-dimensional linear transport equation
• Conservative Discontinuous Galerkin Schemes for Nonlinear Fokker-Planck Collision Operators
• Block-Structured Adaptive Mesh Refinement Algorithms for Vlasov Simulation
• A high-order / low-order (HOLO) algorithm for preserving conservation in time-dependent low-rank transport calculations
• High Order Asymptotic Preserving Nodal Discontinuous Galerkin IMEX Schemes for the BGK Equation
• Conservative discontinuous Galerkin scheme of a gyro-averaged Dougherty collision operator
• Realizability-Preserving DG-IMEX Method for the Two-Moment Model of Fermion Transport
• Micro-macro decomposition based asymptotic-preserving numerical schemes and numerical moments conservation for collisional nonlinear kinetic equations
• Fluid & Gyrokinetic turbulence in open field-line, helical plasmas
• Discontinuous Galerkin schemes for a class of Hamiltonian evolution equations with applications to plasma fluid and kinetic problems