A new three-dimensional, non-linear magnetohydrodynamics (MHD) model has been extended in MIPS code, incorporating parallel heat diffusivity. The model has been benchmarked against the former MHD model used in MIPS code. A preliminary study of the core density collapse event (CDC) observed in the large helical device (LHD) plasma has been performed using the developed model. The equilibrium has been constructed using HINT code for a typical super dense core discharge in LHD, with vacuum magnetic axis configuration and magnetic axis beta plasma. This configuration corresponds to a plasma with a steep pressure gradient and strong Shafranov shift, which makes the plasma potentially unstable in the LHD. The model shows preliminary characteristics of the CDC event. The plasma is destabilized by high-n ballooning modes in the low-field side region during the linear regime, eventually leading to the collapse of the pressure and density profiles, together with the stochastization of the magnetic field and a shift to low-n modes centered at the core of the plasma after the non-linear coupling at the relaxation regime.
This paper presents a new 3D non-linear MHD model that can simulate the core density collapse (CDC) event observed in the Large Helical Device (LHD) plasma. The model shows that the plasma is destabilized by high-n ballooning modes, leading to the collapse of pressure and density profiles, and a shift to low-n modes in the core.