Global and Kinetic MHD Simulation by the Gpic-MHD Code

In order to implement large-scale and high-beta tokamak simulation, a new algorithm of the electromagnetic gyrokinetic PIC (particle-in-cell) code was proposed and installed on the Gpic-MHD code [Gyrokinetic PIC code for magnetohydrodynamic (MHD) simulation]. In the new algorithm, the vorticity equation and the generalized Ohm's law along the magnetic field are derived from the basic equations of the gyrokinetic Vlasov, Poisson, and Ampere system and are used to describe the spatio-temporal evolution of the field quantities of the electrostatic potential phi and the longitudinal component of the vector potential A(z). The basic algorithm is equivalent to solving the reduced-MHD-type equations with kinetic corrections, in which MHD physics related to Alfven modes are well described. The estimation of perturbed electron pressure from particle dynamics is dominant, while the effects of other moments are negligible. Another advantage of the algorithm is that the longitudinal induced electric field, E-Tz = -partial derivative A(z)/partial derivative t, is explicitly estimated by the generalized Ohm's law and used in the equations of motion. Furthermore, the particle velocities along the magnetic field are used (v(z)-formulation) instead of generalized momentums (p(z)-formulation), hence there is no problem of 'cancellation', which would otherwise appear when A(z) is estimated from the Ampere's law in the p(z)-formulation. The successful simulation of the collisionless internal kink mode by the new Gpic-MHD with realistic values of the large-scale and high-beta tokamaks revealed the usefulness of the new algorithm.