Plasma and Fusion Research

Volume 13, 1203088 (2018)

Rapid Communications

An Ingenious Plasma Simulation Model for Large-Scale System with Self-Consistent Electron Dynamics
Tomonori TAKIZUKA, Kenzo IBANO and Masatoshi YAGI1)
Graduate School of Engineering, Osaka University, Suita 565-0871, Japan
National Institutes for Quantum and Radiological Science and Technology, Rokkasho 039-3212, Japan
(Received 12 April 2018 / Accepted 8 June 2018 / Published 10 July 2018)


An ingenious model for large-scale plasma simulations is proposed. An artificial permittivity ε much grater than ε0 (permittivity of free space) is introduced to the Poisson equation, ∇⋅(ε∇φ) = −ρ (φ is electrostatic potential and ρ is charge density), and the Debye length is artificially elongated. A large-scale system with a reasonable number of spatial grids can be treated stably even for including the self-consistent electron dynamics. A dielectric tensor is adopted for the three-dimensional global simulation of tokamak plasmas in the cylindrical coordinates (R, θ, Z). A diagonal element εθθ for toroidal axis θ is set much larger than those for the poloidal plane, εRR = εZZ = ε. Number of toroidal meshes Nθ can be reduced enough smaller than that of poloidal-plane meshes, NR and NZ, in comparison with Nθ > NR, NZ for the usual case of isotropic ε.


simulation model, global plasma, electron dynamics, Poisson equation

DOI: 10.1585/pfr.13.1203088


  • [1] T. Takizuka, Plasma Sci. Technol. 13, 316 (2011).
  • [2] T. Takizuka, Plasma Phys. Control. Fusion 59, 034008 (2017).
  • [3] M. Hosokawa, A. Loarte, G.T.A. Huijsmans, T. Takizuka and N. Hayashi, Plasma Fusion Res. 11, 1403104 (2016).
  • [4] M. Hosokawa et al., 16th Int. Workshop on Plasma Edge Theory in Fusion Devices (Marseille, France, 2017).
  • [5] T. Takizuka, M. Hosokawa, A. Loarte and N. Hayashi, Proc. 15th BPSI meeting (Kasuga, 2017),
  • [6] R.W. Hockney and J.W. Eastwood, Computer Simulation Using Particles (New York, McGraw-Hill, 1981).
  • [7] J.W. Connor and J.B. Taylor, Nucl. Fusion 17, 1047 (1977).
  • [8] T. Takizuka and H. Abe, J. Comput. Phys. 25, 205 (1977).