Search

We present Magnetospheric Multiscale observations showing large numbers of slow electron holes with speeds clustered near the local minimum of double-humped velocity distribution functions of background ions. Theoretical computations show that slow electron holes can avoid the acceleration that otherwise prevents their remaining slow only under these same circumstances. Although the origin of the slow electron holes is still elusive, the agreement between observation and theory about the conditions for their existence is remarkable.

We present Magnetospheric Multiscale observations showing large numbers of slow electron holes with speeds clustered near the local minimum of double-humped velocity distribution functions of background ions. Theoretical computations show that slow electron holes can avoid the acceleration that otherwise prevents their remaining slow only under these same circumstances. Although the origin of the slow electron holes is still elusive, the agreement between observation and theory about the conditions for their existence is remarkable.

Analytic treatment is presented of the electrostatic instability of an initially planar electron hole in a plasma of effectively infinite particle magnetization. It is shown that there is an unstable mode consisting of a rigid shift of the hole in the trapping direction. Its low frequency is determined by the real part of the force balance between the Maxwell stress arising from the transverse wavenumber and the kinematic jetting from the hole's acceleration. The very low growth rate arises from a delicate balance in the imaginary part of the force between the passing-particle jetting, which is destabilizing, and the resonant response of the trapped particles, which is stabilizing. Nearly universal scalings of the complex frequency and with hole depth are derived. Particle in cell simulations show that the slow-growing instabilities previously investigated as coupled hole-wave phenomena occur at the predicted frequency, but with growth rates 2 to 4 times greater than the analytic prediction. This higher rate may be caused by a reduced resonant stabilization because of numerical phase-space diffusion in the simulations. ©2019 Keywords: plasma instabilities; plasma nonlinear phenomena; space plasma physics, NASA (Grant NNX16AG82G)

The mitigation of divertor heat fluxes is an active topic of investigation on existing tokamaks. One approach uses radiation, both inside and outside the last closed flux surface (LCFS), to convert plasma thermal energy, usually directed towards dedicated plasma facing components, to soft X-ray and ultraviolet radiation, spread over a much larger surface area. Recent enhanced D-α H-mode experiments on Alcator C-Mod varied the ICRF input power and radiative power losses via impurity seeding to demonstrate that normalized energy confinement depends strongly on the difference between input power and the radiated power inside the LCFS. These investigations also show that when seeded with either Ne or N2, a factor of two and higher reduction in outer divertor heat flux is achieved while maintaining H[subscript 98,y2] ∼ 1.0. Conversely, when seeding with Ar, confinement is limited to H[subscript 98,y2] ∼ 0.8 for a similar level of exhaust power., United States. Dept. of Energy (DOE Contract Number DEFC0299ER54512)

We present inboard (HFS) and outboard (LFS) radial electric field (E[subscript r]) and impurity temperature (T[subscript z]) measurements in the I-mode and H-mode pedestal of Alcator C-Mod. These measurements reveal strong Er wells at the HFS and the LFS midplane in both regimes and clear pedestals in T[subscript z], which are of similar shape and height for the HFS and LFS. While the H-mode E[subscript r] well has a radially symmetric structure, the E[subscript r] well in I-mode is asymmetric, with a stronger ExB shear layer at the outer edge of the E[subscript r] well, near the separatrix. Comparison of HFS and LFS profiles indicates that impurity temperature and plasma potential are not simultaneously flux functions. Uncertainties in radial alignment after mapping HFS measurements along flux surfaces to the LFS do not, however, allow direct determination as to which quantity varies poloidally and to what extent. Radially aligning HFS and LFS measurements based on the T[subscript z] profiles would result in substantial inboard-outboard variations of plasma potential and electron density. Aligning HFS and LFS E[subscript r] wells instead also approximately aligns the impurity poloidal flow profiles, while resulting in a LFS impurity temperature exceeding the HFS values in the region of steepest gradients by up to 70%. Considerations based on a simplified form of total parallel momentum balance and estimates of parallel and perpendicular heat transport time scales seem to favor an approximate alignment of the E[subscript r] wells and a substantial poloidal asymmetry in impurity temperature., United States. Dept. of Energy (Cooperative Agreement DE-FC02-99ER54512), Swiss National Science Foundation

The ion drag force arising from plasma flow past an embedded spherical grain is calculated self-consistently and non-linearly using particle in cell codes, accounting for ion-neutral collisions. Using ion velocity distribution appropriate for ion drift driven by a force field gives wake potential and force greatly different from a shifted Maxwellian distribution, regardless of collisionality. The low-collisionality forces are shown to be consistent with estimates based upon cross-sections for scattering in a Yukawa (shielded) grain field, but only if non-linear shielding length is used. Finite collisionality initially enhances the drag force, but only by up to a factor of 2. Larger collisionality eventually reduces the drag force. In the collisional regime, the drift distribution gives larger drag than the shift distribution even at velocities where their collisionless drags are equal. Comprehensive practical analytic formulas for force that fit the calculations are provided., United States. Dept. of Energy (National Science Foundation (U.S.). Grant DE-FG02-06ER54982), United States. Dept. of Energy. Office of Science (Graduate Research Fellowship Program)

The observation of parallel ion velocity in the near-lunar wake approximately equal to external solar wind velocity can be explained within uncertainties by an analytic electrostatic expansion model. The one-dimensional model frequently used is inadequate because it does not account for the moon's spherical shape. However, application of a more recent generalization to three dimensions of the solution along characteristics predicts higher velocities and is probably sufficient to account for the SARA observations on the Chandrayaan-1 spacecraft.

The plasma force on grains of specified charge and height in a collisional DC plasma sheath is calculated using the multidimensional particle in cell code COPTIC. The background ion velocity distribution functions for the unperturbed sheath vary substantially with collisionality. The grain force is found to agree quite well with a combination of background electric-field force plus ion drag force. However, the drag force must take account of the non-Maxwellian (and spatially varying) ion distribution function, and the collisional drag enhancement. It is shown how to translate the dimensionless results into practical equilibrium including other forces such as gravity., United States. Dept. of Energy (National Science Foundation (U.S.). Grant DE-FG02-06ER54982)

The wake behind a large object (such as the moon) moving rapidly through a plasma (such as the solar wind) contains a region of depleted density, into which the plasma expands along the magnetic field, transverse to the flow. It is shown here that (in addition to any ion instability) a bump-on-tail which is unstable appears on the electrons' parallel velocity distribution function because of the convective non-conservation of parallel energy (drift-energization). It arises regardless of any non-thermal features on the external electron velocity distribution. The detailed electron distribution function throughout the wake is calculated by integration along orbits; and the substantial energy level of resulting electron plasma (Langmuir) turbulence is evaluated quasi-linearly. It peaks near the wake axis. If the mass of the electrons is artificially enhanced, for example in order to make numerical simulation feasible, then much more unstable electron distributions arise; but these are caused by the unphysical mass ratio., National Science Foundation (U.S.), United States. Dept. of Energy (Grant DE-FG02-06ER54982)

Elementary but general statistical analyses determine the uncertainty arising from photon statistics in measuring a line shift and width. Account is taken of a background as well as the required signal.

Ion cyclotron resonance heating is the primary auxiliary heating on the Alcator C-Mod tokamak and is commonly used on other devices, and is planned for use on ITER. The RF-power density on C-Mod is above 5 MW m−3 providing for a unique opportunity to study wave–particle effects in the high RF power per particle regime. Minority heating produces a highly energetic tail in the minority distribution function which is measured using a compact neutral particle analyser. In this paper, we present the measurements of the fast-ion spectrum between 200 and 2 MeV, compiled over an entire experimental campaign. We also estimate the effective tail temperatures for the fast-ion distribution. We find that the fast-ion distribution is less energetic and less dense with increasing electron density; is more energetic with increasing plasma current; and is more dense but has no measurable change in energy with increasing RF power. Some possible explanations for these findings are discussed., United States. Dept. of Energy (Award DE-FC02-99ER54512)

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[subscript 0] < 1, rotation increments in the countercurrent direction are observed. If the LH-driven current is sufficient to suppress sawteeth and increase q[subscript 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., United States. Dept. of Energy (Contract DE-FC02-99ER54512), United States. Dept. of Energy (Fusion Research Postdoctoral Research Program)

Boundary layer flows in the Alcator C-Mod tokamak are systematically examined as magnetic topology (upper versus lower-null) and plasma density are changed. Utilizing a unique set of scanning Langmuir–Mach probes, including one on the high-field side (HFS) midplane, the poloidal variation of plasma flow components in the parallel, diamagnetic and radial directions are resolved in detail. It is found that the plasma flow pattern can be decomposed into two principal parts: (1) a drift-driven component, which lies within a magnetic flux surface and is divergence-free and (2) a transport-driven component, which gives rise to near-sonic parallel flows on the HFS scrape-off layer (SOL). Toroidal rotation, Pfirsch–Schlüter and transport-driven contributions are unambiguously identified. Transport-driven parallel flows are found to dominate the HFS particle fluxes; the total poloidal-directed flow accounts for ~1/3 to all of the ion flux arriving on the inner divertor. As a result, heat convection is found to be an important player in this region, consistent with the observation of divertor asymmetries that depend on the direction of B × ∇B relative to the active x-point. In contrast, the poloidal projection of parallel flow in the low-field SOL largely cancels with E[subscript r] × B flow; toroidal rotation is the dominant plasma motion there. The magnitude of the transport-driven poloidal flow is found to be quantitatively consistent with fluctuation-induced radial particle fluxes on the low-field side (LFS), identifying this as the primary drive mechanism. Fluctuation-induced fluxes on the HFS are found to be essentially zero, excluding turbulent inward transport as the mechanism that closes the circulation loop in this region., United States. Dept. of Energy (Cooperative Agreement DE-FC02-99ER54512)

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. Dopplertomography 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., United States. Dept. of Energy (Contract DE-FC02-99ER54512), United States. Dept. of Energy (Contract DE-AC02-09CH11466), United States. Dept. of Energy (Fusion Research Postdoctoral Research Program)

In the Alcator C-Mod tokamak, strong, steady-state variations of molybdenum density within a flux surface are routinely observed in plasmas using hydrogen minority ion cyclotron resonant heating. In/out asymmetries, up to a factor of 2, occur with either inboard or outboard accumulation depending on the major radius of the minority resonance layer. These poloidal variations can be attributed to the impurity's high charge and large mass in the neoclassical parallel force balance. The large mass enhances the centrifugal force, causing outboard accumulation while the high charge enhances ion-impurity friction and makes impurities sensitive to small poloidal variations in the plasma potential. Quantitative comparisons between existing parallel high-Z impurity transport theories and experimental results for r/a < 0.7 show good agreement when the resonance layer is on the high-field side of the tokamak but disagree substantially for low-field side heating. Ion-impurity friction is insufficient to explain the experimental results, and the accumulation of impurity density on the inboard side of flux surface is shown to be driven by a poloidal potential variation due to magnetic trapping of non-thermal, cyclotron heated minority ions. Parallel impurity transport theory is extended to account for cyclotron effects and shown to agree with experimentally measured impurity density asymmetries., United States. Dept. of Energy (Agreement DE-FC02-99ER54512), United States. Dept. of Energy. Office of Fusion Energy Sciences (Postdoctoral Research Program)

The forces acting on a spherical conducting particle in a transversely flowing magnetized plasma are calculated in the entire range of magnetization and Debye length, using the particle code SCEPTIC3D (Patacchini and Hutchinson 2010 Plasma Phys. Control. Fusion 52 035005, 2011 Plasma Phys. Control. Fusion 53 025005). In short Debye length (i.e. high density) plasmas, both the ion-drag and Lorentz force arising from currents circulating inside the dust show strong components antiparallel to the convective electric field, suggesting that a free dust particle should gyrate faster than what predicted by its Larmor frequency. In intermediate to large Debye length conditions, by a downstream depletion effect already reported in unmagnetized strongly collisional regimes, the ion-drag in the direction of transverse flow can become negative. The internal Lorentz force, however, remains in the flow direction, and large enough in magnitude so that no spontaneous dust motion should occur., National Science Foundation (U.S.), United States. Dept. of Energy (grant DE-FG02-06ER54891)

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[subscript e]/T[subscript 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., National Science Foundation (U.S.). (United States. Dept. of Energy.) (Grant Number DE-FG02- 06ER54982)

Recent computational and theoretical progress in understanding and calculating ion collection by negatively charged absorbing objects in a flowing plasma is outlined. The results are placed in the context of key theoretical achievements of prior research. Despite the topic's long history, and past profound insights, fully rigorous quantitative solution of the non-linear, multidimensional, self-consistent, kinetic-theory problem has not until recently been feasible. Now we are able to establish the adequacy or inadequacy of approximate treatments, and provide critical quantitative results. In the process, some qualitative surprises have also emerged.