A corrected method for Coulomb scattering in arbitrarily weighted particle-in-cell plasma simulations

• Published in: Journal of computational physics 2020-07, Vol.413 (C), p.109450, Article 109450
• Main Authors: Higginson, Drew Pitney, Holod, Ihor, Link, Anthony
• Format: Article
• Language: English
• Subjects: 70 PLASMA PHYSICS AND FUSION TECHNOLOGY; Collisional plasma; Computational physics; Computer simulation; Coulomb collisions; Energy conservation; Macroparticles; Monte-Carlo methods; Numerical methods; Particle in cell technique; Scattering; Simulation; Weighted particles
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2020.109450

•An accurate method for scattering with arbitrary particle weighting.•Validated through theory and numerical tests.•Tested within and between different particle groups over large weight ratios.•Errors in previous work Nanbu and Yonemura [J. Comput. Phys. 145 (1998) 639] identified. A physically and statistically accurate method to perform binary, pairwise Coulomb collisions for arbitrarily-weighted macroparticles in particle-in-cell simulations is described. The method is shown to conserve global energy and momentum on average and to produce physically accurate scattering frequencies. This is justified using the statistically expected values of the method and through a series of numerical tests showing that the method is valid independently of macroparticle weight. This novel method is required due to an error in previous work on binary scattering, initially by Nanbu and Yonemura (1998) [11] and applied to subsequent work, that causes inaccuracies in the scattering frequencies for unequally weighted simulations; a problem that occurs both for scattering within the group and between different groups.