Search

Kinetic-magnetohydrodynamic hybrid simulations were performed to investigate the linear growth and the nonlinear evolution of off-axis fishbone mode (OFM) destabilized by trapped energetic ions in tokamak plasmas. The spatial profile of OFM is mainly composed of m/n = 2/1 mode inside the q = 2 magnetic flux surface while the m/n = 3/1 mode is predominant outside the q = 2 surface, where m and n are the poloidal and toroidal mode numbers, respectively, and q is the safety factor. The spatial profile of the OFM is a strongly shearing shape on the poloidal plane, suggesting the nonperturbative effect of the interaction with energetic ions. The frequency of the OFM in the linear growth phase is in good agreement with the precession drift frequency of trapped energetic ions, and the frequency chirps down in the nonlinear phase. Two types of resonance conditions between trapped energetic ions and OFM are found. For the first type of resonance, the precession drift frequency matches the OFM frequency, while for the second type, the sum of the precession drift frequency and the bounce frequency matches the OFM frequency. The first type of resonance is the primary resonance for the destabilization of OFM. The resonance frequency which is defined based on precession drift frequency and bounce frequency of the nonlinear orbit for each resonant particle is analyzed to understand the frequency chirping. The resonance frequency of the particles that transfer energy to the OFM chirps down, which may result in the chirping down of the OFM frequency. A detailed analysis of the energetic ion distribution function in phase space shows that the gradient of the distribution function along the E′ = const. line drives or stabilizes the instability, where E′ is a combination of energy and toroidal canonical momentum and conserved during the wave–particle interaction. The distribution function is flattened along the E′ = const. line in the nonlinear phase leading to the saturation of the instability.

Energetic particle transport by a magnetohydrodynamic (MHD) instability driven by helically trapped energetic particles is studied for a high-performance Large Helical Device plasma with kinetic-MHD hybrid simulations. It is observed in the simulation that an MHD mode with poloidal/toroidal mode numbers ????/????=2/1driven by helically trapped energetic particles causes a significant redistribution of perpendicular energetic particle pressure profile. The frequency of the MHD mode decreases rapidly at the saturation of the instability and changes sign, which indicates a reversal of the mode propagation direction. It is found that the helically trapped energetic particles interacting strongly with the MHD mode change the precession drift direction at the same time as the reversal of the MHD mode propagation direction. The helically trapped energetic particles with the precession drift reversal are transported rapidly in the radially outward direction before the original precession drift direction is recovered. The precession drift reversal and the outward transport are caused by interaction with the electric field of the MHD mode. The vast majority of trapped energetic particles which interact strongly with the MHD mode experience precession drift reversal, leading to a significant redistribution of the perpendicular energetic particle pressure profile.Numerical computations were performed on the Plasma Simulator (NEC SX-Aurora TSUBASA) of NIFS with the support and under the auspices of the NIFS Collaboration Research program (No. NIFS20KNST169), the JFRS-1 of the International Fusion Energy Research Centre, and the Supercomputer Fugaku of the RIKEN Center for Computational Science (Project ID: hp210178). This work was supported by MEXT as Program for Promoting Research studies on the Supercomputer Fugaku (Exploration of burning plasma confinement physics).

A new kinetic-magnetohydrodynamic hybrid simulation model where the gyrokinetic particle-in-cell simulation is applied to both thermal ions and energetic particles is presented. Toroidal Alfvén eigenmodes (TAEs) destabilized by energetic ions in tokamak plasmas are simulated with the new simulation model. Energy channeling from energetic ions to thermal ions through Alfvén eigenmodes (AEs) is demonstrated by the simulation. The distribution function fluctuations and the resonance condition are analyzed for both thermal ions and energetic ions. The strong energy transfer between the particles and the AE and the strong particle transport occur when the following conditions are satisfied at the resonance location in phase space: (1) the poloidal resonance number is close to the poloidal mode number of the AE, (2) the AE has a substantial amplitude, (3) the distribution function has a substantial gradient along the E ′ = const. line, where E ′ is a conserved variable for the wave-particle interaction. While the distribution function fluctuations for energetic ions are consistent with the resonance condition with the TAEs, the distribution function fluctuations for thermal ions do not satisfy the resonance condition when the bulk plasma beta is 1%. This indicates that the resonance does not play an important role in the interaction between thermal ions and the TAE for the relatively low bulk plasma temperature. On the other hand, when the bulk plasma beta is 4%, the resonance between thermal ions and the TAEs become important leading to Landau damping.

Energetic electron effects on an energetic-ion driven toroidal Alfvén eigenmode (TAE) are investigated via hybrid simulations of an MHD fluid interacting with energetic particles. Both energetic electrons and energetic ions described by drift-kinetic equations are included in the present work. It is found that the TAE can be effectively stabilized by off-axis peaked energetic electrons which are located near the mode center, while the centrally peaked energetic electrons fail to stabilize the mode. It is confirmed that the spatially localized pressure profile of energetic electrons causes the stabilization of TAE. The stabilized TAE has a more localized mode structure accompanied by a significant reduction in the energetic ion driving rate. The small change of mode frequency and dissipation rate indicate the stabilization mechanism is different from the so-called pressure gradient stabilization that drives the TAE into continuum. The results suggest that the strong plasma non-uniformity induced by the energetic electron beta profile may be responsible for the change of mode structure. It is also found that this stabilizing effect is more effective for a high-n TAE. Moreover, it is numerically verified that the positive (negative) pressure gradient at the TAE center will increase (decrease) the mode frequency. The wave-particle interactions are also analysed for a case with energetic electrons peaked at the inner side of the TAE center. It is found that the power transfer to a resonant barely trapped energetic electron, which taps energy from the wave, can be comparable to the power transfer from a resonant energetic ion. This suggests that if a sufficient number of resonant barely trapped electrons are present, they might stabilize energetic-ion driven TAE through the wave-particle interaction.

In gyrokinetic simulations of turbulent impurity transport, trace impurity species are often treated as passive species, in the sense that they are not included in Maxwell equations. This is consistent with the assumption that impurities with low enough concentrations are impacted by turbulence generated by electrons and main ions, but do not impact it significantly in return. In this work, we relax this assumption, and investigate the active impacts of impurity on impurity transport as a function of its concentration, in the presence of trapped-particle-driven turbulence. We focus on $W^{40+}$ tungsten, which is relevant for modern tokamaks, and adopt a reduced gyrokinetic bounce-averaged model for trapped particles in a simplified tokamak geometry. The impacts depend on the relationship between equilibrium density gradient and temperature gradient. When these gradients are equal, we observe that tungsten can be treated as a passive species for concentrations below $2\times 10^{-4}$. Above this concentration, the impurity significantly impacts both density and heat transport, essentially quenching them for concentrations above $10^{-3}$. This quenching occurs as electric potential fluctuations become in phase with impurity density fluctuations.