Search

In the present magnetically confined plasmas, the prediction of particle loading on material surfaces is a primary concern in view of the protection of plasma facing components for next step devices. Thus, an understanding of filament dynamics is needed. In this context, this work aims to develop an automatic detector for filaments arising in the MAST-U plasma. The identification of the filaments has been done starting from 2D images acquired with a fast visible camera. Therefore, it can be faced as an image object recognition problem. Currently, the object recognition is a key output of deep learning and machine learning algorithms. In this paper, a database of several thousands of images generated by a synthetic diagnostic, which reproduces the statistical properties of experimental filaments in terms of position, size and intensity has been used. The synthetic images are pre-processed by mapping them onto the toroidal midplane of the machine. Then a Faster R-CNN is customized to the problem of identifying the filaments. In particular, in order to enhance the performance of the detector, a suitable definition of the target-boxes defining the filament positions and sizes is adopted with good results.

A new tomographic inversion technique is presented for the identification of plasma filaments in wide-angle visible camera data. The technique works on the assumption that background subtracted images of filaments can be represented as a superposition of uniformly emitting magnetic equilibrium field lines. A large collection of equilibrium magnetic field lines is traced and projected onto the camera field of view and combined to form a geometry matrix describing the coordinate transformation from magnetic field aligned coordinates to image pixel coordinates. Inverting this matrix enables the reprojection of the emission in the camera images onto a field aligned basis, from which filaments are readily identifiable. The inversion is a poorly conditioned problem which is overcome using a least-squares approach with Laplacian regularization. Blobs are identified using the "watershed" algorithm and 2D Gaussians are fitted to get the positions, widths, and amplitudes of the filaments. A synthetic camera diagnostic generating images containing experimentally representative filaments is utilized to rigorously benchmark the accuracy and reliability of the technique. 74% of synthetic filaments above the detection amplitude threshold are successfully detected, with 98.8% of detected filaments being true positives. The accuracy with which filament properties and their probability density functions are recovered is discussed, along with sources of error and methods to minimize them.

Intrinsic Scrape Off Layer (SOL) instabilities are studied using flute approximation and incorporating the appropriate sheath boundary conditions at the target. The linear growth rate and the structure of the modes are obtained. The associated diffusion is estimated using a γ / k ⊥ 2 approach for the fastest growing modes. The model used includes curvature and sheath drives, finite Larmor radius effects and partial line tying at the target. The magnetic geometry is obtained using current carrying wires, representing the plasma current and the divertor coils, and naturally generates X-point geometry and magnetic shear effects. The calculation is performed for ITER relevant parameters and scans in SOL width and distance from the separatrix are presented. In addition to a standard Lower Single Null, Super-X and Snowflake configurations are examined in order to assess the importance of the geometry on the stability of the boundary plasma.

Using a kinetic model for the ions and adiabatic electrons, we solve a steady state, electron-repelling magnetic presheath in which a uniform magnetic field makes a small angle $\alpha \ll 1$ (in radians) with the wall. The presheath characteristic thickness is the typical ion gyroradius $\rho_{\text{i}}$. The Debye length $\lambda_{\text{D}}$ and the collisional mean free path of an ion $\lambda_{\text{mfp}}$ satisfy the ordering $\lambda_{\text{D}} \ll \rho_{\text{i}} \ll \alpha \lambda_{\text{mfp}}$, so a quasineutral and collisionless model is used. We assume that the electrostatic potential is a function only of distance from the wall, and it varies over the scale $\rho_{\text{i}}$. Using the expansion in $\alpha \ll 1$, we derive an analytical expression for the ion density that only depends on the ion distribution function at the entrance of the magnetic presheath and the electrostatic potential profile. Importantly, we have added the crucial contribution of the orbits in the region near the wall. By imposing the quasineutrality equation, we derive a condition that the ion distribution function must satisfy at the magnetic presheath entrance --- the kinetic equivalent of the Chodura condition. Using an ion distribution function at the entrance of the magnetic presheath that satisfies the kinetic Chodura condition, we find numerical solutions for the self-consistent electrostatic potential, ion density and flow across the magnetic presheath for several values of $\alpha$. Our numerical results also include the distribution of ion velocities at the Debye sheath entrance. We find that at small values of $\alpha$ there are substantially fewer ions travelling with a large normal component of the velocity into the wall., Comment: 69 pages, 20 figures

Four numerical codes are employed to investigate the dynamics of scrape-off layer filaments in tokamak relevant conditions. Experimental measurements were taken in the MAST device using visual camera imaging, which allows the evaluation of the perpendicular size and velocity of the filaments, as well as the combination of density and temperature associated with the perturbation. A new algorithm based on the light emission integrated along the field lines associated with the position of the filament is developed to ensure that it is properly detected and tracked. The filaments are found to have velocities of the order of 1 km s(-1), a perpendicular diameter of around 2-3 cm and a density amplitude 2-3.5 times the background plasma. 3D and 2D numerical codes (the STORM module of BOUT++, GBS, HESEL and TOKAM3X) are used to reproduce the motion of the observed filaments with the purpose of validating the codes and of better understanding the experimental data. Good agreement is found between the 3D codes. The seeded filament simulations are also able to reproduce the dynamics observed in experiments with accuracy up to the experimental errorbar levels. In addition, the numerical results showed that filaments characterised by similar size and light emission intensity can have quite different dynamics if the pressure perturbation is distributed differently between density and temperature components. As an additional benefit, several observations on the dynamics of the filaments in the presence of evolving temperature fields were made and led to a better understanding of the behaviour of these coherent structures.

New diagnostic, modelling and plant capability on the Mega Ampère Spherical Tokamak (MAST) have delivered important results in key areas for ITER/DEMO and the upcoming MAST Upgrade, a step towards future ST devices on the path to fusion currently under procurement. Micro-stability analysis of the pedestal highlights the potential roles of micro-tearing modes and kinetic ballooning modes for the pedestal formation. Mitigation of edge localized modes (ELM) using resonant magnetic perturbation has been demonstrated for toroidal mode numbers n = 3, 4, 6 with an ELM frequency increase by up to a factor of 9, compatible with pellet fuelling. The peak heat flux of mitigated and natural ELMs follows the same linear trend with ELM energy loss and the first ELM-resolved Ti measurements in the divertor region are shown. Measurements of flow shear and turbulence dynamics during L-H transitions show filaments erupting from the plasma edge whilst the full flow shear is still present. Off-axis neutral beam injection helps to strongly reduce the redistribution of fast-ions due to fishbone modes when compared to on-axis injection. Low-k ion-scale turbulence has been measured in L-mode and compared to global gyro-kinetic simulations. A statistical analysis of principal turbulence time scales shows them to be of comparable magnitude and reasonably correlated with turbulence decorrelation time. Te inside the island of a neoclassical tearing mode allow the analysis of the island evolution without assuming specific models for the heat flux. Other results include the discrepancy of the current profile evolution during the current ramp-up with solutions of the poloidal field diffusion equation, studies of the anomalous Doppler resonance compressional Alfvén eigenmodes, disruption mitigation studies and modelling of the new divertor design for MAST Upgrade. The novel 3D electron Bernstein synthetic imaging shows promising first data sensitive to the edge current profile and flows. © 2013 IAEA, Vienna.

The interaction between turbulence and a nonlinear tearing mode is investigated numerically using a 2D electrostatic model. Turbulence is found to cause transitions between the different roots for the propagation velocity of the mode. The transitions take the mode towards roots with slower propagation that are characterized by a locally flattened density profile. For sufficiently large islands the transition reduces the drive for the tearing mode.

Magnetic reconnection may occur as bursts of nonlinear plasma dynamics on the electron collisionless skin length scale de=c∕ωpe, during which a large fraction of the magnetic energy is converted to plasma thermal energy and plasma flow energy. An example of such a bursty energy release event is given with a simple set of electron Hall equations. The energization mechanism is the cross-field compression of the electron gas between interacting magnetic islands. The electron energization appears to be consistent with the rapid electron energy flux changes measured by the Cluster spacecraft crossing thin current sheets at −17RE in the geotail. The analysis is relevant to planned spacecraft missions for measuring electron scale magnetic reconnection events in the geomagnetic tail.

The influence of the ion polarization current on the stability of a constant-ψ magnetic island in a tokamak plasma is investigated numerically using a reduced two-fluid model in two-dimensional slab geometry. The polarization current is found to be negligibly small when the island is either too narrow or too wide. However, under certain circumstances, there exists an intermediate regime in which the polarization current is appreciable, and has a stabilizing influence on the island. This effect may account for the metastable nature of neoclassical tearing modes in tokamak plasmas.

The magnetic presheath is a boundary layer occurring when magnetized plasma is in contact with a wall and the angle $\alpha$ between the wall and the magnetic field $\vec{B}$ is oblique. Here, we consider the fusion-relevant case of a shallow-angle, $\alpha \ll 1$, electron-repelling sheath, with the electron density given by a Boltzmann distribution, valid for $\alpha / \sqrt{\tau+1} \gg \sqrt{m_{\text{e}}/m_{\text{i}}}$, where $m_{\text{e}}$ is the electron mass, $m_{\text{i}}$ is the ion mass, $\tau = T_{\text{i}}/ZT_{\text{e}}$, $T_{\text{e}}$ is the electron temperature, $T_{\text{i}}$ is the ion temperature, and $Z$ is the ionic charge state. The thickness of the magnetic presheath is of the order of a few ion sound Larmor radii $\rho_{\text{s}} = \sqrt{m_{\text{i}} \left(ZT_{\text{e}} + T_{\text{i}} \right) } / ZeB$, where $e$ is the proton charge and $B = |\vec{B}|$ is the magnitude of the magnetic field. We study the dependence on $\tau $ of the electrostatic potential and ion distribution function in the magnetic presheath by using a set of prescribed ion distribution functions at the magnetic presheath entrance, parameterized by $\tau$. The kinetic model is shown to be asymptotically equivalent to Chodura's fluid model at small ion temperature, $\tau \ll 1$, for $|\ln \alpha| > 3|\ln \tau | \gg 1$. In this limit, despite the fact that fluid equations give a reasonable approximation to the potential, ion gyro-orbits acquire a spatial extent that occupies a large portion of the magnetic presheath. At large ion temperature, $\tau \gg 1$, relevant because $T_{\text{i}}$ is measured to be a few times larger than $T_{\text{e}}$ near divertor targets of fusion devices, ions reach the Debye sheath entrance (and subsequently the wall) at a shallow angle whose size is given by $\sqrt{\alpha}$ or $1/\sqrt{\tau}$, depending on which is largest., Comment: 44 pages, 8 figures

A complete model of the dynamics of scrape-off layer filaments will be rather complex, including temperature evolution, three dimensional geometry and finite Larmor radius effects. However, the basic mechanism of $\boldsymbol{E}\times\boldsymbol{B}$ advection due to electrostatic potential driven by the diamagnetic current can be captured in a much simpler model; a complete understanding of the physics in the simpler model will then aid interpretation of more complex simulations, by allowing the new effects to be disentangled. Here we consider such a simple model, which assumes cold ions and isothermal electrons and is reduced to two dimensions. We derive the scaling with width and amplitude of the velocity of isolated scrape-off layer filaments, allowing for arbitrary elliptical cross-sections, where previously only circular cross-sections have been considered analytically. We also put the scaling with amplitude in a new and more satisfactory form. The analytical results are extensively validated with two dimensional simulations and also compared, with reasonable agreement, to three dimensional simulations having minimal variation parallel to the magnetic field., Comment: This is an author-created, un-copyedited version of an article accepted for publication/published in Plasma Physics and Controlled Fusion. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. The Version of Record is available online at http://dx.doi.org/10.1088/0741-3335/58/1/014030

This paper presents simulations of isolated 3D filaments in a slab geometry obtained using a 3D reduced fluid code. First, systematic scans were performed to investigate how the dynamics of a filament are affected by its amplitude, perpendicular size and parallel extent. The perpendicular size of the filament was found to have a strong influence on its motions, as it determined the relative importance of parallel currents to polarisation and viscous currents, whilst drift-wave instabilities were observed if the initial amplitude of the blob was increased sufficiently. Next, the 3D simulations were compared to 2D simulations using different parallel closures; namely, the sheath dissipation closure, which neglects parallel gradients, and the vorticity advection closure, which neglects the influence of parallel currents. The vorticity advection closure was found to not replicate the 3D perpendicular dynamics and overestimated the initial radial acceleration of all the filaments studied. In contrast, a more satisfactory comparison with the sheath dissipation closure was obtained, even in the presence of significant parallel gradients, where the closure is no longer valid. Specifically it captured the contrasting dynamics of filaments with different perpen- dicular sizes that were observed in the 3D simulations which the vorticity advection closure failed to replicate. However, neither closure successfully replicated the Boltzmann spinning effects and associated poloidal drift of the blob that was observed in the 3D simulations. Although the sheath dissipation closure was concluded to be more successful in replicating the 3D dynamics, it is emphasised that the vorticity closure may still be relevant for situations where the parallel current is inhibited from closing through the sheath due to effects such as increased resistivity near the targets.