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:
-
- Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
- Rensselaer Polytechnic Inst., Troy, NY (United States)
- 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}
}
Web of Science
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
Monte Carlo evaluation of transport coefficients
journal, January 1981
- Boozer, Allen H.
- Physics of Fluids, Vol. 24, Issue 5
A new bdelta f method for neoclassical transport studies
journal, August 1999
- Wang, W. X.; Nakajima, N.; Okamoto, M.
- Plasma Physics and Controlled Fusion, Vol. 41, Issue 9
Full linearized Fokker–Planck collisions in neoclassical transport simulations
journal, December 2011
- Belli, E. A.; Candy, J.
- Plasma Physics and Controlled Fusion, Vol. 54, Issue 1
A mass, momentum, and energy conserving, fully implicit, scalable algorithm for the multi-dimensional, multi-species Rosenbluth–Fokker–Planck equation
journal, September 2015
- Taitano, W. T.; Chacón, L.; Simakov, A. N.
- Journal of Computational Physics, Vol. 297
Fast elliptic solvers in cylindrical coordinates and the Coulomb collision operator
journal, September 2011
- Pataki, Andras; Greengard, Leslie
- Journal of Computational Physics, Vol. 230, Issue 21
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
- Ku, S.; Chang, C. S.; Diamond, P. H.
- Nuclear Fusion, Vol. 49, Issue 11
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
An efficient explicit numerical scheme for diffusion-type equations with a highly inhomogeneous and highly anisotropic diffusion tensor
journal, April 2007
- Larroche, O.
- Journal of Computational Physics, Vol. 223, Issue 1
Fokker–Planck kinetic modeling of suprathermal α -particles in a fusion plasma
journal, December 2014
- Peigney, B. E.; Larroche, O.; Tikhonchuk, V.
- Journal of Computational Physics, Vol. 278
A code for the solution of the Vlasov-Fokker-Planck equation in 1-D or 2-D
journal, December 1988
- Epperlein, E. M.; Rickard, G. J.; Bell, A. R.
- Computer Physics Communications, Vol. 52, Issue 1
Fokker–Planck modeling of electron transport in laser-produced plasmas
journal, June 1994
- Epperlein, E. M.
- Laser and Particle Beams, Vol. 12, Issue 2
Theory of plasma transport in toroidal confinement systems
journal, April 1976
- Hinton, F. L.; Hazeltine, R. D.
- Reviews of Modern Physics, Vol. 48, Issue 2
Kinetic calculation of neoclassical transport including self-consistent electron and impurity dynamics
journal, July 2008
- Belli, E. A.; Candy, J.
- Plasma Physics and Controlled Fusion, Vol. 50, Issue 9
Scaling to 150K cores: Recent algorithm and performance engineering developments enabling XGC1 to run at scale
journal, July 2009
- Adams, Mark F.; Ku, Seung-Hoe; Worley, Patrick
- Journal of Physics: Conference Series, Vol. 180
Mesh generation for confined fusion plasma simulation
journal, September 2015
- Zhang, Fan; Hager, Robert; Ku, Seung-Hoe
- Engineering with Computers, Vol. 32, Issue 2
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
Gyrokinetic neoclassical study of the bootstrap current in the tokamak edge pedestal with fully non-linear Coulomb collisions
journal, April 2016
- Hager, Robert; Chang, C. S.
- Physics of Plasmas, Vol. 23, Issue 4
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
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
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
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
Development of a Gyrokinetic Particle-in-Cell Code for Whole-Volume Modeling of Stellarators
journal, May 2019
- Moritaka, Toseo; Hager, Robert; Cole, Michael
- Plasma, Vol. 2, Issue 2
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
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
Differential formulation of the gyrokinetic Landau operator
text, January 2017
- Hirvijoki, Eero; Brizard, Alain J.; Pfefferlé, David
- arXiv