Accuracy of momentum and gyrodensity transport in global gyrokinetic particle-in-cell simulations
- Centre for Fusion, Space and Astrophysics, Department of Physics, University of Warwick, CV4 7AL Coventry (United Kingdom)
- Centre de Recherches en Physique des Plasmas, Association Euratom-Confédération Suisse, Ecole Polytechnique Fédérale de Lausanne, PPB, 1015 Lausanne (Switzerland)
Gyrokinetic Particle-In-Cell (PIC) simulations based on conservative Lagrangian formalisms admit transport equations for conserved quantities such as gyrodensity and toroidal momentum, and these can be derived for arbitrary wavelength, even though previous applications have used the long-wavelength approximation. In control-variate PIC simulations, a consequence of the different treatment of the background (f{sub 0}) and perturbed parts (δf), when a splitting f = f{sub 0} + δf is performed, is that analytical transport relations for the relevant fluxes and moments are only reproduced in the large marker number limit. The transport equations for f can be used to write the inconsistency in the perturbed quantities explicitly in terms of the sampling of the background distribution f{sub 0}. This immediately allows estimates of the error in consistency of momentum transport in control-variate PIC simulations. This inconsistency tends to accumulate secularly and is not directly affected by the sources and noise control in the system. Although physical tokamaks often rotate quite strongly, the standard gyrokinetic formalism assumes weak perpendicular flows, comparable to the drift speed. For systems with such weak flows, maintaining acceptably small relative errors requires that a number of markers scale with the fourth power of the linear system size to consistently resolve long-wavelength evolution. To avoid this unfavourable scaling, an algorithm for exact gyrodensity transport has been developed, and this is shown to allow accurate simulations with an order of magnitude fewer markers.
- OSTI ID:
- 22252934
- Journal Information:
- Physics of Plasmas, Vol. 21, Issue 5; Other Information: (c) 2014 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA); ISSN 1070-664X
- Country of Publication:
- United States
- Language:
- English
Similar Records
3D gyrokinetic particle-in-cell simulation of fusion plasma microturbulence on parallel computers
Center for Gyrokinetic/MHD Hybrid Simulation of Energetic Particle Physics in Toroidal Plasmas (CSEPP). Final report