The free energy balance equation applied to gyrokinetic instabilities, the effect of the charge flux constraint, and application to simplified kinetic models

The free energy balance equation for gyrokinetic fluctuations is derived and applied to instabilities. An additional term due to electromagnetic sources is included. This can provide a simpler way to compute the free energy balance in practical applications, and is also conceptually clarifying. The free energy balance, by itself, is not sufficient to determine an eigenfrequency. The preceding results are derived in general geometry. The charge flux constraint in gyrokinetics can provide a necessary additional relation, and the combination of these two can be equivalent to a dispersion relation. The charge flux constraint can prevent the appearance of an unstable eigenmode even though the free energy balance would allow strongly growing fluctuations. The application of these concepts to simplified kinetic models in simplified geometry is also indicated.Comment: 5 page

Geometrical Properties of a "Snow-Flake" Divertor

Using a simple set of poloidal field coils, one can reach the situation where the null of the poloidal magnetic field in the divertor region is of a second order, not of the first order as in the usual X-point divertor. Then, the separatrix in the vicinity of the null-point splits the poloidal plane not into four sectors, but into six sectors, making the whole structure looking like a snow-flake (whence a name, [1]). This arrangement allows one to spread the heat load over much broader area than in the case of a standard divertor. A disadvantage of this configuration is in that it is topologically unstable, and, with the current in the plasma varying with time, it would switch either to the standard X-point mode, or to the mode with two X-points close to each other. To avoid this problem, it is suggested to have a current in the divertor coils by roughly 5% higher than in an 'optimum' regime (the one where a snow-flake separatrix is formed). In this mode, the configuration becomes stable and can be controlled by varying the current in the divertor coils in concert with the plasma current; on the other hand, a strong flaring of the scrape-off layer still remains in force. Geometrical properties of this configurations are analyzed for a simple model. Potential advantages and disadvantages of this scheme are discussed

Multi-Scale Interactions of Microtearing Turbulence in the Tokamak Pedestal

Microtearing turbulence in an idealized pedestal scenario is found to saturate via zonal fields, while also exciting strong zonal flows; a concurrent upshift of the non-linear critical gradient is observed. The zonal flows cause electron-temperature-gradient-driven turbulence to be ameliorated. When applying resonant magnetic perturbations, the prompt charge loss off the flux surface erodes the zonal flow, leading to higher electron-scale fluxes, while leaving microtearing saturation physics unaffected.</p

Gyrokinetic analysis and simulation of pedestals, to identify the culprits for energy losses using fingerprints

Fusion performance in tokamaks hinges critically on the efficacy of the Edge Transport Barrier (ETB) at suppressing energy losses. The new concept of fingerprints is introduced to identify the instabilities that cause the transport losses in the ETB of many of today's experiments, from widely posited candidates. Analysis of the Gyrokinetic-Maxwell equations, and gyrokinetic simulations of experiments, find that each mode type produces characteristic ratios of transport in the various channels: density, heat and impurities. This, together with experimental observations of transport in some channel, or, of the relative size of the driving sources of channels, can identify or determine the dominant modes causing energy transport. In multiple ELMy H-mode cases that are examined, these fingerprints indicate that MHD-like modes are apparently not the dominant agent of energy transport; rather, this role is played by Micro-Tearing Modes (MTM) and Electron Temperature Gradient (ETG) modes, and in addition, possibly Ion Temperature Gradient (ITG)/Trapped Electron Modes (ITG/TEM) on JET. MHD-like modes may dominate the electron particle losses. Fluctuation frequency can also be an important means of identification, and is often closely related to the transport fingerprint. The analytical arguments unify and explain previously disparate experimental observations on multiple devices, including DIII-D, JET and ASDEX-U, and detailed simulations of two DIII-D ETBs also demonstrate and corroborate this

Negative-Energy Perturbations in Circularly Cylindrical Equilibria within the Framework of Maxwell-Drift Kinetic Theory

The conditions for the existence of negative-energy perturbations (which could be nonlinearly unstable and cause anomalous transport) are investigated in the framework of linearized collisionless Maxwell-drift kinetic theory for the case of equilibria of magnetically confined, circularly cylindrical plasmas and vanishing initial field perturbations. For wave vectors with a non-vanishing component parallel to the magnetic field, the plane equilibrium conditions (derived by Throumoulopoulos and Pfirsch [Phys Rev. E {\bf 49}, 3290 (1994)]) are shown to remain valid, while the condition for perpendicular perturbations (which are found to be the most important modes) is modified. Consequently, besides the tokamak equilibrium regime in which the existence of negative-energy perturbations is related to the threshold value of 2/3 of the quantity $\eta_\nu = \frac {\partial \ln T_\nu} {\partial \ln N_\nu}$, a new regime appears, not present in plane equilibria, in which negative-energy perturbations exist for {\em any} value of $\eta_\nu$. For various analytic cold-ion tokamak equilibria a substantial fraction of thermal electrons are associated with negative-energy perturbations (active particles). In particular, for linearly stable equilibria of a paramagnetic plasma with flat electron temperature profile ($\eta_e=0$), the entire velocity space is occupied by active electrons. The part of the velocity space occupied by active particles increases from the center to the plasma edge and is larger in a paramagnetic plasma than in a diamagnetic plasma with the same pressure profile. It is also shown that, unlike in plane equilibria, negative-energy perturbations exist in force-free reversed-field pinch equilibria with a substantial fraction of active particles.Comment: 31 pages, late

Gyrokinetic Simulations Compared with Magnetic Fluctuations Diagnosed with a Faraday-Effect Radial Interferometer-Polarimeter in the DIII-D pedestal

Experimental data on electromagnetic fluctuations in DIII-D, made available by the Faraday-effect Radial Interferometer-Polarimeter (RIP) diagnostic, is examined in comparison with detailed gyrokinetic simulations using Gyrokinetic Electromagnetic Numerical Experiment (GENE). The diagnostic has the unique capability of making internal measurements of fluctuating magnetic fields $\frac{\int n_e \delta B_r dR}{\int n_e dR}$. Local linear simulations identify microtearing modes (MTMs) over a substantial range of toroidal mode numbers (peaking at $n=15$) with frequencies in good agreement with the experimental data. Local nonlinear simulations reinforce this result by producing a magnetic frequency spectrum in good agreement with that diagnosed by RIP. Simulated heat fluxes are in the range of experimental expectations. However, magnetic fluctuation amplitudes are substantially lower than the experimental expectations. Possible sources of this discrepancy are discussed, notably the fact that the diagnostics are localized at the mid-plane -- the poloidal location where the simulations predict the fluctuation amplitudes to be smallest. Despite some discrepancies, several connections between simulations and experiments, combined with general criteria discriminating between potential pedestal instabilities, strongly point to MTMs as the source of the observed magnetic fluctuations

Modeling electron temperature profiles in the pedestal with simple formulas for ETG transport

This paper reports on the refinement (building on Ref.~\cite{hatch_22}) and application of simple formulas for electron heat transport from electron temperature gradient (ETG) driven turbulence in the pedestal. The formulas are improved by (1) improving the parameterization for certain key parameters and (2) carefully accounting for the impact of geometry and shaping in the underlying gyrokinetic simulation database. Comparisons with nonlinear gyrokinetic simulations of ETG transport in the MAST pedestal demonstrate the model's applicability to spherical tokamaks in addition to standard aspect ratio tokamaks. We identify bounds for model applicability: the model is accurate in the steep gradient region, where the ETG turbulence is largely slab-like, but accuracy decreases as the temperature gradient becomes weaker in the pedestal top and the instabilities become increasingly toroidal in nature. We use the formula to model the electron temperature profile in the pedestal for four experimental scenarios while extensively varying input parameters to represent uncertainties. In all cases, the predicted electron temperature pedestal exhibits extreme sensitivity to separatrix temperature and density, which has implications for core-edge integration. The model reproduces the electron temperature profile for high $\eta_e = L_{ne}/L_{Te}$ scenarios but not for low $\eta_e$ scenarios in which microtearing modes have been identified. We develop a proof-of-concept model for MTM transport and explore the relative roles of ETG and MTM in setting the electron temperature profile. We propose that pedestal scenarios predicted for future devices should be tested for compatibility with ETG transport