Since 18 of December 2019 conferences.iaea.org 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

Linear Analysis of Cross-field Dynamics with Feedback Instability on Detached Divertor Plasmas

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

Nice, France

Regular Poster Magnetic Fusion Theory and Modelling P4 Posters 4

Speaker

Hiroki Hasegawa (National Institute for Fusion Science)

Description

The theoretical model of the feedback instability is proposed to explain the mechanism of the correlation between the detachment and the cross-field plasma transport. It is shown that (1) the feedback instability on the detached divertor plasma can be induced in a certain condition in which the recombination frequency $\nu_{\textrm{rec}}$ is larger than the ion cyclotron frequency $\Omega_{\textrm{ci}}$ in the recombination region and the density gradient and the electric field in the direction perpendicular to the magnetic flux surface are not zero, and that (2) the feedback instability can provide the cross-field plasma transport in the boundary layer of magnetic fusion torus devices.

The correlation between the detachment and the cross-field plasma transport in the boundary layer has been reported in various magnetic confinement devices, that is, tokamak$^{1}$, helical$^{2}$, and linear$^{3}$ devices. Such a correlation is expected to expand the width of the heat flux to the divertor target, i.e., $\lambda_{\textrm{q}}$. However, the physical dynamics of the correlation has not been revealed. In this study, we investigate the cross-field dynamics in the detached plasma state with the coupling model between magnetized plasmas characterized by different current mechanisms. In the recombination region in front of a divertor target, $\nu_{\textrm{rec}}$ can be larger than $\Omega_{\textrm{ci}}$ because of the high density and the low temperature. In such a situation, the cross-field motion of ions is mainly in the direction of the electric field, while that of electrons is almost in the direction of the $E \times B$ drift. Thus, the difference in the direction of motion may provide the cross-field current in the recombination region. On the other hand, the cross-field current can be generated by only the polarization, the grad-$B$, and the diamagnetic drifts in the upstream plasma. We have considered whether such a difference between the current mechanisms in each region induces the cross-field plasma transport.

In this study, we have derived the linear dispersion relation from the continuity equations,
$\displaystyle \frac{\partial n^{\text{P}}}{\partial t} + \mathbf{\nabla}_{\perp} \cdot (n^{\textrm{P}} \mathbf{v}_{s \perp}^{\textrm{P}}) + \frac{\Gamma_{s \parallel}^{\textrm{P}}(z = L_{z}) - \Gamma_{s \parallel}^{\textrm{B}}}{L_{z}} = 0 \; \; \; \; \; \; \; (1)$
and
$\displaystyle \frac{\partial n^{\textrm{R}}}{\partial t} + \mathbf{\nabla}_{\perp} \cdot (n^{\textrm{R}} \mathbf{v}_{s \perp}^{\textrm{R}}) + \frac{\Gamma_{s \parallel}^{\textrm{B}} - \Gamma_{s \parallel}^{\textrm{R}}(z = -h)}{L_{z}} = -\alpha \left[(n^{\textrm{R}})^{2}-(n_{0}^{\textrm{R}})^{2} \right], \; \; \; \; \; \; \; (2) $
and the charge conservation equations,
$\displaystyle \mathbf{\nabla}_{\perp} \cdot \mathbf{j}_{s \perp}^{\textrm{P}} + \frac{j_{\parallel}^{\textrm{P}}(z = L_{z}) - j_{\parallel}^{\textrm{B}}}{L_{z}} = 0 \; \; \; \; \; \; \; (3) $
and
$\displaystyle \mathbf{\nabla_{\perp}} \cdot \mathbf{j}_{s \perp}^{\textrm{R}} + \frac{j_{\parallel}^{\textrm{B}} - j_{\parallel}^{\textrm{R}}(z = -h)}{L_{z}} = 0, \; \; \; \; \; \; \; (4) $
in the upstream plasma and the recombination region in the simple configuration as shown in Fig. 1. In this configuration, the magnetic field is parallel to the $z$ axis and the $x$ and $y$ directions correspond to the direction perpendicular to the magnetic flux surface and the toroidal direction in torus devices, respectively. Here, $n$ is the plasma density, $\mathbf{v}_{s \perp}$ is the flow velocity perpendicular to the magnetic field, $\Gamma_{s \parallel}$ is the parallel flux, $\alpha$ is the recombination coefficient, $n_{0}$ is the plasma density at the equilibrium state, $\mathbf{j}_{\perp}$ and $\mathbf{j}_{\parallel}$ are the perpendicular and the parallel currents, the superscripts P, R, and B indicate the quantities in the upstream plasma, in the recombination region, and at the boundary between those regions, respectively, and the subscript $s$ represents the particle species. The upstream plasma flow velocity $\mathbf{v}_{s \perp}^{\text{P}}$ is composed of the $E \times B$, the polarization, the grad-$B$, and the diamagnetic drifts, while the recombination region flow velocity $\mathbf{v}_{s \perp}^{\text{R}}$ includes each drift with the Hall mobility and the motion in the direction of the perpendicular electric field with the Pedersen mobility. Linearizing Eqs. (1)-(4), as a result, we obtain the cubic equation regarding the frequency $\omega$ as the dispersion relation. It is found that the one mode of them has a positive growth rate under a certain condition. Figure 2 shows the dependence of the growth rate $\gamma$ of the unstable mode, i.e., the feedback instability mode, on the wave number $k$ and the propagation direction $\theta$. In Fig. 3 we present the dependence of the group velocity $v_{\textrm{g}} = \partial \omega / \partial k$ of the unstable mode on $k$ and $\theta$. Here, the typical parameters for fusion torus devices are assumed as follows: $B = 5 \; \textrm{T}$, $\partial B / \partial x = -1 \; \textrm{T}/\textrm{m}$, $n_{0}^{\textrm{P}} = 5 \times 10^{19} \; \textrm{m}^{-3}$, $\partial n_{0}^{\textrm{P}} / \partial x = -1.67 \times 10^{21} \; \textrm{m}^{-4}$, the initial electric fields $E_{x0}^{\textrm{P}} = E_{x0}^{\textrm{R}} = -100 \; \textrm{V}/\textrm{m}$, the electron and ion temperatures $T_{\textrm{e}}^{\textrm{P}} = T_{\textrm{i}}^{\textrm{P}} = 50 \; \textrm{eV}$, $T_{\textrm{e}}^{\textrm{R}} = T_{\textrm{i}}^{\textrm{R}} = 0.3 \; \textrm{eV}$, $\nu_{\textrm{rec}} / \Omega_{\textrm{ci}} = 10$, $L_{z} = 10 \; \textrm{m}$, $h = 0.3 \; \textrm{m}$, the ion-to-electron mass ratio $m_{\textrm{i}} / m_{\textrm{e}} = 3.67 \times 10^{3}$, and the ion-to-electron charge ratio $q_{\textrm{i}} / |q_{\textrm{e}}| = 1$. In those figures, the area inside the red curve designates the unstable region in which the feedback instability can be induced. Thus, Figs. 2 and 3 indicate that the waves $k \rho_{\textrm{s}}^{\textrm{P}} > 0.8$ and $ \theta \sim 3\pi/4$ can transport the plasma lump with the speed $\sim 0.002 \; c_{\textrm{s}}^{\textrm{P}}$. The simple estimation shows that the maximum of heat flux density is reduced to $1 - (v_{{\textrm{g}}x}/c_{\textrm{s}}^{\textrm{R}})(\tilde{n}/n_{0}^{\textrm{R}})(h/\lambda_{\textrm{q}}^{\textrm{B}}) \approx 82$% of the initial value and that the heat flux width is expanded to $1 + (h/\lambda_{\textrm{q}}^{\textrm{B}})(v_{\textrm{g}x}/c_{\textrm{s}}^{\textrm{R}}) \approx 280$% of the initial width if $\tilde{n}/n_{0}^{\textrm{R}} \sim 0.1$ and $\lambda_{\textrm{q}}^{\textrm{B}} \sim 3 \; \textrm{mm}$. Here, $v_{{\textrm{g}}x}$ is the $x$ component of $v_{\textrm{g}}$ and $\tilde{n}$ is the time averaged density of the transported plasma lump.

Furthermore, to verify the feedback instability model, the spiraling plasma ejection observed around the recombination front under the detached divertor condition in the NAGDIS-II linear device experiment$^{3}$ is analyzed, in which the NAGDIS-II contributes to the establishments of the detachment and the cross-field transport mechanisms for future fusion reactors such as ITER and DEMO. The radial speed $v_{r} \sim 80 \; \textrm{m}/\textrm{s}$ at $r \sim 20 \; \textrm{mm}$ and azimuthal speed $v_{\theta} \sim 200 \; \textrm{m}/\textrm{s}$ at $r \sim 5 \; \textrm{mm}$ obtained in the experiment are in good agreement with $v_{{\textrm{g}}x} \sim 800 \; \textrm{m}/\textrm{s}$ and $v_{{\textrm{g}}y} \sim 200 \; \textrm{m}/\textrm{s}$ estimated by the theoretical model with the NAGDIS-II parameters if $v_{r}$ is reduced as $r$ increases.

$^{1}$Potzel S et al. 2013 J. Nucl. Mater. 438 S285.
$^{2}$Tanaka H. et al. 2010 Phys. Plasmas 17 102509.
$^{3}$Tanaka H. et al. 2018 Plasma Phys. Control. Fusion 60 075013.

Configuration of the system for the theoretical model. The wave vector $\mathbf{k}$ is assumed to be perpendicular to the magnetic field $\mathbf{B}$ as $\mathbf{k} = (k \cos \theta, k \sin \theta, 0)$.
Dependence of the growth rate $\gamma$ of the unstable mode on the wave number $k$ and the propagation direction $\theta$ whose definitions are shown in Fig. 1. The area inside the red contour line which represents $\gamma = 0$ designates the unstable region. Here, $\rho_{\textrm{s}}^{\textrm{P}} = c_{\textrm{s}}^{\textrm{P}} / \Omega_{\textrm{ci}} $ and $c_{\textrm{s}}^{\textrm{P}}$ is the ion acoustic speed defined by $c_{\textrm{s}}^{\textrm{P}} = (T_{\textrm{e}}^{\textrm{P}}/m_{\textrm{i}})^{1/2}$.
Dependence of the group velocity $\partial \omega / \partial k$  of the unstable mode on $k$ and $\theta$ whose definitions are shown in Fig. 1. The area inside the red curve which represents $\gamma = 0$ designates the unstable region.

Affiliation National Institute for Fusion Science
Country or International Organization Japan

Primary author

Hiroki Hasegawa (National Institute for Fusion Science)

Co-authors

Dr Hirohiko Tanaka (Graduate School of Engineering, Nagoya University) Prof. Seiji Ishiguro (National Institute for Fusion Science )

Presentation Materials

There are no materials yet.