Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC

Mollén, Albert; Adams, M. F.; Knepley, M. G.; Hager, R.; Chang, C. S.

doi:10.1017/s0022377821000441

Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC

Journal Article · Tue May 04 00:00:00 EDT 2021 · Journal of Plasma Physics

DOI:https://doi.org/10.1017/s0022377821000441· OSTI ID:1804610

^[1]; ^[2]; ^[3]; ^[1]; ^[1]

Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
State Univ. of New York (SUNY), Buffalo, NY (United States)

The global total-f gyrokinetic particle-in-cell code XGC, used to study transport in magnetic fusion plasmas or to couple with a core gyrokinetic code while functioning as an edge gyrokinetic code, implements a five-dimensional continuum grid to perform the dissipative operations, such as plasma collisions, or to exchange the particle distribution function information with a core code. To transfer the distribution function between marker particles and a rectangular two-dimensional velocity-space grid, XGC employs a bilinear mapping. The conservation of particle density and momentum is accurate enough in this bilinear operation, but the error in the particle energy conservation can become undesirably large and cause non-negligible numerical heating in a steep edge pedestal. In the present work we update XGC to use a novel mapping technique, based on the calculation of a pseudo-inverse, to exactly preserve moments up to the order of the discretization space. Here we describe the details of the implementation and we demonstrate the reduced interpolation error for a tokamak test plasma using first- and second-order elements with the pseudo-inverse method and comparing with the bilinear mapping.

View Accepted Manuscript (DOE)

Research Organization:: Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States). National Energy Research Scientific Computing Center (NERSC); Princeton Plasma Physics Laboratory (PPPL), Princeton, NJ (United States)

Sponsoring Organization:: USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR); USDOE Office of Science (SC), Fusion Energy Sciences (FES)

Grant/Contract Number:: AC02-05CH11231; AC02-09CH11466

OSTI ID:: 1804610

Journal Information:: Journal of Plasma Physics, Journal Name: Journal of Plasma Physics Journal Issue: 2 Vol. 87; ISSN 0022-3778

Publisher:: Cambridge University PressCopyright Statement

Country of Publication:: United States

Language:: English

References (15)

The Semi-Lagrangian Method for the Numerical Resolution of the Vlasov Equation Sonnendrücker, Eric; Roche, Jean; Bertrand, Pierre Journal of Computational Physics, Vol. 149, Issue 2 https://doi.org/10.1006/jcph.1998.6148	journal	March 1999
Conservative semi-Lagrangian schemes for Vlasov equations Crouseilles, Nicolas; Mehrenberger, Michel; Sonnendrücker, Eric Journal of Computational Physics, Vol. 229, Issue 6 https://doi.org/10.1016/j.jcp.2009.11.007	journal	March 2010
Development of semi-Lagrangian gyrokinetic code for full-f turbulence simulation in general tokamak geometry Kwon, Jae-Min; Yi, Dokkyun; Piao, Xiangfan Journal of Computational Physics, Vol. 283 https://doi.org/10.1016/j.jcp.2014.12.017	journal	February 2015
A new hybrid-Lagrangian numerical scheme for gyrokinetic simulation of tokamak edge plasma Ku, S.; Hager, R.; Chang, C. S. Journal of Computational Physics, Vol. 315 https://doi.org/10.1016/j.jcp.2016.03.062	journal	June 2016
A fully non-linear multi-species Fokker–Planck–Landau collision operator for simulation of fusion plasma Hager, Robert; Yoon, E. S.; Ku, S. Journal of Computational Physics, Vol. 315 https://doi.org/10.1016/j.jcp.2016.03.064	journal	June 2016
Moment preserving constrained resampling with applications to particle-in-cell methods Faghihi, D.; Carey, V.; Michoski, C. Journal of Computational Physics, Vol. 409 https://doi.org/10.1016/j.jcp.2020.109317	journal	May 2020
A generalized inverse for matrices Penrose, R. Mathematical Proceedings of the Cambridge Philosophical Society, Vol. 51, Issue 3 https://doi.org/10.1017/S0305004100030401	journal	July 1955
Compressed ion temperature gradient turbulence in diverted tokamak edge Chang, C. S.; Ku, S.; Diamond, P. H. Physics of Plasmas, Vol. 16, Issue 5 https://doi.org/10.1063/1.3099329	journal	May 2009
Conservative discretization of the Landau collision integral Hirvijoki, E.; Adams, M. F. Physics of Plasmas, Vol. 24, Issue 3 https://doi.org/10.1063/1.4979122	journal	March 2017
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
A tight-coupling scheme sharing minimum information across a spatial interface between gyrokinetic turbulence codes Dominski, J.; Ku, S.; Chang, C. -S. Physics of Plasmas, Vol. 25, Issue 7 https://doi.org/10.1063/1.5044707	journal	July 2018
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
Spatial coupling of gyrokinetic simulations, a generalized scheme based on first-principles Dominski, J.; Cheng, J.; Merlo, G. Physics of Plasmas, Vol. 28, Issue 2 https://doi.org/10.1063/5.0027160	journal	February 2021
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
Numerical heating in particle-in-cell simulations with Monte Carlo binary collisions Alves, E. P.; Mori, W. B.; Fiuza, F. Physical Review E, Vol. 103, Issue 1 https://doi.org/10.1103/PhysRevE.103.013306	journal	January 2021

Similar Records

Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC

Dataset · Thu Apr 01 00:00:00 EDT 2021 · OSTI ID:1814952

First coupled GENE–XGC microturbulence simulations

Journal Article · Thu Dec 31 19:00:00 EST 2020 · Physics of Plasmas · OSTI ID:1818685

Related Subjects

70 PLASMA PHYSICS AND FUSION TECHNOLOGY
fusion plasma
plasma simulation

Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC

Citation Formats

References (15)

Similar Records

Related Subjects