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. The non-linear single-species Fokker–Planck–Landau collision operator developed by Yoon and Chang (2014) [9] is generalized to include multiple particle species. Moreover, 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. As a result, the collision operator's good weak and strong scaling behavior are shown.

Authors:
 [1];  [2];  [1];  [3];  [3];  [1]
  1. Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
  2. Rensselaer Polytechnic Inst., Troy, NY (United States)
  3. Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
Publication Date:
Research Org.:
Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States). Oak Ridge Leadership Computing Facility (OLCF); Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
Sponsoring Org.:
USDOE Office of Science (SC)
OSTI Identifier:
1254653
Alternate Identifier(s):
OSTI ID: 1295136; OSTI ID: 1325292
Report Number(s):
PPPL-5253
Journal ID: ISSN 0021-9991; PII: S0021999116300298
Grant/Contract Number:  
AC02-09CH11466; AC05-00OR22725; SC0008449; AC02-06CH11357; AC02-05CH11231
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 315; Journal Issue: C; Journal ID: ISSN 0021-9991
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
70 PLASMA PHYSICS AND FUSION TECHNOLOGY; plasma; fusion; collision operator; XGC; particle-in-cell; 97 MATHEMATICS AND COMPUTING

Citation Formats

Hager, Robert, Yoon, E. S., Ku, S., D'Azevedo, E. F., Worley, P. H., and Chang, C. S. 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., Ku, S., D'Azevedo, E. F., Worley, P. H., & Chang, C. S. A fully non-linear multi-species Fokker–Planck–Landau collision operator for simulation of fusion plasma. United States. https://doi.org/10.1016/j.jcp.2016.03.064
Hager, Robert, Yoon, E. S., Ku, S., D'Azevedo, E. F., Worley, P. H., and Chang, C. S. 2016. "A fully non-linear multi-species Fokker–Planck–Landau collision operator for simulation of fusion plasma". United States. https://doi.org/10.1016/j.jcp.2016.03.064. https://www.osti.gov/servlets/purl/1254653.
@article{osti_1254653,
title = {A fully non-linear multi-species Fokker–Planck–Landau collision operator for simulation of fusion plasma},
author = {Hager, Robert and Yoon, E. S. and Ku, S. and D'Azevedo, E. F. and Worley, P. H. and Chang, C. S.},
abstractNote = {Fusion edge plasmas can be far from thermal equilibrium and require the use of a non-linear collision operator for accurate numerical simulations. The non-linear single-species Fokker–Planck–Landau collision operator developed by Yoon and Chang (2014) [9] is generalized to include multiple particle species. Moreover, 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. As a result, the collision operator's good weak and strong scaling behavior are shown.},
doi = {10.1016/j.jcp.2016.03.064},
url = {https://www.osti.gov/biblio/1254653}, journal = {Journal of Computational Physics},
issn = {0021-9991},
number = C,
volume = 315,
place = {United States},
year = {Mon Apr 04 00:00:00 EDT 2016},
month = {Mon Apr 04 00:00:00 EDT 2016}
}

Journal Article:

Citation Metrics:
Cited by: 44 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Approximate Fokker–Planck collision operator for transport theory applications
journal, January 1976


Monte Carlo evaluation of transport coefficients
journal, January 1981


A new bdelta f method for neoclassical transport studies
journal, August 1999


Full linearized Fokker–Planck collisions in neoclassical transport simulations
journal, December 2011


Fast elliptic solvers in cylindrical coordinates and the Coulomb collision operator
journal, September 2011


Fokker-Planck Equation for an Inverse-Square Force
journal, July 1957


Fokker–Planck kinetic modeling of suprathermal α -particles in a fusion plasma
journal, December 2014


A code for the solution of the Vlasov-Fokker-Planck equation in 1-D or 2-D
journal, December 1988


Fokker–Planck modeling of electron transport in laser-produced plasmas
journal, June 1994


Theory of plasma transport in toroidal confinement systems
journal, April 1976


Kinetic calculation of neoclassical transport including self-consistent electron and impurity dynamics
journal, July 2008


Scaling to 150K cores: Recent algorithm and performance engineering developments enabling XGC1 to run at scale
journal, July 2009


Mesh generation for confined fusion plasma simulation
journal, September 2015


Works referencing / citing this record:

Neutral recycling effects on ITG turbulence
journal, July 2017


Development and testing of an unstructured mesh method for whole plasma gyrokinetic simulations in realistic tokamak geometry
journal, December 2019


XGC
software, June 2018


Verification of the global gyrokinetic stellarator code XGC-S for linear ion temperature gradient driven modes
journal, August 2019


Full- f version of GENE for turbulence in open-field-line systems
journal, June 2018


Comparative collisionless alpha particle confinement in stellarator reactors with the XGC gyrokinetic code
journal, March 2019


Development of a Gyrokinetic Particle-in-Cell Code for Whole-Volume Modeling of Stellarators
journal, May 2019


A fast low-to-high confinement mode bifurcation dynamics in the boundary-plasma gyrokinetic code XGC1
journal, May 2018


Cross-verification of neoclassical transport solutions from XGCa against NEO
journal, October 2019


Differential formulation of the gyrokinetic Landau operator
text, January 2017