Sandpile modeling of pellet pacing in fusion plasmas.

Sandpile models have been used to provide simple phenomenological models without incorporating the detailed features of a fully featured model. The Chapman sandpile model [Chapman et al., Phys. Rev. Lett. 86, 2814 (2001)PRLTAO0031-900710.1103/PhysRevLett.86.2814] has been used as an analog for the behavior of a plasma edge, with mass loss events being used as analogs for edge-localized modes (ELMs). In this work we modify the Chapman sandpile model by providing for both increased and intermittent driving. We show that the behavior of the sandpile, when continuously fuelled at very high driving, can be determined analytically by a simple algorithm. We observe that the size of the largest avalanches is better reduced by increasing constant driving than by the intermittent introduction of "pellets" of sand. Using the sandpile model as a reduced model of ELMing behavior, we conject that ELM control in a fusion plasma may similarly prove more effective with increased total fuelling than with pellet addition.

Energetic Geodesic Acoustic Modes Associated with Two-Stream-like Instabilities in Tokamak Plasmas.

An unstable branch of the energetic geodesic acoustic mode (EGAM) is found using fluid theory with fast ions characterized by their narrow width in energy distribution and collective transit along field lines. This mode, with a frequency much lower than the thermal GAM frequency ω_{GAM}, is now confirmed as a new type of unstable EGAM: a reactive instability similar to the two-stream instability. The mode can have a very small fast ion density threshold when the fast ion transit frequency is smaller than ω_{GAM}, consistent with the onset of the mode right after the turn-on of the beam in DIII-D experiments. The transition of this reactive EGAM to the velocity gradient driven EGAM is also discussed.

A multichannel magnetic probe system for analysing magnetic fluctuations in helical axis plasmas.

The need to understand the structure of magnetic fluctuations in H-1NF heliac [S. Hamberger et al., Fusion Technol. 17, 123 (1990)] plasmas has motivated the installation of a sixteen former, tri-axis helical magnetic probe Mirnov array (HMA). The new array complements two existing poloidal Mirnov arrays by providing polarisation information, higher frequency response, and improved toroidal resolution. The helical placement is ideal for helical axis plasmas because it positions the array as close as possible to the plasma in regions of varying degrees of favourable curvature in the magnetohydrodynamic sense, but almost constant magnetic angle. This makes phase variation with probe position near linear, greatly simplifying the analysis of the data. Several of the issues involved in the design, installation, data analysis, and calibration of this unique array are presented including probe coil design, frequency response measurements, mode number identification, orientation calculations, and mapping probe coil positions to magnetic coordinates. Details of specially designed digitally programmable pre-amplifiers, which allow gains and filters to be changed as part of the data acquisition initialisation sequence and stored with the probe signals, are also presented. The low shear heliac geometry [R. Jiménez-Gómez et al., Nucl. Fusion 51, 033001 (2011)], flexibility of the H-1NF heliac, and wealth of information provided by the HMA create a unique opportunity for detailed study of Alfvén eigenmodes, which could be a serious issue for future fusion reactors.

Minimally constrained model of self-organized helical states in reversed-field pinches.

We show that the self-organized single-helical-axis (SHAx) and double-axis (DAx) states in reversed field pinches can be reproduced in a minimally constrained equilibrium model using only five parameters. This is a significant reduction on previous representations of the SHAx which have required an infinite number of constraints. The DAx state, which has a nontrivial topology, has not previously been reproduced using an equilibrium model that preserves this topological structure. We show that both states are a consequence of transport barrier formation in the plasma core, in agreement with experimental results. We take the limit of zero pressure in this work, although the model is also valid for finite pressure.

Using Bayesian analysis and Gaussian processes to infer electron temperature and density profiles on the Mega-Ampere Spherical Tokamak experiment.

A unified, Bayesian inference of midplane electron temperature and density profiles using both Thomson scattering (TS) and interferometric data is presented. Beyond the Bayesian nature of the analysis, novel features of the inference are the use of a Gaussian process prior to infer a mollification length-scale of inferred profiles and the use of Gauss-Laguerre quadratures to directly calculate the depolarisation term associated with the TS forward model. Results are presented from an application of the method to data from the high resolution TS system on the Mega-Ampere Spherical Tokamak, along with a comparison to profiles coming from the standard analysis carried out on that system.

The use of bayesian inversion to resolve plasma equilibrium.

Recently, bayesian probability theory has been used at a number of experiments to fold uncertainties and interdependencies in the diagnostic data and forward models, together with prior knowledge of the state of the plasma, to increase accuracy of inferred physics variables. A new probabilistic framework, MINERVA, based on bayesian graphical models, has been used at JET and W7-AS to yield predictions of internal magnetic structure. A feature of the framework is the bayesian inversion for poloidal magnetic flux without the need for an explicit equilibrium assumption. Building on this, we discuss results from a new project to develop bayesian inversion tools that aim to (1) distinguish between competing equilibrium theories, which capture different physics, using the MAST spherical tokamak, and (2) test the predictions of MHD theory, particularly mode structure, using the H-1 Heliac. Specifically, we report on correction of the motional Stark effect, pickup coils, flux-loop constrained bayesian inferred equilibrium for varying toroidal flux.

A high resolution Mirnov array for the Mega Ampere Spherical Tokamak.

Over the past two decades, the increase in neutral-beam heating and alpha particle production in magnetically confined fusion plasmas has led to an increase in energetic particle driven mode activity, much of which has an electromagnetic signature which can be detected by the use of external Mirnov coils. Typically, the frequency and spatial wave number band of such oscillations increase with increasing injection energy, offering new challenges for diagnostic design. In particular, as the frequency approaches the megahertz range, care must be taken to model the stray capacitance of the coil, which limits the resonant frequency of the probe; model transmission line effects in the system, which if unchecked can produce system resonances; and minimize coil conductive shielding, so as to minimize skin currents which limit the frequency response of the coil. As well as optimizing the frequency response, the coils should also be positioned to confidently identify oscillations over a wide wave number band. This work, which draws on new techniques in stray capacitance modeling and coil positioning, is a case study of the outboard Mirnov array for high-frequency acquisition in the Mega Ampere Spherical Tokamak, and is intended as a roadmap for the design of high frequency, weak field strength magnetic diagnostics.

Laboratory evidence for stochastic plasma-wave growth.

The first laboratory confirmation of stochastic growth theory is reported. Floating potential fluctuations are measured in a vacuum arc centrifuge using a Langmuir probe. Statistical analysis of the energy density reveals a lognormal distribution over roughly 2 orders of magnitude, with a high-field nonlinear cutoff whose spatial dependence is consistent with the predicted eigenmode profile. These results are consistent with stochastic growth and nonlinear saturation of a spatially extended eigenmode, the first evidence for stochastic growth of an extended structure.

Statistics of polarization and Stokes parameters of stochastic waves.

Several theories now exist to describe the probability distribution functions (PDFs) for the electric field strength, intensity, and power of signals. In this work, a model is developed for the PDFs of the polarization properties of the superposition of multiple transverse wave populations. The polarization of each transverse wave population is described by a polarization ellipse with fixed axial ratio and polarization angle, and PDFs for the field strength and phase. Wave populations are vectorially added, and expressions found for the Stokes parameters I , U , Q , and V , as well as the degrees of linear and circular polarization, and integral expressions for their statistics. In this work, lognormal distributions are chosen for the electric field, corresponding to stochastic growth, and polarization PDFs are numerically calculated for the superposition of orthonormal mode populations, which might represent the natural modes emitted by a source. Examples are provided of the superposition of linear, circular, and elliptically polarized wave populations in cases where the component field strength PDFs are the same, and where one field strength PDF is dominant.

Plasma instability of a vacuum arc centrifuge.

Ever since conception of the vacuum arc centrifuge in 1980, periodic fluctuations in the ion saturation current and floating potential have been observed in Langmuir probe measurements in the rotation region of a vacuum arc centrifuge. In this work we develop a linearized theoretical model to describe a range of instabilities in the vacuum arc centrifuge plasma column, and then test the validity of the description through comparison with experiment. We conclude that the observed instability is a "universal" instability, driven by the density gradient, in a plasma with finite conductivity.