Toroidal rotation is critical for fusion in tokamaks, since it stabilizes instabilities that can otherwise cause disruptions or degrade confinement. Unlike present-day devices, ITER might not have enough neutral-beam torque to easily avoid these instabilities. We must therefore understand how the plasma rotates intrinsically, that is, without applied torque. Experimentally, torque-free plasmas indeed rotate, with profiles that are often non-flat and even non-monotonic. The rotation depends on many plasma parameters including collisionality and plasma current, and exhibits sudden bifurcations (rotation reversals) at critical parameter values.Since toroidal angular momentum is conserved in axisymmetric systems, and since experimentally inferred momentum transport is much too large to be neoclassical, theoretical work has focused on rotation drive by nondiffusive turbulent momentum fluxes. In the edge, intrinsic rotation relaxes to a steady state in which the total momentum outflux from the plasma vanishes. Ion drift orbits, scrape-off-layer flows, separatrix geometry, and turbulence intensity gradient all play a role. In the core, nondiffusive and viscous momentum fluxes balance to set the rotation gradient at each flux surface. Although many mechanisms have been proposed for the nondiffusive fluxes, most are treated in one of two distinct but related gyrokinetic formulations. In a radially local fluxtube, appropriate for rho star <<1, the lowest-order gyrokinetic formulations exhibit a symmetry that prohibits nondiffusive momentum flux for nonrotating plasmas in an up- down symmetric magnetic geometry with no ExB shear. Many symmetry-breaking mechanisms have been identified, but none have yet been conclusively demonstrated to drive a strong enough flux to explain commonly observed experimental rotation profiles. Radially global gyrokinetic simulations naturally include many symmetry-breaking mechanisms, and have shown cases with experimentally relevant levels of nondiffusive flux. These promising early results motivate further work to analyze, verify, and validate.This article provides a pedagogical introduction to intrinsic rotation in axisymmetric devices. Intended for both newcomers to the topic and experienced practitioners, the article reviews a broad range of topics including experimental and theoretical results for both edge and core rotation, while maintaining a focus on the underlying concepts.
An optimization approach that incorporates the predictive transport code TRANSP is proposed for tokamak scenario development. Optimization methods are often employed to develop open-loop control strategies to aid access to high performance tokamak scenarios. In general, the optimization approaches use control-oriented models, i.e. models that are reduced in complexity and prediction accuracy as compared to physics-oriented transport codes such as TRANSP. In the presented approach, an optimization procedure using the TRANSP code to simulate the tokamak plasma is considered for improved predictive capabilities. As a test case, the neutral beam injection (NBI) power is optimized to develop a control strategy that maximizes the non-inductive current fraction during the ramp-up phase for NSTX-U. Simulation studies towards the achievement of non-inductive ramp up in NSTX-U have already been carried out with the TRANSP code. The optimization-based approach proposed in this work is used to maximize the non-inductive current fraction during ramp-up in NSTX-U, demonstrating that the scenario development task can be automated. An additional test case considers optimization of the current ramp rate in DIII-D for obtaining a stationary plasma characterized by
a flat loop voltage profile in the flattop phase.
Zweben SJ, Myra JR, Diallo A, Russell DA, Scotti F, Stotler DP
Transient small-scale structures were identified in the wake of blobs movingpoloidally through the SOL of high-powered H-mode plasmas in NSTX, using the gaspuff imaging (GPI) diagnostic. These blob wakes had a poloidal wavelength in therange 3.5 cm, which is significantly smaller than the average blob scale of~12 cm, and the wakes had a poloidal velocity of 1.5 km/sec in theelectron diamagnetic direction, which is opposite to the blob poloidal velocity inthese shots. These wakes were radially localized 0-4 cm outside the separatrix andoccurred within ~50 microsec after the passage of a blob through the GPI field of view.The clearest wakes were seen when the GPI viewing angle was well aligned with thelocal B field line, as expected for such small-scale structures given the diagnosticgeometry. A plausible theoretical interpretation of the wakes is discussed: theobserved wakes share some features of drift waves and/or drift-Alfven waves whichcould be excited
The Magnetospheric Multiscale (MMS) mission has given us unprecedented access to high cadence particle and field data of magnetic reconnection at Earth's magnetopause. MMS first passed very near an X-line on 16 October 2015, the Burch event, and has since observed multiple X-line crossings. Subsequent 3D particle-in-cell (PIC) modeling efforts of and comparison with the Burch event have revealed a host of novel physical insights concerning magnetic reconnection, turbulence induced particle mixing, and secondary instabilities. In this study, we employ the Gkeyll simulation framework to study the Burch event with different classes of extended, multi-fluid magnetohydrodynamics (MHD), including models that incorporate important kinetic effects, such as the electron pressure tensor, with physics-based closure relations designed to capture linear Landau damping. Such fluid modeling approaches are able to capture different levels of kinetic physics in global simulations and are generally less costly than fully kinetic PIC. We focus on the additional physics one can capture with increasing levels of fluid closure refinement via comparison with MMS data and existing PIC simulations. In particular, we find that the ten-moment model well captures the agyrotropic structure of the pressure tensor in the vicinity of the X-line and the magnitude of anisotropic electron heating observed in MMS and PIC simulations. However, the ten-moment model has difficulty resolving the lower hybrid drift instability, which has been observed to plays a fundamental role in heating and mixing electrons in the current layer.
Force-driven parallel shear flow in a spatially periodic domain is shown to be linearly unstable
with respect to both the Reynolds number and the domain aspect ratio. This finding is confirmed
by computer simulations, and a simple expression is derived to determine stable flow conditions.
Periodic extensions of Couette and Poiseuille flows are unstable at Reynolds numbers two orders
of magnitude smaller than their aperiodic equivalents because the periodic boundaries impose
fundamentally different constraints. This instability has important implications for designing computational models of nonlinear dynamic processes with periodicity.
The data are 4554 light curves derived from images taken of the globular cluster M4 by the Kepler space telescope during the K2 portion of its mission, specifically during Campaign 2 of that mission, which occurred in 2014. A total of 3856 images were taken over approximately three months at a cadence of approximately half an hour. The purpose of these observations was to find stars and other objects that vary in brightness over time --- variable stars. Also included is a table with associated information for each of the 4554 objects and their light curves.