Since 18 of December 2019 uses Nucleus credentials. Visit our help pages for information on how to Register and Sign-in using Nucleus.
May 10 – 15, 2021
Virtual Event
Europe/Vienna timezone
The Conference will be held virtually from 10-15 May 2021

Study of the ECR-heating influence on the anomalous transport of tungsten ions in T-10 plasma

May 12, 2021, 2:00 PM
4h 45m
Virtual Event

Virtual Event

Regular Poster Magnetic Fusion Experiments P4 Posters 4


Mr Ivan Zemtsov (National Research Centre "Kurchatov Institute")


Tungsten removal scenario in discharge 70358: $I_{pl} = 220$ kA, $P_{EC} = 0.75$ MW (on-axis, co-ECCD)

Calculations results for discharge 70358:(a) predictions of experimental W removal (b) diffusion coefficient, (c) pinching rate.

An increasing number of installations for magnetic plasma confinement implement tungsten as a material for the plasma-facing components. For example, tungsten plates are used in the ITER divertor. Tungsten is the heaviest plasma impurity and despite the fact that its content in plasma $c_{\rm{W}} = \sum{n_\rm{W}}/n_e$ is relatively low ($c_{\rm{W}}\approx10^{-4}–10^{-5}$), due to its high emissivity W ions can radiate from the plasma a significant amount of electron energy and thereby greatly affect the plasma parameters. Besides that, most installations conduct experiments with additional heating, which introduces the energy of fast atoms and high-frequency waves into the plasma. This dramatically affects the transport processes in the plasma. The neoclassical transport of impurities is defined theoretically, but their anomalous transport is not. The study of anomalous transport characteristics in regimes with additional heating is a highly relevant task.

The T-10 tokamak is an installation with a circle cross section and a limiter configuration (major radius $R = 1.5$ m, minor radius $a = 0.3$ m) where the carbon (poloidal) limiter was replaced by a tungsten one in 2015. T-10 is equipped with a gyrotron complex with a total capacity of up to 2.5 MW for heating electrons at the second harmonic of electron cyclotron frequency (ECR-heating). Three gyrotrons are used: “A” (1 MW, 140 GHz, co-injection, capable of rotation along the poloidal angle), “B” (0.5 MW, 129 GHz, with injection perpendicular to the plasma current) and “C” (1 MW, 140 GHz, capable of rotation along the toroidal angle, with co-/counter-injection). The electron temperature $T_e(r)$ is measured using ECE-diagnostics and the slope of the spectrum in the SXR-region. CXRS-diagnostics [1, 2] is used to measure the ion temperature $T_i(r)$ and the density of light impurity nuclei $n_Z(r)$. The electron density $n_e(r)$ is measured with a 16-channel interferometer.

The tungsten inflow into the discharge is estimated from the brightness of the WI atom line $\lambda= 400.1$ nm at the W-limiter region. The density of tungsten ions is determined from the integrated radiation losses power recorded by AXUV-diagnostics. In [3], it is shown that the AXUV signal in the central region of the column (within $a/2$) is almost exclusively generated by tungsten radiation. In [4] the preservation of coronal equilibrium is determined on ASDEX Upgrade. It is confirmed for T-10 conditions (in ohmic and ECR regimes) [5]. This makes it possible to determine the profile of the total tungsten ions density $\sum{n_i}(r)$ from the expression:
$$P_{\rm{W}}^{rad}(r) = {n_e}(r)L_{\rm{W}}^{eff}(r)\sum\limits_i {{n_i}(r)},$$ where $P_{\rm{W}}^{rad}$ is power of radiation losses in the central region, registered by AXUV detectors, $L_{\rm{W}}^{eff}$ is tungsten cooling factor [6] calculated in the coronal equilibrium approximation. In T-10 experiments with ECR-heating, tungsten is removed from the center of the plasma column [7, 8]. In discharges where the removal is exponential, the description of the AXUV-signal experimental dynamics in the region from 0 to $a/2$ in the transport model allows to estimate the transport coefficients of W ions. For this purpose, a system of two continuity equations (dynamic and stationary) is solved: $$\left\{ \begin{array}{l} \frac{{\partial {n_\rm{W}}(r,t)}}{{\partial t}} + {\rm{div}}\left[ { - D(r) \cdot \nabla {n_\rm{W}}(r,t) + V(r) \cdot {n_\rm{W}}(r,t)} \right] = {Q_\rm{W}}(r,t),\\ {\rm{div}}\left[ { - D(r) \cdot \nabla {n_\rm{W}}(r) + V(r) \cdot {n_\rm{W}}(r)} \right] = {Q_\rm{W}}(r), \end{array} \right.\ (1)$$ where $n_\rm{W}$ is the full density of tungsten ions, $D$ and $V$ are the desired diffusion coefficient and pinching velocity, $Q_\rm{W}$ is the sum of sources and sinks ($Q_\rm{W}\approx0$ for the central region of the plasma). The solution is carried out by choosing $D(r)$ and $V(r)$ that allow us to describe the temporal behavior of the tungsten radiation profile recorded by the AXUV detector in the best way. The scenario of such a discharge is shown in Fig. 1: at the ECRH stage, there is a quasistationary part (600 – 640 ms) where the electron density and temperature are nearly constant, but the decay of the AXUV-signal in the central region is pronounced. It indicates the removal of tungsten from the center of the plasma. The red line on last plot is an exponential approximation of the AXUV signal to find the experimental decay time $\tau_{exp}$. The system (1) is solved using the STRAHL code integrated in the ASTRA code. Transport coefficients are considered as the sum of neoclassical and anomalous components: $D = D_{neo} + D_{an}$ and $V = V_{neo} + V_{an}$. Neoclassic is calculated using the NEOART code. To conduct the complete study, the selection of discharges with different ECR heating parameters (on-axis / off-axis heating, presence/ absence of sawtooth oscillations, active / suppressed MHD mode $m = 2$) suitable for processing is carried out. The dynamics of tungsten removal from the central region of the cord is simulated for time intervals in which the exponential decay of AXUV-signals occurs without changes in $n_e$ and $T_e$. An example of an AXUV decay description in a discharge with central heating with a power $P_{EC} = 0.75$ MW is shown in Fig. 1 (red line). The calculation is made with the coefficients $D$ and $V$ shown in Fig. 2. The characteristic decay time is $\tau_{exp}\approx14$ ms. This time can be obtained by adding diffusion to the 1.5 – 2 m$^2$/s over the ohmic level. The pinch term in this case turns out to be practically zero at radii up to the $a/2$. Thus, the introduction of on-axis ECR-heating leads to the tungsten removal from the plasma center mainly due to an increase in its anomalous diffusion.


  1. V.A. Krupin Physics of Atomic Nuclei, 2015, Vol. 78, No. 10, pp. 1164–1173
  2. Sushkov A. High-resolution multiwire proportional soft x-ray diagnostic measurements on TCV // Rev. Sci. Instrum. 2008. Vol. 79, No 2 Pt 1. P. 023506
  3. V.A. Krupin Nucl. Fusion 57 (2017) 066041 (9pp)
  4. Pütterich T. Modelling of measured tungsten spectra from ASDEX Upgrade and predictions for ITER // Plasma Phys. Controlled Fusion. IOP Publishing, 2008. Vol. 50, No 8. P. 085016
  5. Zemtsov I. et al. — In: 44th EPS Conf. on Plasma Physics, 2017, p. 4.175.
  6. Demura A.V. Tungsten Ions in Plasmas: Statistical Theory of Radiative-Collisional Processes // Atoms Peace Int. J. Multidisciplinary Digital Publishing Institute, 2015. Vol. 3, No 2. P. 162–181
  7. V.A. Krupin Nucl. Fusion 57 (2017) 066041 (9pp)
  8. Sertoli M. Local effects of ECRH on argon transport in L-mode discharges at ASDEX Upgrade // Plasma Phys. Controlled Fusion. 2011. Vol. 53, No 3. P. 035024
Country or International Organization Russia
Affiliation NRC "Kurchatov Institute"

Primary authors

Mr Ivan Zemtsov (National Research Centre "Kurchatov Institute") Mr Vadim Krupin (NRC "Kurchatov Institute") Mr Maxim Nurgaliev (National Research Centre "Kurchatov Institute") Mr Leonid Klyuchnikov (NRC "Kurchatov Institute") Mr Anton Nemets (NRC "Kurchatov Institute") Alexey Dnestrovskiy (NRC Kurchatov Institute)


Mrs Tatiana Myalton (NRC "Kurchatov Institute") Mr Dmitrii Sarychev (NRC "Kurchatov Institute") Georgy Subbotin (MIPT) Mr Dmitrii Sergeev (NRC "Kurchatov Institute") Mr Nikita Solovev (NRC "Kurchatov Institute") Mr Alexander Borschegovskiy (NRC "Kurchatov Institute")

Presentation materials