Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userFlux driven pedestal formation in tokamaks: turbulence simulations validated against the isotope effect
52
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The spontaneous pedestal formation above a power threshold at the edge of magnetically confined plasma is modelled for the first time in flux driven three-dimensional fluid simulations of electromagnetic turbulence with the code EMEDGE3D. The role of the collisional friction between trapped and passing particles is shown to be the key ingredient for shearing the radial electric field, hence stabilizing the turbulence, rather than the Maxwell and Reynolds stresses. The isotope effect, observed in many tokamaks worldwide, is recovered in EMEDGE3D simulations: the power threshold for pedestal formation is lower for Tritium than for Deuterium. The turbulence auto-correlation time is found to increase with the ion mass easing the radial electric shear stabilization, hence the pedestal formation.
rate research
Read More
In some conditions, I-mode plasmas can feature pedestal relaxation events (PREs) that transiently enhance the energy reaching the divertor target plates. To shed light into their appearance, characteristics and energy reaching the divertor targets, a comparative study between two tokamaks $-$ Alcator C-Mod and ASDEX Upgrade $-$ is carried out. It is found that PREs appear only in a subset of I-mode discharges, mainly when the plasma is close to the H-mode transition. Also, the nature of the triggering instability is discussed by comparing measurements close to the separatrix in both devices. The PRE relative energy loss from the confined region increases with decreasing pedestal top collisionality $ u_{mathrm{ped}}^*$. In addition, the relative electron temperature drop at the pedestal top, which is related to the conductive energy loss, rises with decreasing $ u_{mathrm{ped}}^*$. Finally, the peak parallel energy fluence due to the PRE measured on the divertor in both devices is compared to the model introduced in [1] for type-I ELMs. The model is shown to provide an upper boundary for PRE energy fluence data, while a lower boundary is found by dividing the model by three. These two boundaries are used to make projections to future devices such as DEMO and ARC.
Cross-field neoclassical transport of heat, particles and momentum is studied in sharp density pedestals, with a focus on isotope and profile effects, using a radially global approach. Global effects -- which tend to reduce the peak ion heat flux, and shift it outward -- increase with isotope mass for fixed profiles. The heat flux reduction exhibits a saturation with a favorable isotopic trend. A significant part of the heat flux can be convective even in pure plasmas, unlike in the plasma core, and it is sensitive to how momentum sources are distributed between the various species. In particular, if only ion momentum sources are allowed, in global simulations of pure plasmas the ion particle flux remains close to its local value, while this may not be the case for simulations with isotope mixtures or electron momentum sources. The radial angular momentum transport that is a finite orbit width effect, is found to be strongly correlated with heat sources.
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.
A long-lasting debate in space plasma physics concerns the nature of subproton-scale fluctuations in solar wind (SW) turbulence. Over the past decade, a series of theoretical and observational studies were presented in favor of either kinetic Alfven wave (KAW) or whistler turbulence. Here, we investigate numerically the nature of the subproton-scale turbulent cascade for typical SW parameters by means of unprecedented high-resolution simulations of forced hybrid-kinetic turbulence in two real-space and three velocity-space dimensions. Our analysis suggests that small-scale turbulence in this model is dominated by KAWs at $betagtrsim1$ and by magnetosonic/whistler fluctuations at lower $beta$. The spectral properties of the turbulence appear to be in good agreement with theoretical predictions. A tentative interpretation of this result in terms of relative changes in the damping rates of the different waves is also presented. Overall, the results raise interesting new questions about the properties and variability of subproton-scale turbulence in the SW, including its possible dependence on the plasma $beta$, and call for detailed and extensive parametric explorations of driven kinetic turbulence in three dimensions.
We are developing a new continuum gyrokinetic code, Gkeyll, for use in edge plasma simulations, and here present initial simulations of turbulence on open field lines with model sheath boundary conditions. The code implements an energy conserving discontinuous Galerkin scheme, applicable to a general class of Hamiltonian equations. Several applications to test problems have been done, including a calculation of the parallel heat-flux on divertor plates resulting from an ELM crash in JET, for a 1x/1v SOL scenario explored previously, where the ELM is modeled as a time-dependent intense upstream source. Here we present initial simulations of turbulence on open field lines in the LAPD linear plasma device. We have also done simulations in a helical open-field-line geometry. While various simplifications have been made at present, this still includes some of the key physics of SOL turbulence, such as bad-curvature drive for instabilities and rapid parallel losses with sheath boundary conditions. This is useful for demonstrating the overall feasibility of this approach and for initial physics studies of SOL turbulence. We developed a novel version of DG that uses Maxwellian-weighted basis functions while still preserving exact particle and energy conservation. The Maxwellian-weighted DG method achieves the same error with 4 times less computational cost in 1v, or 16 times lower cost in the 2 velocity dimensions of gyrokinetics (assuming memory bandwidth is the limiting factor).
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
