DOE PAGES 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:
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. Mon . "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},
journal = {Journal of Computational Physics},
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

  • Hirshman, S. P.; Sigmar, D. J.
  • Physics of Fluids, Vol. 19, Issue 10
  • DOI: 10.1063/1.861356

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


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

  • Rosenbluth, Marshall N.; MacDonald, William M.; Judd, David L.
  • Physical Review, Vol. 107, Issue 1
  • DOI: 10.1103/PhysRev.107.1

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

  • Lu, Z. X.; Lauber, Ph.; Hayward-Schneider, T.
  • Physics of Plasmas, Vol. 26, Issue 12
  • DOI: 10.1063/1.5124376

XGC
software, June 2018


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

  • Cole, M. D. J.; Hager, R.; Moritaka, T.
  • Physics of Plasmas, Vol. 26, Issue 8
  • DOI: 10.1063/1.5109259

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

  • Pan, Q.; Told, D.; Shi, E. L.
  • Physics of Plasmas, Vol. 25, Issue 6
  • DOI: 10.1063/1.5008895

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

  • Cole, M. D. J.; Hager, R.; Moritaka, T.
  • Physics of Plasmas, Vol. 26, Issue 3
  • DOI: 10.1063/1.5085349

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

  • Ku, S.; Chang, C. S.; Hager, R.
  • Physics of Plasmas, Vol. 25, Issue 5
  • DOI: 10.1063/1.5020792

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

  • Hager, R.; Dominski, J.; Chang, C. S.
  • Physics of Plasmas, Vol. 26, Issue 10
  • DOI: 10.1063/1.5121308

Differential formulation of the gyrokinetic Landau operator
text, January 2017