Search

A decade-long operation of the Korean Superconducting Tokamak Advanced Research (KSTAR) has contributed significantly to the operation of superconducting tokamak devices and the advancement of tokamak physics which will be beneficial for the ITER and K-DEMO programs. Even with limited heating capability, various conventional as well as new operating regimes have been explored and have achieved improved performance. As examples, a long pulse high-confinement mode operation with and without an edge-localized mode (ELM) crash was well over 70 and 30 s, respectively. The unique capabilities of KSTAR allowed it to improve the capability of controlling harmful instabilities, and they have been instrumental in uncovering much new physics. The highlights are that the L/H transition threshold power is sensitive to the resonant magnetic perturbation (RMP) and insensitive to non-resonant magnetic perturbation. Co-Ip offset rotation dominated by an electron channel predicted by general neoclassical toroidal viscosity theory was confirmed. Improved heat dispersal in a divertor system using three rows of rotating RMP was demonstrated and predictive control of the ELM-crash with a priori modeling was successfully tested. In magnetohydrodynamic physics, validation of the full reconnection model (i.e. q0  >  1 right after the sawtooth crash) and self-consistent validation of the anisotropic distribution of turbulence amplitude and flow in the presence of the 2/1 island with theoretical models were achieved. The turbulence amplitude induced by RMP was linearly increased with the slow RMP coil current ramp-up time (i.e. the magnetic diffusion time scale). The Dα spikes (i.e. ELM-crash amplitude) was linearly decreased with the turbulence amplitude and not correlated with the perpendicular electron flow. In the turbulence area, a non-diffusive ‘avalanche’ transport event and the role of a quiescent coherent mode in confinement were studied. To accommodate the anticipation of a higher performance of the KSTAR plasmas with the increased heating powers, a new divertor/internal interface with a full active cooling system will be implemented after a full test of the new heating (neutral beam injection II and electron cyclotron heating) and current drive (CD) (Helicon and lower hybrid CD) systems. An upgrade plan for the internal hardware, heating systems and efficient CD system may allow for a long pulse operation of higher performance plasmas at βN  >  3.0 with f bs ~ 0.5 and Ti  >  10 keV. [ABSTRACT FROM AUTHOR]

Minor and major disruptions by explosive MHD instabilities were observed with the novel quasi 3D electron cyclotron emission imaging (ECEI) system in the KSTAR plasma. The fine electron temperature (Te) fluctuation images revealed two types of minor disruptions: a small minor disruption is a localized fast transport event due to a single m/n  =  2/1 magnetic island growth, while a large minor disruption is partial collapse of the region with two successive fast heat transport events by the correlated m/n  =  2/1 and m/n  =  1/1 instabilities. The m/n  =  2/1 magnetic island growth during the minor disruption is normally limited below the saturation width. However, as the additional interchange-like perturbation grows near the inner separatrix of the 2/1 island, the 2/1 island can expand beyond the limit through coupling with the cold bubble formed by the interchange-like perturbation. [ABSTRACT FROM AUTHOR]

Fast sawtooth crash of m  =  1 mode (tens of µs) is a common MHD phenomenon in the core of tokamak plasma. In KSTAR plasmas with localized heating and current drive close to the sawtooth inversion radius (ρinv), very slow sawtooth crashes are observed with time scales of a few ms. Scan of heating deposition along with power modulation revealed two different mechanisms of the slow crash with the help of detailed 2D imaging diagnostics. The slow crash appears with slow heat transport of the remained m  =  1 mode after the partial crash of outer multiple m/n  =  1 modes, or with the growth of an island formed on ρinv when the heating deposition locates just inside ρinv or just outside ρinv. In the heating power modulation experiment, the slow crash with the island was sustained for the redistribution time scale of the driven current by a very short blip of heating. [ABSTRACT FROM AUTHOR]

Temporal and spatial modulation of the edge localized mode (ELM) structure has been observed during the inter-ELM-crash period by toroidally-separated two electron cyclotron emission imaging systems. The observed modulation is interpreted as a beat wave of two modes with adjacent toroidal mode number. An additional assumption is that each mode has to have a different poloidal rotation speed. In nonlinear simulation, the low-n mode can be driven by locking between the dominant modes. The modulation is reconstructed using beat waves not the locking of modes. [ABSTRACT FROM AUTHOR]

Long-pulse plasma operation at high normalized beta, βN, above the n = 1 ideal MHD no-wall stability limit in KSTAR is presently limited by tearing instabilities rather than resistive wall modes. H-mode plasma operation during the recent KSTAR device campaign produced discharges having strong m/n = 2/1 tearing instabilities at βN lower than the ideal MHD no-wall beta limit. The unstable tearing mode consequently reduced plasma confinement and toroidal plasma rotation significantly. The experiment confirmed that an extended duration of electron cyclotron heating (ECH) at the initial phase of the discharge plays a critical role in mode destabilization. To study destabilizing mechanisms that affect the mode growth, the stability of the observed tearing modes from plasmas with significantly different βN is computed by using the resistive DCON code and the M3D-C1 code employing different physics. The computed tearing stability index, Δ′, differs between the mode that is destabilized by the early ECH at lower βN, and the mode that is destabilized at higher βN with observed mode triggering activity. Equilibrium reconstructions that include constraints from internal profile diagnostics are used as input for reliable computation of stability. The modified Rutherford equation (MRE) describing the evolution of the neoclassical tearing mode (NTM) island width has been constructed for KSTAR plasmas by using plasma parameters computed by the TRANSP code. In preparation for long-pulse plasma operation at higher beta utilizing increased plasma heating power, a resistive wall mode (RWM) active feedback control algorithm that includes magnetic sensor compensation of the prompt applied field and the field from the induced current on the passive conductors has been completed and enabled on KSTAR. Use of multiple toroidal sensor arrays is enabled for increased control performance by including the effect of varied mode helicities in the outboard region where the mode measurement is made. This analysis on beta-limiting instabilities and active mode control provides the required foundation for high confinement plasma operation on KSTAR without disruption. [ABSTRACT FROM AUTHOR]

The bifurcation of perpendicular mode velocity near the pedestal top () at the onset of RMP-driven edge-localized mode (ELM)-crash suppression has been directly measured by using an electron cyclotron emission imaging diagnostic on KSTAR. The ELM crashes are suppressed along with a sudden reduction of , which synchronizes with the transition into and out of the ELM-crash suppression. The change of flow appears mainly responsible for the change of , whose magnitude is the smallest near the normalized flux surface during the ELM-crash suppression. The plasma response to the RMP field was most enhanced in the vicinity of during the ELM-crash suppression and showed a hysteresis behavior with respect to resonant field strength. With these changes, a strong nonlinear coupling between turbulent eddies was observed in the ELM-crash suppression, which is an important feature of the ELM-crash suppression. [ABSTRACT FROM AUTHOR]

Using a high-speed and broadband radio frequency (RF) (0.1–1 GHz) spectrum analyzer developed on the KSTAR tokamak, it is found that several distinct stages of RF emission appear at the pedestal collapse in high confinement discharges. Comparison with 2D electron cyclotron emission (ECE) images has revealed that each stage is related to the instantaneous condition at the outboard mid-plane edge. First, high-harmonic ion cyclotron emissions (ICE) are intensified with the appearance of a non-modal filamentary perturbation in the edge within several tens of microseconds before the collapse. Then, the RF emission becomes broad toward high-frequency range (<500 MHz) at the burst onset of the non-modal filament. During the pedestal collapse initiated by the filament burst, rapid chirping (1–3 μs) appear with additional filament bursts. The strong correlation between the RF spectra and the perturbation structure provides important clues on the stability of edge-localized modes and on the ion dynamics in the plasma boundary. [ABSTRACT FROM AUTHOR]

The central safety factor (q0) during sawtooth oscillation has been measured with a great accuracy with the motional Stark effect (MSE) system on KSTAR and the measured value was However, this measurement alone cannot validate the disputed full and partial reconnection models definitively due to non-trivial off-set error (~0.05). Supplemental experiment of the excited m  =  2, m  =  3 modes that are extremely sensitive to the background q0 and core magnetic shear definitively validates the ‘full reconnection model’. The radial position of the excited modes right after the crash and time evolution into the 1/1 kink mode before the crash in a sawtoothing plasma suggests that in the MHD quiescent period after the crash and before the crash. Additional measurement of the long lived m  =  3, m  =  5 modes in a non-sawtoothing discharge (presumably ) further validates the ‘full reconnection model’. [ABSTRACT FROM AUTHOR]

Broadband density fluctuations with peak frequency ranging from 150 to 400 kHz were measured using a multichannel microwave imaging reflectometer in core region of the low-density L-mode plasmas heated by neutral beam injection on KSTAR. These fluctuations have been studied by comparing the dominant mode scales estimated from the measurement with those predicted from linear gyrokinetic simulation. The measured poloidal wavenumbers are qualitatively comparable to those of the ‘fastest growing modes’ from simulations, whereas they are larger than those of the ‘transport-dominant modes’ by about a factor of three. The agreement on wavenumbers between the measurement and linear simulation (for the fastest growing modes) is probably due to sufficiently weak flow shear compared to the maximum linear growth rate. Meanwhile, the transport-dominant modes seem to be related to the fluctuations in lower frequencies (∼80–150 kHz) observed in some of the measurement. [ABSTRACT FROM AUTHOR]

Multiscale interaction between the magnetic island and turbulence has been demonstrated through simultaneous two-dimensional measurements of turbulence and temperature and flow profiles. The magnetic island and turbulence can mutually interact via coupling between the electron temperature (Te) gradient, the Te turbulence, and the poloidal flow. The Te gradient altered by the magnetic island steepens outside and flattens inside the island. The Te turbulence can appear in increased Te gradient regions. The combined effects of the Te gradient and the poloidal flow shear determines the two-dimensional distribution of the Te turbulence. When the poloidal vortex flow forms, it can maintain the steepest Te gradient and the magnetic island acts more like an electron heat transport barrier. Interestingly, when the Te gradient, the Te turbulence, and the vortex flow shear increase beyond critical levels, the magnetic island turns into a fast electron heat transport channel, which directly leads to the minor disruption. [ABSTRACT FROM AUTHOR]

In this contribution we analyze modulation electron cyclotron resonance heating (MECH) experiment and discuss higher harmonic frequency dependence of transport coefficients. We use the bidirectional heat pulse propagation method, in which both inward propagating heat pulse and outward propagating heat pulse are analyzed at a radial range, in order to distinguish frequency dependence of transport coefficients due to hysteresis from that due to other reasons, such as radially dependent transport coefficients, a finite damping term, or boundary effects. The method is applied to MECH experiments performed in various helical and tokamak devices, i.e. Large Helical Device (LHD), TJ-II, Korea Superconducting Tokamak Advanced Research (KSTAR), and Doublet III-D (DIII-D) with different plasma conditions. The frequency dependence of transport coefficients are clearly observed, showing a possibility of existence of transport hysteresis in flux-gradient relation. [ABSTRACT FROM AUTHOR]

Multiple (two or more) flux tubes are commonly observed inside and/or near the q = 1 flux surface in KSTAR tokamak plasmas with localized electron cyclotron resonance heating and current drive (ECH/CD). Detailed 2D and quasi-3D images of the flux tubes obtained by an advanced imaging diagnostic system showed that the flux tubes are m/n = 1/1 field-aligned structures co-rotating around the magnetic axis. The flux tubes typically merge together and become like the internal kink mode of the usual sawtooth, which then collapses like a usual sawtooth crash. A systematic scan of ECH/CD beam position showed a strong correlation with the number of flux tubes. In the presence of multiple flux tubes close to the q = 1 surface, the radially outward heat transport was enhanced, which explains naturally temporal changes of electron temperature. We emphasize that the multiple flux tubes are a universal feature distinct from the internal kink instability and play a critical role in the control of sawteeth using ECH/CD. [ABSTRACT FROM AUTHOR]

A detailed study of edge-localized mode (ELM) dynamics in the KSTAR tokamak is performed using a two-dimensional (2D) electron cyclotron emission imaging (ECEI) diagnostic system. Highly coherent mode structures rotating in the poloidal view plane are routinely observed in the inter-ELM pedestal region where the optical thickness for ECE rapidly changes and the interpretation of emission intensity is complicated. To have confidence on the measurements, the observed images are compared with synthetic images of the ELM structure deduced from three-field BOUT++ simulations. The synthetic process considers instrumental effects of the ECEI diagnostic, intrinsic broadening of the ECE and background noise. The synthetic 2D images highly resemble the observed structure, providing confidence that the ELM dynamics can be visualized by ECEI. [ABSTRACT FROM AUTHOR]