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

Journal Article · · Journal of Computational Physics
 [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)

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.

Research Organization:
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 Organization:
USDOE Office of Science (SC)
Grant/Contract Number:
AC02-09CH11466; AC05-00OR22725; SC0008449; AC02-06CH11357; AC02-05CH11231
OSTI ID:
1254653
Alternate ID(s):
OSTI ID: 1295136; OSTI ID: 1325292
Report Number(s):
PPPL-5253; PII: S0021999116300298
Journal Information:
Journal of Computational Physics, Vol. 315, Issue C; ISSN 0021-9991
Publisher:
ElsevierCopyright Statement
Country of Publication:
United States
Language:
English
Citation Metrics:
Cited by: 44 works
Citation information provided by
Web of Science

References (16)

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
A mass, momentum, and energy conserving, fully implicit, scalable algorithm for the multi-dimensional, multi-species Rosenbluth–Fokker–Planck equation journal September 2015
Fast elliptic solvers in cylindrical coordinates and the Coulomb collision operator journal September 2011
Full-f gyrokinetic particle simulation of centrally heated global ITG turbulence from magnetic axis to edge pedestal top in a realistic tokamak geometry journal September 2009
Fokker-Planck Equation for an Inverse-Square Force journal July 1957
An efficient explicit numerical scheme for diffusion-type equations with a highly inhomogeneous and highly anisotropic diffusion tensor journal April 2007
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

Cited By (11)

Neutral recycling effects on ITG turbulence journal July 2017
Gyrokinetic neoclassical study of the bootstrap current in the tokamak edge pedestal with fully non-linear Coulomb collisions journal April 2016
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