Abstract
The steady-state behavior of recycling neutral atoms in a tokamak edge region has been analyzed through a two-dimensional Monte Carlo simulation. A particle tracking algorithm used in earlier research on the neutral particle transport is applied to this Monte Carlo simulation in order to perform more accurate calculations with the EDGETRAN code which was previously developed for a two-dimensional edge plasma transport in the authors' laboratory. The physical model of neutral recycling includes charge-exchange and ionization interactions between plasmas and neutral atoms. The reflection processes of incident particles on the device wall are described by empirical formulas. Calculations for density, energy, and velocity distributions of neutral deuterium-tritium atoms have been carried out for a medium-sized tokamak with a double-null configuration based on the KT-2 conceptual design. The input plasma parameters such as plasma density, ion and electron temperatures, and ion fluid velocity are provided from the EDGETRAN calculations. As a result of the present numerical analysis, it is noticed that a significant drop of the neutral atom density appears in the region of high plasma density and that the similar distribution of neutral energy to that of plasma ions is present as frequently reported in other studies. Relations between edge plasma conditions and the neutral recycling behavior are discussed from the numerical results obtained herein.