Topic/Type: 2.2 Kinetic methods, Particle-In-Cell and Vlasov, Poster
Plasma and Space Science Center, National Cheng Kung University, Tainan 70101, Taiwan
A new algorithm is developed to calculate ion beam trajectories in linear magnetic confinement devices by applying perturbation method . The equation of motion (the newton equation) is solved including the Lorentz force term and Coulomb collisional relaxation term . The frictional force is regarded as a small term, . The basic concept is as follows. Normalizing the equation of motion by ion cyclotron frequency,
for the lowest order equation and
for the first order equation in . The solution then is given by the summation , . The lowest order solution is periodic when is axis-symmetric (Noether\'s theorem) and accordingly the first order solution. The crux in Eq.(5) and (6) are the changes in particle velocity () and particle\'s displacement () both induced by the small friction force (the term). By storing the first periodic motion of both the lowest and the higher order solution, the algorithm can predict further periodic motion and thus can reduce computation time. The perturbation method is also useful since we only need to change the constant when the the plasma parameters change e.g. background densities and temperatures (and do not need to recalculate the whole trajectories). We focus on the high energy particles whose temperature is much larger than that of the background thermal plasma. The scheme is useful in studying ion beam heating and macroscopic stability of plasma columns such as magnetic mirrors and theta pinch plasmas. In general, the algorithm can be applied to periodic motion under perturbative frictional forces, such as guiding center trajectories or satellite motion. This work is partially supported by NCKU Top University Project.
 Y.Nishimura, Diploma thesis, Osaka University (1991).
 K.Miyamoto, Plasma Physics for Nuclear Fusion, (Iwanami, Tokyo, 1986), p.76.
 T.Ishimura, Phys. Fluids 27, 2139 (1984).
 Y.Nishimura, Contributions to Plasma Physics 48, 224 (2008).