Multispecies structure-preserving particle discretization of the Landau collision operator
arXiv:2212.00663 · doi:10.1063/5.0105182
Abstract
This paper proposes a novel numerical integrator for modeling multispecies Coulomb collisions in kinetic plasmas. The proposed scheme provides an energy-, momentum-, and positivity-preserving particle discretization of the nonlinear Landau collision operator, extending the works of J.A. Carrillo et al., Journal of Computational Physics, 7, 100066 (2020) and E. Hirvijoki, Plasma Physics and Controlled Fusion, 63, 044003 (2021). The discrete-time conservation properties are analyzed both algebraically and numerically, and an efficient, GPU-parallelized implementation is validated against inhomogeneous temperature relaxation, isotropization and thermalization examples. The results agree with analytical estimates, confirming the method capable of reproducing physics.
arXiv comment23 pages, 14 figures
References in corpus (17)
- GEMPIC: Geometric ElectroMagnetic Particle-In-Cell Methods
- Variational formulation of particle algorithms for kinetic plasma simulations
- Geometric integration of the Vlasov-Maxwell system with a variational particle-in-cell scheme
- Explicit high-order non-canonical symplectic particle-in-cell algorithms for Vlasov-Maxwell systems
- Canonical symplectic particle-in-cell method for long-term large-scale simulations of the Vlasov-Maxwell system
- Structure-Preserving Geometric Particle-in-Cell Methods for Vlasov-Maxwell Systems
- Hamiltonian integration methods for Vlasov-Maxwell equations
- Metriplectic Integrators for the Landau Collision Operator
- Explicit high-order noncanonical symplectic algorithms for ideal two-fluid systems
- Conservative discretization of the Landau collision integral
- Explicit Structure-Preserving Geometric Particle-in-Cell Algorithm in Curvilinear Orthogonal Coordinate Systems and Its Applications to Whole-Device 6D Kinetic Simulations of Tokamak Physics
- Structure-preserving marker-particle discretizations of Coulomb collisions for particle-in-cell codes
- Landau Collision Integral Solver with Adaptive Mesh Refinement on Emerging Architectures
- Structure-preserving strategy for conservative simulation of relativistic nonlinear Landau--Fokker--Planck equation
- A fully implicit, scalable, conservative nonlinear relativistic Fokker-Planck 0D-2P solver for runaway electron
- Subcycling of particle orbits in variational, geometric electromagnetic particle-in-cell methods
- Exascale Landau collision operator in the Cuda programming model applied to thermal quench plasmas