No Arabic abstract
The XGC1 edge gyrokinetic code is used for a high fidelity prediction for the width of the heat-flux to divertor plates in attached plasma condition. The simulation results are validated against the empirical scaling $lambda_q propto B_P^{-gamma}$ obtained from present tokamak devices, where $lambda_q$ is the divertor heat-flux width mapped to the outboard midplane and $gamma_q=1.19$ as defined by T. Eich et al. [Nucl. Fusion 53 (2013) 093031], and $B_P$ is the magnitude of the poloidal magnetic field at outboard midplane separatrix surface. This empirical scaling predicts $lambda_q leq 1mm$ when extrapolated to ITER, which would require operation with very high separatrix densities $(n_{sep}/n_{Greenwald} > 0.6)$ in the Q=10 scenario to achieve semi-detached plasma operation and high radiative fractions leading to acceptable divertor power fluxes. XGC1 predicts, however, that $lambda_q$ for ITER is over 5 mm, suggesting that operation in the ITER Q=10 scenario with acceptable divertor power loads could be obtained over a wider range of plasma separatrix densities and radiative fractions. The physics reason behind this difference is, according to the XGC1 results, that while the ion magnetic drift contribution to the divertor heat-flux width is wider in the present tokamaks, the turbulent electron contribution is wider in ITER. A high current C-Mod discharge is found to be in a mixed regime: While the heat-flux width by the ion neoclassical magnetic drift is still wider than the turbulent electron heat-flux width, the heat-flux magnitude is dominated by the narrower electron heat-flux.
The distribution of particles and power to plasma-facing components is of key importance in the design of next-generation fusion devices. Power and particle decay lengths have been measured in a number of MAST L-mode and H-mode discharges in order to determine their parametric dependencies, by fitting power and particle flux profiles measured by divertor Langmuir probes, to a convolution of an exponential decay and a Gaussian function. In all discharges analysed, it is found that exponential decay lengths mapped to the midplane are mostly dependent on separatrix electron density and plasma current (or parallel connection length). The widths of the convolved Gaussian functions have been used to derive an approximate diffusion coefficient, which is found to vary from 1m2/s to 7m2/s, and is systematically lower in H-mode compared with L-mode.
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 present paper deals with the surface heat flux estimation with thermocouples (TC) and fiber Bragg grating (FBG) embedded in the plasma facing components (PFC) of the WEST tokamak. A 2D heat transfer model combined with the conjugate gradient method (CGM) and the adjoint state is used to estimate the plasma heat flux deposited on the PFC. The plasma heat flux is characterized by the time evolution of its amplitude and spatial shape on the target (heat flux decay length $lambda^t_q$, power spreading in the private flux region $S^t$ and the strike point location $x_0$). As a first step, five ohmic pulses have been investigated with different magnetic configuration and divertor X-point height varying from 44 to 68 mm from the surface. Despite an outboard shift, the relative displacements of the outer strike point as well as the heat flux decay length derived from the TC/FBG systems are consistent with the magnetic equilibrium reconstruction.
Particle transport due to Ion Temperature Gradient/Trapped Electron (ITG/TE) mode turbulence is investigated using the gyrokinetic code GENE. Both a reduced quasilinear (QL) treatment and nonlinear (NL) simulations are performed for typical tokamak parameters corresponding to ITG dominated turbulence. A selfconsistent treatment is used, where the stationary local profiles are calculated corresponding to zero particle flux simultaneously for electrons and trace impurities. The scaling of the stationary profiles with magnetic shear, safety factor, electron-to-ion temperature ratio, collisionality, toroidal sheared rotation, triangularity, and elongation is investigated. In addition, the effect of different main ion mass on the zero flux condition is discussed. The electron density gradient can significantly affect the stationary impurity profile scaling. It is therefore expected, that a selfconsistent treatment will yield results more comparable to experimental results for parameter scans where the stationary background density profile is sensitive. This is shown to be the case in scans over magnetic shear, collisionality, elongation, and temperature ratio, for which the simultaneous zero flux electron and impurity profiles are calculated. A slight asymmetry between hydrogen, deuterium and tritium with respect to profile peaking is obtained, in particular for scans in collisionality and temperature ratio.
The free-streaming plus recycling model (FSRM) has recently been developed to understand and predict tungsten gross erosion rates from the divertor during edge localized modes (ELMs). In this work, the FSRM was tested against experimental measurements of W sputtering during ELMs, conducted via fast WI spectroscopy. Good agreement is observed using a variety of controlling techniques, including gas puffing, neutral beam heating, and plasma shaping to modify the pedestal stability boundary and thus the ELM behavior. ELM mitigation by pellet pacing was observed to strongly reduce W sputtering by flushing C impurities from the pedestal and reducing the divertor target electron temperature. No reduction of W sputtering was observed during the application of resonant magnetic perturbations (RMPs), in contrast to the prediction of the FSRM. Potential sources of this discrepancy are discussed. Finally, the framework of the FSRM is utilized to predict intra-ELM W sputtering rates in ITER. It is concluded that W erosion during ELMs in ITER will be caused mainly by free-streaming fuel ions, but free-streaming seeded impurities (N or Ne) may increase the erosion rate significantly if present in the pedestal at even the 1% level. Impurity recycling is not expected to cause significant W erosion in ITER due to the very low target electron temperature.