Energy conservation in a nonlinear gyrokinetic particle-in-cell code for ion-temperature-gradient-driven modes in θ-pinch geometry

Abstract

A global nonlinear simulation code for the time evolution of ion-temperature-gradient-driven modes in θ-pinch geometry as a first approximation to the stellarator Wendelstein 7-X (W7-X) [Grieger et al., Proceedings of the 13th International Conference on Plasma Physics and Controlled Nuclear Fusion Research, Washington, DC, 1990 (International Atomic Energy Agency, Vienna, 1991), Vol. 3, p. 525] has been developed. A δf particle-in-cell (PIC) method is used to solve the coupled system of gyrokinetic equations for the ions, in the electrostatic approximation, and the quasineutrality equation, assuming adiabatically responding electrons. The focus has been on adherence to conservation laws, i.e., particle number and energy conservation. Besides other improvements it has been shown that a well-chosen initial distribution of the markers in reduced phase space makes optimal use of the δf PIC method to reduce the statistical noise for a given number of markers. In a model including all (1351) physically relevant modes, it has been possible to achieve energy conservation beyond the saturation of the instability.

Global particle-in-cell simulations of fast-particle effects on shear Alfvén waves

“…The integrals over the gyroangle Eq. (4) are approximated with an N-point discrete sum [29,31,32]. The cancellation problem [27,33] is solved using the iterative scheme No.…”

First principles based simulations of instabilities and turbulence

“…There is great interest in these phenomena, especially in the context of magnetic confinement fusion research, which has lasted for several decades and is still growing. The understanding of anomalous transport has greatly improved over the past years, thanks in large part to the role of large-scale numerical simulations [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16]. Predictions used for the design of the next generation of experimental fusion reactors have been largely made on extrapolations from a range of existing experiments worldwide and assumed semi-empirical scaling laws.…”

Electromagnetic gyrokinetic PIC simulation with an adjustable control variates method

“…Furthermore, in PIC simulation the gyro-center distribution is represented by a sum over δ functions at the marker positions [15], [16] …”