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.
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.
Woods, B. J. Q.; Duarte, V. N.; Fredrickson, E. D.; Gorelenkov, N. N.; Podestà, M.; Vann, R. G. L.
Abrupt large events in the Alfvenic and sub-Alfvenic frequency bands in tokamaks are typically correlated with increased fast-ion loss. Here, machine learning is used to speed up the laborious process of characterizing the behavior of magnetic perturbations from corresponding frequency spectrograms that are typically identified by humans. The analysis allows for comparison between different mode character (such as quiescent, fixed frequency, and chirping, avalanching) and plasma parameters obtained from the TRANSP code, such as the ratio of the neutral beam injection (NBI) velocity and the Alfven velocity (v_inj./v_A), the q-profile, and the ratio of the neutral beam beta and the total plasma beta (beta_beam,i / beta). In agreement with the previous work by Fredrickson et al., we find a correlation between beta_beam,i and mode character. In addition, previously unknown correlations are found between moments of the spectrograms and mode character. Character transition from quiescent to nonquiescent behavior for magnetic fluctuations in the 50200-kHz frequency band is observed along the boundary v_phi ~ (1/4)(v_inj. - 3v_A), where v_phi is the rotation velocity.
Since 1850 the concentration of atmospheric methane (CH4), a potent greenhouse gas, has more than doubled. Recent studies suggest that emission inventories may be missing sources and underestimating emissions. To investigate whether offshore oil and gas platforms leak CH4 during normal operation, we measured CH4 mole fractions around eight oil and gas production platforms in the North Sea which were neither flaring gas nor off-loading oil. We use the measurements from summer 2017, along with meteorological data, in a Gaussian plume model to estimate CH4 emissions from each platform. We find CH4 mole fractions of between 11 and 370 ppb above background concentrations downwind of the platforms measured, corresponding to a median CH4 emission of 6.8 g CH4 s-1 for each platform, with a range of 2.9 to 22.3 g CH4 s-1. When matched to production records, during our measurements individual platforms lost between 0.04% and 1.4% of gas produced with a median loss of 0.23%. When the measured platforms are considered collectively, (i.e. the sum of platforms’ emission fluxes weighted by the sum of the platforms’ production), we estimate the CH4 loss to be 0.19% of gas production. These estimates are substantially higher than the emissions most recently reported to the National Atmospheric Emission Inventory (NAEI) for total CH4 loss from United Kingdom platforms in the North Sea. The NAEI reports CH4 losses from the offshore oil and gas platforms we measured to be 0.13% of gas production, with most of their emissions coming from gas flaring and offshore oil loading, neither of which were taking place at the time of our measurements. All oil and gas platforms we observed were found to leak CH4 during normal operation and much of this leakage has not been included in UK emission inventories. Further research is required to accurately determine total CH4 leakage from all offshore oil and gas operations and to properly include the leakage in national and international emission inventories.
The Electromagnetic Particle Injector (EPI) concept is advanced through the simulation of ablatant deposition into ITER H-mode discharges with calculations showing penetration past the H-mode pedestal for a range of injection velocities and granule sizes concurrent with the requirements of disruption mitigation. As discharge stored energy increases in future fusion devices such as ITER, control and handling of disruption events becomes a critical issue. An unmitigated disruption could lead to failure of the plasma facing components resulting in financially and politically costly repairs. Methods to facilitate the quench of an unstable high current discharge are required. With the onset warning time for some ITER disruption events estimated to be less than 10 ms, a disruption mitigation system needs to be considered which operates at injection speeds greater than gaseous sound speeds. Such an actuator could then serve as a means to augment presently planned pneumatic injection systems. The EPI uses a rail gun concept whereby a radiative payload is delivered into the discharge by means of the JxB forces generated by an external current pulse, allowing for injection velocities in excess of 1 km/s. The present status of the EPI project is outlined, including the addition of boost magnetic coils. These coils augment the self-generated rail gun magnetic field and thus provide a more efficient acceleration of the payload. The coils and the holder designed to constrain them have been modelled with the ANSYS code to ensure structural integrity through the range of operational coil cu
In homogeneous drift-wave (DW) turbulence, zonal flows (ZFs) can be generated via a modulational instability (MI) that either saturates monotonically or leads to oscillations of the ZF energy at the nonlinear stage. This dynamics is often attributed as the predator-prey oscillations induced by ZF collisional damping; however, similar dynamics is also observed in collisionless ZFs, in which case a different mechanism must be involved. Here, we propose a semi-analytic theory that explains the transition between the oscillations and saturation of collisionless ZFs within the quasilinear Hasegawa-Mima model. By analyzing phase-space trajectories of DW quanta (driftons) within the geometrical-optics (GO) approximation, we argue that the parameter that controls this transition is N ~ \gamma_MI/\omega_DW, where \gamma_MI is the MI growth rate and \omega_DW is the linear DW frequency. We argue that at N << 1, ZFs oscillate due to the presence of so-called passing drifton trajectories, and we derive an approximate formula for the ZF amplitude as a function of time in this regime. We also show that at N >~ 1, the passing trajectories vanish and ZFs saturate monotonically, which can be attributed to phase mixing of higher-order sidebands. A modification of N that accounts for effects beyond the GO limit is also proposed. These analytic results are tested against both quasilinear and fully-nonlinear simulations. They also explain the earlier numerical results by Connaughton et al. [J. Fluid Mech. 654, 207 (2010)] and Gallagher et al. [Phys. Plasmas 19, 122115 (2012)] and offer a revised perspective on what the control parameter is that determines the transition from the oscillations to saturation of collisionless ZFs.
Schwartz, Jacob; Emdee, Eric; Goldston, Robert; Jaworski, Michael
The lithium vapor box divertor is a potential solution for power exhaust in toroidal confinement devices. The divertor plasma interacts with a localized, dense cloud of lithium vapor, leading to volumetric radiation, cooling, recombination, and detachment. To minimize contamination of the core plasma, lithium vapor is condensed on cool (300°C to 400°C) baffles upstream of the detachment point. Before implementing this in a toroidal plasma device with a slot divertor geometry, we consider an experiment with a scaled baffled-pipe geometry in the high-power linear plasma device Magnum-PSI. Three 15 cm-scale open cylinders joined by 6 cm diameter ‘nozzles’ are positioned on the plasma beam axis upstream of a target. The central box may be loaded with several tens of grams of lithium, which can be evaporated at 650°C to produce a vapor predicted, using a simple plasma-neutral interaction model, to be dense enough to cause volumetric detachment in the plasma. The power delivered to the target and box walls as measured by increases in their temperatures after a 10 s plasma pulse can be compared to determine the effectiveness of the vapor in detaching the plasma. Direct Simulation Monte Carlo simulations are performed to estimate the flow rates of lithium vapor between the boxes and to estimate the trapping of H2 delivered by the plasma in the boxes, which could inadvertently lead to detachment. Details of the geometry, simulations, and possible diagnostic techniques are presented.
This work continues a series of papers where we propose an algorithm for quasioptical modeling of electromagnetic beams with and without mode conversion. The general theory was reported in the first paper of this series, where a parabolic partial differential equation was derived for the field envelope that may contain one or multiple modes with close group velocities. Here, we present a corresponding code PARADE (PAraxial RAy DEscription) and its test applications to single-mode beams in vacuum and also in inhomogeneous magnetized plasma. The numerical results are compared, respectively, with analytic formulas from Gaussian-beam optics and also with cold-plasma ray tracing. Quasioptical simulations of mode-converting beams are reported in the next, third paper of this series.
This work continues a series of papers where we propose an algorithm for quasioptical modeling of electromagnetic beams with and without mode conversion. The general theory was reported in the first paper of this series, where a parabolic partial differential equation was derived for the field envelope that may contain one or multiple modes with close group velocities. In the second paper, we presented a corresponding code PARADE (PAraxial RAy DEscription) and its test applications to single-mode beams. Here, we report quasioptical simulations of mode-converting beams for the first time. We also demonstrate that PARADE can model splitting of two-mode beams. The numerical results produced by PARADE show good agreement with those of one-dimensional full-wave simulations and also with conventional ray tracing (to the extent that one-dimensional and ray-tracing simulations are applicable).
A new model of heating, current drive, torque and other effects of neutral beam injection on NSTX-U that uses neural networks has been developed. The model has been trained and tested on the results of the Monte Carlo code NUBEAM for the database of experimental discharges taken during the first operational campaign of NSTX-U. By projecting flux surface quantities onto empirically derived basis functions, the model is able to efficiently and accurately reproduce the behavior of both scalars, like the total neutron rate and shine through, and profiles, like beam current drive and heating. The model has been tested on the NSTX-U real-time computer, demonstrating a rapid execution time orders of magnitude faster than the Monte Carlo code that is well suited for the iterative calculations needed to interpret experimental results, optimization during scenario development activities, and real-time plasma control applications. Simulation results of a proposed design for a nonlinear observer that embeds the neural network calculations to estimate the poloidal flux profile evolution, as well as effective charge and fast ion diffusivity, are presented.