A formulation of the anisotropic pressure magnetohydrodynamic equilibrium problem for three-dimensional plasmas with imposed nested magnetic surfaces is developed based on a bi-Maxwellian model of the distribution function for the energetic particle species. The hot particle distribution function satisfies the constraint . Large parallel and perpendicular anisotropy factors can be explored within the model through the choice of the hot particle perpendicular to parallel temperature ratio T⊥/T‖. A fixed boundary version of the VMEC code has been adapted to numerically compute three-dimensional anisotropic pressure equilibria. Applications to a 10-field period Heliotron device and a 2-field quasiaxisymmetric stellarator demonstrate that the pressures do not vary significantly around the magnetic surfaces when the total parallel pressure p‖ is larger than its perpendicular counterpart p⊥. For off-axis hot particle deposition with p⊥ > p‖, p⊥ concentrates in the region where the energetic particles are generated. On the other hand, p‖ is distributed roughly uniformly around the flux surfaces in the Heliotron but is localized on the low field side in the quasiaxisymmetric machine. The hot particle density structure correlates more closely with the corresponding perpendicular rather than with parallel pressure. The specific form for the definition of β that best correlates with the Shafranov shift is identified.

Reversed-shear Alfvén eigenmodes were observed for the first time in a helical plasma having negative q₀'' (the curvature of the safety factor q at the zero shear layer). The frequency is swept downward and upward sequentially via the time variation in the maximum of q. The eigenmodes calculated by ideal MHD theory are consistent with the experimental data. The frequency sweeping is mainly determined by the effects of energetic ions and the bulk pressure gradient. Coupling of reversed-shear Alfvén eigenmodes with energetic ion driven geodesic acoustic modes generates a multitude of frequency-sweeping modes.

The structure of the radial electric field and heat transport at the magnetic island in the Large Helical Device is investigated by measuring the radial profile of poloidal flow with charge exchange spectroscopy. The convective poloidal flow inside the island is observed when the n/m=1/1 external perturbation field becomes large enough to increase the magnetic island width above a critical value (15-20% of minor radius) in LHD. This convective poloidal flow results in a non-flat space potential inside the magnetic island. The sign of the curvature of the space potential depends on the radial electric field at the boundary of the magnetic island. The heat transport inside the magnetic island is studied with a cold pulse propagation technique. The experimental results show the existence of the radial electric field shear at the boundary of the magnetic island and a reduction of heat transport inside the magnetic island

Characteristics of MHD instabilities and their impacts on plasma confinement are studied in current free plasmas of the Large Helical Device(LHD). Spontaneous L-H transition is often observed in high beta plasmas in the range of 2% averaged beta at low toroidal field (B t ≤ 0.6T). The stored energy rapidly rises by the transition, but quickly saturates due to the growth of m=2/n=3 and m=2/n=2 modes (m and n: poloidal and toroidal mode numbers) excited in the plasma edge region. Even in low beta plasmas, ELM like activities are sometimes induced in high performance plasmas with steep edge pressure gradient, and transiently reduce the stored energy by about 10%. Energetic ion driven MHD modes such as Alfven eigenmodes are studied in the very wide range of characteristic parameters: the averaged beta of energetic ions <β b// > up to 5% and the ratio of energetic ion velocity to the Alfven velocity V b// /V A up to 2.5. In addition to the observation of toroidicity induced Alfven eigenmodes (TAEs), coherent magnetic fluctuations of helicity induced Alfven eigenmodes (HAEs) have been observed for the first time in NBI heated plasmas. Transition of TAE to global Alfven eigenmode(GAE) is also observed in a discharge with temporal evolution of the rotational transform profile, having a similarity to the phenomenon in a reversed shear tokamak. At the low magnetic field, bursting TAEs transiently induce a significant loss of energetic ions, but lead to the transient improvement of bulk plasma confinement in the plasma central region.

Recent experimental results in the Large Helical Device have indicated that a large pressure gradient can be formed beyond the stability criterion for the Mercier (high-n) mode. While the stability against an interchange mode is violated in the inward-shifted configuration due to an enhancement of the magnetic hill, the neoclassical transport and confinement of high-energy particle are, in contrast, improved by this inward shift. Mitigation of the unfavourable effects of MHD instability has led to a significant extension of the operational regime. Achievements of the stored energy of 1 MJ and the volume-averaged beta of 3% are representative results from this finding. A confinement enhancement factor above the international stellarator scaling ISS95 is also maintained around 1.5 towards a volume-averaged beta, β , of