Asymmetric one-dimensional slow electron holes.

Slow solitary positive-potential peaks sustained by trapped electron deficit in a plasma with asymmetric ion velocity distributions are in principle asymmetric, involving a potential change across the hole. It is shown theoretically how to construct such asymmetric electron holes, thus providing fully consistent solutions of the one-dimensional Vlasov-Poisson equation for a wide variety of prescribed background ion velocity distributions. Because of ion reflection forces experienced by the hole, there is generally only one discrete slow hole velocity that is in equilibrium. Moreover the equilibrium is unstable unless there is a local minimum in the ion velocity distribution, in which the hole velocity then resides. For stable equilibria with Maxwellian electrons, the potential drop across the hole is shown to be Δϕ≃2/9f^{'''}T_{e}/e(eψ/m_{i})^{2}, where ψ is the hole peak potential, f^{'''} is the third derivative of the background ion velocity distribution function at the hole velocity, and T_{e} is the electron temperature. Potential asymmetry is small for holes of the amplitudes usually observed, ψ≲0.5T_{e}/e.

How can slow plasma electron holes exist?

One-dimensional analysis is presented of solitary positive potential plasma structures whose velocity lies within the range of ion distribution velocities that are strongly populated: "slow" electron holes. It is shown that to avoid the self-acceleration of the hole velocity away from ion velocities it must lie within a local minimum in the ion velocity distribution. Quantitative criteria for the existence of stable equilibria are obtained. The background ion distributions required are generally stable to ion-ion modes unless the electron temperature is much higher than the ion temperature. Since slow positive potential solitons are shown not to be possible without a significant contribution from trapped electrons, it seems highly likely that such observed slow potential structures are indeed electron holes.

Transverse instability magnetic field thresholds of electron phase-space holes.

A detailed comparison is presented of analytical and particle-in-cell simulation investigation of the transverse instability, in two dimensions, of initially one-dimensional electron phase-space hole equilibria. Good quantitative agreement is found between the shift-mode analysis and the simulations for the magnetic field (B) threshold at which the instability becomes overstable (time oscillatory) and for the real and imaginary parts of the frequency. The simulation B threshold for full stabilization exceeds the predictions of shift-mode analysis by 20-30%, because the mode becomes substantially narrower in spatial extent than a pure shift. This threshold shift is qualitatively explained by the kinematic mechanism of instability.

Kinematic Mechanism of Plasma Electron Hole Transverse Instability.

It is shown through multidimensional particle-in-cell simulations that at least in Maxwellian background plasmas the long-wavelength transverse instability of plasma electron holes is caused not by the previously proposed focusing of trapped particles but instead by kinematic jetting of marginally passing electrons. The mechanism is explained and heuristic analytic estimates obtained which agree with the growth rates and transverse wave numbers observed in the simulations.

Effects of magnetic shear on toroidal rotation in tokamak plasmas with lower hybrid current drive.

Application of lower hybrid (LH) current drive in tokamak plasmas can induce both co- and countercurrent directed changes in toroidal rotation, depending on the core q profile. For discharges with q(0) <1, rotation increments in the countercurrent direction are observed. If the LH-driven current is sufficient to suppress sawteeth and increase q(0) above unity, the core toroidal rotation change is in the cocurrent direction. This change in sign of the rotation increment is consistent with a change in sign of the residual stress (the divergence of which constitutes an intrinsic torque that drives the flow) through its dependence on magnetic shear.

X-ray imaging crystal spectroscopy for use in plasma transport research.

This research describes advancements in the spectral analysis and error propagation techniques associated with x-ray imaging crystal spectroscopy (XICS) that have enabled this diagnostic to be used to accurately constrain particle, momentum, and heat transport studies in a tokamak for the first time. Doppler tomography techniques have been extended to include propagation of statistical uncertainty due to photon noise, the effect of non-uniform instrumental broadening as well as flux surface variations in impurity density. These methods have been deployed as a suite of modeling and analysis tools, written in interactive data language (IDL) and designed for general use on tokamaks. Its application to the Alcator C-Mod XICS is discussed, along with novel spectral and spatial calibration techniques. Example ion temperature and radial electric field profiles from recent I-mode plasmas are shown, and the impact of poloidally asymmetric impurity density and natural line broadening is discussed in the context of the planned ITER x-ray crystal spectrometer.

Intergrain forces in low-Mach-number plasma wakes.

Large-scale particle-in-cell calculations of the plasma wake interactions of two negatively charged grains smaller than the Debye length are carried out using the coptic code over a wide range of subsonic plasma flow velocities. In plasmas with the temperature ratio T(e)/T(i)=100, it is found that a single grain's oscillatory wake disappears for flow Mach numbers M less than approximately 0.3, which is the parameter regime where Landau damping is expected to be strong. Neutral collisions suppress potential oscillations above M=0.3, but not the trailing attractive potential peak caused by ion focusing. The transverse (grain-aligning) force on a downstream particle in the wake of another is obtained rigorously from the code in three-dimensional simulations. It shows general agreement with the force that would be deduced from the single-grain wake potential gradient. Except for relatively large grains in the nonlinear collisional regime, the grain-aligning force is very small for slow flow.

Forces on a small grain in the nonlinear plasma wake of another.

The transverse force on a spherical charged grain lying in the plasma wake of another grain is analyzed to assess the importance of ion-drag perturbation, in addition to the wake-potential-gradient. The ion-drag perturbation is intrinsically one order smaller than the wake-potential force in the ratio of grain size (r(p)) to Debye length (λ(De)). So ion-drag perturbation is important only in nonlinear wakes. Rigorous particle-in-cell calculations of the force are performed in the nonlinear regime with two interacting grains. It is found that even for quite large grains, r(p)/λ(De)=0.1, the force is dominated by the wake-potential gradient. The wake-potential structure can then help explain the preferred alignment of floating dust grains.

Observation of self-generated flows in tokamak plasmas with lower-hybrid-driven current.

In Alcator C-Mod discharges lower hybrid waves have been shown to induce a countercurrent change in toroidal rotation of up to 60 km/s in the central region of the plasma (r/a approximately <0.4). This modification of the toroidal rotation profile develops on a time scale comparable to the current redistribution time (approximately 100 ms) but longer than the energy and momentum confinement times (approximately 20 ms). A comparison of the co- and countercurrent injected waves indicates that current drive (as opposed to heating) is responsible for the rotation profile modifications. Furthermore, the changes in central rotation velocity induced by lower hybrid current drive (LHCD) are well correlated with changes in normalized internal inductance. The application of LHCD has been shown to generate sheared rotation profiles and a negative increment in the radial electric field profile consistent with a fast electron pinch.

Two dimensional radiated power diagnostics on Alcator C-Mod.

The radiated power diagnostics for the Alcator C-Mod tokamak have been upgraded to measure two dimensional structure of the photon emissivity profile in order to investigate poloidal asymmetries in the core radiation. Commonly utilized unbiased absolute extreme ultraviolet (AXUV) diode arrays view the plasma along five different horizontal planes. The layout of the diagnostic set is shown and the results from calibrations and recent experiments are discussed. Data showing a significant, 30%-40%, inboard/outboard emissivity asymmetry during ELM-free H-mode are presented. The ability to use AXUV diode arrays to measure absolute radiated power is explored by comparing diode and resistive bolometer-based emissivity profiles for highly radiative L-mode plasmas seeded with argon. Emissivity profiles match in the core but disagree radially outward resulting in an underprediction of P(rad) of nearly 50% by the diodes compared to P(rad) determined using resistive bolometers.

Ion collection by oblique surfaces of an object in a transversely flowing strongly magnetized plasma.

The equations governing a collisionless obliquely flowing plasma around an ion-absorbing object in a strong magnetic field are shown to have an exact analytic solution even for an arbitrary (two-dimensional) object shape, when temperature is uniform, and diffusive transport can be ignored. The solution has an extremely simple geometric embodiment. It shows that the ion collection flux density to a convex body's surface depends only upon the orientation of the surface and provides the theoretical justification and calibration of oblique "Mach probes." The exponential form of this exact solution helps explain the approximate fit of this function to previous numerical solutions.

Observation of anomalous momentum transport in tokamak plasmas with no momentum input.

Anomalous momentum transport has been observed in Alcator C-Mod tokamak plasmas through analysis of the time evolution of core impurity toroidal rotation velocity profiles. Following the L-mode to EDA (enhanced D(alpha)) H-mode transition, the ensuing cocurrent toroidal rotation velocity, which is generated in the absence of any external momentum source, is observed to propagate in from the edge plasma to the core. The steady state toroidal rotation velocity profiles are relatively flat and the momentum transport can be simulated with a simple diffusion model. Velocity profiles during edge localized mode free (ELM-free) H-modes are centrally peaked, which suggests the addition of inward momentum convection. In all operating regimes the observed momentum diffusivities are much larger than the neoclassical values.