A fully non-linear multi-species Fokker–Planck–Landau collision operator for simulation of fusion plasma
- 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)
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
Web of Science
Similar Records
Encoder–decoder neural network for solving the nonlinear Fokker–Planck–Landau collision operator in XGC
A Study of Turbulence-Driven Inward Momentum Flux & Construction of a Difference Scheme for Fokker-Planck-Landau Operator