skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A fully non-linear multi-species 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 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.

Authors:
 [1]
  1. 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 0021-9991
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; ALGORITHMS; COLLISIONS; COMPUTERIZED SIMULATION; FOKKER-PLANCK EQUATION; IMPLEMENTATION; NONLINEAR PROBLEMS; PERFORMANCE; PLASMA; THERMAL EQUILIBRIUM

Citation Formats

Hager, Robert, Yoon, E.S., E-mail: yoone@rpi.edu, Ku, S., E-mail: sku@pppl.gov, D'Azevedo, E.F., E-mail: dazevedoef@ornl.gov, Worley, P.H., E-mail: worleyph@ornl.gov, and Chang, C.S., E-mail: cschang@pppl.gov. A fully non-linear multi-species 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., E-mail: yoone@rpi.edu, Ku, S., E-mail: sku@pppl.gov, D'Azevedo, E.F., E-mail: dazevedoef@ornl.gov, Worley, P.H., E-mail: worleyph@ornl.gov, & Chang, C.S., E-mail: cschang@pppl.gov. A fully non-linear multi-species 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., E-mail: yoone@rpi.edu, Ku, S., E-mail: sku@pppl.gov, D'Azevedo, E.F., E-mail: dazevedoef@ornl.gov, Worley, P.H., E-mail: worleyph@ornl.gov, and Chang, C.S., E-mail: cschang@pppl.gov. Wed . "A fully non-linear multi-species 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 non-linear multi-species Fokker–Planck–Landau collision operator for simulation of fusion plasma},
author = {Hager, Robert and Yoon, E.S., E-mail: yoone@rpi.edu and Ku, S., E-mail: sku@pppl.gov and D'Azevedo, E.F., E-mail: dazevedoef@ornl.gov and Worley, P.H., E-mail: worleyph@ornl.gov and Chang, C.S., E-mail: cschang@pppl.gov},
abstractNote = {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.},
doi = {10.1016/J.JCP.2016.03.064},
journal = {Journal of Computational Physics},
issn = {0021-9991},
number = ,
volume = 315,
place = {United States},
year = {2016},
month = {6}
}

Works referencing / citing this record:

Neutral recycling effects on ITG turbulence
journal, July 2017