Simulation of density fluctuations before the L-H transition for Hydrogen and Deuterium plasmas in the DIII-D tokamak using the BOUT++ code

Abstract

A six-field two-fluid model has been used to simulate density fluctuations. The equilibrium is generated by experimental measurements for both Deuterium (D) and Hydrogen (H) plasmas at the lowest densities of DIII-D low to high confinement (L-H) transition experiments. In linear simulations, the unstable modes are found to be resistive ballooning modes with the most unstable mode number n = 30 or .The ion diamagnetic drift and convection flow are balanced when the radial electric field (Er) is calculated from the pressure profile without net flow. The curvature drift plays an important role in this stage. Two poloidally counter propagating modes are found in the nonlinear simulation of the D plasma at electron density m−3 near the separatrix while a single ion mode is found in the H plasma at the similar lower density, which are consistent with the experimental results measured by the beam emission spectroscopy (BES) diagnostic on the DIII-D tokamak. The frequency of the electron modes and the ion modes are about 40 kHz and 10 kHz respectively. The poloidal wave number is about 0.2 cm−1 () for both ion and electron modes.The particle flux, ion and electron heat fluxes are ∼3.5–6 times larger for the H plasma than the D plasma, which makes it harder to achieve H-mode for the same heating power. The change of the atomic mass number A from 2 to 1 using D plasma equilibrium make little difference on the flux. Increase the electric field will suppress the density fluctuation. The electric field scan and ion mass scan results show that the dual-mode results primarily from differences in the profiles rather than the ion mass.