No Arabic abstract
The physical processes taking place at the edge region are crucial for the operation of tokamaks as they govern the interaction of hot plasma with the vessel walls. Numerical modeling of the edge with state-of-the-art codes attempts to elucidate interactions between neoclassical drifts, turbulence, poloidal and parallel flows that control the physical set-up of the SOL region. Here, we present post-processing analysis of simulations from the gyrokinetic code XGC1, comparing edge turbulence characteristics from a simulation of DIII-D against one of C-Mod. We find that the equilibrium $E times B$ flux across the separatrix has a similar poloidal pattern in both discharges which can be explained by magnetic drifts and trapped ion excursions. However, collisionality is noted to play a major role in that it prevents local charge accumulations from having global effects in C-Mod. In both cases, turbulent electron heat flux is higher than the ion one. This seems to be a universal characteristic of the tokamak edge. We identify turbulent frequencies and growth rates of the dominant mode in both simulations. In C-Mod, these numbers point to the presence of a drift wave. In DIII-D, linear simulations with Gene reveal a trapped electron mode. Furthermore, we present the amplitude and size distributions of the blobs from both simulations. Amplitude distributions are in qualitative agreement with experimental observations while size distributions are consistent with the fact that most blobs are not connecting to the divertor plates and suggest that they are generated by the shearing of the turbulent modes.
Understanding the multi-scale neoclassical and turbulence physics in the edge region (pedestal + scrape-off layer) is required in order to reliably predict performance in future fusion devices. We explore turbulent characteristics in the edge region from a multiscale neoclassical and turbulent XGC1 gyrokinetic simulation in a DIII-D like tokamak geometry, here excluding neutrals and collisions. For an H-mode type plasma with steep pedestal, it is found that the electron density fluctuations increase towards the separatrix, and stay high well into the SOL, reaching a maximum value of $delta n_e / bar{n}_e sim 0.18$. Blobs are observed, born around the magnetic separatrix surface and propagate radially outward with velocities generally less than 1 km/s. Strong poloidal motion of the blobs is also present, near 20 km/s, consistent with $E times B$ rotation. The electron density fluctuations show a negative skewness in the closed field line pedestal regions, consistent with the presence of holes, followed by a transition to strong positive skewness across the separatrix and into the SOL. These simulations indicate that not only neoclassical phenomena, but also turbulence, including the blob-generation mechanism, can remain important in the steep H-mode pedestal and SOL. Qualitative comparisons will be made to experimental observations.
Edge plasma density fluctuations are shown to have a significant effect on the electron cyclotron resonance heating (ECRH) beam in the DIII-D tokamak. Experimental measurements of the ECRH deposition profile have been taken in three operating scenarios: L-mode, H-mode and negative triangularity. Each scenario corresponds to distinct turbulence characteristics in the edge region through which the beam must propagate. The measured ECRH deposition profile is significantly broadened by comparison to the profile predicted by the ray tracing code TORAY-GA and has been shown to scale with the severity of edge turbulence. Conventional ray tracing does not include the effects of turbulence and therefore a 3D full-wave cold plasma finite difference time domain code EMIT-3D is presented and used for the simulations. The turbulence is generated through the Hermes model in the BOUT++ framework which takes as input the measured time averaged electron density, temperature and magnetic field profiles for the specific shot in question. The simulated turbulence is constrained to match the experimentally measured (by use of the BES and DBS systems) correlation length and normalised fluctuation levels. The predictions of the beam broadening from the simulations are found to agree very well with the experimentally-observed broadening in all cases: L-mode, H-mode and negative triangularity. Due to the large gradients within the H-mode edge, the resolution uncertainty and error in the measurement from Thomson scattering and BES diagnostics result in a spread in the simulated turbulence amplitude. In light of this a parameter scan through the range in experimental diagnostic measurement uncertainty has been conducted to explore the impact on beam broadening predictions.
The guiding-center kinetic neoclassical transport code, XGC0, [C.S. Chang et. al, Phys. Plasmas 11, 2649 (2004)] is used to compute the heat fluxes and the heat-load width in the outer divertor plates of Alcator C-Mod and DIII-D tokamaks. The dependence of the width of heat-load fluxes on neoclassical effects, neutral collisions and anomalous transport is investigated using the XGC0 code. The XGC0 code includes realistic X-point geometry, a neutral source model, the effects of collisions, and a diffusion model for anomalous transport. It is observed that width of the XGC0 neoclassical heat-load is approximately inversely proportional to the total plasma current $I_{rm p}$. The scaling of the width of the divertor heat-load with plasma current is examined for an Alcator C-Mod discharge and four DIII-D discharges. The scaling of the divertor heat-load width with plasma current is found to be weaker in the Alcator C-Mod discharge compared to scaling found in the DIII-D discharges. The effect of neutral collisions on the $1/I_{rm p}$ scaling of heat-load width is shown not to be significant. Although inclusion of poloidally uniform anomalous transport results in a deviation from the $1/I_{rm p}$ scaling, the inclusion of the anomalous transport that is driven by ballooning-type instabilities results in recovering the neoclassical $1/I_{rm p}$ scaling. The Bohm or Gyro-Bohm scalings of anomalous transport does not strongly affect the dependence of the heat-load width on plasma current. The inclusion of anomalous transport, in general, results in widening the width of neoclassical divertor heat-load and enhances the neoclassical heat-load fluxes on the divertor plates. Understanding heat transport in the tokamak scrape-off layer plasmas is important for strengthening the basis for predicting divertor conditions in ITER.
The calculation presented in A neoclassical calculation of rotation profiles and comparison with DIII-D measurements by Stacey, Johnson, and Mandrekas, [Physics of Plasmas, 13, (2006)], contains several errors, including the neglect of the toroidal electric field, an unphysical expression for the electrostatic potential, and an unevaluated relation among its parameters. An alternative formulation is discussed.
Bursty fluctuations in the scrape-off layer (SOL) of Alcator C-Mod have been analyzed using gas puff imaging data. This reveals many of the same fluctuation properties as Langmuir probe measurements, including normal distributed fluctuations in the near SOL region while the far SOL plasma is dominated by large amplitude bursts due to radial motion of blob-like structures. Conditional averaging reveals burst wave forms with a fast rise and slow decay and exponentially distributed waiting times. Based on this, a stochastic model of burst dynamics is constructed. The model predicts that fluctuation amplitudes should follow a Gamma distribution. This is shown to be a good description of the gas puff imaging data, validating this aspect of the model.