Numerical Study of the Two-Species Vlasov-Ampère System: Energy-Conserving Schemes and the Current-Driven Ion-Acoustic Instability

In this paper, we propose energy-conserving Eulerian solvers for the two-species Vlasov-Ampère (VA) system and apply the methods to simulate current-driven ion-acoustic instability. The algorithm is generalized from our previous work for the single-species VA system [9] and Vlasov-Maxwell (VM) system [8]. The main feature of the schemes is their ability to preserve the total particle number and total energy on the fully discrete level regardless of mesh size. Those are desired properties of numerical schemes especially for long time simulations with under-resolved mesh. The conser-vation is realized by explicit and implicit energy-conserving temporal discretizations, and the discontinuous Galerkin (DG) spatial discretizations. We benchmarked our algorithms on a test example to check the one-species limit, and the current-driven ion-acoustic instability. To simulate the current-driven ion-acoustic instability, a slight modification for the implicit method is necessary to fully decouple the split equations. This is achieved by a Gauss-Seidel type iteration technique. Numerical results verified the conservation and performance of our methods