Plasma and Fusion Research
Volume 13, 1203088 (2018)
Rapid Communications
- Graduate School of Engineering, Osaka University, Suita 565-0871, Japan
- 1)
- National Institutes for Quantum and Radiological Science and Technology, Rokkasho 039-3212, Japan
Abstract
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 ε∗.
Keywords
simulation model, global plasma, electron dynamics, Poisson equation
Full Text
References
- [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), http://www.riam.kyushu-u.ac.jp/sosei/bpsi/2017/BPSI2017_proc.pdf
- [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).