Search

This paper summarizes recent progress in modeling the interaction between ion cyclotron range of frequency (ICRF) waves and edge plasma with application to ASDEX Upgrade. The basic theories, the development of ICRF and edge plasma codes, the integrated modeling methods and some key results are reviewed. In particular, the following physical aspects are discussed: (1) ICRF power coupling; (2) slow wave propagation; (3) ICRF-rectified sheath; (4) ICRF-induced convection; (5) ICRF-edge turbulence interaction. Moreover, comprehensive integrated modeling strategies by including all necessary codes in one package and solving multiple physical issues self-consistently are discussed.

Major developments on the Mega Amp Spherical Tokamak (MAST) have enabled important advances in support of ITER and the physics basis of a spherical tokamak (ST) based component test facility (CTF), as well as providing new insight into underlying tokamak physics. For example, L–H transition studies benefit from high spatial and temporal resolution measurements of pedestal profile evolution (temperature, density and radial electric field) and in support of pedestal stability studies the edge current density profile has been inferred from motional Stark effect measurements. The influence of the q-profile and E × B flow shear on transport has been studied in MAST and equilibrium flow shear has been included in gyro-kinetic codes, improving comparisons with the experimental data. H-modes exhibit a weaker q and stronger collisionality dependence of heat diffusivity than implied by IPB98(y,2) scaling, which may have important implications for the design of an ST-based CTF. ELM mitigation, an important issue for ITER, has been demonstrated by applying resonant magnetic perturbations (RMPs) using both internal and external coils, but full stabilization of type-I ELMs has not been observed. Modelling shows the importance of including the plasma response to the RMP fields. MAST plasmas with q > 1 and weak central magnetic shear regularly exhibit a long-lived saturated ideal internal mode. Measured plasma braking in the presence of this mode compares well with neo-classical toroidal viscosity theory. In support of basic physics understanding, high resolution Thomson scattering measurements are providing new insight into sawtooth crash dynamics and neo-classical tearing mode critical island widths. Retarding field analyser measurements show elevated ion temperatures in the scrape-off layer of L-mode plasmas and, in the presence of type-I ELMs, ions with energy greater than 500 eV are detected 20 cm outside the separatrix. Disruption mitigation by massive gas injection has reduced divertor heat loads by up to 70%.

The National Spherical Torus Experiment (NSTX) is currently being upgraded to operate at twice the toroidal field and plasma current (up to 1 T and 2 MA), with a second, more tangentially aimed neutral beam (NB) for current and rotation control, allowing for pulse lengths up to 5 s. Recent NSTX physics analyses have addressed topics that will allow NSTX-Upgrade to achieve the research goals critical to a Fusion Nuclear Science Facility. These include producing stable, 100% non-inductive operation in high-performance plasmas, assessing plasma–material interface (PMI) solutions to handle the high heat loads expected in the next-step devices and exploring the unique spherical torus (ST) parameter regimes to advance predictive capability. Non-inductive operation and current profile control in NSTX-U will be facilitated by co-axial helicity injection (CHI) as well as radio frequency (RF) and NB heating. CHI studies using NIMROD indicate that the reconnection process is consistent with the 2D Sweet–Parker theory. Full-wave AORSA simulations show that RF power losses in the scrape-off layer (SOL) increase significantly for both NSTX and NSTX-U when the launched waves propagate in the SOL. Toroidal Alfven eigenmode avalanches and higher frequency Alfven eigenmodes can affect NB-driven current through energy loss and redistribution of fast ions. The inclusion of rotation and kinetic resonances, which depend on collisionality, is necessary for predicting experimental stability thresholds of fast growing ideal wall and resistive wall modes. Neutral beams and neoclassical toroidal viscosity generated from applied 3D fields can be used as actuators to produce rotation profiles optimized for global stability. DEGAS-2 has been used to study the dependence of gas penetration on SOL temperatures and densities for the MGI system being implemented on the Upgrade for disruption mitigation. PMI studies have focused on the effect of ELMs and 3D fields on plasma detachment and heat flux handling. Simulations indicate that snowflake and impurity seeded radiative divertors are candidates for heat flux mitigation in NSTX-U. Studies of lithium evaporation on graphite surfaces indicate that lithium increases oxygen surface concentrations on graphite, and deuterium–oxygen affinity, which increases deuterium pumping and reduces recycling. In situ and test-stand experiments of lithiated graphite and molybdenum indicate temperature-enhanced sputtering, although that test-stand studies also show the potential for heat flux reduction through lithium vapour shielding. Non-linear gyro kinetic simulations have indicated that ion transport can be enhanced by a shear-flow instability, and that non-local effects are necessary to explain the observed rapid changes in plasma turbulence. Predictive simulations have shown agreement between a microtearing-based reduced transport model and the measured electron temperatures in a microtearing unstable regime. Two Alfven eigenmode-driven fast ion transport models have been developed and successfully benchmarked against NSTX data. Upgrade construction is moving on schedule with initial physics research operation of NSTX-U planned for mid-2015.

In a magnetized plasma, such as in fusion devices or the Earth's magnetosphere, several different kinds of waves can simultaneously exist, having very different physical properties. Under the right conditions one wave can quite suddenly convert to another type. Depending on the case, this can be either a great benefit or a problem for the use of waves to heat and control fusion plasmas. Understanding and accurately modeling such behavior is a major computational challenge.

The stability of plasma blobs which have both density and temperature higher than the surrounding plasma, and can transport heat as well as particles, is considered. It is shown that the internal blob temperature profile Te(r) can drive azimuthal rotation or spin vθ(r) about the blob axis, which produces a robust m=2 rotational instability in the interchange limit (k∥=0). The theory includes the effects of the centrifugal and Coriolis forces, the sheared velocity vθ(r), and the axial sheath boundary condition. Estimates show that finite-Larmor-radius stabilization is ineffective, but the sheath conductivity can be strongly stabilizing. The blob rotational instability has only a small direct impact on the particle and energy transport, but it serves as a useful diagnostic for the underlying blob spin, which is an important variable in determining the blob’s radial velocity. A separate branch of temperature-gradient-driven sheath instabilities, predicted in the eikonal limit, is not observed for low mode nu...

The resistive X point mode is shown to be the dominant mode in boundary plasmas in X point divertor geometry. The poloidal fluctuation phase velocity from the simulation results of the resistive X point turbulence shows a structure across the separatrix that is experimentally measured in many fusion devices. The fluctuation phase velocity is larger than the E × B velocity in both L and H mode phases. It is also demonstrated that there is a strong poloidal asymmetry of particle flux in the proximity of the separatrix. Turbulence suppression in the L-H transition results when sources of energy and particles drive sufficient gradients, as in the experiments.

Author(s): Menard, JE; Allain, JP; Battaglia, DJ; Bedoya, F; Bell, RE; Belova, E; Berkery, JW; Boyer, MD; Crocker, N; Diallo, A; Ebrahimi, F; Ferraro, N; Fredrickson, E; Frerichs, H; Gerhardt, S; Gorelenkov, N; Guttenfelder, W; Heidbrink, W; Kaita, R; Kaye, SM; Kriete, DM; Kubota, S; Leblanc, BP; Liu, D; Lunsford, R; Mueller, D; Myers, CE; Ono, M; Park, JK; Podesta, M; Raman, R; Reinke, M; Ren, Y; Sabbagh, SA; Schmitz, O; Scotti, F; Sechrest, Y; Skinner, CH; Smith, DR; Soukhanovskii, V; Stoltzfus-Dueck, T; Yuh, H; Wang, Z; Waters, I; Ahn, JW; Andre, R; Barchfeld, R; Beiersdorfer, P; Bertelli, N; Bhattacharjee, A; Brennan, D; Buttery, R; Capece, A; Canal, G; Canik, J; Chang, CS; Darrow, D; Delgado-Aparicio, L; Domier, C; Ethier, S; Evans, T; Ferron, J; Finkenthal, M; Fonck, R; Gan, K; Gates, D; Goumiri, I; Gray, T; Hosea, J; Humphreys, D; Jarboe, T; Jardin, S; Jaworski, MA; Koel, B; Kolemen, E; Ku, S; La Haye, RJ; Levinton, F; Luhmann, N; Maingi, R; Maqueda, R; McKee, G; Meier, E; Myra, J; Perkins, R | Abstract: 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, the H-mode accessed, magnetic error fields identified and mitigated, and the first physics research campaign carried out. During ten run weeks of operation, NSTX-U surpassed NSTX record pulse-durations and toroidal fields (TF), and high-performance ∼1 MA 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 H98y,2 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 code. Lastly, the materials analysis and particle probe 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.

Results are presented from the 3-D non-local electromagnetic turbulence code BOUT and the linearized shooting code BAL for studies of turbulence in tokamak boundary plasmas and its relationship to the L-H transition, in a realistic divertor plasma geometry. The key results include: (1)?the identification of the dominant resistive X?point mode in divertor geometry; (2)?turbulence suppression in the L-H transition by shear in the E ? B drift velocity, ion diamagnetism and finite polarization; (3) the fact that calculated values of the thermal diffusivities such as ?i are comparable to those inferred from the experimentally measured plasma density and temperature profiles. On the basis of simulation results, a parameterization of the transport is given that includes the dependence on the relevant physical parameters.

The present work is motivated by the importance of r.f. sheaths in determining the antenna-plasma interaction and the sensitivity of the sheaths to the complicated three-dimensional structure of modern ion cyclotron range of frequency (ICRF) antennas. To analyze r.f. sheaths on the plasma facing regions of the launcher, we first calculate the contact points of the tokamak magnetic field lines on the surface of the antenna Faraday screen and nearby limiters for realistic three-dimensional magnetic flux surface and antenna geometries. Next, the r.f. voltage that can drive sheaths at the contact points is determined and used to assess the resulting sheath power dissipation, r.f.-driven sputtering, and r.f.-induced convective cells (which produce edge profile modification). The calculations are embodied in a computer code, ansat (Antenna Sheath Analysis Tool), and sample ansat runs are shown to highlight the physics and geometry-dependent characteristics of the r.f. sheaths and their relationship to the antenna design. One use of ansat is therefore as a design tool, to assess the strengths and weaknesses of a given design with respect to critical voltage handling and edge plasma interaction issues. Additionally, examples are presented where ansat has been useful in the analysis and interpretation of ICRF experiments

Here, the structure and motion of discrete plasma blobs (a.k.a. filaments) in the edge and scrape-off layer of NSTX is studied for representative Ohmic and H-mode discharges. Individual blobs were tracked in the 2D radial versus poloidal plane using data from the gas puff imaging diagnostic taken at 400 000 frames s-1. A database of blob amplitude, size, ellipticity, tilt, and velocity was obtained for ~45 000 individual blobs. Empirical relationships between various properties are described, e.g. blob speed versus amplitude and blob tilt versus ellipticity. The blob velocities are also compared with analytic models.

Progress is reviewed on the simulation of wave-particle interactions in the ion cyclotron range of frequencies (ICRF). Two important aspects of this problem are described. First, mode conversion from a long wavelength fast magnetosonic wave to short wavelength ion Bernstein waves (IBW) and ion cyclotron waves (ICW) is simulated and validation tests of the simulations against experiment are presented. Second, simulations of the quasilinear evolution of nonthermal ion tails during the minority heating are reviewed and experimental validation tests are also discussed. In this paper we describe how access to teraflop computing capability has made it possible to advance the state of the art in this area. We also discuss two aspects of the wave-particle interaction where future work is needed and where in particular access to sub-petaflop and petaflop computing capability would be highly desirable. This work involves the interaction of ICRF waves with energetic neutral beam ions at high ion cyclotron harmonic number and addresses the inclusion of finite ion drift orbit effects in the nonthermal ion tail evolution and the inclusion of nonlinear effects such as RF sheaths in the antenna – edge plasma coupling.

Edge turbulence images have been made using an ultra-high speed CCD camera on both NSTX and Alcator C-Mod. In both cases, the D-alpha or HeI (587.6 nm) line emission from localized deuterium or helium gas puffs was viewed along a local magnetic field line near the outer midplane. Fluctuations in this line emission reflect fluctuations in electron density and/or electron temperature through the atomic excitation rates, which can be modeled using the DEGAS-2 code. The 2-D structure of the measured turbulence can be compared with theoretical simulations based on 3-D fluid models.

In the course of this grant, Lodestar has collaborated with the TFTR ICRF team to carry out an extensive investigation of rf-edge physics relevant to the understanding and operation of the TFTR fast wave (FW) and ion Bernstein wave (IBW) antennas. Their FW work covered four main areas: rf sheath geometry on TFTR, ICRF convective cells and edge profile modifications, impurity production, and plasma glow and Fermi acceleration. Their IBW work, carried out over the last two years, was in support of the IBW transport barrier and poloidal flow generation experiments. A theoretical survey of rf-edge physics expected to be critical for these antennas was carried out. While the physics that was known to be important for FW systems was considered (especially sheath power dissipation, edge density modifications, and impurities), additional physics enters for IBW waves because of the shorter wavelength, slower group velocity, larger electric fields and the different polarization. In particular, ponderomotive forces, while normally negligible in the FW context, can be substantial for the IBW, but the jitter excursion distance of an electron is too large for validity of the usual theory. A new model for the resulting electron heating and density expulsion was developed. A newmore » model of antenna coupling to the electron plasma wave, mode transformation to the IBW, and the effect of reflected waves at the lower hybrid layer was also developed.« less

Reduced model simulations of turbulence in the edge and scrape-off-layer (SOL) region of a spherical torus or tokamak plasma are employed to address the physics of the scrape-off-layer heat-flux width. The simulation model is an electrostatic two-dimensional fluid turbulence model, applied in the plane perpendicular to the magnetic field at the outboard midplane of the torus. The model contains curvature-driven-interchange modes, sheath losses, and both perpendicular turbulent diffusive and convective (blob) transport. These transport processes compete with classical parallel transport to set the SOL width. Midplane SOL profiles of density, temperature, and parallel heat flux are obtained from the simulation and compared with experimental results from the National Spherical Torus Experiment [S. M. Kaye et al., Phys. Plasmas 8, 1977 (2001)] to study the scaling of the heat-flux width with power and plasma current. It is concluded that midplane turbulence is the main contributor to the SOL heat-flux width f...

The relevance of rf-sheath-plasma interactions to ICRF operation in TFTR has been studied. Screen-screen and screen-limiter rf sheaths have been identified and analyzed for the old and new Bay M geometries on TFIR, including sheath voltages driven by both feeders and the central antenna conductor. Calculations have been carried out to determine the effects of ICRF on edge transport (profile flattening by E {times} B convection), impurity influx ({delta}Z{sub eff}) and edge electron heating (FS glow phenomena) as a function of edge plasma parameters, rf power and antenna phasing. It was found that many of the experimentally observed ICRF-edge plasma interactions on TFIR are consistent with the predictions of rf sheath theory.

Work on ICRF interaction with the edge plasma is reported. ICRF generated convective cells have been established as an important mechanism for influencing edge transport and interaction with the H-mode, and for controlling profiles in the tokamak scrape-off-layer. Power dissipation by rf sheaths has been shown to be significant for some misaligned ICRF and IIBW antenna systems. Near-field antenna sheath work has been extended to the far-field case, important for experiments with low single pass absorption. Impurity modeling and Faraday screen design support has been provided for the ICRF community. In the area of core-ICRF physics, the kinetic theory of heating by applied ICRF waves has been extended to retain important geometrical effects relevant to modeling minority heated tokamak plasmas, thereby improving on the physics base that is standard in presently employed codes. Both the quasilinear theory of ion heating, and the plasma response function important in wave codes have been addressed. In separate studies, it has been shown that highly anisotropic minority heated plasmas can give rise to unstable field fluctuations in some situations. A completely separate series of studies have contributed to the understanding of tokamak confinement physics. Additionally, a diffraction formalism has been produced which will bemore » used to access the focusability of lower hybrid, ECH, and gyrotron scattering antennas in dynamic plasma configurations.« less