Home > Accuracy of momentum and gyrodensity transport in global gyrokinetic particle-in-cell simulations > HTML MARC |

000201458 001__ 201458 000201458 005__ 20190228220028.0 000201458 0247_ $$2doi$$a10.1063/1.4873387 000201458 022__ $$a1070-664X 000201458 02470 $$2ISI$$a000337107200029 000201458 037__ $$aARTICLE 000201458 245__ $$aAccuracy of momentum and gyrodensity transport in global gyrokinetic particle-in-cell simulations 000201458 260__ $$aMelville$$bAmer Inst Physics$$c2014 000201458 269__ $$a2014 000201458 300__ $$a10 000201458 336__ $$aJournal Articles 000201458 520__ $$aGyrokinetic 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(0)) and perturbed parts (delta f), when a splitting f = f(0) + delta 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(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. (C) 2014 AIP Publishing LLC. 000201458 700__ $$aMcmillan, B. F. 000201458 700__ $$aVillard, L. 000201458 773__ $$j21$$k5$$q052501$$tPhysics Of Plasmas 000201458 909C0 $$pCRPP 000201458 909C0 $$0252028$$pSPC$$xU12272$$xU12268$$xU10558$$xU10635$$xU12266$$xU10636$$xU10137$$xU12270$$xU10557$$xU12273$$xU10559$$xU12271$$xU12269$$xU12267$$xU10136 000201458 909CO $$ooai:infoscience.tind.io:201458$$pSB$$particle 000201458 917Z8 $$x112823 000201458 917Z8 $$x105317 000201458 937__ $$aEPFL-ARTICLE-201458 000201458 973__ $$aEPFL$$rREVIEWED$$sPUBLISHED 000201458 980__ $$aARTICLE