National Spherical Torus Experiment Research Articles

The ion cyclotron parametric instabilities and the anomalous heating of ions in the scrape-off layer of tokamak plasma in the high harmonic fast wave heating regime

The electrostatic parametric instabilities of a plasma, driven by the high harmonic fast wave (HHFW) with frequency at ion cyclotron (IC) harmonics of order 30–50 are investigated numerically. The derived numerical results are consistent with parametric decay of HHFW into the high harmonic IC (Bernstein) wave and an unobserved IC quasimode under conditions of the negligible small electron Landau damping. The detected instability develops in the finite interval of the HHFW wavelength along the toroidal magnetic field. The development of this ion kinetic quasimode decay instability is accompanied by the anomalous anisotropic heating of ions, resulted from the interaction of ions with IC parametric turbulence. It was found that the anomalous heating rate of ions across the magnetic field is much larger than the heating rate of ions along the magnetic field. The anisotropic heating of the scrape-off layer ions was observed on the National Spherical Torus Experiment experiments with HHFW heating and current drive at IC harmonics of order 10 [G. Taylor et al. Phys. Plasmas 17, 056114 (2010), and J. Hosea et al. Phys. Plasmas 15, 056104 (2008)]. The derived results predict that these experimental results will be reproduced qualitatively in the experiments with HHFW heating and current drive using HHFW at IC harmonics of order 30–50.

Energy transport analysis of NSTX plasmas with the TGLF turbulent and NEO neoclassical transport models

This work presents a study of plasma transport at low aspect ratio on the National Spherical Torus Experiment tokamak, where the turbulent and neoclassical energy fluxes calculated by the quasilinear Trapped Gyro Landau Fluid (TGLF) model and the multi species drift-kinetic Neoclassical solver (NEO) are validated against experimental data. The turbulent energy transport of two plasma discharges, one in the L-mode confinement regime and another in the H-mode regime, is dominated by electrostatic drift-wave instabilities, while the ion heat transport has a significant neoclassical contribution. The data analysis workflow is described in detail to understand how the variations of mapping and fitting of experimental data affect the power balance solution and subsequent flux-matching plasma profile predictions with the TGYRO solver. On average, the predicted plasma profiles are consistent with experimental data. However, the solutions are sensitive to various input parameters, including boundary conditions, and the electron-ion coupling. Linear gyrokinetic stability analysis demonstrates close agreement of the real frequencies of unstable modes between TGLF and CGYRO gyrokinetic simulations, but higher growth rates are predicted by TGLF, especially for the H-mode case. Estimates of the low-k modes’ contributions to the total flux are consistent with linear stability analysis and the E × B suppression of turbulence in TGLF simulations with the SAT1 saturation model, while the SAT2 saturation model over-predicts the low-k modes’ contribution in the H-mode case. Moreover, the results with SAT1 model are consistent with power balance analysis, which indicates only neoclassical ion energy fluxes inside ρ < 0.4 in the L-mode case and in the H-mode case. The presence of multi-scale turbulence and ion-scale driven zonal flow mixing effects are also observed in TGLF scans of the electron turbulent heat flux over a range of temperature gradients and the electron-ion temperature ratio, which could explain the strong model sensitivity to variations of input parameters.

The role of fast ions in spontaneous neoclassical tearing mode instabilities in NSTX

Spontaneous neoclassical tearing modes (spontaneous NTMs) have been observed in the National Spherical Torus Experiment (NSTX). Recently, a TRANSP based model for the modified Rutherford equation analysis was developed to accurately predict the effect of fast ions on the growth of magnetic islands. The goal of this paper is to utilize the TRANSP NTM model to understand the role of fast ions in NSTX spontaneous NTMs. Our analysis shows that when the passing fast ion driven kinetic neoclassical polarization current term is larger than 1% of the bootstrap current term in the modified Rutherford equation, fast ions play a decisive role in the early growth phase of spontaneous NTM. This result is applicable to any tokamak plasmas with a fast ion drift orbit width, or the fast ion poloidal Larmor radius, comparable to the magnetic island width.

COMPASS-U Global Heat Balance Calculations

COMPASS-U is a medium size, high magnetic field experimental tokamak ( <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$R =0.9$ </tex-math></inline-formula> m, <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$B_{t} = 5$ </tex-math></inline-formula> T, and <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$I_{p} = 2$ </tex-math></inline-formula> MA), built at the Institute of Plasma Physics, Czech Academy of Sciences (IPP-CAS). This global heat balance calculation for COMPASS-U was done at Princeton Plasma Physics Laboratory (PPPL). Based on our previous experience of building global thermal model for National Spherical Torus Experiment-Upgrade (NSTX-U) at PPPL, this 2-D global thermal model geometry represents a typical cross section of COMPASS-U machine. The model includes thermal radiation, conduction, and convection among components and also between the machine and outer environment. Helium gas heating/cooling was modeled with fluid element and surface element. This model was used to calculate component temperatures and heat distribution during heat up, cryogenic cool down, normal operation, and fault operation scenarios. Definition of 14 main thermal scenarios was provided by IPP-CAS. First nine are normal operation scenarios. Scenarios #10–#14 are fault scenarios. Some results of thermal scenario #5 will be given and discussed in this article, including peak temperatures, temperature ratcheting, energy distribution, and required cooling power.

Novel angular velocity estimation technique for plasma filaments.

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.

Bus Design for the Poloidal Field Coils of the NSTX-Upgrade Fusion Device

The National Spherical Torus eXperiment (NSTX) has undergone a major upgrade to NSTX-U at Princeton Plasma Physics Laboratory (PPPL). NSTX upgrade (NSTX-U) will double the toroidal field, plasma current, and neutral beam injection heating power, as well as significantly increase the pulse duration. NSTX-U uses three poloidal field (PF) coils at the vessel top and three at the bottom near the divertor areas to control the local plasma shape there. These shaping coils operate at a maximum terminal voltage of 2 kV, corresponding to a maximum current about 20 kA. To supply the power to the PF coils, hard copper bus bars are typically used near the vacuum vessel, while water cooled flexible cables are used away from the vessel. The detailed design of the hard copper bus bars will be covered in this article. During operation, these hard bus bars are subject to high electromagnetic (EM) forces, thermal displacement loads, as well as plasma halo loads due to disruption. The EM, thermal, and structural analysis were performed, and the results revealed that, with the worst loads combined, the new design meets the NSTX-U thermal, structural, and fatigue cycle requirements. The manufacturing and installation process of the bus work will be discussed at the end of this article.

Design Considerations for the Implementation of a High-Field-Side Transient CHI System on QUEST

Transient coaxial helicity injection (T-CHI), a method first developed on the small helicity injected torus (HIT-II) experiment and then validated on the much larger National Spherical Torus Experiment (NSTX) device, is a method to initiate an inductive-like tokamak plasma discharge without reliance on the central solenoid. A CHI discharge is initiated by driving current along magnetic flux that connects the inner and outer divertor plates on one end of the tokamak. To permit this, on both HIT-II and NSTX, toroidal ceramic insulators were used to electrically separate the inner and outer vessel components. The use of such large toroidal vacuum insulators may not be easy to implement in reactors. To address this issue, the Q-shu University experiments with steady-state spherical tokamak (ST) (QUEST) is developing a reactor-relevant CHI configuration in which one of the divertor plates is electrically insulated from the rest of the vessel. The first application of T-CHI on QUEST biased the CHI electrode to the outer vessel. While the CHI discharges could be easily generated, it was found that as the discharge filled the vessel, the separation distance between the injector magnetic flux footprints widened, a condition that is not favorable for the generation of closed flux surfaces. Biasing the electrode to the inner wall is a configuration similar to that used on NSTX and HIT-II, but initial testing in this configuration has proved to be challenging. The design described here overcomes the present limitation by locating the CHI electrode much closer to the CHI injector flux coil and using an NSTX-like gas injection manifold to enable high-field-side T-CHI startup on QUEST. The concepts described in this article should also benefit the future implementation of T-CHI systems in other tokamaks and spherical tokamaks.

NSTX-U Machine Core Vacuum Seals Upgrade Design

The vacuum in the National Spherical Torus Experiment Upgrade (NSTX-U) is in the range of <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$2.66\,\times \,10 ^{-6}$ </tex-math></inline-formula> Pa ( <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$2\,\times \,10 ^{-8}$ </tex-math></inline-formula> torr). To enhance vacuum seals and lower gas permeation in the NSTX-U and therefore ensuring operation safety and quality, the machine core vacuum seal design has been upgraded using double elastomeric seals and pumped interspaces (guard channels). In NSTX-U machine core, there are six major vacuum seals, including upper and lower vacuum vessel main flanges, upper and lower bellows flanges, and two in the ceramic isolation ring assembly that are mounted to the upper vacuum vessel main flange. The elastomeric seals are commercially available and customized fluorocarbon (Viton) O-rings. To ensure appropriate O-ring compression ratios for the vacuum seal, the O-ring grooves are customized with appropriate tolerances considering specific situations for each seal. In this upgrade, the vacuum in each interspace is designed to reach 13.33 Pa (0.1 torr), thus the overall pressure difference crossing each double O-ring seal is expected to be reduced by at least three orders of magnitude, and accordingly the overall rates of leakage and permeation will be lowered by three orders of magnitude. The corresponding pump-down times for each seal, at room temperature, are estimated at about 60 s. In addition, the vacuum pumped interspaces will be able to function as a real-time leak monitor.

Global gyrokinetic study of shaping effects on electromagnetic modes at NSTX aspect ratio with ad hoc parallel magnetic perturbation effects

Plasma shaping may have a stronger effect on global turbulence in tight-aspect-ratio tokamaks than in conventional-aspect-ratio tokamaks due to the higher toroidicity and more acute poloidal asymmetry in the magnetic field. In addition, previous local gyrokinetic studies have shown that it is necessary to include parallel magnetic field perturbations in order to accurately compute growth rates of electromagnetic modes in tight-aspect-ratio tokamaks. In this work, the effects of elongation and triangularity on global, ion-scale, linear electromagnetic modes are studied at National Spherical Torus Experiment (NSTX) aspect ratio and high plasma β using the global gyrokinetic particle-in-cell code XGC. The effects of compressional magnetic perturbations are approximated via a well-known modification to the particle drifts that was developed for flux-tube simulations [Joiner et al., Phys. Plasmas 17, 072104 (2010)], without proof of its validity in a global simulation, with the gyrokinetic codes GENE and GEM being used for local verification and global cross-verification. Magnetic equilibria are re-constructed for each distinct plasma profile that is used. Coulomb collision effects are not considered. Within the limitations imposed by the present study, it is found that linear growth rates of electromagnetic modes (collisionless microtearing modes and kinetic ballooning modes) are significantly reduced in a high-elongation and high-triangularity NSTX-like geometry compared to a circular NSTX-like geometry. For example, growth rates of kinetic ballooning modes at high-β are reduced to the level of that of collisionless trapped electron modes.

The National Spherical Torus Experiment-Upgrade poloidal high-k scattering system pitch angle design modifications.

A 693GHz, eight-channel, poloidal high-k (k refers to wavenumber) collective scattering system is under development for the National Spherical Torus Experiment-Upgrade device. It will replace the previous 280GHz, five-channel, tangential scattering system to study high-k electron density fluctuations, thereby providing a measurement of the kθ-spectrum of both electron temperature gradient and ion temperature gradient modes. A tool is under development to calculate the wavenumber that exists in the presence of strong magnetic pitch angles. We use this tool to motivate a new receiver optical design for significantly improved performance, details of which are presented herein.

Internal rotation of ELM filaments on NSTX

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.

Transient CHI System Design Studies for Pegasus-III

Transient coaxial helicity injection (transient CHI), first developed on the Helicity Injected Torus-II (HIT-II) and later on the National Spherical Torus Experiment (NSTX) for implementing solenoid-free plasma current startup capability in a spherical tokamak (ST), is now planned to be tested on the Pegasus-III ST using a novel double-biased configuration. Such a configuration is likely needed for transient CHI deployment in a reactor. The transient CHI system optimization will be studied on Pegasus-III to enable startup toroidal persisting currents at the limits permitted by the external poloidal field coils. A transient CHI discharge is generated by driving injector current along magnetic field lines that connect the inner and outer divertor plates on one end of the ST. Simulations using the Tokamak Simulation Code are used to assess the transient CHI toroidal current generation potential and electrode gap location on the Pegasus-III. While past transient CHI systems have used high-voltage, oil-filled capacitors for driving the injector current, for improved safety, Pegasus-III will use a high-current capacitor bank based on low-voltage electrolytic capacitors. The designed and fabricated system is capable of over 32 kA. The modular design features permit the system to be upgraded to higher currents, as needed, to meet the future needs of the Pegasus-III facility.

Application of simulations and theory of sub-cyclotron frequency modes to DIII-D

Theoretical and numerical study of Alfvén eigenmodes (AEs) in the sub-cyclotron frequency range explains the observed frequency scaling with beam parameters in experiments across different devices. Global Alfvén (GAEs) and compressional Alfvén eigenmodes (CAEs) are frequently excited during neutral beam injection (NBI) in the National Spherical Torus Experiment (NSTX/NSTX-U), as well as other beam-heated devices such as MAST and DIII-D. These modes are driven unstable through the Doppler shifted cyclotron resonance with the NBI ions and can be excited in ITER due to super-Alfvénic velocities and anisotropy of the fast ion distribution. Numerical models and theory for sub-cyclotron frequency modes, previously developed and used to study the excitation of GAEs/CAEs in the NSTX(-U), have been successfully applied to explain the DIII-D observations. Simulations for DIII-D demonstrate that the modes with ω/ω ci ∼ 0.6, previously mis-identified as CAEs, have shear polarization and should be identified as the GAEs. Simulation results match the observed frequencies and estimated toroidal mode numbers in DIII-D experiments. Good agreement is found between scaling of the experimentally observed GAE frequencies with NBI parameters for NSTX, NSTX-U and DIII-D and scaling predicted by the theory and simulations. The DIII-D simulation results also demonstrate that GAEs can be excited for larger aspect-ratio devices and smaller (sub-Alfvénic) beam injection velocities, therefore they can be unstable in ITER.

AC compensation of 3D magnetic diagnostic signals in DIII-D and National Spherical Torus Experiment-Upgrade (NSTX-U) for real-time application.

A time domain algorithm has been developed to remove the vacuum pickup generated by both coil current (DC) and induced vessel current (AC) in real time from three dimensional (3D) magnetic diagnostic signals in the National Spherical Torus Experiment-Upgrade (NSTX-U) and DIII-D tokamaks. The possibility of detecting 3D plasma perturbations in real time is essential in modern and future tokamaks to avoid and control MHD instabilities. The presence of vacuum field pickup, due to toroidally asymmetric (3D) coils or to misalignment between sensors and axisymmetric (2D) coils, pollutes the measured plasma 3D field, making the detection of the magnetic field produced by the plasma challenging. Although the DC coupling between coils and sensors can be easily calculated and removed, the AC part is more difficult. An algorithm based on a layered low-pass filter approach for the AC compensation and its application for DIII-D and NSTX-U data is presented, showing that this method reduces the vacuum pickup to the noise level. Comparison of plasma response measurements with and without vacuum compensation shows that accurate mode locking detection and plasma response identification require precise AC and DC compensations.

HTS Cable Conductor for Compact Fusion Tokamak Solenoids

Significant progress has been made recently in the US fusion community to develop a strategic plan to enable engineering design and construction of a fusion pilot plant (FPP). Princeton Plasma Physics Laboratory (PPPL) is working on developing high current density HTS conductors for next fusion experiments. Partnering with the US industry, we are evaluating feasibility and affordability of Conductor on Round Core (CORC <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"><tex-math notation="LaTeX">${}^\circledR$</tex-math></inline-formula> ) developed by Advanced Conductor Technologies (ACT) for the next compact fusion tokamak facility. High current density achieved by a CORC <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"><tex-math notation="LaTeX">${}^\circledR$</tex-math></inline-formula> cable based on a four-layer model coil recently tested at National High Magnetic Field Laboratory (NHMFL) motivated its consideration for low cost, reduced size fusion magnet application. This is of interest to PPPL because of its scalability to tokamak central solenoid (CS) coils in terms of required flux swings for plasma startup operations. Partnering with ACT, we designed and built a two-layer model coil solenoid directly wound with CORC to demonstrate its applicability for compact solenoids such as that used in the national spherical torus experiment (NSTX), NSTX-upgrade (NSTX-U) and the US sustained high-power density test facility (SHPD). The ∼160 mm diameter solenoid wound by a two-layer CORC is being tested in early 2022 under electromagnetic cyclic loading at NHMFL in a unique 160 mm bore, 14 T background field magnet facility.

On the origin of high harmonic fast wave edge losses in NSTX

Many mechanisms have been proposed to explain the high harmonic fast wave edge losses in the National Spherical Torus Experiment (NSTX) (Perkins et al 2012 Phys. Rev. Lett. 109 045001). Here, we discuss a newly discovered edge loss mechanism (Tierens et al 2021 Plasma Phys. Control. Fusion 64 035001) which involves the resonant excitation of surface waves on filaments. We show that this edge loss mechanism can explain the main parameter dependencies observed in NSTX experiments.

Neural net modeling of equilibria in NSTX-U

Neural networks (NNs) offer a path towards synthesizing and interpreting data on faster timescales than traditional physics-informed computational models. In this work we develop two NNs relevant to equilibrium and shape control modeling, which are part of a suite of tools being developed for the National Spherical Torus Experiment-Upgrade for fast prediction, optimization, and visualization of plasma scenarios. The networks include Eqnet, a free-boundary equilibrium solver trained on the EFIT01 (Equilibrium FITtting 01) reconstruction algorithm, and Pertnet, which is trained on the Gspert code and predicts the non-rigid plasma response, a nonlinear term that arises in shape control modeling. The NNs are trained with different combinations of inputs and outputs in order to offer flexibility in use cases. In particular, Eqnet can use magnetic diagnostics as inputs and act as an EFIT-like reconstruction algorithm, or, by using pressure and current profile information the NN can act as a forward Grad–Shafranov equilibrium solver. This forward-mode version is envisioned to be implemented in the suite of tools for simulation of plasma scenarios. The reconstruction-mode version gives some performance improvements compared to the online reconstruction code real-time EFIT, especially when vessel eddy currents are significant. We report strong performance for all NNs indicating that the models could reliably be used within closed-loop simulations or other applications. Some limitations are discussed.

Ideal MHD Limited Electron Temperature in Spherical Tokamaks.

It is well documented that the central electron temperature in the national spherical torus experiment (NSTX) remains largely unchanged as the external heating power, and hence the normalized volume averaged plasma pressure β increases [Stutman, Phys. Rev. Lett. 102, 115002 (2009)PRLTAO0031-900710.1103/PhysRevLett.102.115002]. Here we present a hypothesis that low n, pressure driven ideal magnetohydrodynamic (MHD) instabilities that are nondisruptive, can break magnetic surfaces in the central region and thereby flatten the electron temperature profiles. We demonstrate this mechanism in a 3D resistive MHD simulation of a NSTX discharge. By varying the toroidal magnetic field strength, and/or the heating power, we show that there is a critical value of β, above which the central temperature profile no longer peaks on axis.

Kinetic modeling of neutral transport for a continuum gyrokinetic code

We present the first-of-its-kind coupling of a continuum full-f gyrokinetic turbulence model with a 6D continuum model for kinetic neutrals, carried out using the Gkeyll code. Our objective is to improve the first-principle understanding of the role of neutrals in plasma fueling, detachment, and their interaction with edge plasma profiles and turbulence statistics. Our model includes only atomic hydrogen and incorporates electron-impact ionization, charge exchange, and wall recycling. These features have been successfully verified with analytical predictions and benchmarked with the DEGAS2 Monte Carlo neutral code. We carry out simulations for a scrape-off layer (SOL) with simplified geometry and National Spherical Torus Experiment parameters. We compare these results to a baseline simulation without neutrals and find that neutral interactions reduce the normalized density fluctuation levels and associated skewness and kurtosis, while increasing auto-correlation times. A flatter density profile is also observed, similar to the SOL density shoulder formation in experimental scenarios with high fueling.

Turbulent broadening of electron heat-flux width in electromagnetic gyrokinetic simulations of a helical scrape-off layer model

We demonstrate that cross field transport in the scrape-off layer (SOL) can be moderately increased by electromagnetic effects in high-beta regimes, resulting in broadening of the electron heat-flux width on the endplates. This conclusion is taken from full-f electromagnetic gyrokinetic simulations of a helical SOL model that roughly approximates the SOL of the National Spherical Torus Experiment. The simulations have been performed with the Gkeyll code, which recently became the first code to demonstrate the capability to simulate electromagnetic gyrokinetic turbulence on open magnetic field lines with sheath boundary conditions. We scan the source rate and thus β, so that the normalized pressure gradient (the MHD ballooning parameter α∝∂β/∂r∝β/Lp) is scanned over an experimentally relevant range, α=0.3−1.5. While there is little change in the pressure gradient scale length Lp near the midplane as beta is increased, a 10% increase in cross field transport near the midplane results in an increase in the electron heat-flux width λq and a 25% reduction of the peak electron heat flux to the endplates.