Convergence of an Eulerian scheme for the Vlasov-Poisson-BGK model

TL;DR

This paper provides the first convergence analysis for a finite-difference Eulerian scheme applied to the Vlasov-Poisson-BGK model, addressing stability challenges and establishing error estimates.

Contribution

It introduces a novel convergence theory for a non-splitting Eulerian scheme on the full phase-space grid for the VPBGK model, overcoming stability obstacles.

Findings

Established error estimates for the distribution function in weighted L∞ norm.

Proved convergence of the scheme under a truncated velocity domain with Neumann boundary conditions.

Developed a modified lower bound estimate for ionized systems to ensure stability.

Abstract

The Vlasov-Poisson-BGK (VPBGK) model is a kinetic model for describing the dynamics of collisional plasmas. Although various numerical schemes have been developed for it, a corresponding convergence theory has been absent. This paper fills this gap by presenting the first convergence analysis for a non-splitting, finite-difference Eulerian scheme discretized on the full phase-space grid. A major theoretical obstacle is the mixing of velocity indices induced by the electric field, which hinders the derivation of a uniform lower bound for the discrete solution. To overcome this stability challenge, we propose a modified lower bound estimate suitable for ionized systems that incorporates the step-wise degradation. Under a truncated velocity domain with a Neumann boundary condition, we establish error estimates for the distribution function in a weighted $L^{\infty}$ norm and for the electric field in a $L^{\infty}$ norm, respectively.