ANTIC: A code for calculation of neutral transport in cylindrical plasmas

A scheme for computing the steady state transport of neutral atoms cylindrical plasmas is described. The physical model used represents atoms emerging from charge exchange collisions by an isotropic source of neutrals with energy equal to 3/2 times the local ion temperature. The process of charge exchange and impact ionization by electrons and ions are included. The transport can be described by an integral equation for the neutral source density, and this equation approximated by a set of algebraic equations for the zone average source densities. A derivation of the appropriate kernel and a technique for computing it are presented. Comparison of computed results with predictions of other, more exact, codes exhibit quite satisfactory agreement for temperatures as high as 10 keV.