No Arabic abstract
A procedure is presented to map from the spatial correlation parameters of a turbulent density field (the radial and binormal correlation lengths and wavenumbers, and the fluctuation amplitude) to correlation parameters that would be measured by a Beam Emission Spectroscopy (BES) diagnostic. The inverse mapping is also derived, which results in resolution criteria for recovering correct correlation parameters, depending on the spatial response of the instrument quantified in terms of Point-Spread Functions (PSFs). Thus, a procedure is presented that allows for a systematic comparison between theoretical predictions and experimental observations. This procedure is illustrated using the MAST BES system and the validity of the underlying assumptions is tested on fluctuating density fields generated by direct numerical simulations using the gyrokinetic code GS2. The measurement of the correlation time, by means of the cross-correlation time-delay (CCTD) method, is also investigated and is shown to be sensitive to the fluctuating radial component of velocity, as well as to small variations in the spatial properties of the PSFs.
The Beam Emission Spectroscopy (BES) turbulence diagnostic on MAST is to be upgraded in June 2010 from a 1D trial system to a 2D imaging system (8 radial times 4 poloidal channels) based on a newly developed APD array camera. The spatial resolution of the new system is calculated in terms of the point spread function (PSF) to account for the effects of field-line curvature, observation geometry, the finite lifetime of the excited state of the beam atoms, and beam attenuation and divergence. It is found that the radial spatial resolution is ~ 2-3 cm and the poloidal spatial resolution ~ 1-5 cm depending on the radial viewing location. The absolute number of detected photons is also calculated, hence the photon noise level can be determined
Numerical simulation of plasma turbulence in the Large Plasma Device (LAPD) [Gekelman et al, Rev. Sci. Inst., 62, 2875, 1991] is presented. The model, implemented in the BOUndary Turbulence (BOUT) code [M. Umansky et al, Contrib. Plasma Phys. 180, 887 (2009)], includes 3-D collisional fluid equations for plasma density, electron parallel momentum, and current continuity, and also includes the effects of ion-neutral collisions. In nonlinear simulations using measured LAPD density profiles but assuming constant temperature profile for simplicity, self-consistent evolution of instabilities and nonlinearly-generated zonal flows results in a saturated turbulent state. Comparisons of these simulations with measurements in LAPD plasmas reveal good qualitative and reasonable quantitative agreement, in particular in frequency spectrum, spatial correlation and amplitude probability distribution function of density fluctuations. For comparison with LAPD measurements, the plasma density profile in simulations is maintained either by direct azimuthal averaging on each time step, or by adding particle source/sink function. The inferred source/sink values are consistent with the estimated ionization source and parallel losses in LAPD. These simulations lay the groundwork for more a comprehensive effort to test fluid turbulence simulation against LAPD data.
We analyze plasma heating in weakly collisional kinetic Alfven wave (KAW) turbulence using high resolution gyrokinetic simulations spanning the range of scales between the ion and the electron gyroradii. Real space structures that have a higher than average heating rate are shown not to be confined to current sheets. This novel result is at odds with previous studies, which use the electromagnetic work in the local electron fluid frame, i.e. $mathbf{J} !cdot! (mathbf{E} + mathbf{v}_etimesmathbf{B})$, as a proxy for turbulent dissipation to argue that heating follows the intermittent spatial structure of the electric current. Furthermore, we show that electrons are dominated by parallel heating while the ions prefer the perpendicular heating route. We comment on the implications of the results presented here.
Observations of ion-scale (k_y*rho_i <= 1) density turbulence of relative amplitude dn_e/n_e <= 0.2% are available on the Mega Amp Spherical Tokamak (MAST) using a 2D (8 radial x 4 poloidal channel) imaging Beam Emission Spectroscopy (BES) diagnostic. Spatial and temporal characteristics of this turbulence, i.e., amplitudes, correlation times, radial and perpendicular correlation lengths and apparent phase velocities of the density contours, are determined by means of correlation analysis. For a low-density, L-mode discharge with strong equilibrium flow shear exhibiting an internal transport barrier (ITB) in the ion channel, the observed turbulence characteristics are compared with synthetic density turbulence data generated from global, non-linear, gyro-kinetic simulations using the particle-in-cell (PIC) code NEMORB. This validation exercise highlights the need to include increasingly sophisticated physics, e.g., kinetic treatment of trapped electrons, equilibrium flow shear and collisions, to reproduce most of the characteristics of the observed turbulence. Even so, significant discrepancies remain: an underprediction by the simulations of the turbulence amplituide and heat flux at plasma periphery and the finding that the correlation times of the numerically simulated turbulence are typically two orders of magnitude longer than those measured in MAST. Comparison of these correlation times with various linear timescales suggests that, while the measured turbulence is strong and may be `critically balanced, the simulated turbulence is weak.
Turbulence is a major factor limiting the achievement of better tokamak performance as it enhances the transport of particles, momentum and heat which hinders the foremost objective of tokamaks. Hence, understanding and possibly being able to control turbulence in tokamaks is of paramount importance, not to mention our intellectual curiosity of it.