No Arabic abstract
Long range frequency chirping of Bernstein-Greene-Kruskal modes, whose existence is determined by the fast particles, is investigated in cases where these particles do not move freely and their motion is bounded to restricted orbits. An equilibrium oscillating potential, which creates different orbit topologies of energetic particles, is included into the bump-on-tail instability problem of a plasma wave. With respect to fast particles dynamics, the extended model captures the range of particles motion (trapped/passing) with energy and thus represents a more realistic 1D picture of the long range sweeping events observed for weakly damped modes, e.g. global Alfven eigenmodes, in tokamaks. The Poisson equation is solved numerically along with bounce averaging the Vlasov equation in the adiabatic regime. We demonstrate that the shape and the saturation amplitude of the nonlinear mode structure depends not only on the amount of deviation from the initial eigenfrequency but also on the initial energy of the resonant electrons in the equilibrium potential. Similarly, the results reveal that the resonant electrons following different equilibrium orbits in the electrostatic potential lead to different rates of frequency evolution. As compared to the previous model [Breizman B.N. 2010 Nucl. Fusion 50 084014], it is shown that the frequency sweeps with lower rates. The additional physics included in the model enables a more complete 1D description of the range of phenomena observed in experiments.
In recent years, a strong reduction of plasma turbulence in the presence of energetic particles has been reported in a number of magnetic confinement experiments and corresponding gyrokinetic simulations. While highly relevant to performance predictions for burning plasmas, an explanation for this primarily nonlinear effect has remained elusive so far. A thorough analysis finds that linearly marginally stable energetic particle driven modes are excited nonlinearly, depleting the energy content of the turbulence and acting as an additional catalyst for energy transfer to zonal modes (the dominant turbulence saturation channel). Respective signatures are found in a number of simulations for different JET and ASDEX Upgrade discharges with reduced transport levels attributed to energetic ion effects.
The nonlinear dynamics of energetic-particle (EP) driven geodesic acoustic modes (EGAM) is investigated here. A numerical analysis with the global gyrokinetic particle-in-cell code ORB5 is performed, and the results are interpreted with the analytical theory, in close comparison with the theory of the beam-plasma instability. Only axisymmetric modes are considered, with a nonlinear dynamics determined by wave-particle interaction. Quadratic scalings of the saturated electric field with respect to the linear growth rate are found for the case of interest. The EP bounce frequency is calculated as a function of the EGAM frequency, and shown not to depend on the value of the bulk temperature. Near the saturation, we observe a transition from adiabatic to non-adiabatic dynamics, i.e., the frequency chirping rate becomes comparable to the resonant EP bounce frequency. The numerical analysis is performed here with electrostatic simulations with circular flux surfaces, and kinetic effects of the electrons are neglected.
Turbulence in tokamaks generates radially sheared zonal flows. Their oscillatory counterparts, geodesic acoustic modes (GAMs), appear due to the action of the magnetic field curvature. The GAMs can be driven unstable by an anisotropic energetic particle (EP) population leading to the formation of global radial structures, called EGAMs. The EGAMs can redistribute EP energy to the bulk plasma through collisionless wave-particle interaction. In such a way, the EGAMs might contribute to the plasma heating. Thus, investigation of EGAM properties, especially in the velocity space, is necessary for precise understanding of the transport phenomena in tokamak plasmas. In this work, the nonlinear dynamics of EGAMs without considering the mode interaction with the turbulence is investigated with the help of a Mode-Particle-Resonance (MPR) diagnostic implemented in the global gyrokinetic particle-in-cell code ORB5. An ASDEX Upgrade discharge is chosen as a reference case for this investigation due to its rich EP nonlinear dynamics. An experimentally relevant magnetic field configuration, thermal species profiles and an EP density profile are taken for EGAM chirping modelling and its comparison with available empirical data. The same magnetic configuration is used to explore energy transfer by the mode from the energetic particles to the thermal plasma including kinetic electron effects. For a given EGAM level the plasma heating by the mode can be significantly enhanced by varying the EP parameters. Electron dynamics decreases the EGAM saturation amplitude and consequently reduces the plasma heating, even though the mode transfers its energy to thermal ions much more than to electrons.
This paper presents a study of the interaction between Alfven modes and zonal structures, considering a realistic ASDEX Upgrade equilibrium. The results of gyrokinetic simulations with the global, electromagnetic, particle-in-cell code ORB5 are presented, where the modes are driven unstable by energetic particles with a bump-on-tail equilibrium distribution function, with radial density gradient. Two regimes have been observed: at low energetic particles concentration, the Alfven mode saturates at much higher level in presence of zonal structures; on the other hand at high energetic particles concentration the difference is less pronounced. The former regime is characterized by the zonal structure (identified as an energetic particle driven geodesic acoustic mode), being more unstable than the Alfven mode. In the latter regime the Alfven mode is more unstable than the zonal structure. The theoretical explanation is given in terms of a 3-wave coupling of the energetic particle driven geodesic acoustic mode and Alfven mode, mediated by the curvature-pressure coupling term of the energetic particles.
The aim of the present study is to analyze the stability of the pressure gradient driven modes (PM) and Alfven eigenmodes (AE) in the Large Helical Device (LHD) plasma if the rotational transform profile is modified by the current drive of the tangential neutral beam injectors (NBI). This study forms a basic search for optimized operation scenarios with reduced mode activity. The analysis is performed using the code FAR3d which solves the reduced MHD equations describing the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations for density and parallel velocity moments of the energetic particle (EP) species, including the effect of the acoustic modes. The Landau damping and resonant destabilization effects are added via the closure relation. On-axis and off-axis NBI current drive modifies the rotational transform which becomes strongly distorted as the intensity of the neutral beam current drive (NBCD) increases, leading to wider continuum gaps and modifying the magnetic shear. The simulations with on-axis NBI injection show that a counter (ctr-) NBCD in inward shifted and default configurations leads to a lower growth rate of the PM, although strong n=1 and 2 AEs can be destabilized. For the outward shifted configurations, a co-NBCD improves the AEs stability but the PM are further destabilized if the co-NBCD intensity is 30 kA/T. If the NBI injection is off-axis, the plasma stability is not significantly improved due to the further destabilization of the AE and energetic particle modes (EPM) in the middle and outer plasma region.