Numerical scheme for treatment of Uehling-Uhlenbeck equation for two-particle interactions in relativistic plasma

JOURNAL OF COMPUTATIONAL PHYSICS

卷:373

Article

Prakapenia, M. A.^1,2  Siutsou, I. A.¹  Vereshchagin, G. V.^1,3 

[1] Natl Acad Sci Belarus, Inst Phys, ICRANet Minsk, Nezaleznasci Av 68-2, Minsk 220072, BELARUS

[2] Belarusian State Univ, Dept Theoret Phys & Astrophys, Nezaleznasci Av 4, Minsk 220030, BELARUS

[3] ICRANet, Piazza Repubbl 10, I-65122 Pescara, Italy

关键词: Uehling-Uhlenbeck equations;    Collision integral;    Binary interactions;    Relativistic plasma;   

DOI  :  10.1016/j.jcp.2018.07.010

来源: Elsevier

【 摘 要 】

We present a new efficient method to compute Uehling-Uhlenbeck collision integral for all two-particle interactions in relativistic plasma with drastic improvement in computation time with respect to existing methods. Plasma is assumed isotropic in momentum space. The set of reactions consists of: Moeller and Bhabha scattering, Compton scattering, two-photon pair annihilation, and two-photon pair production, which are described by QED matrix elements. In our method exact energy and particle number conservation laws are fulfilled. Reaction rates are compared, where possible, with the corresponding analytical expressions and convergence of numerical rates is demonstrated. (C) 2018 Published by Elsevier Inc.

【 授权许可】

Free