Suprathermal electron studies in Tokamak plasmas by means of diagnostic measurements and modeling

To achieve reactor-relevant conditions in a tokamak plasma, auxiliary heating systems are required and can be realized by waves injected in the plasma that heat ions or electrons under certain conditions. Electron cyclotron resonant heating (ECRH) is a very flexible and robust technique featuring localized power deposition and current drive (CD) capabilities. Its fundamental principles such as damping on the cyclotron resonance are well understood and the application of ECRH is a proven and established tool; electron cyclotron current drive (ECCD) is regularly used to develop advanced scenarios and control magnetohydrodynamics (MHD) instabilities in the plasma by tailoring the current profile. There remain important open questions, such as the phase space dynamics, the observed radial broadening of the suprathermal electron distribution function (e.d.f.) and discrepancies in predicted and experimental CD efficiency. These are addressed in this thesis. One of its main goals is indeed to improve the understanding of wave-particle interaction in plasmas and current drive mechanisms. This was accomplished by combined experimental and numerical studies, strongly based on the conjunction of hard X-ray (HXR) bremsstrahlung measurements and Fokker-Planck modeling, characterizing the suprathermal electron population. The hard X-ray tomographic spectrometer (HXRS) diagnostic was purposely developed to perform these studies, in particular by investigating spatial HXR emission asymmetries in the co- and counter-current directions and within the poloidal plane. The system uses cadmium-telluride (CdTe) detectors and digital acquisition to store the complete time history of incoming photon pulses. An extensive study of digital pulse processing algorithms was performed and its consequent application allows the HXRS to handle high count rates in a noisy tokamak environment. Numerous other numerical tools were developed in the course of this thesis, among others to improve the time resolution by conditional averaging and to obtain local information with the general tomographic inversion (GTI) package. The interfaces of the comparatively new LUKE code and well-established CQL3D Fokker-Planck (F-P) code to the tokamak à configuration variable (TCV) data were refurbished and a detailed benchmarking of these two codes was performed for the first time. Indeed, the theory-predicted toroidal and poloidal emission asymmetries could be consistently verified by experiment and modeling in many cases, including scans of a variety of plasma and wave parameters. The effects of supra-thermal electron diffusion and radio frequency (RF) wave scattering, both resulting in a radial broadening of the HXR emission, were separated by a poloidal deposition location angle scan. Furthermore, previous results on anomalous diffusion and CD efficiency were reproduced with increased confidence arising from enhanced diagnostic specifications. The plasma response to electron cyclotron (EC) absorption and the role of quasi-linear effects were investigated using the coherent averaging capabilities of the HXRS. Several MHD instabilities can occur in the plasma center and better understanding of these modes and events is indispensable for their mitigation in order to prevent their negative effects on confinement and stability. Sawtooth crashes are such a major instability and localized at the q=1 surface. They can be described as the evolution of an internal m=1 kink mode leading to [...]

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

Pedestal Characteristics and MHD Stability of H-Mode Plasmas in TCV

Andreas Pitzschke

The tokamak à configuration variable (TCV) is unique in its ability to create a variety of plasma shapes and to heat the electron population in high density regimes using microwave power at the third harmonic of the electron cyclotron frequency. In the frame of this thesis, the impact of plasma shaping and heating on the properties of the edge transport barrier (ETB) in the high confinement mode (H-mode) was studied. This mode of operation is foreseen as one of the reference scenarios for ITER, the International Tokamak Experimental Reactor, which is being built to demonstrate the feasibility of thermonuclear fusion using magnetic confinement. A feature of H-mode regime operation are edge localized modes (ELMs), instabilities driven by the steep pressure gradients that form in the plasma edge region due to a transport barrier. During an ELM event, energy and particles are expelled from the plasma in a short burst. This will cause serious problems with respect to the heat load on plasma facing components in a tokamak of the size of ITER. Understanding of the phenomena associated with ELMs is thus required and dedicated investigations of their theory and experimental observations are carried out in many laboratories worldwide. This thesis presents several experimental and numerical investigations of tokamak behavior for configurations where the plasma edge plays an important role. From the experimental viewpoint, studies of transport barriers are challenging, as plasma parameters change strongly within a narrow spatial region. As part of the work presented here, the TCV Thomson scattering system was upgraded to meet the requirements for diagnosing electron temperature and density with high spatial resolution in the region of internal and external transport barriers. Simultaneously, the data analysis was significantly improved to cope with statistical uncertainties and alleviate eventual systematic errors. For measurements of the time evolution of density and temperature profile during the ELM cycle, the low repetition rate of the lasers used for Thomson scattering is a limiting. Although the system on TCV comprises 3 laser units that may be triggered in sequence with time separations down to 1 ms, time evolution over longer periods can only be reconstructed from repetitive events. In this context, an adjustment of the laser trigger to improve the synchronization with the ELM event is an advantage. A method was developed and implemented to generate a synchronizing trigger sequence, by a real-time monitoring of the D-alpha emission, which provides a marker for the ELM event. Recently, a "snowflake" (SF) divertor configuration, proposed as a possible solution to reduce the plasma-wall interaction by changing the divertor's poloidal magnetic field topology, was generated, for the first time, in TCV. A numerical code (KINX), based on a magnetohydrodynamic model (ideal MHD), was used to investigate the stability limits of this configuration under H-mode conditions and compare them with a similar standard single-null equilibrium. In a series of experiments, improved energy confinement was found and explained by improved stability of the edge region in the SF configuration. The influence of the pedestal structure in ELMy H-mode plasmas on the energy confinement and on ELM energy losses was investigated. The different ELM regimes found in TCV were analyzed, in particular the transition between type-III to type-I ELMs. The operational boundary of each ELM regime was characterized and verified by ideal MHD stability simulations for the ETB region. Recent studies on the scaling of the pedestal width with normalized poloidal pressure were confirmed. Using the capabilities of TCV, the influence of plasma shaping on pedestal parameters and MHD stability limits was investigated. In the past, models were developed to describe the onset of type-I ELMs, which are associated with modes in the ETB region arising from a coupling of pressure- and current-driven instabilities (coupled kink-ballooning modes). Experimental studies were performed to trace the temporal evolution of pedestal parameters characterizing the ETB during an ELM cycle. The results of these experiments were analyzed using information from MHD stability calculations. It is concluded that these models are capable of predicting limits as necessary conditions for ELM activity, but are not sufficient to fully explain ELM triggering.

EPFL2011

Extremely Shaped Plasmas to Improve the Tokamak Concept

Francesco Piras

Energy is essential for human existence and our future depends on plentiful and accessible sources of energy. The world population is fast growing and the average energy used per capita increases. One of the greatest challenges for human beings is that of meeting the growing demand for energy in a responsible, equitable and sustainable way. The possibility to obtain energy by "fusing" light atoms addresses these needs. Nuclear fusion reactions are clean, safe and the amount of fuel present on Earth (hydrogen isotopes) is practically inexhaustible and well distributed. Nuclear fusion is a natural process that occurs in all active stars like our Sun. Since the first demonstration of a deuterium fusion reaction (Rutherford 1933), researchers worldwide have tried to replicate this process on Earth by building a thermonuclear fusion reactor. Nevertheless, the challenge posed by the construction of a nuclear fusion reactor is greater than the one presented earlier by the development of a fission reactor. During the IAEA Conference in Geneva in the early 1958, L.A.Artsimovich declared: "Plasma physics is very difficult. Worldwide collaboration is needed for progress" and E.Teller, at the same conference: "Fusion technology is very complex. It is almost impossible to build a fusion reactor in this century". They were right. The extremely high temperature and density necessary to fuse hydrogen isotopes makes it difficult indeed to create a successful fusion reactor. Even though the physics of the fusion reaction appears clear, we are still facing problems on the road towards bulding the "box" that can efficiently confine the hot gas in the state of plasma. The best results so far have been obtained confining a plasma with strong magnetic fields in a toroidal configuration ("tokamak"). The Centre de Recherches en Physique des Plasmas in Switzerland actively studies this promising configuration towards the development of a nuclear fusion reactor. The experimental activity of the Tokamak à Configuration Variable (TCV) mainly focuses on the research of optimized plasma shapes capable of improving the global performance and solve the technological challenges of a tokamak reactor. Several theoretical and experimental results show the importance of the plasma shape in tokamaks. The maximum value of β (an indicator of the confinement efficiency) is for example related to the ratio between the height and the width of the plasma. The plasma shape can also affect the power necessary to access improved confinement regimes, as well as the plasma stability. This thesis reports on a contribution towards the optimization of the tokamak plasma shape. In particular, it describes the theoretical and experimental studies carried out in the TCV tokamak on two innovative plasma shapes: the doublet shaped plasma and the snowflake divertor. Doublet shaped plasmas have been studied in the past by the General Atomics group. Since then, the development of new plasma diagnostics and the discovery of new confinement regimes have given new reasons for interest in this unusual configuration. TCV is the only tokamak worldwide theoretically able to establish and control this configuration. This thesis illustrates new motivations for creating doublet plasmas. The vertical stability of the configuration is studied using a rigid model and the results are compared with those obtained with the KINX MHD stability code. The best strategy for controlling a doublet on TCV is also investigated, and a possible setup of the TCV control system is suggested for the doublet configuration. Analyzing the possible scenarios for doublet creation, the most promising scenario consists of the creation of two independent plasmas, which are subsequently merged to establish a doublet. For this reason, particular attention needs to be devoted to the problem of the plasma start-up. In this thesis, a general analysis of the TCV ohmic and assisted with ECH plasma start-up is presented, and recent attempts to create a doublet plasma are reported. Since the magnetic field reconstruction at the breakdown time is important to better diagnose these plasmas, the entire magnetic system of TCV has been calibrated with an original technique, also described in the manuscript. The last part of this thesis is devoted to the snowflake divertor configuration. This innovative plasma shape has been proposed and theoretically studied by Dr. D.D.Ryutov from the Lawrence Livermore National Laboratory. In Ryutov's articles, this configuration was proposed to alleviate the problems of the plasma-wall interaction and possibly affect the plasma edge stability. The TCV tokamak was the first to report the creation and control of a snowflake configuration, and the candidate was the principal investigator of this work. These results are accordingly discussed in this thesis. Details are provided in particular on the strategy used to establish the configuration. An edge-localized mode (ELM) H-mode regime, supported by electron cyclotron heating, has been successfully established in a snowflake. This regime exhibits 2 to 3 times lower ELM frequency but only a 20%-30% increase in normalized ELM energy (ΔWELM/WP ) compared to an identically-shaped, conventional, single-null, diverted H-mode. Enhanced stability of mid- to high-toroidal-mode-number ideal modes is consistent with the different snowflake ELM phenomenology. Finally, the capability of the snowflake to redistribute the edge power on the additional strike points has been confirmed experimentally and is also reported in this thesis.

EPFL2011

Pedestal Characteristics and MHD Stability of H-Mode Plasmas in TCV

Andreas Pitzschke

EPFL2011

Electron Cyclotron Heating and Suprathermal Electron Dynamics in the TCV Tokamak

Silvano Gnesin

EPFL2011

Extremely Shaped Plasmas to Improve the Tokamak Concept

Francesco Piras

EPFL2011