J.R. Ferron

Equilibrium analysis of current profiles in tokamaks

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.

Edge localized modes and the pedestal: A model based on coupled peeling–ballooning modes

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.

Real time equilibrium reconstruction for tokamak discharge control

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 real time equilibrium reconstruction to discharge shape control is described.

Sustained rotational stabilization of DIII-D plasmas above the no-wall beta limit

Sustained stabilization of the n=1 kink mode by plasma rotation at beta approaching twice the stability limit calculated without a wall has been achieved in DIII-D by a combination of error field reduction and sufficient rotation drive. Previous experiments have transiently exceeded the no-wall beta limit. However, demonstration of sustained rotational stabilization has remained elusive because the rotation has been found to decay whenever the plasma is wall stabilized. Recent theory [Boozer, Phys. Rev. Lett. 86, 5059 (2001)] predicts a resonant response to error fields in a plasma approaching marginal stability to a low-n kink mode. Enhancement of magnetic nonaxisymmetry in the plasma leads to strong damping of the toroidal rotation, precisely in the high-beta regime where it is needed for stabilization. This resonant response, or “error field amplification” is demonstrated in DIII-D experiments: applied n=1 radial fields cause enhanced plasma response and strong rotation damping at beta above the no wal...

ELMs and constraints on the H-mode pedestal: peeling-ballooning stability calculation and comparison with experiment

We review and test the peeling–ballooning model for edge localized modes (ELMs) and pedestal constraints, a model based upon theoretical analysis of magnetohydrodynamic (MHD) instabilities that can limit the pedestal height and drive ELMs. A highly efficient MHD stability code, ELITE, is used to calculate quantitative stability constraints on the pedestal, including constraints on the pedestal height. Because of the impact of collisionality on the bootstrap current, these pedestal constraints are dependent on the density and temperature separately, rather than simply on the pressure. ELITE stability calculations are directly compared with experimental data for a series of plasmas in which the density is varied and ELM characteristics change. In addition, a technique is developed whereby peeling–ballooning pedestal constraints are calculated as a function of key equilibrium parameters via ELITE calculations using series of model equilibria. This technique is used to successfully compare the expected pedestal height as a function of density, triangularity and plasma current with experimental data. Furthermore, the technique can be applied for parameter ranges beyond the purview of present experiments, and we present a brief projection of peeling–ballooning pedestal constraints for burning plasma tokamak designs.

Wall stabilization of high beta plasmas in DIII-D

Detailed analysis of recent high beta discharges in the DIII-D tokamak demonstrates that the resistive vacuum vessel can provide stabilization of low n magnetohydrodynamic (MHD) modes. The experimental beta values reaching up to {beta}{sub T} = 12.6% are more than 30% larger than the maximum stable beta calculated with no wall stabilization. Plasma rotation is essential for stabilization. When the plasma rotation slows sufficiently, unstable modes with the characteristics of the predicted {open_quotes}resistive wall{close_quotes} mode are observed. Through slowing of the plasma rotation between the q = 2 and q = 3 surfaces with the application of a non-axisymmetric field, the authors have determined that the rotation at the outer rational surfaces is most important, and that the critical rotation frequency is of the order of {Omega}/2{pi} = 1 kHz.

MHD Equilibrium Reconstruction in the DIII-D Tokamak

Abstract Physics elements and advances crucial for the development of axisymmetric magnetohydrodynamic equilibrium reconstruction to support plasma operation and data analysis in the DIII-D tokamak are reviewed. A response function formalism and a Picard linearization scheme are used to efficiently combine the equilibrium and the fitting iterations and search for the optimum solution vector. Algorithms to incorporate internal current and pressure profile measurements, topological constraints, and toroidal plasma rotation into the equilibrium reconstruction are described. Choice of basis functions and boundary conditions essential for accurate reconstruction of L- and H-mode equilibrium plasma boundary and current and pressure profiles is discussed. The computational structure used to efficiently integrate these elements into the equilibrium reconstruction code EFIT is summarized.

Design and simulation of extremum-seeking open-loop optimal control of current profile in the DIII-D tokamak

In a magnetic fusion reactor, the achievement of a certain type of plasma current profiles, which are compatible with magnetohydrodynamic stability at high plasma pressure, is key to enable high fusion gain and non-inductive sustainment of the plasma current for steady-state operation. The approach taken toward establishing such plasma current profiles at the DIII-D tokamak is to create the desired profile during the plasma current ramp-up and early flattop phases. The evolution in time of the current profile is related to the evolution of the poloidal flux, which is modeled in normalized cylindrical coordinates using a partial differential equation usually referred to as the magnetic diffusion equation. The control problem is formulated as an open-loop, finite-time, optimal control problem for a nonlinear distributed parameter system, and is approached using extremum seeking. Simulation results, which demonstrate the accuracy of the considered model and the efficiency of the proposed controller, are presented.

Progress toward fully noninductive, high beta conditions in DIII-D

The DIII-D Advanced Tokamak (AT) program in the DIII-D tokamak [J. L. Luxon, Plasma Physics and Controlled Fusion Research, 1986, Vol. I (International Atomic Energy Agency, Vienna, 1987), p. 159] is aimed at developing a scientific basis for steady-state, high-performance operation in future devices. This requires simultaneously achieving 100% noninductive operation with high self-driven bootstrap current fraction and toroidal beta. Recent progress in this area includes demonstration of 100% noninductive conditions with toroidal beta, βT=3.6%, normalized beta, βN=3.5, and confinement factor, H89=2.4 with the plasma current driven completely by bootstrap, neutral beam current drive, and electron cyclotron current drive (ECCD). The equilibrium reconstructions indicate that the noninductive current profile is well aligned, with little inductively driven current remaining anywhere in the plasma. The current balance calculation improved with beam ion redistribution that was supported by recent fast ion diagno...

Modification of high mode pedestal instabilities in the DIII-D tokamak

The amplitude and frequency of modes driven in the edge region of tokamak high mode (H-mode) discharges [type I edge-localized modes (ELMs)] are shown to depend on the discharge shape. The measured pressure gradient threshold for instability and its scaling with discharge shape are compared with predictions from ideal magnetohydrodynamic theory for low toroidal mode number (n) instabilities driven by pressure gradient and current density and good agreement is found. Reductions in mode amplitude are observed in discharge shapes with either high squareness or low triangularity where the stability threshold in the edge pressure gradient is predicted to be reduced and the most unstable mode is expected to have higher values of n. The importance of access to the ballooning mode second stability regime is demonstrated through the changes in the ELM character that occur when second regime access is not available. An edge stability model is presented that predicts that there is a threshold value of n for second r...