Plasma and Fusion Research
Volume 5, S1051 (2010)
Regular Articles
- Graduate School of Engineering, Hokkaido University, Sapporo 060-8628, Japan
Abstract
The binary interaction approximation (BIA) to N-body problems is proposed. The BIA conserves total linear momenta in principle. Other invariants, such as the total angular momentum and total energy, are conserved to at least 12 effective digits for a two-dimensional hydrogen plasma of T = 10 keV and n = 1020 m−3. For such a plasma, the total CPU time of the BIA is found to scale as approximately N1.9, while the conventional direct integration method scales as approximately N3.
Keywords
N-body problem, algebraic approximation, binary interaction approximation, variable step size, parallel computation.
Full Text
References
- [1] A.W. Appel, An Efficient Program for Many-Body Simulation, SIAM J. Sci. Stat. Comput. 6, 85 (1985).
- [2] J.E. Barnes and P. Hutt, A hierarchical O(NlogN) forcecalculation algorithm, Nature 324, 446 (1986).
- [3] P.P. Brieu, F.J. Summers and J.P. Ostriker, Cosmological Simulations Using Special Purpose Computers: Implementing P3M on GRAPE, APJ 453, 566 (1995).
- [4] J. Makino, M. Taiji, T. Ebisuzaki, and D. Sugimoto, GRAPE-4: A Massively Parallel Special-Purpose Computer for Collisional N-Body Simulations, APJ 480, 432 (1997).
- [5] S. Oikawa and H. Funasaka, Algebraic analysis approach for multibody problems, Plasma Fusion Res. 36, S1073 (2008).
- [6] S. Oikawa, K. Higashi, H. Funasaka and Y. Kitagawa, Algebraic analysis approach for multibody problems II: Variance of velocity changes, ITC-18, P2-58(2008); also submitted to PFR.
- [7] 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).
- [8] E. Fehlberg, Klassische Runge-Kutta-Formeln vierter und niedrigerer Ordnung mit Schrittweiten-Kontrolle und ihre Anwendung auf Warmeleitungsprobleme, Computing (Arch. Elektron. Rechnen), vol. 6, 61 (1970).
This paper may be cited as follows:
Shun-ichi OIKAWA and Hideo FUNASAKA, Plasma Fusion Res. 5, S1051 (2010).