Tokamak elongation - how much is too much? Part 2. Numerical results

The analytic theory presented in Paper I is converted into a form convenient for numerical analysis. A fast and accurate code has been written using this numerical formulation. The results are presented by first defining a reference set of physical parameters based on experimental data from high performance discharges. Scaling relations of maximum achievable elongation ( ) versus inverse aspect ratio ( ) are obtained numerically for various values of poloidal beta ( ), wall radius and feedback capability parameter ( ) in ranges near the reference values. It is also shown that each value of occurs at a corresponding value of optimized triangularity ( ), whose scaling is also determined as a function of . The results show that the theoretical predictions of are slightly higher than experimental observations for high performance discharges, as measured by high average pressure. The theoretical values are noticeably lower. We suggest that the explanation is associated with the observation that high performance involves not only MHD stability, but also MHD modes described by in the Troyon limit and transport as characterized by . Operation away from the MHD optimum may still lead to higher performance if there are more than compensatory gains in and . Unfortunately, while the empirical scaling of and with the elongation ( ) has been determined, the dependence on has still not been quantified. This information is needed in order to perform more accurate overall optimizations in future experimental designs.