Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userGlobal Theory to Understand Toroidal Drift Waves in Steep Gradient
63
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Toroidal drift waves with unconventional mode structures and non-ground eigenstates, which differ from typical ballooning structure mode, are found to be important recently by large scale global gyrokinetic simulations and especially become dominant at strong gradient edge plasmas [cf., Xie and Xiao, Phys. Plasmas, 22, 090703 (2015)]. The global stability and mode structures of drift wave in this steep edge density and temperature gradients are examined by both direct numerical solutions of a model two-dimensional eigen equation and analytical theory employing WKB-ballooning approach. Theory agrees with numerical solutions quite well. Our results indicate that (i) non-ground eigenstates and unconventional mode structures generally exist and can be roughly described by two parameters `quantum number $l$ and ballooning angle $vartheta_k$, (ii) local model can overestimate the growth rate largely, say, $>50%$, and (iii) the narrow steep equilibrium profile leads to twisting (triangle-like) radial mode structures. With velocity space integral, semi-local theory predicts that the critical jump gradient of the most unstable ion temperature gradient mode from ground state $l=0$ to non-ground state $l=1$ is $L_T^{-1}Rsim50$. These features can have important consequences to turbulent transport.
rate research
Read More
First principle gyrokinetic simulation of the edge turbulent transport in toroidal plasmas finds a reverse trend in the turbulent transport coefficients under strong gradients. It is found that there exist both linear and nonlinear critical gradients for the nonmonotonicity of transport characteristics. The discontinuity of transport flux slope around the turning gradient shows features of second order phase transition. Under strong gradient the most unstable modes are in non-ground eigenstates with unconventional mode structures, which significantly reduces the effective correlation length and thus reverse the transport trend. Our results suggest a completely new mechanism for the low to high confinement mode transition without invoking shear flow or zonal flow.
The gyrokinetic toroidal code (GTC) has been upgraded for global simulations by coupling the core and scrape-off layer (SOL) regions across the separatrix with field-aligned particle-grid interpolations. A fully kinetic particle pusher for high frequency waves (ion cyclotron frequency and beyond) and a guiding center pusher for low frequency waves have been implemented using cylindrical coordinates in a global toroidal geometry. The two integrators correctly capture the particle orbits and agree well with each other, conserving energy and canonical angular momentum. As a verification and application of this new capability, ion guiding center simulations have been carried out to study ion orbit losses at the edge of the DIII-D tokamak for single null magnetic separatrix discharges. The ion loss conditions are examined as a function of the pitch angle for cases without and with a radial electric field. The simulations show good agreement with past theoretical results and with experimentally observed feature in which high energy ions flow out along the ion drift orbits and then hit the divertor plates. A measure of the ion direct orbit loss fraction shows that the loss fraction increases with the ion energy for DIII-D in the initial velocity space. Finally, as a further verification of the capability of the new code, self-consistent simulations of zonal flows in the core region of the DIII-D tokamak were carried out. All DIII-D simulations were perfomed in the absence of turbulence.
This paper presents the current state of the global gyrokinetic code ORB5 as an update of the previous reference [Jolliet et al., Comp. Phys. Commun. 177 409 (2007)]. The ORB5 code solves the electromagnetic Vlasov-Maxwell system of equations using a PIC scheme and also includes collisions and strong flows. The code assumes multiple gyrokinetic ion species at all wavelengths for the polarization density and drift-kinetic electrons. Variants of the physical model can be selected for electrons such as assuming an adiabatic response or a ``hybrid model in which passing electrons are assumed adiabatic and trapped electrons are drift-kinetic. A Fourier filter as well as various control variates and noise reduction techniques enable simulations with good signal-to-noise ratios at a limited numerical cost. They are completed with different momentum and zonal flow-conserving heat sources allowing for temperature-gradient and flux-driven simulations. The code, which runs on both CPUs and GPUs, is well benchmarked against other similar codes and analytical predictions, and shows good scalability up to thousands of nodes.
We show that zonal flow can be preferentially excited by intermediate-scale toroidal electron temperature gradient (ETG) turbulence in tokamak plasmas. Previous theoretical studies that yielded an opposite conclusion assumed a fluid approximation for ETG modes. Here, we carry out a gyrokinetic analysis which ultimately yields a nonlinear Schr{o}dinger equation for the ETG dynamics with a Navier-Stokes type nonlinearity. For typical tokamak parameters, it is found that zonal flow generation plays an important role in the intermediate-scale ETG turbulence. This finding offers an explanation for recent multi-scale gyrokinetic simulations.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
