Multispecies structure-preserving particle discretization of the Landau collision operator

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.

23 pages, 14 figures

References in corpus (4)

• Explicit high-order noncanonical symplectic algorithms for ideal two-fluid systems
• Structure-preserving marker-particle discretizations of Coulomb collisions for particle-in-cell codes
• Structure-preserving strategy for conservative simulation of relativistic nonlinear Landau--Fokker--Planck equation
• Exascale Landau collision operator in the Cuda programming model applied to thermal quench plasmas

Cited by in corpus (1)

• Energy-momentum-conserving stochastic differential equations and algorithms for nonlinear Landau-Fokker-Planck equation