Skip navigation
Please use this identifier to cite or link to this item:
Title: Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC
Contributors: Mollen Albert
Adams Mark F.
Knepley Matthew G.
Hager Robert
Chang C. S.
U. S. Department of Energy
Keywords: Fusion Plasma
Plasma Simulation
Plasma Transport
Issue Date: Mar-2021
Publisher: Princeton Plasma Physics Laboratory, Princeton University
Related Publication: Journal of Plasma Physics
Abstract: 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 5-dimensional (5D) 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 2D 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. We describe the details of the implementation and we demonstrate the reduced interpolation error for a tokamak test plasma by using 1st- and 2nd-order elements with the pseudo-inverse method and comparing to the bilinear mapping.
Appears in Collections:Theory and Computation

Files in This Item:
File Description SizeFormat 
README.txt646 BTextView/Download
ARK_DATA.zip814.65 MBUnknownView/Download

Items in Dataspace are protected by copyright, with all rights reserved, unless otherwise indicated.