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

JOURNAL OF COMPUTATIONAL PHYSICS

卷:413

Article

Higginson, Drew Pitney¹  Holod, Ihor¹  Link, Anthony¹ 

[1] Lawrence Livermore Natl Lab, Livermore, CA 94551 USA

关键词: Coulomb collisions;    Weighted particles;    Numerical methods;    Monte-Carlo methods;    Collisional plasma;   

DOI  :  10.1016/j.jcp.2020.109450

来源: Elsevier

【 摘 要 】

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. (C) 2020 Elsevier Inc. All rights reserved.

【 授权许可】

Free