Real-time Equilibrium Reconstruction with Kinetic Constraints on DIII-D
Many discharge relevant parameters can not be measured directly. However, they can be evaluated indirectly from available data from real-time diagnostics. A least squares solution of the diagnostic data to the Grad-Shafranov model allows this evaluation (rtEFIT). The Grad-Shafranov model describes the force balance of the tokamak equilibrium, and allows for a distributed current source. rtEFIT determines the distribution of poloidal flux and toroidal current density in the R,Z-plane.
If kinetic constraints are absent, the pressure profiles and current profiles resulting from equilibrium reconstruction do not exhibit the features that one would typically expect (no pedestal in pressure, no bootstrap peak in current). Our work focuses on using real-time kinetic profile data from diagnostics like Thomson Scattering (TS) and Charge Exchange Recombination Spectroscopy (CER) to generate pressure profile fits (and current profile fits) which are sent to rtEFIT, to be used to constrain the solution of equilibrium reconstruction. See the Figure below for an example of the qualitative improvements to the features of the reconstructed current and pressure profile when the kinetically constrained version is compared against equilibrium reconstruction using only magnetics:
Advanced ECE and ECCD Applications
Electron Cyclotron Emission (ECE) is a crucial diagnostic for measuring the electron temperature in tokamaks as a function of time and radius. Electron Cyclotron Current Drive (ECCD) is a technique used to inject power into the plasma and can be used to suppress dangerous instabilities. We have long been involved in the development of these system on DIII-D, which has some of the most advanced diagnostic systems in the world. For example, we developed a set of real-time steerable mirrors for the DIII-D system and designed algorithms to reduce the power requirements and time needed to suppress 3/2 and 2/1 neoclassical tearing modes (NTMs). This allows operators to achieve more stable plasmas and avoid disruptions.
Recently, we have been looking into ways to extend the information gathered by ECE to improve equilibrium reconstructions on future machines. We have developed a new algorithm that can generate a fully-constrained kinetic equilibrium using only a ECE diagnostic and an image of the plasma boundary. This will be useful in future machines with intense neutron heat fluxes, since both ECE and plasma imaging can be robustly designed in harsh neutron environments, unlike many other diagnostics commonly employed today. This is made possible by measurement of several different rational q surfaces simultaneously using only a single diagnostic, which allows us to transform the traditional constraints of pressure and current into constrains on electron temperature and q, as depicted here:
Pedestal Transport Studies
The H-mode pedestal is an extremely important and complicated region of tokamak plasmas. It lies near the plasma edge and is characterized by steep pressure gradients, which are responsible not only for improved plasma performance but also for dangerous transient instabilities called ELMs. In order to better understand these phenomena we are focusing on two separate physics influences: that of neutral ionization and that of micro-instability transport through the pedestal.
To study the effect of neutrals, we are comparing the pedestal structure on similar discharges fueled by pellet injection and gas puffing, which deposit neutrals closer to the core and edge, respectively. We have found that the deposition of neutrals can play a large role in setting the structure of the H-mode pedestal and thus should be considered in predictive models of the pedestal height and width.
In a separate research thrust, we are using novel experimental tools to perturb the edge current in the pedestal region. This allows us to decouple the instability drive and resonant location for micro-tearing modes (MTMs), which are likely responsible for a large portion of the electron heat flux through the pedestal region. This investigation has so far led to the first time-dependent experimental identification of MTMs in an H-mode pedestal, paving the way for future improvements to predictive pedestal models. An example experimental identification of a m/n=16/3 MTM in the H-mode edge is shown here:
Automated equilibria and SOL modeling on DIII-D
As the understanding of plasma physics progresses, more and more detailed modeling is required to solve the outstanding problems in the field. Unfortunately, much of this modeling requires time-consuming manual oversight, which severely limits our abilities to conduct database studies of large numbers of discharges. On DIII-D, we have developed two automated codes to help with this issue. First, the CAKE code is now operational and regularly provides high-quality automatically-generated kinetic equilibria for use in physics studies. Second, we have recently developed the autoUEDGE algorithm, with self-consistently solves transport equations in the scrape-off-layer (SOL) in order to facilitate database studies of edge plasma transport and to provide improved initial conditions for expert simulations. An example of a converged autoUEDGE transport solution is shown here: