The ability of an injected lithium granule to promptly trigger an edge localized mode (ELM) has been established in multiple experiments. By horizontally injecting granules ranging in diameter from 200 microns to 1mm in diameter into the low field side of EAST H-mode discharges we have determined that granules with diameter > 600 microns are successful in triggering ELMs more than 95% of the time. It was also demonstrated that below 600 microns the triggering efficiency decreased roughly with granule size. Granules were radially injected from the outer midplane with velocities ~ 80 m/s into EAST upper single null discharges with an ITER like tungsten monoblock divertor. These granules were individually tracked throughout their injection cycle in order to determine their efficacy at triggering an ELM. For those granules of sufficient size, ELM triggering was a prompt response to granule injection. By simulating the granule injection with an experimentally benchmarked neutral gas shielding (NGS) model, the ablatant mass deposition required to promptly trigger an ELM is calculated and the fractional mass deposition is determined.
Stotler, D.P.; Battaglia, D.J.; Hager, R.; Kim, K.; Koskela, T.; Park, G.; Reinke, M.L.
Modifications of the drift-kinetic transport code XGC0 to include the
transport, ionization, and recombination of individual charge states,
as well as the associated radiation, are described. The code is first
applied to a simulation of an NSTX H-mode discharge with carbon
impurity to demonstrate the approach to coronal equilibrium. The
effects of neoclassical phenomena on the radiated power profile are
examined sequentially through the activation of individual physics
modules in the code. Orbit squeezing and the neoclassical inward
pinch result in increased radiation for temperatures above a few
hundred eV and changes to the ratios of charge state emissions at a
given electron temperature. Analogous simulations with a neon
impurity yield qualitatively similar results.
Non-axisymmetric control coils and the so-called snowflake divertor configuration are two potential solutions proposed to solve two separate outstanding issues on the path towards self-sustained burning plasma operations, namely the transient energy bursts caused by edge localized modes and the steady state heat exhaust problem. In a reactor, these two proposed solutions would have to operate simultaneously and it is, therefore, important to investigate their compatibility and to identify possible conflicts that could prevent them from operating simultaneously. In this work, single- and two-fluid resistive magnetohydrodynamic calculations are used to investigate the effect of externally applied magnetic perturbations on the snowflake divertor configuration. The calculations are based on simulated NSTX-U plasmas and the results show that additional and longer magnetic lobes are created in the null-point region of the snowflake configuration, compared to those in the conventional single-null. The intersection of these longer and additional lobes with the divertor plates are expected to cause more striations in the particle and heat flux target profiles. In addition, the results indicate that the size of the magnetic lobes, in both single-null and snowflake configurations, are more sensitive to resonant magnetic perturbations than to non-resonant magnetic perturbations. The results also suggest that lower values of current in non-axisymmetric control coils would be required to suppress edge localized modes in plasmas with the snowflake configuration.
In this paper we present initial simulations of pedestal control by Lithium Granule Injection (LGI) in NSTX. A model for small granule ablation has been implemented in the M3D-C1 code , allowing the simulation of realistic Lithium granule injections. 2D simulations in NSTX L-mode and H-mode plasmas are done and the effect of granule size, injection angle and velocity on the pedestal gradient increase are studied. For H-mode cases, the amplitude of the local pressure perturbation caused by the granules is highly dependent on the solid granule size. In our simulations, reducing the granule injection velocity allows one to inject more particles at the pedestal top.
Lunsford, R.; Bortolon, A.; Roquemore, A.L.; Mansfield, D.K.; Jaworski, M.A.; Kaita, R.; Maingi, R.; Nagy, A.
By employing a neutral gas shielding (NGS) model to characterize impurity granule
injection the pedestal atomic deposition for three different species of granule:
lithium, boron, and carbon are determined. Utilizing the duration of ablation
events recorded on experiments performed at DIII-D to calibrate the NGS model we
are able to quantify the ablation rate and mass deposition location with respect
to the plasma density profile. The species specific granule shielding constant
is then used to model granule ablation within NSTX-U discharges. Simulations of
300, 500 and 700 micron diameter granules injected at 50 m/sec are presented for
NSTX-U L-mode type plasmas as well as H-mode discharges with low natural ELM
frequencies. Additionally, ablation calculations of 500 micron granules of each
species are presented at velocities ranging from 50 � 150 m/sec. In H-mode type
discharges these simulations show that the majority of the injected granule is
ablated within or just past the steep gradient region of the discharge. At this
radial position, the perturbation to the background plasma generated by the ablating
granule can lead to conditions advantageous for the rapid triggering of an ELM crash
Stotler, D.P.; Lang, J.; Chang, C.S.; Churchill, R.M.; Ku, S.-H.
The effects of recycled neutral atoms on tokamak ion temperature
gradient (ITG) driven turbulence have been investigated in a steep
edge pedestal, magnetic separatrix configuration, with the full-f
edge gryokinetic code XGC1. Ion temperature gradient turbulence is
the most fundamental and robust edge plasma instability, having a long
radial correlation length and an ability to impact other forms of
pedestal turbulence. The neutral atoms enhance the ITG turbulence,
first, by increasing the ion temperature gradient in the pedestal via
the cooling effects of charge exchange and, second, by a relative
reduction in the ExB shearing rate.
Results of 3D nonlinear simulations of neutral-beam-driven compressional Alfven eigenmodes (CAEs) in the National Spherical Torus Experiment (NSTX) are presented. Hybrid MHD-particle simulations for the H-mode NSTX discharge (shot 141398) using the HYM code show unstable CAE modes for a range of toroidal mode numbers, n=4-9, and frequencies below the ion cyclotron frequency. It is found that the essential feature of CAEs is their coupling to kinetic Alfven wave (KAW) that occurs on the high-field side at the Alfven resonance location. High-frequency Alfven eigenmodes are frequently observed in beam-heated NSTX plasmas, and have been linked to flattening of the electron temperature profiles at high beam power. Coupling between CAE and KAW suggests an energy channeling mechanism to explain these observations, in which beam-driven CAEs dissipate their energy at the resonance location,
therefore significantly modifying the energy deposition profile. Nonlinear simulations demonstrate that CAEs can channel the energy of the beam ions from the injection region near the magnetic axis to the location of the resonant mode conversion at the edge of the beam density profile. A set of nonlinear simulations show that the CAE instability saturates due to nonlinear particle trapping, and a large fraction of beam energy can be transferred to several unstable CAEs of relatively large amplitudes and absorbed at the resonant location. Absorption rate shows a strong scaling with the beam power.
A scintillator type fast ion loss detector measures the gyroradius and pitch angle distribution of superthermal ions escaping from a magnetically confined fusion plasma at a single location. Described here is a technique for optimizing the angular orientation of such a detector in an axisymmetric tokamak geometry in order to intercept losses over a useful and interesting ranges of pitch angle. The method consists of evaluating the detector acceptance as a function of the fast ion constants of motion, i.e. energy, canonical toroidal momentum, and magnetic moment. The detector acceptance can then be plotted in a plane of constant energy and compared with the relevant orbit class boundaries and fast ion source distributions. Knowledge of expected or interesting mechanisms of loss can further guide selection of the detector orientation. The example of a fast ion loss detector for the National Spherical Torus Experiment-Upgrade (NSTX-U) is considered.
The National Spherical Torus Experiment (NSTX) has undergone a major upgrade, and the NSTX Upgrade (NSTX-U) Project was completed in the summer of 2015. NSTX-U first plasma was subsequently achieved, diagnostic and control systems have been commissioned, H-Mode accessed, magnetic error fields identified and mitigated, and the first physics research campaign carried out. During 10 run weeks of operation, NSTX-U surpassed NSTX-record pulse-durations and toroidal fields, and high-performance ~1MA H-mode plasmas comparable to the best of NSTX have been sustained near and slightly above the n=1 no-wall stability limit and with H-mode confinement multiplier H98y2 above 1. Transport and turbulence studies in L-mode plasmas have identified the coexistence of at least two ion-gyro-scale turbulent micro-instabilities near the same radial location but propagating in opposite (i.e. ion and electron diamagnetic) directions. These modes have the characteristics of ion-temperature gradient and micro-tearing modes, respectively, and the role of these modes in contributing to thermal transport is under active investigation. The new second more tangential neutral beam injection was observed to significantly modify the stability of two types of Alfven Eigenmodes. Improvements in offline disruption forecasting were made in the areas of identification of rotating MHD modes and other macroscopic instabilities using the Disruption Event Characterization and Forecasting (DECAF) code. Lastly, the Materials Analysis and Particle Probe (MAPP) was utilized on NSTX-U for the first time and enabled assessments of the correlation between boronized wall conditions and plasma performance. These and other highlights from the first run campaign of NSTX-U are described.
The WallDYN package has recently been applied to a number of tokamaks to self-consistently model
the evolution of mixed-material plasma facing surfaces. A key component of the WallDYN model is the
concentration-dependent surface sputtering rate, calculated using SDTRIM.SP. This modeled sputtering
rate is strongly influenced by the surface binding energies (SBEs) of the constituent materials, which
are well known for pure elements but often are poorly constrained for mixed-materials. This work examines
the sensitivity of WallDYN surface evolution calculations to different models for mixed-material
SBEs, focusing on the carbon/lithium/oxygen/deuterium system present in NSTX. A realistic plasma background is reconstructed from a high density, H-mode NSTX discharge, featuring an attached outer strike
point with local density and temperature of 4e20 m^-3 and 4 eV, respectively. It is found that various
mixed-material SBE models lead to significant qualitative and quantitative changes in the surface evolution
profile at the outer divertor, with the highest leverage parameter being the C-Li binding model.
Uncertainties of order 50%, appearing on time scales relevant to tokamak experiments, highlight the importance of choosing an appropriate mixed-material sputtering representation when modeling the surface
evolution of plasma facing components. These results are generalized to other fusion-relevant materials
with different ranges of SBEs.
Linear stability analysis of the national spherical torus experiment (NSTX) Li-conditioned
ELM-free H-mode equilibria is carried out in the context of the extended
magneto-hydrodynamic (MHD) model in NIMROD. The purpose is to investigate the physical
cause behind edge localized mode (ELM) suppression in experiment after the Li-coating of
the divertor and the first wall of the NSTX tokamak. Besides ideal MHD modeling, including
finite-Larmor radius effect and two-fluid Hall and electron diamagnetic drift contributions,
a non-ideal resistivity model is employed, taking into account the increase of Z eff after
Li-conditioning in ELM-free H-mode. Unlike an earlier conclusion from an eigenvalue code
analysis of these equilibria, NIMROD results find that after reduced recycling from divertor
plates, profile modification is necessary but insufficient to explain the mechanism behind
complete ELMs suppression in ideal two-fluid MHD. After considering the higher plasma
resistivity due to higher Z eff , the complete stabilization could be explained. A thorough
analysis of both pre-lithium ELMy and with-lithium ELM-free cases using ideal and
non-ideal MHD models is presented, after accurately including a vacuum-like cold halo
region in NIMROD to investigate ELMs.
In this paper we present data from experiments on NSTX-U where it is shown for the first time that small amounts of high pitch-angle beam ions can strongly suppress the counter-propagating Global Alfvén Eigenmodes (GAE). GAE have been implicated in the redistribution of fast ions and modification of the electron power balance in previous experiments on NSTX. The ability to predict the stability of Alfvén modes, and developing methods to control them, is important for fusion reactor like the International Tokamak Experimental Reactor (ITER) which are heated by a large population of non-thermal, super-Alfvénic ions consisting of fusion generated alphas and beam ions injected for current profile control. We present a qualitative interpretation of these observations using an analytic model of the Doppler-shifted ion-cyclotron resonance drive responsible for GAE instability which has an important dependence on k⊥ρL. A quantitative analysis of this data with the HYM stability code predicts both the frequencies and instability of the GAE prior to, and suppression of the GAE after the injection of high pitch-angle beam ions.
Radio-frequency (RF) rectification is an important sheath phenomenon for wave heating of plasma in fusion devices and is proposed to be the mechanism responsible for converting highharmonic fast-wave (HHFW) power in the National Spherical Torus eXperiment (NSTX) into a heat ux to the divertor. RF rectification has two aspects: current rectification and voltage recti- fication, and, while the latter is emphasized in many application, we demonstrate the importance of current rectification in analysis of the NSTX divertor during HHFW heating. When rectified currents are accounted for in first-principle models for the heat ux to the tiles, we predict a sizeable enhancement for the heat ux in the presence of an RF field: for one case studied, the predicted heat ux increases from 0:103 MW=m2 to 0:209 MW=m2. We also demonstrate how this rectification scales with injected HHFW power by tracking probe characteristics during a HHFW power ramp; the rectified current may be clamped at a certain level. This work is important for minimizing SOL losses of HHFW power in NSTX-U but may also have implications for near-field studies of ICRF antennae: ignoring rectified current may lead to underestimated heat uxes and overestimated rectified voltages.
The 2-D radial vs. poloidal cross-correlation functions of edge plasma turbulence were measured near the outer midplane using the gas puff imaging (GPI) diagnostic on NSTX. These correlation functions were evaluated at radii r= 0 cm, ±3 cm, and ±6 cm from the separatrix and poloidal locations p=0 cm and ±7.5 cm from the GPI poloidal center line for 20 different shots. The ellipticity ε and tilt angle φ of the positive cross- correlation regions, and the minimum negative cross-correlation “cmin” and total negative over positive values “neg/pos” were evaluated for each of these cases. The average results over this data set were ε=2.2±0.9, φ=87±34o (i.e. poloidally oriented), cmin= -0.30±0.15, and neg/pos=0.25±0.24. Thus there was significant variation in these correlation results within this database, with dependences on the location within the image, the magnetic geometry, and the plasma parameters. Possible causes for this variation are discussed, including the misalignment of the GPI view with the local B field line, the magnetic shear of field lines in the edge, the poloidal flow shear of the turbulence, blob-hole correlations, and the neutral density ‘shadowing’ effect in GPI.
A compact and multi-view Solid State Neutral Particle Analyzer (SSNPA) diagnostic based on silicon photodiode arrays has been successfully tested on the National Spherical Torus Experiment-Upgrade (NSTX-U). The SSNPA diagnostic provides spatially, temporally, and pitch-angle resolved measurements of fast-ion distribution by detecting fast neutral flux resulting from charge exchange (CX) reactions. The system consists of three 16-channel subsystems: t-SSNPA viewing the plasma mid-radius and neutral beam (NB) line #2 tangentially, r-SSNPA viewing the plasma core and NB line #1 radially and p-SSNPA with no intersection with any NB lines. Due to the setup geometry, the active CX signals of t-SSNPA and r-SSNPA are mainly sensitive to passing and trapped particles respectively. In addition, both t-SSNPA and r-SSNPA utilize three vertically stacked arrays with different filter thickness to obtain coarse energy information. The experimental data show that all channels are operational. The signal to noise ratio is typically larger than 10 and the main noise is x-ray induced signal. The active and passive CX signals are clearly observed on t-SSNPA and r-SSNPA during NB modulation. The SSNPA data also indicate significant losses of passing particles during sawteeth, while trapped particles are weakly affected. Fluctuations up to 120 kHz, have been observed on SSNPA, and they are strongly correlated with magnetohydrodynamics (MHD) instabilities.
Flesch, K.; Kremeyer, T.; Schmitz, O.; Soukhanovskii, V.; Wenzel, U.
Direct measurements of the helium (He) fractional neutral pressure in the neutral gas around fusion devices is challenging because of the small mass difference between the abundant D2 molecules and the He atoms which will be produced by D-T fusion. To study He exhaust, Penning gauges were used to measure total neutral pressure assisted by spectroscopy to resolve the D and He partial pressures. In this contribution, initial results are shown from developing this technique into a miniaturized configuration for direct in-situ measurements in the divertor of fusion devices. The configuration is based off a gauge originally designed for the National Spherical Tokamak Experiment-Upgrade (NSTX-U). The goal of this new miniaturized design it to reduce the space required by the gauge on the device and use of the inherent magnetic field of the machine rather than permanent magnets inside the gauge, enabling it to be adapted into a system that can be extended directly into the divertor region. The feasibility test of the method for NSTX-U and the Wendelstein 7-X (W7-X) stellarator are surveyed. For W7-X, a commercial Penning Gauge has been installed on an outboard vacuum flange as a generic feasibility test in the neutral gas environment of a stellarator. At an integration time of 25s, helium lines can be seen down to 10^-5 mbar and H-alpha lines down to 10^-6 mbar. Successful measurement of the total as well as the fractional neutral pressures of He and H has been shown. A first prototype of the miniature Penning gauge has been tested in Madison and shows a near linear power law scaling between current and pressure: I = C*P^n with n = 1.0 - 1.2. Pressure measurements were achieved starting at 10^-3 mbar and down to 10^-6 mbar. A modular gauge is being assembled, which allows easy interchangeability of the anode to test new anode geometries, in order to improve optical access and increase spectroscopic sensitivity. This shall enable an increase of the time resolution of the spectroscopically assisted fractional neutral pressure measurements to up to 1kHz.
A comprehensive set of spectroscopic diagnostics is planned in the National Spherical Torus Experi- ment Upgrade to connect measurements of molybdenum and tungsten divertor sources to scrape-o↵ layer (SOL) and core impurity transport, supporting the installation of high-Z plasma facing compo- nents which is scheduled to begin with a row of molybdenum tiles. Imaging with narrow-bandpass interference filters and high-resolution spectroscopy will be coupled to estimate divertor impurity influxes. Vacuum ultraviolet and extreme ultraviolet spectrometers will allow connecting high-Z sources to SOL transport and core impurity content. The high-Z diagnostics suite complements the existing measurements for low-Z impurities (carbon and lithium), critical for the characterization of sputtering of high-Z materials.
A sequence of H-mode discharges with increasing levels of pre-discharge lithium evaporation (�dose�) was conducted in high triangularity and elongation boundary shape in NSTX. Energy confinement increased, and recycling decreased with increasing lithium dose, similar to a previous lithium dose scan in medium triangularity and elongation plasmas. Data-constrained SOLPS interpretive modeling quantified the edge transport change: the electron particle diffusivity decreased by 10-30x. The electron thermal diffusivity decreased by 4x just inside the top of the pedestal, but increased by up to 5x very near the separatrix. These results provide a baseline expectation for lithium benefits in NSTX-U, which is optimized for a boundary shape similar to the one in this experiment.
The control of divertor heat loads - both steady state and transient - remains a key challenge for the successful operation of ITER and FNSF. Magnetic perturbations provide a promising technique to control ELMs (transients), but understanding their detailed impact is difficult due to their symmetry breaking nature. One approach for reducing steady state heat loads are so called 'advanced divertors' which aim at optimizing the magnetic field configuration: the snowflake and the (super-)X-divertor. It is likely that both concepts - magnetic perturbations and advanced divertors - will have to work together, and we explore their interaction based on the NSTX-U setup. An overview of different divertor configurations under the impact of magnetic perturbations is presented, and the resulting impact on plasma edge transport is investigated with the EMC3-EIRENE code.
Variations in size of the magnetic footprint of the perturbed separatrix are found, which is related to the level of flux expansion on the divertor target. Non-axisymmetric peaking of the heat flux related to the perturbed separatrix is found at the outer strike point, but only in locations where flux expansion is not too large.
The Far-infrared Tangential Interferometer/Polarimeter (FIReTIP) system has been refurbished and
is being reinstalled on the National Spherical Torus Experiment-Upgrade (NSTX-U) to supply
real-time line-integrated core electron density measurements for use in the NSTX-U plasma control
system (PCS) to facilitate real-time density feedback control of the NSTX-U plasma. Inclusion
of a visible light heterodyne interferometer in the FIReTIP system allows for real-time vibration
compensation due to movement of an internally mounted retroreflector and the FIReTIP front-end
optics. Real-time signal correction is achieved through use of a National Instruments CompactRIO
field-programmable gate array.