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

Isotope effects in internal transport barrier strength on Large Helical Device

13 May 2021, 14:00
4h 45m
Nice, France

Nice, France

Regular Poster Magnetic Fusion Experiments P6 Posters 6


Tatsuya Kobayashi (National Institute for Fusion Science)


Strength of the internal transport barrier (ITB) is quantified in deuterium (D) and hydrogen (H) plasmas using a unique criterion based on the transport nonlinearity on the temperature. Two distinct isotope effects are found: (i) Stronger ITBs in D plasmas and (ii) edge confinement degradation accompanied by the ITB formation emerging in H plasmas. Principal component analysis reveals the important role of the density peaking for the strong ITB formation.
The isotope effect in magnetically confined plasma physics has been a long-standing mystery unsolved. In general, there are a number of experimental case studies on the isotope effect in tokamaks, but less in stellarators/heliotrons. The isotope effect is particularly prominent in transport barrier formation in tokamaks. For systematic understanding of the background physics, the isotope effect in the transport barrier property in stellarators/heliotrons need be assessed. Here, we report a recent progress on isotope effect studies for the ITB in Large Helical Device (LHD).
Unlike the case of tokamak plasmas, there is no generally accepted criterion for the ITB strength in stellarators/heliotrons. Here, a new criterion for the ITB strength is proposed by defining a unique scalar coefficient. The typical L-mode plasmas in LHD are characterized by the dome-shaped temperature profile with the diffusion coefficient being proportional to the temperature to the power of a factor $\alpha$, i.e., $\chi \propto T^\alpha$, where $\alpha=1$ is widely applicable. The reference L-mode profile $T_{\rm L}^{\rm ref}(r)$ is defined by extrapolating the edge temperature profile to the core according to the solution of the thermal diffusion equation with $\chi \propto T^\alpha$. By comparing $T_{\rm L}^{\rm ref}(r)$ with the entire temperature profile, the profile gain factor as the ITB strength is defined as $G_{1.0} = \int n (r) T (r) dV/\int n (r) T_{\rm L}^{\rm ref} (r) dV$, where the subscript 1.0 represents the $\alpha$ value used [A].
Radial profiles of the ion temperature in (a) D and (b) H plasmas, and (c) $\bar{n}_{\rm e}$ dependence of the ITB strength $G_{1.0}$ with inserts showing the electron density profiles and the $G_{1.0}$ dependence of the edge ion temperature, respectively. Stronger ITBs formed in D plasmas are revealed.
The proposed method is applied to both D and H discharges with the line averaged density $\bar{n}_{\rm e}$ scan in the shot-to-shot basis. In LHD, the ITB in the ion temperature profile is typically formed when an intense neutral beam heating is applied to low $\bar{n}_{\rm e}$ plasmas. Figures 1 (a) and (b) compare the high and low $\bar{n}_{\rm e}$ discharges for corresponding pairs of D and H plasmas, where the neoclassical transport is optimized by shifting the magnetic axis inward. Symbols and curves correspond to the measured data and $T_{\rm L}^{\rm ref}(r)$, respectively. Inserts present the electron density $n_{\rm e}$ profiles in the low $\bar{n}_{\rm e}$ discharges, showing that $n_{\rm e}$ tends to peak in D plasmas. However, the peaking of $n_{\rm e}$ is very weak and no ITB in $n_{\rm e}$ profile is present. Stronger ion temperature ITB is formed in the low $\bar{n}_{\rm e}$ D plasma with $G_{1.0} \sim 1.5$ compared to the H plasma with $G_{1.0} \sim 1.3$. The ITB foot-points in D and H plasmas are similar. Unlike the L-H transition, it is hard to define the threshold power for the ITB in LHD. In stead, we discuss $\bar{n}_{\rm e}$ dependence of the ITB strength $G_{1.0}$ in Fig. 1 (c), where the systematic comparison between D and H plasmas is performed. When $\bar{n}_{\rm e}$ is high, $G_{1.0} \sim 1$ that corresponds to the L-mode, and decreasing $\bar{n}_{\rm e}$ leads to a non-monotonic increase in $G_{1.0}$. Larger $G_{1.0}$ is routinely observed in D plasmas in $\bar{n}_{\rm e} < 2 \times 10^{19} {\rm m}^{-3}$. In addition, the edge temperature decreases when the ITB is formed in the H case as shown in Fig. 1 (b). The insert in Fig. 1 (c) shows the $G_{1.0}$ dependence of the edge ion temperature $T_{\rm i,edge}$ averaged in $0.68 < r_{\rm eff}/a_{99} < 0.83$. As $G_{1.0}$ increases, $T_{\rm i,edge}$ clearly decays in H plasmas. This limits the core temperature increment even in the presence of the ITB in H plasmas [B].
The ITB strength as a function of the line averaged density and inverse scale length of the electron density.
It is essential to identify parameters that play a role for determining the ITB strength. To resolve the relation among parameters that change with $\bar{n}_{\rm e}$, the principal component analysis (PCA) is performed. As appeared in literature, the inverse electron density gradient length $L_{n_{\rm e}}^{-1}$, the carbon density $n_{\rm c}$ and its inverse gradient length $L_{n_{\rm c}}^{-1}$ are predicted to impact on the ITB formation. The PCA is applied for the database in the five-dimensional parameter space ($\bar{n}_{\rm e},G_{1.0},L_{n_{\rm e}}^{-1},n_{\rm c},L_{n_{\rm c}}^{-1}$), where data with different magnetic configurations and the isotope mass are taken all at once. Here, the local parameters are obtained at the ITB foot-point of $r_{\rm eff}/a_{99} = 0.6$. Note that the density profiles apparently having no ITBs are treated as the background elements that possibly affect the ITB property. It is found that the primary and secondary principal components are almost equivalent to $\bar{n}_{\rm e}$ and $L_{n_{\rm e}}^{-1}$, respectively, and $\bar{n}_{\rm e}$ and $L_{n_{\rm e}}^{-1}$ have high correlations with $G_{1.0}$. According to the PCA, $G_{1.0}$ is plotted as a function of $\bar{n}_{\rm e}$ and $L_{n_{\rm e}}^{-1}$ for D and H plasmas in Fig. 2. Stronger ITBs are formed when $L_{n_{\rm e}}^{-1}$ is large and $\bar{n}_{\rm e}$ is small. A working hypothesis is drawn from the observation: In D plasmas $n_{\rm e}$ profile tends to peak due to the particle transport or the beam fueling. The peaked $n_{\rm e}$ profile enhances the ITB formation, possibly due to the ion temperature gradient mode stabilization [C]. The radial electric field shear remains approximately unchanged with $\bar{n}_{\rm e}$, showing less impact on the ITB formation.
The isotope effect in the transport barrier formation property is stronger than that in the L-mode plasma confinement in LHD [D-F]. This is qualitatively similar to the tokamak case, where the isotope effect in the transport barrier threshold power is much clearer with respect to that in the confinement scaling exponent. A detailed comparison between LHD and tokamaks stimulates the phenomenological understanding of the isotope effect.
[A] T. Kobayashi et al., Plasma Phys. Control. Fusion 61, 085005 (2019)
[B] T. Kobayashi et al., Sci. Rep. 9, 15913 (2019)
[C] M. Z. Tokar et al., Phys. Rev. Lett. 84 895 (2000)
[D] H. Yamada et al., Phys. Rev. Lett. 123, 185001 (2019)
[E] T. Kobayashi et al., submitted to Nucl. Fusion (2020)
[F] H. Yamada et al., this conference (2020)

Affiliation National Institute for Fusion Science
Country or International Organization Japan

Primary author

Tatsuya Kobayashi (National Institute for Fusion Science)


Hiromi Takahashi (National Institute for Fusion Science) Kenichi Nagaoka (National Institute for Fusion Science) Makoto Sasaki (Research Institute for Applied Mechanics Kyushu University) Motoki Nakata (National Institute for Fusion Science) Masayuki Yokoyama (National Institute for Fusion Science) Ryousuke Seki (National Institute for Fusion Science) Mikirou Yoshinuma (National Institute for Fusion Science) Katsumi Ida (National Institute for Fusion Science)

Presentation Materials

There are no materials yet.