Monte Carlo simulations of tungsten redeposition at the divertor target

Recent modeling of controlled edge-localized modes (ELMs) in ITER with tungsten (W) divertor target plates by the SOLPS code package predicted high electron temperatures (>100 eV) and densities (>1 × 1021 m−3) at the outer target. Under certain scenarios W sputtered during ELMs can penetrate into the core in quantities large enough to cause deterioration of the discharge performance, as was shown by coupled SOLPS5.0/STRAHL/ASTRA runs. The net sputtering yield, however, was expected to be dramatically reduced by the ‘prompt redeposition’ during the first Larmor gyration of W1+ (Fussman et al 1995 Proc. 15th Int. Conf. on Plasma Physics and Controlled Nuclear Fusion Research (Vienna: IAEA) vol 2, p 143). Under high ne/Te conditions at the target during ITER ELMs, prompt redeposition would reduce W sputtering by factor p−2 ∼ 104 (with p ≡ τionωgyro ∼ 0.01). However, this relation does not include the effects of multiple ionizations of sputtered W atoms and the electric field in the magnetic pre-sheath (MPS, or ‘Chodura sheath’) and Debye sheath (DS). Monte Carlo simulations of W redeposition with the inclusion of these effects are described in the paper. It is shown that for p ≪ 1, the inclusion of multiple W ionizations and the electric field in the MPS and DS changes the physics of W redeposition from geometrical effects of circular gyro-orbits hitting the target surface, to mainly energy considerations; the key effect is the electric potential barrier for ions trying to escape into the main plasma. The overwhelming majority of ions are drawn back to the target by a strong attracting electric field. It is also shown that the possibility of a W self-sputtering avalanche by ions circulating in the MPS can be ruled out due to the smallness of the sputtered W neutral energies, which means that they do not penetrate very far into the MPS before ionizing; thus the W ions do not gain a large kinetic energy as they are accelerated back to the surface by the MPS/DS electric field; this leads to modest self-sputtering yields. The results of these simulations are applicable to a wide range of plasma conditions at the target plates that can be encountered in various magnetic confinement fusion devices.

[1]  Jeffrey N. Brooks,et al.  Near‐surface sputtered particle transport for an oblique incidence magnetic field plasma , 1990 .

[2]  J. Roth,et al.  Erosion and redeposition in the ASDEX Upgrade divertor , 1994 .

[3]  P. C. Stangeby,et al.  Obtaining reactor-relevant divertor conditions in tokamaks , 2011 .

[4]  M. Sugihara,et al.  A full tungsten divertor for ITER: Physics issues and design status , 2013 .

[5]  J. Roth,et al.  Threshold energy for sputtering and its dependence on angle of incidence , 1993 .

[6]  David N Ruzic,et al.  Modeling of particle-surface reflections including surface roughness characterized by fractal geometry , 1989 .

[7]  R. Dux,et al.  Main chamber sources and edge transport of tungsten in H-mode plasmas at ASDEX Upgrade , 2011 .

[8]  L. Horton,et al.  Type-I ELM substructure on the divertor target plates in ASDEX Upgrade , 2005 .

[9]  Effect of surface roughness and substrate material on carbon deposition in tokamak TEXTOR , 2008 .

[10]  K. Hoshino,et al.  Multi-dimensional impurity transport code by Monte Carlo method including gyro-orbit effects , 2003 .

[11]  A. Loarte,et al.  Progress on the application of ELM control schemes to ITER scenarios from the non-active phase to DT operation , 2014 .

[12]  P. C. Stangeby,et al.  The Chodura sheath for angles of a few degrees between the magnetic field and the surface of divertor targets and limiters , 2012 .

[13]  M. W. Thompson II. The energy spectrum of ejected atoms during the high energy sputtering of gold , 1968 .

[14]  F L Tabarés,et al.  Tritium inventory control during ITER operation under carbon plasma-facing components by nitrogen-based plasma chemistry: a review , 2013 .

[15]  D. C. Griffin,et al.  Electron-impact ionization of atomic ions in the W isonuclear sequence , 2005 .

[16]  M. Groth,et al.  1D kinetic modelling of the JET SOL with tungsten divertor plates , 2013, Journal of nuclear materials. Journal des materiaux nucleaires.