A fully nonlinear multispecies Fokker–Planck–Landau collision operator for simulation of fusion plasma
Abstract
Fusion edge plasmas can be far from thermal equilibrium and require the use of a nonlinear collision operator for accurate numerical simulations. In this article, the nonlinear singlespecies 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 nonlinear Fokker–Planck–Landau operator in the gyrokinetic particleincell codes XGC1 and XGCa is described and results of a verification study are discussed. Finally, the numerical techniques that make our nonlinear collision operator viable on highperformance 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.
 Authors:

 Princeton Plasma Physics Laboratory, P.O. Box 451, Princeton, NJ 08543 (United States)
 Publication Date:
 OSTI Identifier:
 22572321
 Resource Type:
 Journal Article
 Journal Name:
 Journal of Computational Physics
 Additional Journal Information:
 Journal Volume: 315; Other Information: Copyright (c) 2016 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA); Journal ID: ISSN 00219991
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; ALGORITHMS; COLLISIONS; COMPUTERIZED SIMULATION; FOKKERPLANCK EQUATION; IMPLEMENTATION; NONLINEAR PROBLEMS; PERFORMANCE; PLASMA; THERMAL EQUILIBRIUM
Citation Formats
Hager, Robert, Yoon, E.S., Email: yoone@rpi.edu, Ku, S., Email: sku@pppl.gov, D'Azevedo, E.F., Email: dazevedoef@ornl.gov, Worley, P.H., Email: worleyph@ornl.gov, and Chang, C.S., Email: cschang@pppl.gov. A fully nonlinear multispecies Fokker–Planck–Landau collision operator for simulation of fusion plasma. United States: N. p., 2016.
Web. doi:10.1016/J.JCP.2016.03.064.
Hager, Robert, Yoon, E.S., Email: yoone@rpi.edu, Ku, S., Email: sku@pppl.gov, D'Azevedo, E.F., Email: dazevedoef@ornl.gov, Worley, P.H., Email: worleyph@ornl.gov, & Chang, C.S., Email: cschang@pppl.gov. A fully nonlinear multispecies Fokker–Planck–Landau collision operator for simulation of fusion plasma. United States. doi:10.1016/J.JCP.2016.03.064.
Hager, Robert, Yoon, E.S., Email: yoone@rpi.edu, Ku, S., Email: sku@pppl.gov, D'Azevedo, E.F., Email: dazevedoef@ornl.gov, Worley, P.H., Email: worleyph@ornl.gov, and Chang, C.S., Email: cschang@pppl.gov. Wed .
"A fully nonlinear multispecies Fokker–Planck–Landau collision operator for simulation of fusion plasma". United States. doi:10.1016/J.JCP.2016.03.064.
@article{osti_22572321,
title = {A fully nonlinear multispecies Fokker–Planck–Landau collision operator for simulation of fusion plasma},
author = {Hager, Robert and Yoon, E.S., Email: yoone@rpi.edu and Ku, S., Email: sku@pppl.gov and D'Azevedo, E.F., Email: dazevedoef@ornl.gov and Worley, P.H., Email: worleyph@ornl.gov and Chang, C.S., Email: cschang@pppl.gov},
abstractNote = {Fusion edge plasmas can be far from thermal equilibrium and require the use of a nonlinear collision operator for accurate numerical simulations. In this article, the nonlinear singlespecies 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 nonlinear Fokker–Planck–Landau operator in the gyrokinetic particleincell codes XGC1 and XGCa is described and results of a verification study are discussed. Finally, the numerical techniques that make our nonlinear collision operator viable on highperformance 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.},
doi = {10.1016/J.JCP.2016.03.064},
journal = {Journal of Computational Physics},
issn = {00219991},
number = ,
volume = 315,
place = {United States},
year = {2016},
month = {6}
}
Works referencing / citing this record:
Neutral recycling effects on ITG turbulence
journal, July 2017
 Stotler, D. P.; Lang, J.; Chang, C. S.
 Nuclear Fusion, Vol. 57, Issue 8