No Arabic abstract
The gyrokinetic turbulence code GS2 was used to investigate the effects of plasma beta on linear, collisionless ion temperature gradient (ITG) modes and trapped electron modes (TEM) in National Compact Stellarator Experiment (NCSX) geometry. Plasma beta affects stability in two ways: through the equilibrium and through magnetic fluctuations. The first was studied here by comparing ITG and TEM stability in two NCSX equilibria of differing beta values, revealing that the high beta equilibrium was marginally more stable than the low beta equilibrium in the adiabatic-electron ITG mode case. However, the high beta case had a lower kinetic-electron ITG mode critical gradient. Electrostatic and electromagnetic ITG and TEM mode growth rate dependencies on temperature gradient and density gradient were qualitatively similar. The second beta effect is demonstrated via electromagnetic ITG growth rates dependency on GS2s beta input parameter. A linear benchmark with gyrokinetic codes GENE and GKV-X is also presented.
An analytic equilibrium, the Toroidal Bessel Function Model, is used in conjunction with the gyrokinetic code GYRO to investigate the nature of microinstabilities in a reversed field pinch (RFP) plasma. The effect of the normalized electron plasma pressure ({beta}) on the characteristics of the microinstabilities is studied. A transition between an ion temperature gradient (ITG) driven mode and a microtearing mode as the dominant instability is found to occur at a {beta} value of approximately 4.5%. Suppression of the ITG mode occurs as in the tokamak, through coupling to shear Alfven waves, with a critical {beta} for stability higher than its tokamak equivalent due to a shorter parallel connection length. There is a steep dependence of the microtearing growth rate on temperature gradient suggesting high profile stiffness. There is evidence for a collisionless microtearing mode. The properties of this mode are investigated, and it is found that curvature drift plays an important role in the instability.
The linear gyrokinetic stability properties of magnetically confined electron-positron plasmas are investigated in the parameter regime most likely to be relevant for the first laboratory experiments involving such plasmas, where the density is small enough that collisions can be ignored and the Debye length substantially exceeds the gyroradius. Although the plasma beta is very small, electromagnetic effects are retained, but magnetic compressibility can be neglected. The work of a previous publication (Helander, 2014) is thus extended to include electromagnetic instabilities, which are of importance in closed-field-line configurations, where such instabilities can occur at arbitrarily low pressure. It is found that gyrokinetic instabilities are completely absent if the magnetic field is homogeneous: any instability must involve magnetic curvature or shear. Furthermore, in dipole magnetic fields, the stability threshold for interchange modes with wavelengths exceeding the Debye radius coincides with that in ideal MHD. Above this threshold, the quasilinear particle flux is directed inward if the temperature gradient is sufficiently large, leading to spontaneous peaking of the density profile.
The Dupree-Weinstock renormalization is used to prove that a reactive closure exists for drift wave turbulence in magnetized plasmas. The result is used to explain recent results in gyrokinetic simulations and is also related to the Mattor-Parker closure. The level of closure is found in terms of applied external sources.
The plasma beta, ratio of kinetic to magnetic pressure, inside a tokamak should stay below the Troyon limit to avoid major plasma instabilities. However, this paper argues that Troyon limit occurs only when current profiles cannot sustain high beta equilibria. This paper shows that plasmas with peak beta on the order of unity are stable to ideal MHD modes. While the theoretical (ideal-MHD) stability of unity beta equilibrium has been resolved, the experimental existence of high beta equilibria in tokamaks remains an unsolved problem of plasma physics.
To help reveal the complete picture of linear kinetic drift modes, four independent numerical approaches, based on integral equation, Euler initial value simulation, Euler matrix eigenvalue solution and Lagrangian particle simulation, respectively, are used to solve the linear gyrokinetic electrostatic drift modes equation in Z-pinch with slab simplification and in tokamak with ballooning space coordinate. We identify that these approaches can yield the same solution with the difference smaller than 1%, and the discrepancies mainly come from the numerical convergence, which is the first detailed benchmark of four independent numerical approaches for gyrokinetic linear drift modes. Using these approaches, we find that the entropy mode and interchange mode are on the same branch in Z-pinch, and the entropy mode can have both electron and ion branches. And, at strong gradient, more than one eigenstate of the ion temperature gradient mode (ITG) can be unstable and the most unstable one can be on non-ground eigenstates. The propagation of ITGs from ion to electron diamagnetic direction at strong gradient is also observed, which implies that the propagation direction is not a decisive criterion for the experimental diagnosis of turbulent mode at the edge plasmas.