Investigation of the density shoulder formation by using self-consistent simulations of plasma turbulence and neutral kinetic dynamics

Simulations of plasma turbulence in a lower single-null magnetic configurations are presented. The plasma dynamics is modelled by the drift-reduced two-fluid Braginskii equations that are coupled with a kinetic model for a single neutral species that considers ionization, charge-exchange, recombination and elastic collisions. The effect of increased core fuelling on the plasma scrape-off layer (SOL) density and temperature profile is investigated. The increase in core fuelling leads to an increase of the e-folding length in the near SOL and, in the far SOL, to an increase of the plasma density. These results are in agreement with experimental measurements, and in particular with the observations of the formation of a density shoulder in high-fuelling scenarios. The physical mechanisms underlying the increase of the far SOL density are analysed comparing parallel and perpendicular fluxes in the SOL and considering also simulations with similar parameters but without neutrals. Despite the increase of the blob size, the increase of the far SOL density observed at high-fuelling rates is found to be mainly caused by the strong decrease of parallel transport, due to the cooling of electrons resulting from ionization events.

Electron Cyclotron Heating and Suprathermal Electron Dynamics in the TCV Tokamak

Silvano Gnesin

This thesis is concerned with the physics of suprathermal electrons in thermonuclear, magnetically confined plasmas. Under a variety of conditions, in laboratory as well as space plasmas, the electron velocity distribution function is not in thermodynamic equilibrium owing to internal or external drives. Accordingly, the distribution function departs from the equilibrium Maxwellian, and in particular generally develops a high-energy tail. In tokamak plasmas, this occurs especially as a result of injection of high-power electromagnetic waves, used for heating and current drive, as well as a result of internal magnetohydrodynamic (MHD) instabilities. The physics of these phenomena is intimately tied to the properties and dynamics of this suprathermal electron population. This motivates the development of instrumental apparatus to measure its properties as well as of numerical codes to simulate their dynamics. Both aspects are reflected in this thesis work, which features advanced instrumental development and experimental measurements as well as numerical modeling. The instrumental development consisted of the complete design of a spectroscopic and tomographic system of four multi-detector hard X-ray (HXR) cameras for the TCV tokamak. The goal is to measure bremsstrahlung emission from suprathermal electrons with energies in the 10-300 keV range, with the ultimate aim of providing the first full tomographic reconstruction at these energies in a noncircular plasma. In particular, suprathermal electrons are generated in TCV by a high-power electron cyclotron heating (ECH) system and are also observed in the presence of MHD events, such as sawtooth oscillations and disruptive instabilities. This diagnostic employs state-of-the-art solid-state detectors and is optimized for the tight space requirements of the TCV ports. It features a novel collimator concept that combines compactness and flexibility as well as full digital acquisition of the photon pulses, greatly enhancing its potential for full spectral analysis in high-fluency scenarios. Additional flexibility is afforded by the possibility to rotate the orientation of two of the cameras, permitting the crucial comparison of radiation emitted perpendicular and parallel to the primary magnetic field. The design of the HXR system was optimized through an extensive iterative simulation process with the aid of tomographic reconstruction codes as well as quasilinear Fokker-Planck modeling of ECH-driven TCV plasmas. In parallel, the selection of the detectors for this system was performed through comprehensive laboratory testing of several candidate detectors available on the market. While the design was completed in the course of the thesis work, commissioning of the system has only commenced recently with one of the four cameras installed on TCV. The first preliminary results, discussed in the last part of this thesis, include basic parameter scans of ECH wave-plasma interaction and the investigation of the dynamic response of suprathermal electrons to modulated ECH. In addition, the cameras possess the novel ability to discriminate against very high-energy γ-ray radiation that cannot be collimated and must thus be excluded from spatial distribution analysis. A basic study of the conditions for γ-ray suppression was conducted in preparation for future experiments. The Fokker-Planck modeling tool used in this diagnostic development was acquired through a collaboration with CEA-Cadarache, initially with the primary motivation of studying the simultaneous plasma heating by 2nd and 3rd harmonic electron cyclotron waves that is uniquely possible on TCV. This motivated a dedicated study, both theoretical and experimental, of one particular instance of this combined heating, which became a second primary subject of this thesis work. The particular scenario studied here is one in which a single ECH frequency is resonant at both harmonics in the same plasma. The primary objective of this study was to determine whether a synergy effect existed, permitting an enhancement of the intrinsically weak 3rd harmonic absorption by the suprathermal electrons generated at the 2nd harmonic resonance. An associated question was whether this effect, if it existed, was experimentally measurable or was in fact observed in TCV. The simulations performed in the course of this study indeed predict the existence of such a synergy, although the answer to the second question was ultimately negative, at least within the current technical limitations. This study has proven nevertheless highly valuable in providing new insight into the complex velocity-space dynamics that govern ECH wave-particle interaction and suprathermal electron dynamics.

EPFL2011

Impact of neutral atoms on plasma turbulence in the tokamak edge region

Rogério Manuel Cabete De Jesus Jorge, Federico David Halpern, Paola Paruta, Paolo Ricci, Fabio Riva, Christoph Wersal

A novel first-principles self-consistent model that couples plasma and neutral atom physics suitable for the simulation of turbulent plasma behaviour in the tokamak edge region has been developed [1] and implemented in the GBS code [2]. While the plasma is modelled by the drift-reduced two fluid Braginskii equations, a kinetic model is used for the neutrals, valid in short and in long mean free path scenarios. The model includes ionization, charge-exchange, recombination, and elastic collisional processes. The model was used to study the transition from the sheath to the conduction limited regime by increasing the plasma density in the system. We compared the simulation results with the predictions of an expanded and refined two-point model, estimating the drop of electron and ion temperature along the magnetic field lines in the SOL. The model is now being applied to investigate the impact of neutrals on turbulent edge features, such as the broad shoulder observed in the far SOL at high plasma density. References [1] C. Wersal and P. Ricci 2015 Nucl. Fusion 55 123014 [2] P Ricci et al 2012 Plasma Phys. Control. Fusion 54 124047

2016

Self-consistent simulation of multi-component plasma turbulence and neutral dynamics in the tokamak boundary

André Calado Coroado

The interaction between neutral particles and the plasma plays a key role in determining the dynamics of the tokamak boundary that, in turn, significantly impact the overall performance of the device. Leveraging the work in Wersal and Ricci [Nucl. Fusion 55, 123014 (2015)], the present thesis describes the development and implementation in the GBS code of a mass-conserving multi-component self-consistent model to simulate the interplay between neutrals and plasma in the tokamak boundary. The simulation results are shown to be a useful tool to disentangle the physics at play in the tokamak boundary, and in particular for the interpretation of gas puff imaging (GPI).Developed in the last decade, the GBS code [Ricci et al, Plasma Phys. Control. Fusion 54, 124047 (2012)] allows for self-consistent three-dimensional numerical simulations of the turbulent plasma and neutral dynamics in the tokamak boundary. In GBS, a set of Braginskii equations in the drift limit describes the plasma time evolution, and the neutrals are modelled by solving a kinetic equation using the method of characteristics. While GBS enables simulations in arbitrary magnetic configurations, here we focus on limited plasmas. We first describe the geometrical operators and proper boundary conditions to ensure that mass conservation is satisfied. In comparison to the previous non-mass-conserving model of GBS, the mass-conserving simulations capture more accurately the sharp transition of the plasma and neutral quantities between the edge and scrape-off layer regions. In addition, we show that mass conserving simulations allow for reliable quantitative studies of particle fluxes in the tokamak boundary.A multi-component model of the neutral-plasma interaction is then developed by extending the single-component model to the description of a deuterium plasma that includes electrons, D+ ions, D atoms, D2 molecules and D2+ ions. The molecular dynamics is introduced through a set of drift-reduced Braginskii equations for the D2+ species and considering a kinetic equation for D2 molecules, in addition to the kinetic equation for D atoms, thus resulting in a coupled system of kinetic equations for the atomic and molecular neutral distribution functions. The first multi-component GBS simulations show that, in the sheath-limited regime under consideration, most of the D2 molecules cross the last closed flux surface (LCFS) and are dissociated or ionized in the edge region, thus giving rise to sources of D atoms inside the LCFS. This leads to an inward radial shift of the peak of the plasma source due to ionization of D atoms with respect to the single-component simulations.The multi-component model is applied to the simulation of GPI diagnostics, where the presence of a molecular gas puff is simulated self-consistently with the plasma and neutral dynamics of the tokamak boundary. The injected molecules interact with the boundary plasma, resulting in the emission of light in the D alpha wavelength that can be measured to infer the turbulent properties of the plasma. The simulated mechanisms underlying the light emission, which include the excitation of D atoms and dissociation of both D2 and D2+, provide a reliable tool for the interpretation of GPI experimental measurements. The impact of neutral fluctuations on the D alpha emission rate is investigated, as well as the correlation between the D alpha emission and the plasma and neutral quantities.