Plasma and Fusion Research
Volume 7, 2401101 (2012)
Regular Articles
- Faculty of Engineering, Hokkaido University, N-13, W-8, Sapporo 060-8628, Japan
- 1)
- Graduate School of Engineering, Hokkaido University, N-13, W-8, Sapporo 060-8628, Japan
Abstract
The BIA (Binary Interaction Approximation) formulation in the presence of neutralizing immovable background ion is presented for analysis of multiple electron motion. Such a BIA scheme is applied to electrons in plasmas. A test calculation shows that 1) the plasma oscillation and its frequency are successfully detected, 2) the CPU time for the BIA are less than 1.5 sec and 1 hour for two and three dimensional analysis, while 127 sec and 13 hours for the direct integration method (DIM) by using a Runge-Kutta-Fehlberg integrator with an absolute error tolerance of 10−16, and 3) the number of time steps for the DIM, in such a case, are as many as 5.8 × 104 and 3.6 × 106, while those for the BIA are only 256 and 512.
Keywords
N-body problem, binary interaction approximation, plasma oscillation, background ion, Fourier analysis
Full Text
References
- [1] S. Oikawa and H. Funasaka, Plasma Fusion Res. 5, S1051 (2010).
- [2] S. Oikawa, K. Higashi and H. Funasaka, Plasma Fusion Res. 5, S1048 (2010).
- [3] S. Oikawa and H. Funasaka, Plasma Fusion Res. 3, S1073 (2008).
- [4] A.W. Appel, SIAM J. Sci. Stat. Comput. 6, 85 (1985).
- [5] J.E. Barnes and P. Hutt, Nature 324, 446 (1986).
- [6] P.P. Brieu, F.J. Summers and J.P. Ostriker, APJ 453, 566 (1995).
- [7] J. Makino, M. Taiji, T. Ebisuzaki and D. Sugimoto, APJ 480, 432 (1997).
- [8] E. Fehlberg, “Low-order classical Runge-Kutta formulas with step size control and their application to some heat transfer problems,” NASA Technical Report 315 (1969).
This paper may be cited as follows:
Shun-ichi OIKAWA, Takanori KAMEI and Masaki GOTO, Plasma Fusion Res. 7, 2401101 (2012).