Stewart J. Zweben

The objectives of NSTX-U research are to reinforce the advantages of STs while addressing the challenges. To extend confinement physics of low-A, high beta plasmas to lower collisionality levels, understanding of the transport mechanisms that set confinement performance and pedestal profiles is being advanced through gyrokinetic simulations, reduced model development, and comparison to NSTX experiment, as well as improved simulation of RF heating. To develop stable non-inductive scenarios needed for steady-state operation, various performance-limiting modes of instability were studied, including MHD, tearing modes, and energetic particle instabilities. Predictive tools were developed, covering disruptions, runaway electrons, equilibrium reconstruction, and control tools. To develop power and particle handling techniques to optimize plasma exhaust in high performance scenarios, innovative lithium-based solutions are being developed to handle the very high heat flux levels that the increased heating power and compact geometry of NSTX-U will produce, and will be seen in future STs. Predictive capabilities accounting for plasma phenomena, like edge harmonic oscillations, ELMs, and blobs, are being tested and improved. In these ways, NSTX-U researchers are advancing the physics understanding of ST plasmas to maximize the benefit that will be gained from further NSTX-U experiments and to increase confidence in projections to future devices.

Magnetic field aligned filaments such as blobs and edge localized mode filaments carry significant amounts of heat and particles to the plasma facing components and they decrease their lifetime. The dynamics of these filaments determine at least a part of the heat and particle loads. These dynamics can be characterized by their translation and rotation. In this paper, we present an analysis method novel for fusion plasmas, which can estimate the angular velocity of the filaments on frame-by-frame time resolution. After pre-processing, the frames are two-dimensional (2D) Fourier-transformed, then the resulting 2D Fourier magnitude spectra are transformed to log-polar coordinates, and finally the 2D cross-correlation coefficient function (CCCF) is calculated between the consecutive frames. The displacement of the CCCF’s peak along the angular coordinate estimates the angle of rotation of the most intense structure in the frame. The proposed angular velocity estimation method is tested and validated for its accuracy and robustness by applying it to rotating Gaussian-structures. The method is also applied to gas-puff imaging measurements of filaments in National Spherical Torus Experiment plasmas.

Edge localized modes (ELMs) are a threat to tokamaks due to their high heat and particle loads on plasma facing components. A significant portion of this energy is carried and deposited by the emerging ELM filaments, whose dynamics are directly connected to their impact. Therefore, understanding their underlying physics is important for the operation of future fusion reactors. Our paper extends our knowledge of ELM filaments by reporting on their internal rotation (spinning) around the magnetic field lines along which they are extended. Our analysis of gas-puff imaging data on National Spherical Torus Experiment shows that ELM filaments are characterized by internal rotation in the direction of the ion-gyromotion with ω=15.2 krad/s median angular velocity, which is approximately three times faster than the blob rotation in the background turbulence. The characteristic size of the ELM filament was also assessed and found to be similar to the blobs. A nearly linear trend was found between the angular velocity and the radial velocity of the ELM filament. The angular velocity was found to be linearly dependent on the distance of the filament from the separatrix, as well. An analytical model called the shear-induced rotation model was identified as a candidate for explaining the physics of the observations. Our results show that the modeled mechanism could significantly influence the rotation of the ELM filament; however, it cannot be a sole contributor.

The time dependence of the blob pulse shape and the waiting time between blobs was found using data from the gas puff imaging diagnostic in NSTX. The database used was of 103 shots from 2010 as described in a previous paper (Zweben et al., Phys. Plasmas 29, 012505 (2022).]. The blob pulse shape was well fit by an exponential rise and fall where the average rise time was τr = 9.0 ± 2.7 μs and the average fall time was τf = 16.6 ± 5.8 μs. The waiting times between blob pulses above a threshold of three times the mean had a broad distribution with an average of τw = 1.2 ± 0.85 ms over the database. The blob intermittency parameter γb = τd/τw, where the blob pulse duration was τd = τr + τf, ranged from γb ∼1% to 5% for shots in this database and increased almost linearly with the blob fraction.

Previously, edge-localized mode (ELM)-free and inter-ELM divertor peak heat flux reduction induced by an edge harmonic oscillation (EHO) was observed in NSTX [Gan et al., Nucl. Fusion 57, 126053 (2017)]. This paper introduces new analysis of the EHO impact on the divertor heat flux. It was found that enhanced edge turbulence significantly altered the divertor heat flux footprint in NSTX H-modes, relative to ELM-free discharges. When the background edge fluctuation level was low, the EHO significantly increased the heat flux width and decreased the divertor peak heat flux. When the background edge fluctuation level was high, the EHO actually increased the divertor peak heat flux. Furthermore, it was also found that the heat flux width increased with the increasing frequency of the EHO.

This paper presents a simple analysis of the fraction of blobs observed in the scrape-off layer of NSTX and calculates the correlation between this blob fraction and various plasma parameters. The measurements were made using the gas puff imaging (GPI) diagnostic on a database of 103 shots from 2010 with neutral beam power P_NBI = 0–6 MW. The blob fraction ranged from f_b = 0.1% to 4.8%, where f_b is defined as the fraction of time the GPI signal was larger than three times its average value within 0–7.5 cm outside the separatrix. The blob fraction generally decreased with increasing neutral beam injection power and was lower on average in H-mode plasmas than Ohmic and L-mode plasmas. The blob fraction had its highest correlation with the measured poloidal turbulence velocity and its radial gradient just inside the separatrix. This is in part consistent with a model in which the blob generation rate increases with the poloidal velocity shear.

Anomalous plasma transport in the boundary region of a tokamak plasma is commonly associated with the formation and evolution of coherent density structures known as blobs. Recently, a theory for a universal mechanism of plasma blob formation has been put forward. It is based on a breaking process of a radially elongated streamer due to poloidal and radial velocity shears. The theory is well supported by two-dimensional and three-dimensional numerical simulation results but lacks experimental validation. In this work, we report the first ever experimental validation of this universal criterion by testing it against NSTX data on blobs obtained using the gas-puff imaging diagnostic. It is found that the criterion is widely satisfied in most L-mode discharges and may explain the significantly larger number of blob events. We also validate the theoretical criterion against ADITYA Langmuir probe data taken in the scrape-off layer region.

Dynamics of fast transient events are challenging to be analyzed with high time resolution. Such events can occur in fusion plasmas such as the filaments during edge-localized modes (ELMs). In this paper, we present a robust method—the spatial displacement estimation—for estimating the displacements of structures with fast dynamics from high spatial and time resolution imaging diagnostics [e.g., gas-puff imaging (GPI)] with sampling time temporal resolution. First, a background suppression method is shown, which suppresses the slowly time-evolving and spatially non-uniform background in the signal. In the second step, a two-dimensional polynomial trend subtraction method is presented to tackle the remaining polynomial order trend in the signal. After performing these pre-processing steps, the spatial displacement of the propagating structure is estimated from the two-dimensional spatial cross-correlation coefficient function calculated between consecutive frames. The method is tested for its robustness and accuracy by simulated Gaussian events and spatially displaced random noise. An example application of the method is presented on propagating ELM filaments measured by the GPI system on the National Spherical Torus Experiment spherical tokamak.