A fully non-linear multi-species Fokker-Planck-Landau collision operator for simulation of fusion plasma

JOURNAL OF COMPUTATIONAL PHYSICS

卷:315

Article

Hager, Robert¹  Yoon, E. S.²  Ku, S.¹  D'Azevedo, E. F.³  Worley, P. H.³  Chang, C. S.¹ 

[1] Princeton Plasma Phys Lab, POB 451, Princeton, NJ 08543 USA

[2] Rensselaer Polytech Inst, 110 8th St, Troy, NY 12180 USA

[3] Oak Ridge Natl Lab, POB 2008, Oak Ridge, TN 37831 USA

关键词: Plasma;    Fusion;    Collision operator;    XGC;    Particle-in-cell;   

DOI  :  10.1016/j.jcp.2016.03.064

来源: Elsevier

【 摘 要 】

Fusion edge plasmas can be far from thermal equilibrium and require the use of a non-linear collision operator for accurate numerical simulations. In this article, the non-linear single-species Fokker-Planck-Landau collision operator developed by Yoon and Chang (2014) [9] is generalized to include multiple particle species. The finite volume discretization used in this work naturally yields exact conservation of mass, momentum, and energy. The implementation of this new non-linear Fokker-Planck-Landau operator in the gyrokinetic particle-in-cell codes XGC1 and XGCa is described and results of a verification study are discussed. Finally, the numerical techniques that make our non-linear collision operator viable on high-performance computing systems are described, including specialized load balancing algorithms and nested OpenMP parallelization. The collision operator's good weak and strong scaling behavior are shown. (C) 2016 Elsevier Inc. All rights reserved.

【 授权许可】

Free