No Arabic abstract
Self-consistent simulations of neoclassical and electrostatic turbulent transport in a DIII-D H-mode edge plasma under resonant magnetic perturbations (RMPs) have been performed using the global total-f gyrokinetic particle-in-cell code XGC, in order to study density-pump out and electron heat confinement. The RMP field is imported from the extended magneto-hydrodynamics (MHD) code M3D-C1, taking into account the linear two-fluid plasma response. With both neoclassical and turbulence physics considered together, the XGC simulation reproduces two key features of experimentally observed edge transport under RMPs: increased radial particle transport in the pedestal region that is sufficient to account for the experimental pump-out rate, and suppression of the electron heat flux in the steepest part of the edge pedestal. In the simulation, the density fluctuation amplitude of modes moving in the electron diamagnetic direction increases due to interaction with RMPs in the pedestal shoulder and outward, while the electron temperature fluctuation amplitude decreases.
Type-I Edge Localised Modes (ELMs) have been mitigated in MAST through the application of n = 3, 4 and 6 resonant magnetic perturbations (RMPs). For each toroidal mode number of the non-axisymmetric applied fields, the frequency of the ELMs has been increased significantly, and the peak heat flux on the divertor plates reduced commensurately. This increase in ELM frequency occurs despite a significant drop in the edge pressure gradient, which would be expected to stabilise the peeling-ballooning modes thought to be responsible for type-I ELMs. Various mechanisms which could cause a destabilisation of the peeling-ballooning modes are presented, including pedestal widening, plasma rotation braking, three dimensional corrugation of the plasma boundary and the existence of radially extended lobe structures near to the X-point. This leads to a model aimed at resolving the apparent dichotomy of ELM control, that is to say ELM suppression occurring due to the pedestal pressure reduction below the peeling-ballooning stability boundary, whilst the reduction in pressure can also lead to ELM mitigation, which is ostensibly a destabilisation of peeling-ballooning modes. In the case of ELM mitigation, the pedestal broadening, 3d corrugation or lobes near the X-point degrade ballooning stability so much that the pedestal recovers rapidly to cross the new stability boundary at lower pressure more frequently, whilst in the case of suppression, the plasma parameters are such that the particle transport reduces the edge pressure below the stability boundary which is only mildly affected by negligible rotation braking, small edge corrugation or short, broad lobe structures.
Global electromagnetic gyrokinetic simulations show the existence of near threshold conditions for both a high-$n$ kinetic ballooning mode (KBM) and an intermediate-$n$ kinetic version of peeling-ballooning mode (KPBM) in the edge pedestal of two DIII-D H-mode discharges. When the magnetic shear is reduced in a narrow region of steep pressure gradient, the KPBM is significantly stabilized, while the KBM is weakly destabilized and hence becomes the most-unstable mode. Collisions decrease the KBMs critical $beta$ and increase the growth rate.
Sustained ELM mitigation has been achieved using RMPs with a toroidal mode number of n=4 and n=6 in lower single null and with n=3 in connected double null plasmas on MAST. The ELM frequency increases by up to a factor of eight with a similar reduction in ELM energy loss. A threshold current for ELM mitigation is observed above which the ELM frequency increases approximately linearly with current in the coils. A comparison of the filament structures observed during the ELMs in the natural and mitigated stages shows that the mitigated ELMs have the characteristics of type I ELMs even though their frequency is higher, their energy loss is reduced and the pedestal pressure gradient is decreased. During the ELM mitigated stage clear lobe structures are observed in visible-light imaging of the X-point region. The size of these lobes is correlated with the increase in ELM frequency observed. The RMPs produce a clear 3D distortion to the plasma and it is likely that these distortions explain why ELMs are destabilised and hence why ELM mitigation occurs.
Understanding the multi-scale neoclassical and turbulence physics in the edge region (pedestal + scrape-off layer) is required in order to reliably predict performance in future fusion devices. We explore turbulent characteristics in the edge region from a multiscale neoclassical and turbulent XGC1 gyrokinetic simulation in a DIII-D like tokamak geometry, here excluding neutrals and collisions. For an H-mode type plasma with steep pedestal, it is found that the electron density fluctuations increase towards the separatrix, and stay high well into the SOL, reaching a maximum value of $delta n_e / bar{n}_e sim 0.18$. Blobs are observed, born around the magnetic separatrix surface and propagate radially outward with velocities generally less than 1 km/s. Strong poloidal motion of the blobs is also present, near 20 km/s, consistent with $E times B$ rotation. The electron density fluctuations show a negative skewness in the closed field line pedestal regions, consistent with the presence of holes, followed by a transition to strong positive skewness across the separatrix and into the SOL. These simulations indicate that not only neoclassical phenomena, but also turbulence, including the blob-generation mechanism, can remain important in the steep H-mode pedestal and SOL. Qualitative comparisons will be made to experimental observations.
The quasilinear particle flux arising from gyrokinetic instabilities is calculated in the electrostatic and collisionless approximation, keeping the geometry of the magnetic field arbitrary. In particular, the flux of electrons and heavy impurity ions is studied in the limit where the former move quickly, and the latter slowly, along the field compared with the mode frequency. Conclusions are drawn about how the particle fluxes of these species depend on the magnetic-field geometry, mode structure and frequency of the instability. Under some conditions, such as everywhere favourable or unfavourable magnetic curvature and modest temperature gradients, it is possible to make general statements about the fluxes independently of the details of the instability. In quasi-isodynamic stellarators with favourable bounce-averaged curvature for most particles, the particle flux is always outward if the temperature gradient is not too large, suggesting that it might be difficult to fuel such devices with gas puffing from the wall. In devices with predominantly unfavourable magnetic curvature, the particle flux can be inward, resulting in spontaneous density peaking in the centre of the plasma. In the limit of highly charged impurities, ordinary diffusion (proportional to the density gradient) dominates over other transport channels and the diffusion coefficient becomes independent of mass and charge. An estimate for the level of transport caused by magnetic-field fluctuations arising from ion-temperature-gradient instabilities is also given and is shown to be small compared with the electrostatic component.