Monte Carlo estimation of neoclassical transport for the TJ-II stellarator

Abstract

The neoclassical transport properties of TJ-II stellarator [C. Alejaldre et al., Fusion Technol. 13, 521 (1988)] are studied with the monoenergetic Monte Carlo technique. A compromise between the number of modes and particles and the required computing time to obtain reliable estimates, from the computational point of view, of the monoenergetic diffusion coefficients is shown to be of one thousand particles and one hundred harmonics, because of the rich magnetic-field structure of TJ-II. Although, these requirements are probably too demanding in making the transport estimations. The data base containing the normalized monoenergetic diffusion coefficient for several radial positions, radial electric fields and collisionalities have been fitted using a neural network. This fit reduces the number of points necessary in the data base and allows a smooth interpolation and extrapolation to perform the convolutions of the monoenergetic coefficients with the Maxwellian. For two different typical TJ-II discharges the ambipolar radial electric field, and the neoclassical particle and heat fluxes are presented both showing rather large positive radial electric fields at the plasma core and small negative fields at the edge. The neoclassical particle and energy confinement time are in surprisingly good agreement with the experimental energy balance analysis and the international stellarator scaling. Although no satisfactory explanation is available yet the large neoclassical diffusion caused by the complex ripple structure of TJ-II magnetic field may be an important ingredient.