U.S. Department of EnergyOffice of Scientific and Technical Information
Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC
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.
- Authors:
-
- Princeton University (PPPL)
- Publication Date:
- DOE Contract Number:
- AC02-09CH11466
- Research Org.:
- Princeton Plasma Physics Laboratory (PPPL), Princeton, NJ (United States)
- Sponsoring Org.:
- USDOE Office of Science (SC)
- Subject:
- Fusion Plasma; Particle-In-Cell; Plasma Simulation; Plasma Transport; XGC
- OSTI Identifier:
- 1814952
- DOI:
- https://doi.org/10.11578/1814952
Citation Formats
Albert, Mollen, F, Adams Mark, G, Knepley Matthew, Robert, Hager, and S, Chang C. Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC. United States: N. p., 2021.
Web. doi:10.11578/1814952.
Albert, Mollen, F, Adams Mark, G, Knepley Matthew, Robert, Hager, & S, Chang C. Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC. United States. doi:https://doi.org/10.11578/1814952
Albert, Mollen, F, Adams Mark, G, Knepley Matthew, Robert, Hager, and S, Chang C. 2021.
"Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC". United States. doi:https://doi.org/10.11578/1814952. https://www.osti.gov/servlets/purl/1814952. Pub date:Thu Apr 01 04:00:00 UTC 2021
@article{osti_1814952,
title = {Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC},
author = {Albert, Mollen and F, Adams Mark and G, Knepley Matthew and Robert, Hager and S, Chang C},
abstractNote = {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.},
doi = {10.11578/1814952},
journal = {},
number = ,
volume = ,
place = {United States},
year = {Thu Apr 01 04:00:00 UTC 2021},
month = {Thu Apr 01 04:00:00 UTC 2021}
}