Numerical calculation of axisymmetric non-neutral plasma equilibria

Abstract

Efficient techniques for computing axisymmetric non-neutral plasma equilibria are described. These equilibria may be obtained either by requiring global thermal equilibrium, by specifying the midplane radial density profile, or by specifying the radial profile of ∫n dz. Both splines and finite-differences are used, and the accuracy of the two is compared by using a new characterization of the thermal equilibrium density profile which gives a simple formula for estimating the radial and axial gradient scale lengths of thermal equilibria. It is found that for global thermal equilibrium 1% accuracy is achieved with splines if the distance between neighboring splines is about two Debye lengths while finite differences require a grid spacing of about one-half Debye length to achieve the same accuracy.