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

Non-resonant global mode in LHD partial collapse with net toroidal current

May 12, 2021, 8:30 AM
Virtual Event

Virtual Event

Regular Poster Magnetic Fusion Theory and Modelling P3 Posters 3


Katsuji ICHIGUCHI (National Institute for Fusion Science)


A transition from an interchange mode with high mode numbers $(m,n)$ for $m$ (poloidal) and $n$ (toroidal) to a non-resonant $(m,n)=(1,1)$ mode is found in the nonlinear magnetohydrodynamic (MHD) simulation for a Large Helical Device (LHD) plasma with net toroidal current. This transition occurs when the rotational transform is closed to unity in the core region. Because partial collapses caused by (1,1) modes are observed in the LHD experiments with net toroidal current, this transition to the non-resonant mode is a candidate for explaining the partial collapses.

In the magnetic confinement of plasmas, it is crucial that the plasmas are stable against MHD instabilities. Therefore, the stability property is extensively studied in the LHD experiments. In such experiments, partial collapse phenomena are observed when net toroidal current is driven by the neutral beam injection so that the rotational transform $\iota / 2 \pi$ is increased [1]. These collapses are always caused by (1,1) modes. These modes are considered to be pressure driven modes because the equilibria are strongly Mercier unstable. However, according to the theory of the pressure driven modes, the linear growth rate is larger for higher mode numbers. Therefore, it has been required to explain why the partial collapses are caused by the (1,1) modes in LHD.

This problem is investigated by means of three-dimensional nonlinear MHD simulations for LHD plasmas with net toroidal current. In the simulations, the HINT[2] and the MIPS[3] codes are utilized for the equilibrium and the nonlinear dynamics calculations, respectively. In the equilibrium calculation, the vacuum magnetic configuration corresponding to the experiment in Ref.[1] is used. The pressure profile $P_{eq}$ and the axis beta $\beta_0$ are employed as $ P_{eq} = P_0 ( 1 - 0. 68 \rho^2 - 0. 32 \rho^4 ) $ and $\beta_0 =1.6\%$, respectively. Here $\rho $ denotes the square root of the normalized toroidal magnetic flux. The current density of the net toroidal current is assumed as $ J_{eq} = J_0 ( 1 - \rho^2 )^4 $ and the total current $I=80$kA for $B_0 = 3$T are chosen. The values of $\beta_0$ and $I$ are close to those of the experiment in Ref.[1]. In this equilibrium, $\iota / 2\pi $ is close to unity and a no-shear point appears in the core region as shown in Fig.1 (a) and (b). Similar rotational transform profiles are often observed in the experiments.

Profiles of the rotational transform (dashed line) and dominant components of the perturbed pressure (solid lines) at (a)  $t/\tau_A=500$ and (b) $t/\tau_A=800$ in the nonlinear evolution for $I=80$kA.

In the nonlinear dynamics calculation, the values of $S=5.6\times 10^6$, $\nu = \chi_\perp = 2.5 $m$^2$/s are employed for the magnetic Reynolds number, the viscosity and the perpendicular heat conductivity, respectively. These values are relevant to the experiments. Also, $\chi_\parallel = 10^{3} \chi_\perp $ is assumed for the parallel heat conductivity. At $t/\tau_A = 500$ in the early nonlinear phase, the (3,3) component is dominant as shown in Fig.1 (a), where $\tau_A$ denotes the Alfv$\acute{\mbox{e}}$n time. This mode is a typical interchange mode resonant at the $\iota /2 \pi = 1$ surface. These mode numbers are higher than those in the observations. However, at $t/\tau_A =800$ in the further nonlinear phase, the (1,1) component appears and becomes dominant as shown in Fig.1 (b), which corresponds to the observations. This mode is not resonant at the $\iota / 2 \pi = 1 $ surface but is localized around the no-shear region. Therefore, a transition from the resonant (3,3) mode to the non-resonant (1,1) mode occurs. Figure 2 (a) shows that the magnetic field becomes locally stochastic in the unstable region due to the $m=1$ global convection. This convection also causes the partial collapse of the total pressure from one poloidal direction as shown in Fig.2 (b). This collapse tendency is similar to the experimental observation as shown in Fig.2 (c).

The case with no net toroidal current ($I=0$) is also examined as a reference. Here, $\beta_0 =3.0\%$ is employed. As shown in Fig.3 (a), the rotational transform has a monotonically increasing profile. In early nonlinear phase, the (3,2) component is dominant and the (1,1) component is negligibly small. This perturbation is a typical interchange mode resonant at the $\iota / 2 \pi = 2/3$ surface. The total pressure profile decays with a triangular deformation as shown in Fig.3 (b). In the further time evolution, other sideband modes are enhanced and the total pressure decays in multiple poloidal directions. These mode numbers and decay property do not match the observed partial collapse shown in Fig. 2(c). Therefore, it is necessary for the transition that the rotational transform is changed by the net toroidal current so as to have no-shear region where the value is close to unity. Such a profile in the safety factor can appear also in tokamak sawtooth collapses. Therefore, extension of the analysis for flat rotational transform profiles is expected to provide common knowledge for avoiding collapses in current carrying plasmas.

(a) Puncture plots of the magnetic field lines and patterns of positive (red) and negative (blue) pressure perturbations, (b)  bird's eye view of the total pressure and (c) the pressure profile along $z=0$ line at $t/\tau_A=800$ for  $I=80$kA. A typical profile obtained in [1] is also plotted in (c).

(a) Profiles of the rotational transform (thick dashed line) and dominant components of the perturbed pressure (solid lines) and (b) bird's eye view of the total pressure at an early phase of the nonlinear evolution for $I=0$.

[1] S. Sakakibara, et al., 2015 Nucl. Fusion, 55 083020.
[2] Y. Suzuki, et al, Nucl. Fusion 46 (2006) L19.
[3] Y. Todo, et al, Plasma and Fusion Res. 5 (2010) S2062.

Affiliation National Institute for Fusion Science
Country or International Organization Japan

Primary author

Katsuji ICHIGUCHI (National Institute for Fusion Science)


Yasuhiro Suzuki (National Institute for Fusion Science) Yasushi Todo (National Institute for Fusion Science) Satoru Sakakibara (National Institute for Fusion Science) Katsumi Ida (National Institute for Fusion Science) Yuki Takemura (National Institute for Fusion Science) Masahiko Sato (National Institute for Fusion Science) Satoshi Ohdachi (National Institute for Fusion Science) Yoshiro NARUSHIMA (National Institute for Fusion Science) Linda Sugiyama (M.I.T.) Benjamin CARRERAS (BACV Solutions Inc.)

Presentation materials