Abstract
Compact toroidal configuration is of simpler construction than the conventional tokamak and has important advantages due to the novel physics properties of low aspect ratio. In this paper we are developing a numerical program to study the magnetic dynamo or relaxation of CT's characterized by arbitrary tight aspect ratio. It is shown that the numerical method (Collocation Method), used here, works quite well to calculate numerically the lowest zero flux eigenvalues mu of Taylor's relaxed plasma state equation (del) over right arrow x (B) over right arrow = mu(B) over right arrow for an axisymmetric tokamaks of circular cross section. An excellent fulfillment of the toroidal flux vanishing boundary condition RRB integral integral B(phi)drdz = 0 along the whole boundary for such tokamaks are achieved. Dependence of mu on the aspect ratio is also obtained. Several runs of the program for various wave numbers k showed that mu is very insensitive to the choice of k. Besides, the poloidal magnetic field topologies inside the tokamak are well represented.