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

Mollén, Albert ; Adams, Mark F. ; Knepley, Matthew G. ; Hager, Robert ; Chang, C. S.
Issue date: 2021
Rights:
Creative Commons Attribution 4.0 International (CC BY)
Cite as:
Mollén, Albert, Adams, Mark F., Knepley, Matthew G., Hager, Robert, & Chang, C. S. (2021). Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC [Data set]. Princeton Plasma Physics Laboratory, Princeton University. https://doi.org/10.11578/1814952
@electronic{molln_albert_2021,
  author      = {Mollén, Albert and
                Adams, Mark F. and
                Knepley, Matthew G. and
                Hager, Robert and
                Chang, C. S.},
  title       = {{Implementation of higher-order velocity
                mapping between marker particles and gri
                d in the particle-in-cell code XGC}},
  publisher   = {{Princeton Plasma Physics Laboratory, Pri
                nceton University}},
  year        = 2021,
  url         = {https://doi.org/10.11578/1814952}
}
Description:

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.

Show More
# Filename Filesize
1 README.txt 646 Bytes
2 ARK_DATA.zip 834 MB