A multi-dimensional, energy- and charge-conserving, nonlinearly implicit, electromagnetic Vlasov-Darwin particle-in-cell algorithm

Abstract

For decades, the Vlasov-Darwin model has been recognized to be attractive for particle-in-cell (PIC) kinetic plasma simulations in non-radiative electromagnetic regimes, to avoid radiative noise issues and gain computational efficiency. However, the Darwin model results in an elliptic set of field equations that renders conventional explicit time integration unconditionally unstable. Here, we explore a fully implicit PIC algorithm for the Vlasov-Darwin model in multiple dimensions, which overcomes many difficulties of traditional semi-implicit Darwin PIC algorithms. The finite-difference scheme for Darwin field equations and particle equations of motion is space-time-centered, employing particle sub-cycling and orbit-averaging. The algorithm conserves total energy, local charge, canonical-momentum in the ignorable direction, and preserves the Coulomb gauge exactly. An asymptotically well-posed fluid preconditioner allows efficient use of large time steps and cell sizes, which are determined by accuracy considerations, not stability, and can be orders of magnitude larger than required in a standard explicit electromagnetic PIC simulation. We demonstrate the accuracy and efficiency properties of the algorithm with various numerical experiments in 2D-3V.

35 pages, 6 figures

References in corpus (1)

• Darwin-Vlasov Simulations of magnetized Plasmas

Cited by in corpus (27)

• Exactly Energy Conserving Semi-Implicit Particle in Cell Formulation
• Asymptotic-Preserving methods and multiscale models for plasma physics
• Electromagnetic full-$f$ gyrokinetics in the tokamak edge with discontinuous Galerkin methods
• Energy-conserving time propagation for a geometric particle-in-cell Vlasov--Maxwell solver
• Sparse grid techniques for particle-in-cell schemes
• Gauss's Law Satisfying Energy-Conserving Semi-Implicit Particle-in-Cell Method
• An energy-conserving and asymptotic-preserving charged-particle orbit implicit time integrator for arbitrary electromagnetic fields
• A semi-implicit, energy- and charge-conserving particle-in-cell algorithm for the relativistic Vlasov-Maxwell equations
• A fully implicit, conservative, non-linear, electromagnetic hybrid particle-ion/fluid-electron algorithm
• Quadratic conservative scheme for relativistic Vlasov--Maxwell system
• Enabling particle applications for exascale computing platforms
• ViDA: a Vlasov-DArwin solver for plasma physics at electron scales
• Verification of long wavelength electromagnetic modes with a gyrokinetic-fluid hybrid model in the XGC code
• Modeling of very high frequency large-electrode capacitively coupled plasmas with a fully electromagnetic particle-in-cell code
• An implicit, conservative and asymptotic-preserving electrostatic particle-in-cell algorithm for arbitrarily magnetized plasmas in uniform magnetic fields
• Verification of a Fully Implicit Particle-in-Cell Method for the $v_{\parallel }$ Formalism of Electromagnetic Gyrokinetics in the XGC Code
• An efficient, conservative, time-implicit solver for the fully kinetic arbitrary-species 1D-2V Vlasov-Ampère system
• An unsupervised machine-learning checkpoint-restart algorithm using Gaussian mixtures for particle-in-cell simulations
• Arbitrary-order time-accurate semi-Lagrangian spectral approximations of the Vlasov-Poisson system
• A Semi-Lagrangian Spectral Method for the Vlasov-Poisson System based on Fourier, Legendre and Hermite Polynomials
• A conservative implicit-PIC scheme for the hybrid kinetic-ion fluid-electron plasma model on curvilinear meshes
• Magnetic Fluctuations in Gyrokinetic Simulations of Tokamak Scrape-Off Layer Turbulence
• An implicit, conservative electrostatic particle-in-cell algorithm for paraxial magnetic nozzles
• Hamiltonian reduction of Vlasov-Maxwell to a dark slow manifold
• Kinetic Simulation of Collisional Magnetized Plasmas with Semi-Implicit Time Integration
• Finite spatial-grid effects in energy-conserving particle-in-cell algorithms
• Bridging kinetic plasma descriptions and single fluid models