A fully non-linear multi-species Fokker–Planck–Landau collision operator for simulation of fusion plasma

Hager, Robert; Yoon, E. S.; Ku, S.; D'Azevedo, E. F.; Worley, P. H.; Chang, C. S.

doi:10.1016/j.jcp.2016.03.064

A fully non-linear multi-species Fokker–Planck–Landau collision operator for simulation of fusion plasma

Journal Article · Mon Apr 04 00:00:00 EDT 2016 · Journal of Computational Physics

DOI:https://doi.org/10.1016/j.jcp.2016.03.064· OSTI ID:1254653

Hager, Robert ^[1]; Yoon, E. S. ^[2]; Ku, S. ^[1]; D'Azevedo, E. F. ^[3]; Worley, P. H. ^[3]; Chang, C. S. ^[1]

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.

View Accepted Manuscript (DOE)

Research Organization:: Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States). Oak Ridge Leadership Computing Facility (OLCF); Princeton Plasma Physics Laboratory (PPPL), Princeton, NJ (United States)

Sponsoring Organization:: DOE Office of Science; USDOE

Grant/Contract Number:: AC02-09CH11466; AC05-00OR22725; SC0008449; AC02-06CH11357; AC05-00OR22725; AC02-05CH11231

OSTI ID:: 1254653

Alternate ID(s):: OSTI ID: 1325292
OSTI ID: 22572321
OSTI ID: 1295136

Report Number(s):: PPPL--5253; PII: S0021999116300298

Journal Information:: Journal of Computational Physics, Journal Name: Journal of Computational Physics Journal Issue: C Vol. 315; ISSN 0021-9991

Publisher:: ElsevierCopyright Statement

Country of Publication:: United States

Language:: English

References (16)

Mesh generation for confined fusion plasma simulation Zhang, Fan; Hager, Robert; Ku, Seung-Hoe Engineering with Computers, Vol. 32, Issue 2 https://doi.org/10.1007/s00366-015-0417-y	journal	September 2015
A code for the solution of the Vlasov-Fokker-Planck equation in 1-D or 2-D Epperlein, E. M.; Rickard, G. J.; Bell, A. R. Computer Physics Communications, Vol. 52, Issue 1 https://doi.org/10.1016/0010-4655(88)90165-8	journal	December 1988
An efficient explicit numerical scheme for diffusion-type equations with a highly inhomogeneous and highly anisotropic diffusion tensor Larroche, O. Journal of Computational Physics, Vol. 223, Issue 1 https://doi.org/10.1016/j.jcp.2006.09.016	journal	April 2007
Fast elliptic solvers in cylindrical coordinates and the Coulomb collision operator Pataki, Andras; Greengard, Leslie Journal of Computational Physics, Vol. 230, Issue 21 https://doi.org/10.1016/j.jcp.2011.07.005	journal	September 2011
Fokker–Planck kinetic modeling of suprathermal α -particles in a fusion plasma Peigney, B. E.; Larroche, O.; Tikhonchuk, V. Journal of Computational Physics, Vol. 278 https://doi.org/10.1016/j.jcp.2014.08.033	journal	December 2014
A mass, momentum, and energy conserving, fully implicit, scalable algorithm for the multi-dimensional, multi-species Rosenbluth–Fokker–Planck equation Taitano, W. T.; Chacón, L.; Simakov, A. N. Journal of Computational Physics, Vol. 297 https://doi.org/10.1016/j.jcp.2015.05.025	journal	September 2015
Fokker–Planck modeling of electron transport in laser-produced plasmas Epperlein, E. M. Laser and Particle Beams, Vol. 12, Issue 2 https://doi.org/10.1017/S0263034600007722	journal	June 1994
Approximate Fokker–Planck collision operator for transport theory applications Hirshman, S. P.; Sigmar, D. J. Physics of Fluids, Vol. 19, Issue 10 https://doi.org/10.1063/1.861356	journal	January 1976
Monte Carlo evaluation of transport coefficients Boozer, Allen H. Physics of Fluids, Vol. 24, Issue 5 https://doi.org/10.1063/1.863445	journal	January 1981
Full-f gyrokinetic particle simulation of centrally heated global ITG turbulence from magnetic axis to edge pedestal top in a realistic tokamak geometry Ku, S.; Chang, C. S.; Diamond, P. H. Nuclear Fusion, Vol. 49, Issue 11 https://doi.org/10.1088/0029-5515/49/11/115021	journal	September 2009
A new bdelta f method for neoclassical transport studies Wang, W. X.; Nakajima, N.; Okamoto, M. Plasma Physics and Controlled Fusion, Vol. 41, Issue 9 https://doi.org/10.1088/0741-3335/41/9/303	journal	August 1999
Kinetic calculation of neoclassical transport including self-consistent electron and impurity dynamics Belli, E. A.; Candy, J. Plasma Physics and Controlled Fusion, Vol. 50, Issue 9 https://doi.org/10.1088/0741-3335/50/9/095010	journal	July 2008
Full linearized Fokker–Planck collisions in neoclassical transport simulations Belli, E. A.; Candy, J. Plasma Physics and Controlled Fusion, Vol. 54, Issue 1 https://doi.org/10.1088/0741-3335/54/1/015015	journal	December 2011
Scaling to 150K cores: Recent algorithm and performance engineering developments enabling XGC1 to run at scale Adams, Mark F.; Ku, Seung-Hoe; Worley, Patrick Journal of Physics: Conference Series, Vol. 180 https://doi.org/10.1088/1742-6596/180/1/012036	journal	July 2009
Fokker-Planck Equation for an Inverse-Square Force Rosenbluth, Marshall N.; MacDonald, William M.; Judd, David L. Physical Review, Vol. 107, Issue 1 https://doi.org/10.1103/PhysRev.107.1	journal	July 1957
Theory of plasma transport in toroidal confinement systems Hinton, F. L.; Hazeltine, R. D. Reviews of Modern Physics, Vol. 48, Issue 2 https://doi.org/10.1103/RevModPhys.48.239	journal	April 1976

Cited By (15)

XGC Koh, Sehoon; Ku, Seung-Hoe; Yoon, Eisung DOE CODE https://doi.org/10.11578/dc.20180627.11	software	June 2018
Analyzing Performance of Selected NESAP Applications on the Cori HPC System Kurth, Thorsten; Arndt, William; Barnes, Taylor High Performance Computing https://doi.org/10.1007/978-3-319-67630-2_25	book	January 2017
Differential formulation of the gyrokinetic Landau operator Hirvijoki, Eero; Brizard, Alain J.; Pfefferlé, David Journal of Plasma Physics, Vol. 83, Issue 1 https://doi.org/10.1017/s0022377816001203	journal	January 2017
Study of up–down poloidal density asymmetry of high- impurities with the new impurity version of XGCa Dominski, J.; Chang, C. S.; Hager, R. Journal of Plasma Physics, Vol. 85, Issue 5 https://doi.org/10.1017/s0022377819000722	journal	October 2019
Gyrokinetic neoclassical study of the bootstrap current in the tokamak edge pedestal with fully non-linear Coulomb collisions Hager, Robert; Chang, C. S. Physics of Plasmas, Vol. 23, Issue 4 https://doi.org/10.1063/1.4945615	journal	April 2016
Full- f version of GENE for turbulence in open-field-line systems Pan, Q.; Told, D.; Shi, E. L. Physics of Plasmas, Vol. 25, Issue 6 https://doi.org/10.1063/1.5008895	journal	June 2018
A fast low-to-high confinement mode bifurcation dynamics in the boundary-plasma gyrokinetic code XGC1 Ku, S.; Chang, C. S.; Hager, R. Physics of Plasmas, Vol. 25, Issue 5 https://doi.org/10.1063/1.5020792	journal	May 2018
Comparative collisionless alpha particle confinement in stellarator reactors with the XGC gyrokinetic code Cole, M. D. J.; Hager, R.; Moritaka, T. Physics of Plasmas, Vol. 26, Issue 3 https://doi.org/10.1063/1.5085349	journal	March 2019
Verification of the global gyrokinetic stellarator code XGC-S for linear ion temperature gradient driven modes Cole, M. D. J.; Hager, R.; Moritaka, T. Physics of Plasmas, Vol. 26, Issue 8 https://doi.org/10.1063/1.5109259	journal	August 2019
Cross-verification of neoclassical transport solutions from XGCa against NEO Hager, R.; Dominski, J.; Chang, C. S. Physics of Plasmas, Vol. 26, Issue 10 https://doi.org/10.1063/1.5121308	journal	October 2019
Development and testing of an unstructured mesh method for whole plasma gyrokinetic simulations in realistic tokamak geometry Lu, Z. X.; Lauber, Ph.; Hayward-Schneider, T. Physics of Plasmas, Vol. 26, Issue 12 https://doi.org/10.1063/1.5124376	journal	December 2019
Neutral recycling effects on ITG turbulence Stotler, D. P.; Lang, J.; Chang, C. S. Nuclear Fusion, Vol. 57, Issue 8 https://doi.org/10.1088/1741-4326/aa7807	journal	July 2017
Fast Low-to-High Confinement Mode Bifurcation Dynamics in a Tokamak Edge Plasma Gyrokinetic Simulation Chang, C. S.; Ku, S.; Tynan, G. R. Physical Review Letters, Vol. 118, Issue 17 https://doi.org/10.1103/physrevlett.118.175001	journal	April 2017
Development of a Gyrokinetic Particle-in-Cell Code for Whole-Volume Modeling of Stellarators Moritaka, Toseo; Hager, Robert; Cole, Michael Plasma, Vol. 2, Issue 2 https://doi.org/10.3390/plasma2020014	journal	May 2019
Differential formulation of the gyrokinetic Landau operator Hirvijoki, Eero; Brizard, Alain J.; Pfefferlé, David arXiv https://doi.org/10.48550/arxiv.1707.05901	text	January 2017

Similar Records

A fully non-linear multi-species Fokker–Planck–Landau collision operator for simulation of fusion plasma

Journal Article · Wed Jun 15 00:00:00 EDT 2016 · Journal of Computational Physics · OSTI ID:22572321

Gyrokinetic Fokker-Planck Collision Operator

Journal Article · Fri May 13 00:00:00 EDT 2011 · Physical Review Letters · OSTI ID:21538335

Related Subjects

70 PLASMA PHYSICS AND FUSION TECHNOLOGY
XGC
collision operator
fusion
particle-in-cell
plasma

A fully non-linear multi-species Fokker–Planck–Landau collision operator for simulation of fusion plasma

Citation Formats

References (16)

Cited By (15)

Similar Records

Related Subjects