Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userGyrokinetic modelling of stationary electron and impurity profiles in tokamaks
386
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Particle transport due to Ion Temperature Gradient/Trapped Electron (ITG/TE) mode turbulence is investigated using the gyrokinetic code GENE. Both a reduced quasilinear (QL) treatment and nonlinear (NL) simulations are performed for typical tokamak parameters corresponding to ITG dominated turbulence. A selfconsistent treatment is used, where the stationary local profiles are calculated corresponding to zero particle flux simultaneously for electrons and trace impurities. The scaling of the stationary profiles with magnetic shear, safety factor, electron-to-ion temperature ratio, collisionality, toroidal sheared rotation, triangularity, and elongation is investigated. In addition, the effect of different main ion mass on the zero flux condition is discussed. The electron density gradient can significantly affect the stationary impurity profile scaling. It is therefore expected, that a selfconsistent treatment will yield results more comparable to experimental results for parameter scans where the stationary background density profile is sensitive. This is shown to be the case in scans over magnetic shear, collisionality, elongation, and temperature ratio, for which the simultaneous zero flux electron and impurity profiles are calculated. A slight asymmetry between hydrogen, deuterium and tritium with respect to profile peaking is obtained, in particular for scans in collisionality and temperature ratio.
rate research
Read More
Gyrokinetic simulations of ion temperature gradient mode and trapped electron mode driven impurity transport in a realistic tokamak geometry are presented and compared with results using simplified geometries. The gyrokinetic results, obtained with the GENE code in both linear and non-linear modes are compared with data and analysis for a dedicated impurity injection discharge at JET. The impact of several factors on heat and particle transport is discussed, lending special focus to tokamak geometry and rotational shear. To this end, results using s-alpha and concentric circular equilibria are compared with results with magnetic geometry from a JET experiment. To further approach experimental conditions, non-linear gyrokinetic simulations are performed with collisions and a carbon background included. The impurity peaking factors, computed by finding local density gradients corresponding to zero particle flux, are discussed. The impurity peaking factors are seen to be reduced by a factor of ~2 in realistic geometry compared with the simplified geometries, due to a reduction of the convective pinch. It is also seen that collisions reduce the peaking factor for low-Z impurities, while increasing it for high charge numbers, which is attributed to a shift in the transport spectra towards higher wavenumbers with the addition of collisions. With the addition of roto-diffusion, an overall reduction of the peaking factors is observed, but this decrease is not sufficient to explain the flat carbon profiles seen at JET.
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).
The XGC1 edge gyrokinetic code is used for a high fidelity prediction for the width of the heat-flux to divertor plates in attached plasma condition. The simulation results are validated against the empirical scaling $lambda_q propto B_P^{-gamma}$ obtained from present tokamak devices, where $lambda_q$ is the divertor heat-flux width mapped to the outboard midplane and $gamma_q=1.19$ as defined by T. Eich et al. [Nucl. Fusion 53 (2013) 093031], and $B_P$ is the magnitude of the poloidal magnetic field at outboard midplane separatrix surface. This empirical scaling predicts $lambda_q leq 1mm$ when extrapolated to ITER, which would require operation with very high separatrix densities $(n_{sep}/n_{Greenwald} > 0.6)$ in the Q=10 scenario to achieve semi-detached plasma operation and high radiative fractions leading to acceptable divertor power fluxes. XGC1 predicts, however, that $lambda_q$ for ITER is over 5 mm, suggesting that operation in the ITER Q=10 scenario with acceptable divertor power loads could be obtained over a wider range of plasma separatrix densities and radiative fractions. The physics reason behind this difference is, according to the XGC1 results, that while the ion magnetic drift contribution to the divertor heat-flux width is wider in the present tokamaks, the turbulent electron contribution is wider in ITER. A high current C-Mod discharge is found to be in a mixed regime: While the heat-flux width by the ion neoclassical magnetic drift is still wider than the turbulent electron heat-flux width, the heat-flux magnitude is dominated by the narrower electron heat-flux.
Linear and nonlinear modelling of Alfvenic instabilities, most notably toroidal Alfven eigenmodes (TAEs), obtained by using the global nonlinear electromagnetic gyrokinetic model of the code ORB5 are presented for the 15 MA scenario of the ITER tokamak. Linear simulations show that elliptic Alfven eigenmodes and odd-parity TAEs are only weakly damped but not excited by alpha particles, whose drive favours even-parity TAEs. Low mode number TAEs are found to be global, requiring global treatment. Nonlinearly, even with double the nominal EP density, single mode simulations lead to saturation with negligible EP transport however multi-mode simulations predict that with double the nominal EP density, enhanced saturation and significant EP redistribution will occur.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
