D. Shiraki

The high beta tokamak-extended pulse magnetohydrodynamic mode control research program

The high beta tokamak-extended pulse (HBT-EP) magnetohydrodynamic (MHD) mode control research program is studying ITER relevant internal modular feedback control coil configurations and their impact on kink mode rigidity, advanced digital control algorithms and the effects of plasma rotation and three-dimensional magnetic fields on MHD mode stability. A new segmented adjustable conducting wall has been installed on the HBT-EP and is made up of 20 independent, movable, wall shell segments instrumented with three distinct sets of 40 saddle coils, totaling 120 in-vessel modular feedback control coils. Each internal coil set has been designed with varying toroidal angular coil coverage of 5, 10 and 15°, spanning the toroidal angle range of an ITER port plug based internal coil to test resistive wall mode (RWM) interaction and multimode MHD plasma response to such highly localized control fields. In addition, we have implemented 336 new poloidal and radial magnetic sensors to quantify the applied three-dimensional fields of our control coils along with the observed plasma response. This paper describes the design and implementation of the new control shell incorporating these control and sensor coils on the HBT-EP, and the research program plan on the upgraded HBT-EP to understand how best to optimize the use of modular feedback coils to control instability growth near the ideal wall stabilization limit, answer critical questions about the role of plasma rotation in active control of the RWM and the ferritic resistive wall mode, and to improve the performance of MHD control systems used in fusion experiments and future burning plasma systems.

Error field detection in DIII-D by magnetic steering of locked modes

Optimal correction coil currents for the n?=?1 intrinsic error field of the DIII-D tokamak are inferred by applying a rotating external magnetic perturbation to steer the phase of a saturated locked mode with poloidal/toroidal mode number m/n?=?2/1. The error field is detected non-disruptively in a single discharge, based on the toroidal torque balance of the resonant surface, which is assumed to be dominated by the balance of resonant electromagnetic torques. This is equivalent to the island being locked at all times to the resonant 2/1 component of the total of the applied and intrinsic error fields, such that the deviation of the locked mode phase from the applied field phase depends on the existing error field. The optimal set of correction coil currents is determined to be those currents which best cancels the torque from the error field, based on fitting of the torque balance model. The toroidal electromagnetic torques are calculated from experimental data using a simplified approach incorporating realistic DIII-D geometry, and including the effect of the plasma response on island torque balance based on the ideal plasma response to external fields. This method of error field detection is demonstrated in DIII-D discharges, and the results are compared with those based on the onset of low-density locked modes in ohmic plasmas. This magnetic steering technique presents an efficient approach to error field detection and is a promising method for ITER, particularly during initial operation when the lack of auxiliary heating systems makes established techniques based on rotation or plasma amplification unsuitable.

A Kalman filter for feedback control of rotating external kink instabilities in the presence of noise

The simulation and experimental optimization of a Kalman filter feedback control algorithm for n=1 tokamak external kink modes are reported. In order to achieve the highest plasma pressure limits in ITER, resistive wall mode stabilization is required [T. C. Hender et al., Nucl. Fusion 47, S128 (2007)] and feedback algorithms will need to distinguish the mode from noise due to other magnetohydrodynamic activity. The Kalman filter contains an internal model that captures the dynamics of a rotating, growing n=1 mode. This model is actively compared with real-time measurements to produce an optimal estimate for the mode’s amplitude and phase. On the High Beta Tokamak-Extended Pulse experiment [T. H. Ivers et al., Phys. Plasmas 3, 1926 (1996)], the Kalman filter algorithm is implemented using a set of digital, field-programmable gate array controllers with 10 μs latencies. Signals from an array of 20 poloidal sensor coils are used to measure the n=1 mode, and the feedback control is applied using 40 poloidally...

Feedback suppression of rotating external kink instabilities in the presence of noise

The authors report on the first experimental demonstration of active feedback suppression of rotating external kink modes near the ideal wall limit in a tokamak using Kalman filtering to discriminate the n=1 kink mode from background noise. The Kalman filter contains an internal model that captures the dynamics of a rotating, growing n=1 mode. Suppression of the external kink mode is demonstrated over a broad range of phase angles between the sensed mode and applied control field, and performance is robust at noise levels that render proportional gain feedback ineffective. Suppression of the kink mode is accomplished without excitation of higher frequencies as was observed in previous experiments using lead-lag loop compensation [A. J. Klein et al., Phys Plasmas 12, 040703 (2005)].

Spatiotemporal Evolution of Runaway Electron Momentum Distributions in Tokamaks.

Novel spatial, temporal, and energetically resolved measurements of bremsstrahlung hard-x-ray (HXR) emission from runaway electron (RE) populations in tokamaks reveal nonmonotonic RE distribution functions whose properties depend on the interplay of electric field acceleration with collisional and synchrotron damping. Measurements are consistent with theoretical predictions of momentum-space attractors that accumulate runaway electrons. RE distribution functions are measured to shift to a higher energy when the synchrotron force is reduced by decreasing the toroidal magnetic field strength. Increasing the collisional damping by increasing the electron density (at a fixed magnetic and electric field) reduces the energy of the nonmonotonic feature and reduces the HXR growth rate at all energies. Higher-energy HXR growth rates extrapolate to zero at the expected threshold electric field for RE sustainment, while low-energy REs are anomalously lost. The compilation of HXR emission from different sight lines into the plasma yields energy and pitch-angle-resolved RE distributions and demonstrates increasing pitch-angle and radial gradients with energy.

Measurements of the toroidal torque balance of error field penetration locked modes

Detailed measurements from the DIII-D tokamak of the toroidal dynamics of error field penetration locked modes under the influence of slowly evolving external fields, enable study of the toroidal torques on the mode, including interaction with the intrinsic error field. The error field in these low density Ohmic discharges is well known based on the mode penetration threshold, allowing resonant and non-resonant torque effects to be distinguished. These m/n = 2/1 locked modes are found to be well described by a toroidal torque balance between the resonant interaction with n = 1 error fields, and a viscous torque in the electron diamagnetic drift direction which is observed to scale as the square of the perturbed field due to the island. Fitting to this empirical torque balance allows a time-resolved measurement of the intrinsic error field of the device, providing evidence for a time-dependent error field in DIII-D due to ramping of the Ohmic coil current.

Measurement of 3D plasma response to external magnetic perturbations in the presence of a rotating external kink

The detailed measurements of the 3D plasma response to applied external magnetic perturbations in the presence of a rotating external kink are presented, and compared with the predictions of a single-helicity linear model of kink mode dynamics. The modular control coils of the High Beta Tokamak-Extended Pulse (HBT-EP) device are used to apply resonant m/n = 3/1 magnetic perturbations to wall-stabilized tokamak plasmas with a pre-existing rotating 3/1 kink mode. The plasma response is measured in high-resolution with the extensive magnetic diagnostic set of the HBT-EP device. The spatial structures of both the naturally rotating kink mode and the externally driven response are independently measured and observed to be identical, while the temporal dynamics are consistent with the independent evolution and superposition of the two modes. This leads to the observation of a characteristic change in 3D field dynamics as a function of the applied field amplitude. This amplitude dependence is found to be different for poloidal and radial fields. The measured 3D response is compared to and shown to be consistent with the predictions of the linear single-helicity model in the “high-dissipation” regime, as reported previously [M. E. Mauel et al., Nucl. Fusion 45, 285 (2005)].

High resolution detection and excitation of resonant magnetic perturbations in a wall-stabilized tokamaka)

We report high-resolution detection of the 3D plasma magnetic response of wall-stabilized tokamak discharges in the High Beta Tokamak-Extended Pulse [T. H. Ivers et al., Phys. Plasmas 3, 1926 (1996)] device. A new adjustable conducting wall has been installed on HBT-EP made up of 20 independent, movable, wall segments instrumented with three distinct sets of 40 modular coils that can be independently driven to generate a wide variety of magnetic perturbations. High-resolution detection of the plasma response is made with 216 poloidal and radial magnetic sensors that have been located and calibrated with high-accuracy. Static and dynamic plasma responses to resonant and non-resonant magnetic perturbations are observed through measurement of the step-response following a rapid change in the toroidal phase of the applied perturbations. Biorthogonal decomposition of the full set of magnetic sensors clearly defines the structures of naturally occurring external kinks as being composed of independent m/n = 3/1 and 6/2 modes. Resonant magnetic perturbations were applied to discharges with pre-existing, saturated m/n = 3/1 external kink mode activity. This m/n = 3/1 kink mode was observed to lock to the applied perturbation field. During this kink mode locked period, the plasma resonant response is characterized by a linear, a saturated, and a disruptive plasma regime dependent on the magnitude of the applied field and value of the edge safety factor and plasma rotation.We report high-resolution detection of the 3D plasma magnetic response of wall-stabilized tokamak discharges in the High Beta Tokamak-Extended Pulse [T. H. Ivers et al., Phys. Plasmas 3, 1926 (1996)] device. A new adjustable conducting wall has been installed on HBT-EP made up of 20 independent, movable, wall segments instrumented with three distinct sets of 40 modular coils that can be independently driven to generate a wide variety of magnetic perturbations. High-resolution detection of the plasma response is made with 216 poloidal and radial magnetic sensors that have been located and calibrated with high-accuracy. Static and dynamic plasma responses to resonant and non-resonant magnetic perturbations are observed through measurement of the step-response following a rapid change in the toroidal phase of the applied perturbations. Biorthogonal decomposition of the full set of magnetic sensors clearly defines the structures of naturally occurring external kinks as being composed of independent m/n = 3/1 ...

Array magnetics modal analysis for the DIII-D tokamak based on localized time-series modelling

Time-series analysis of magnetics data in tokamaks is typically done using block-based fast Fourier transform methods. This work presents the development and deployment of a new set of algorithms for magnetic probe array analysis. The method is based on an estimation technique known as stochastic subspace identification (SSI). Compared with the standard coherence approach or the direct singular value decomposition approach, the new technique exhibits several beneficial properties. For example, the SSI method does not require that frequencies are orthogonal with respect to the timeframe used in the analysis. Frequencies are obtained directly as parameters of localized time-series models. The parameters are extracted by solving small-scale eigenvalue problems. Applications include maximum-likelihood regularized eigenmode pattern estimation, detection of neoclassical tearing modes, including locked mode precursors, and automatic clustering of modes, and magnetics-pattern characterization of sawtooth pre- and postcursors, edge harmonic oscillations and fishbones.

Resolving runaway electron distributions in space, time, and energy

Areas of agreement and disagreement with present-day models of runaway electron (RE) evolution are revealed by measuring MeV-level bremsstrahlung radiation from runaway electrons (REs) with a pinhole camera. Spatially resolved measurements localize the RE beam, reveal energy-dependent RE transport, and can be used to perform full two-dimensional (energy and pitch-angle) inversions of the RE phase-space distribution. Energy-resolved measurements find qualitative agreement with modeling on the role of collisional and synchrotron damping in modifying the RE distribution shape. Measurements are consistent with predictions of phase-space attractors that accumulate REs, with non-monotonic features observed in the distribution. Temporally resolved measurements find qualitative agreement with modeling on the impact of collisional and synchrotron damping in varying the RE growth and decay rate. Anomalous RE loss is observed and found to be largest at low energy. Possible roles for kinetic instability or spatial tr...