High-frequency energetic particle driven geodesic acoustic modes (EGAM) observed in the large helical device plasmas are investigated using a hybrid simulation code for energetic particles and magnetohydrodynamics (MHD). Energetic particle inertia is incorporated in the MHD momentum equation for the simulation where the beam ion density is comparable to the bulk plasma density. Bump-on-tail type beam ion velocity distribution created by slowing down and charge exchange is considered. It is demonstrated that EGAMs have frequencies higher than the geodesic acoustic modes and the dependence on bulk plasma temperature is weak if (1) energetic particle density is comparable to the bulk plasma density and (2) charge exchange time (τcx) is sufficiently shorter than the slowing down time (τs) to create a bump-on-tail type distribution. The frequency of high-frequency EGAM rises as the energetic particle pressure increases under the condition of high energetic particle pressure. The frequency also increases as the energetic particle pitch angle distribution shifts to higher transit frequency. It is found that there are two kinds of particles resonant with EGAM: (1) trapped particles and (2) passing particles with transit frequency close to the mode frequency. The EGAMs investigated in this work are destabilized primarily by the passing particles whose transit frequencies are close to the EGAM frequency.

A nonlinear simulation of the energetic particle driven instabilities in the Chinese First Quasi-Axisymmetric Stellarator (CFQS) has been conducted for the first time. MEGA, a hybrid simulation code for energetic particles interacting with a magneto-hydrodynamic (MHD) fluid, was used in the present work. Both the m/n=3/1 energetic-particle-mode (EPM) like mode and the m/n=5/2 toroidal Alfvén eigenmode (TAE) were found, where m is the poloidal mode number and n is the toroidal mode number. Four important results were obtained as follows. First, the instability in the CFQS in three-dimensional form was shown for the first time. Second, strong toroidal mode coupling was found for the spatial profiles, and it is consistent with the theoretical prediction. Third, the resonant condition caused by the absence of axial symmetry in CFQS was demonstrated for the first time. The general resonant condition is fmode=Nfφ-Lfθ, where fmode, fφ, and fθ are mode frequency, particle toroidal transit frequency, and particle poloidal transit frequency, respectively; N and L are arbitrary integers, represent toroidal and poloidal resonance numbers. For EPM-like mode, the dominant and subdominant resonant conditions are fmode=3fφ-7fθ and fmode=fφ-fθ, respectively. For TAE, the dominant and subdominant resonant conditions are fmode=4fφ-9fθ and fmode=2fφ-3fθ, respectively. On the one hand, the toroidal resonance numbers are different from the toroidal mode numbers by 2. This indicates that the 2-fold rotational symmetry affects the resonance condition. On the other hand, the subdominant resonances satisfy N=n, which is expected for the axisymmetric plasmas and most of the toroidal plasmas including stellarators. Fourth, the nonlinear frequency chirpings in CFQS were demonstrated for the first time. Hole and clump structures were formed in the pitch angle and energy phase space, and the particles comprising the hole and clump were kept resonant with the modes during the mode frequencies chirping.

We have investigated the waveform distortion of energetic particle driven off-axis fishbone mode (OFM) in tokamak plasmas with kinetic magnetohydrodynamic (MHD) hybrid simulations. We extended our previous simulations (Li et al 2022 Nucl. Fusion 62 026013) by considering higher-n harmonics in the MHD fluid, where n is toroidal mode number. The waveform distortion is successfully reproduced in the simulation for both magnetic fluctuations and temperature fluctuations. It is clarified that the waveform distortion arises from the superposition of the n = 2 harmonics on the fundamental n = 1 harmonics of OFM, where the n = 2 harmonics are generated by the MHD nonlinearity from the n = 1 OFM. Two types of waveform distortion can occur depending on the phase relationship between the n = 1 and n = 2 harmonics and the relative amplitude of the n = 2 harmonics to the n = 1 harmonics. Lissajous curve analyses indicate that the wave couplings between the n = 1 and n = 2 harmonics with phase-lock ∼ π and ∼ 0 lead to “rising distortion” and “falling distortion”, respectively. The two types of waveform distortion can be attributed to the strong shearing profile of radial MHD velocity with n = 2 around the q = 2 magnetic flux surface. The dependence of waveform distortion on viscosity is investigated. It is found that the viscosity which is needed to reproduce the waveform distortion is larger than that in the experiment.