Abstract
The direct implicit method with a second-order implicit integration scheme is formulated for and applied to the electron parallel nonlinearity in global electrostatic gyrokinetic particle-in-cell simulations of toroidal fusion plasmas. The method shows improved numerical accuracy and stability properties compared to the direct implicit method with a first-order integration scheme. The conservation of total energy and toroidal angular momentum are analyzed by both techniques and the results are presented.
Original language | English |
---|---|
Pages (from-to) | 22-29 |
Journal | Journal of Computational Physics |
Volume | 239 |
DOIs | |
Publication status | Published - 2013 |
MoE publication type | A1 Journal article-refereed |
Keywords
- fusion plasma physics
- gyrokinetic particle simulation
- kinetic electrons