An efficient method is given to reconstruct the current profile parameters, the plasma shape, and a current profile consistent with the magnetohydrodynamic equilibrium constraint from external magnetic measurements, based on a Picard iteration approach which approximately conserves the measurements. Computational efforts are reduced by parametrizing the current profile linearly in terms of a number of physical parameters. Results of detailed comparative calculations and a sensitivity study are described. Illustrative calculations to reconstruct the current profiles and plasma shapes in ohmically and auxiliarily heated Doublet III plasmas are given which show many interesting features of the current profiles.

The predictions of gyrokinetic and gyrofluid simulations of ion-temperature-gradient (ITG) instability and turbulence in tokamak plasmas as well as some tokamak plasma thermal transport models, which have been widely used for predicting the performance of the proposed ITER tokamak, are compared. These comparisons provide information on effects of differences in the physics content of the various models and on the fusion-relevant figures of merit of plasma performance predicted by the models. Many of the comparisons are undertaken for a simplified plasma model and geometry which is an idealization of the plasma conditions and geometry in a DIII-D H-mode experiment. Most of the models show good agreements in their predictions and assumptions for the linear growth rates and frequencies. There are some differences associated with different equilibria. However, there are significant differences in the transport levels between the models. The causes of some of the differences are examined in some detail, with particular attention to numerical convergence in the turbulence simulations (with respect to simulation mesh size, system size and, for particle-based simulations, the particle number). The implications for predictions of fusion plasma performance are also discussed.

A model based on magnetohydrodynamic (MHD) stability of the tokamak plasma edge region is presented, which describes characteristics of edge localized modes (ELMs) and the pedestal. The model emphasizes the dual role played by large bootstrap currents driven by the sharp pressure gradients in the pedestal region. Pedestal currents reduce the edge magnetic shear, stabilizing high toroidal mode number (n) ballooning modes, while at the same time providing drive for intermediate to low n peeling modes. The result is that coupled peeling–ballooning modes at intermediate n (3<n<20) are often the limiting instability which constrains the pedestal and triggers ELMs. These modes are characterized in shaped tokamak equilibria using an efficient new numerical code, and simplified models are developed for pedestal limits and the ELM cycle. Results are compared to several experiments, and nonideal MHD effects are briefly discussed.

The maximum normalized beta achieved in long-pulse tokamak discharges at low collisionality falls significantly below both that observed in short pulse discharges and that predicted by the ideal MHD theory. Recent long-pulse experiments, in particular those simulating the International Thermonuclear Experimental Reactor ͑ITER͒ ͓M. Rosenbluth et al., Plasma Physics and Controlled Nuclear Fusion ͑International Atomic Energy Agency, Vienna, 1995͒, Vol. 2, p. 517͔ scenarios with low collisionality e * , are often limited by low-m/n nonideal magnetohydrodynamic ͑MHD͒ modes. The effect of saturated MHD modes is a reduction of the confinement time by 10%-20%, depending on the island size and location, and can lead to a disruption. Recent theories on neoclassical destabilization of tearing modes, including the effects of a perturbed helical bootstrap current, are successful in explaining the qualitative behavior of the resistive modes and recent results are consistent with the size of the saturated islands. Also, a strong correlation is observed between the onset of these low-m/n modes with sawteeth, edge localized modes ͑ELM͒, or fishbone events, consistent with the seed island required by the theory. We will focus on a quantitative comparison between both the conventional resistive and neoclassical theories, and the experimental results of several machines, which have all observed these low-m/n nonideal modes. This enables us to single out the key issues in projecting the long-pulse beta limits of ITER-size tokamaks and also to discuss possible plasma control methods that can increase the soft ␤ limit, decrease the seed perturbations, and/or diminish the effects on confinement.

An efficient method is given for self-consistent reconstruction of the tokamak current profiles and their associated magnetic topology using the magnetohydrodynamic (MHD) equilibrium constraint from external magnetic measurements, kinetic profile measurements, internal poloidal magnetic field measurements, and topological information from soft X-ray (SXR) measurements. Illustrative examples for beam heated H-mode divertor discharges in the DIII-D tokamak are presented, using the experimentally measured kinetic profile information and external magnetic data from the existing diagnostics. Comparative reconstructions of the current profile using various combinations of diagnostics are given. Also presented is an alternative magnetic analysis method in which the MHD equilibrium is reconstructed using external magnetic data and a constraint on the edge pressure gradient. The results of a sensitivity study are given which show that the axial safety factor q(0) can be more accurately determined when additional information from internal poloidal magnetic measurements is used in conjunction with the external magnetic, kinetic and SXR topological data.

No abstract

No abstract

A practical method for performing a tokamak equilibrium reconstruction in real time for arbitrary time-varying discharge shapes and current profiles is described. An approximate solution to the Grad-Shafranov equilibrium relation is found which best fits the diagnostic measurements. Thus a solution for the spatial distribution of poloidal flux and toroidal current density is available in real time that is consistent with plasma force balance, allowing accurate evaluation of parameters such as discharge shape and safety factor profile. The equilibrium solutions are produced at a rate sufficient for discharge control. This equilibrium reconstruction algorithm has been implemented on the digital plasma control system for the DIII-D tokamak. The first application of a real time equilibrium reconstruction to discharge shape control is described.