Theoretical studies of equilibrium beta limit in LHD plasmas

Abstract

This study deals with the equilibrium beta limit for high-beta plasmas of the Large Helical Device using a 3D equilibrium calculation code, HINT, which is an initial value solver based on the relaxation method without an assumption of nested flux surfaces. For the finite beta equilibrium, due to the nonlinear 3D equilibrium response, the flux surface breaking begins, and the magnetic field line becomes stochastic in the peripheral region. However, although the magnetic axis shifts until the conventional theoretical beta limit, the separatrix limiting the equilibrium does not appear in the plasma core. For the high beta equilibrium, breaking of the force balance begins in the strongly stochastized region. To keep the force balance in that region, the pressure gradient decreases and the fixed profile is relaxed. As a result, the volume averaged beta, which is equivalent to the plasma stored energy, saturates, although the peak beta increases. The beta value, where the force balance begins breaking in the inside of the stochastic region, is proposed as a new index of the equilibrium beta limit for the typical heliotron-type stellarator. According to this new index of the equilibrium beta limit, the achieved beta value in the experiment of the Large Helical Device is still lower than the equilibrium beta limit. This suggests that if the heating power can be increased, the achieved beta value may still increase.