An efficient, conservative, time-implicit solver for the fully kinetic arbitrary-species 1D-2V Vlasov-Ampère system

Abstract

We consider the solution of the fully kinetic (including electrons) Vlasov-Ampère system in a one-dimensional physical space and two-dimensional velocity space (1D-2V) for an arbitrary number of species with a time-implicit Eulerian algorithm. The problem of velocity-space meshing for disparate thermal and bulk velocities is dealt with by an adaptive coordinate transformation of the Vlasov equation for each species, which is then discretized, including the resulting inertial terms. Mass, momentum, and energy are conserved, and Gauss's law is enforced to within the nonlinear convergence tolerance of the iterative solver through a set of nonlinear constraint functions while permitting significant flexibility in choosing discretizations in time, configuration, and velocity space. We mitigate the temporal stiffness introduced by, e.g., the plasma frequency through the use of high-order/low-order (HOLO) acceleration of the iterative implicit solver. We present several numerical results for canonical problems of varying degrees of complexity, including the multiscale ion-acoustic shock wave problem, which demonstrate the efficacy, accuracy, and efficiency of the scheme.

References in corpus (14)

• An Energy- and Charge-conserving, Implicit, Electrostatic Particle-in-Cell Algorithm
• A positivity-preserving high-order semi-Lagrangian discontinuous Galerkin scheme for the Vlasov-Poisson equations
• Structure and structure-preserving algorithms for plasma physics
• A multi-dimensional, energy- and charge-conserving, nonlinearly implicit, electromagnetic Vlasov-Darwin particle-in-cell algorithm
• Energy-conserving discontinuous Galerkin methods for the Vlasov-Ampère system
• Asymptotic-Preserving methods and multiscale models for plasma physics
• Energy-conserving Discontinuous Galerkin Methods for the Vlasov-Maxwell System
• Fluid preconditioning for Newton-Krylov-based, fully implicit, electrostatic particle-in-cell simulations
• Ion Species Stratification Within Strong Shocks in Two-Ion Plasmas
• Gauss's Law Satisfying Energy-Conserving Semi-Implicit Particle-in-Cell Method
• A Rescaling Velocity Method for Dissipative Kinetic Equations - Applications to Granular Media
• Quadratic conservative scheme for relativistic Vlasov--Maxwell system
• An Adaptive, Implicit, Conservative 1D-2V Multi-Species Vlasov-Fokker-Planck Multiscale Solve in Planar Geometry
• A conservative phase-space moving-grid strategy for a 1D-2V Vlasov-Fokker-Planck Equation

Cited by in corpus (3)

• An Eulerian Vlasov-Fokker-Planck Algorithm for Spherical Implosion Simulations of Inertial Confinement Fusion Capsules
• An Energy Conserving Vlasov Solver That Tolerates Coarse Velocity Space Resolutions: Simulation of MMS Reconnection Events
• Fully kinetic simulations of strong steady-state collisional planar plasma shocks