No Arabic abstract
The existence of kinetic ballooning mode (KBM) high order (non-ground) eigenstates for tokamak plasmas with steep gradient is demonstrated via gyrokinetic electromagnetic eigenvalue solutions, which reveals that eigenmode parity transition is an intrinsic property of electromagnetic plasmas. The eigenstates with quantum number $l=0$ for ground state and $l=1,2,3ldots$ for non-ground states are found to coexist and the most unstable one can be the high order states ($l eq0$). The conventional KBM is the $l=0$ state. It is shown that the $l=1$ KBM has the same mode structure parity as the micro-tearing mode (MTM). In contrast to the MTM, the $l=1$ KBM can be driven by pressure gradient even without collisions and electron temperature gradient. The relevance between various eigenstates of KBM under steep gradient and edge plasma physics is discussed.
A unified ballooning theory, constructed on the basis of two special theories [Y. Z. Zhang, S. M. Mahajan, X. D. Zhang, Phys. Fluids B4, 2729 (1992); Y. Z. Zhang, T. Xie, Nucl. Fusion & Plasma Phys. 33, 193 (2013)], shows that a weak up-down asymmetric mode structure is normally formed in an up-down symmetric equilibrium; the weak up-down asymmetry in mode structure is the manifestation of non-trivial higher order effects beyond the standard ballooning equation. It is shown that the asymmetric mode may have even higher growth rate than symmetric modes. Salient features of the theory are illustrated by investigating a fluid model for the ion temperature gradient (ITG) mode. The two dimensional (2D) analytical form of ITG mode, solved in ballooning representation, is then converted into the radial-poloidal space to provide the natural boundary condition for solving the 2D mathematical local eigenmode problem. We find the analytical expression of mode structure in good agreement with finite difference solution. This sets a reliable framework for quasi-linear computation.
We have performed fully-kinetic simulations of X-B and O-X-B mode conversion in one and two dimensional setups using the PIC code EPOCH. We have recovered the linear dispersion relation for electron Bernstein waves by employing relatively low amplitude incoming waves. The setups presented here can be used to study non-linear regimes of X-B and O-X-B mode conversion.
In the present work the zonal flow (ZF) growth rate in toroidal ion-temperature-gradient (ITG) mode turbulence including the effects of elongation is studied analytically. The scaling of the ZF growth with plasma parameters is examined for typical tokamak parameter values. The physical model used for the toroidal ITG driven mode is based on the ion continuity and ion temperature equations whereas the ZF evolution is described by the vorticity equation. The results indicate that a large ZF growth is found close to marginal stability and for peaked density profiles and these effects may be enhanced by elongation.
The turbulent transport of main ion and trace impurities in a tokamak device in the presence of steep electron density gradients has been studied. The parameters are chosen for trapped electron (TE) mode turbulence, driven primarily by steep electron density gradients relevant to H-mode physics, but with a transition to temperature gradient driven turbulence as the density gradient flattens. Results obtained through non-linear (NL) and quasilinear (QL) gyrokinetic simulations using the GENE code are compared with results obtained from a fluid model. Main ion and impurity transport is studied by examining the balance of convective and diffusive transport, as quantified by the density gradient corresponding to zero particle flux (peaking factor). Scalings are obtained for the impurity peaking with the background electron density gradient and the impurity charge number. It is shown that the impurity peaking factor is weakly dependent on impurity charge and significantly smaller than the driving electron density gradient.
Spherical tokamak plasmas are typically overdense and thus inaccessible to externally-injected microwaves in the electron cyclotron range. The electrostatic electron Bernstein wave (EBW), however, provides a method to access the plasma core for heating and diagnostic purposes. Understanding the details of the coupling process to electromagnetic waves is thus important both for the interpretation of microwave diagnostic data and for assessing the feasibility of EBW heating and current drive. While the coupling is reasonably well-understood in the linear regime, nonlinear physics arising from high input power has not been previously quantified. To tackle this problem, we have performed one- and two-dimensional fully kinetic particle-in-cell simulations of the two possible coupling mechanisms, namely X-B and O-X-B mode conversion. We find that the ion dynamics has a profound effect on the field structure in the nonlinear regime, as high amplitude short-scale oscillations of the longitudinal electric field are excited in the region below the high-density cut-off prior to the arrival of the EBW. We identify this effect as the instability of the X wave with respect to resonant scattering into an EBW and a lower-hybrid wave. We calculate the instability rate analytically and find this basic theory to be in reasonable agreement with our simulation results.