Drift-related transport and plasma-neutral interaction in the TCV divertor

Thermonuclear fusion of light atoms is the primary energy source of stars, such as our Sun, that led to the emergence of life on Earth. However, its economic exploitation as a virtually unlimited and clean energy source is yet to be developed. One of the most promising designs for a nuclear fusion reactor is the tokamak, a toroidal device that uses strong magnetic fields to contain matter so it may be heated to fusion-relevant temperatures of $\sim 10^8\;^\circ$ C. A substantial fraction of the heating power, required to maintain fusion-relevant plasma conditions, is transported across magnetic flux surfaces and channeled through a thin region of open magnetic field lines that surround the confined zone, the scrape-off layer, towards the plasma wall interface at the divertor targets. If left unmitigated, anticipated target heat fluxes in current reactor designs exceed the material limit and raise the need for an accurate predictive model of heat transport and power dissipation in the plasma edge to achieve a viable reactor configuration.

This doctoral thesis is devoted to further understanding of divertor power exhaust, a key challenge on the way towards fusion energy, and contributes to the validation of such a predictive model, the SOLPS-ITER code. This work employs stringent comparisons between simulation results and experiments on the TCV tokamak to validate this model, where possible, and/or indicate remaining problems. This thesis presents the first scrape-off layer simulations that fully account for drifts, currents, carbon impurities and kinetic neutrals, for the TCV tokamak, and thereby provides unprecedented insight into drift-driven transport and plasma-neutral interaction within its divertor.

In an initial stage, the SOLPS-ITER code was used to predict the effectiveness of divertor gas baffles to guide the first baffled TCV campaign. The experimental assessment in 2019 confirmed the essential simulation predictions. Drift simulations identify the diamagnetic current as the dominant cross-field contribution to the divertor charge balance. The resulting characteristics of target current profiles are tested and confirmed in TCV measurements. It is demonstrated, for the first time, that such electric currents give rise to a potential well below the magnetic X-point for the unfavorable magnetic field direction for H-mode access in highly resistive, i.e. low temperature, divertors. The simulation result is supported by reduced analytic models that highlight the underlying physics. The prediction is tested and confirmed in TCV experiments following novel diagnostic capabilities. The presence of such a potential well is shown to substantially reshape the divertor transport in detached divertor conditions. The simulations identify the $E\times B$ -drift as the dominant radial transport channel for particles, heat and momentum in the divertor. A stringent comparison of simulation results to a wide set of edge-relevant diagnostics yields typically excellent qualitative agreement, but quantitative differences remain: the simulations conclude a colder and denser divertor plasma state with an overestimated neutral pressure. The identified shortcomings in this approach will contribute to further development towards a predictive divertor model of present day, and future fusion devices.

Étude du transport d'énergie thermique dans les plasmas du tokamak à configuration variable au moyen de chauffage électronique cyclotronique

Yann Camenen

The development of controlled thermonuclear fusion, a quasi-unlimited energy source suitable for large scale electricity production, is one of the main goals of plasma physics research. Among the directions explored to date, the use of toroidal devices called tokamaks to create and confine hot plasmas using strong magnetic fields is particularly promising. The energy used to heat the plasma must remain well confined in order to achieve plasma temperatures higher than one hundred millon degrees for a sufficiently long period to obtain numerous fusion reactions. In the frame of efficient electricity production, maximising the energy confinement is also essential to achieve the required temperature with the lowest heating power. In tokamak plasmas, energy losses are mainly due to radiation and radial energy transport from the plasma core to the edge. A significant fraction of plasma physics research is therefore dedicated to the study of radial transport in tokamaks and the exploration of new operation regimes characterised by a low transport level. The development and optimisation of diagnostics used to observe plasmas is also part of this work. This thesis work, performed on the Tokamak à Configuration Variable (TCV) in Lausanne, covers the implementation and exploitation of a multi-channel soft X-ray detector with high spatial and temporal resolution, together with the development of the tomographic inversion routines used for data analysis. The detector, comprised of two superposed wire chambers, has been tested and calibrated using an X-ray source and then installed onto the tokamak. The position of the detector was chosen such as to observe the whole plasma cross-section with maximum spatial resolution leading to high quality tomographic inversions. A mobile absorber holder was installed between the plasma and the wire chambers. The energy range of the soft X-ray emission observed by the detector was thus chosen by selecting the appropriate absorber. These various features have made possible the use of the detector for numerous studies and in particular for the spatial and temporal characterisation of the plasma internal transport barrier formation. Plasma shaping abilities covering a wide range of plasma elongations and triangularities, including negative values, are one of the strengths of the TCV tokamak. For instance, plasmas with elongated cross-sections offer higher energy confinement as well as higher plasma current and pressure limits. However, the increase of the plasma vertical instability growth rate with elongation makes the vertical control of elongated plasmas difficult, in particular if the plasma current profile is too peaked. As the current profile is usually peaked for low plasma currents, current profile broadening is required there to achieve high elongation. During this thesis, a current profile broadening method based on temperature profile modification by localised EC heating has been studied in detail. The mechanism of this method has been documented and the optimal conditions for the EC power deposition determined. Using these conditions, the TCV operational space has been extended towards higher elongation at low current. The highest elongation obtained at low current has been increased by over 25% permitting the exploration of the plasma transport properties in this regime. The flexibility of the TCV EC heating system has also been used to investigate radial electron heat transport in L-mode plasmas. For the first time, the normalised temperature gradient has been varied by a factor of four and its influence on electron heat transport has been separated from that of the electron temperature. Electron heat transport increases strongly with the normalised temperature gradient, for values between 6 and 10, and then becomes independent of this parameter. In addition, electron heat transport increases with increasing electron temperature, decreasing density and increasing effective charge. The electron heat transport dependence on these three parameters can be cast as a single dependence on the plasma collisionality. TCV shaping abilities have then been used to test the influence of plasma triangularity. The main variations of the level of electron heat transport are described by a decrease of the electron heat diffusivity towards negative triangularity and high collisionality. At constant collisionality, electron heat transport is two times lower at a negative triangularity of –0.4 than at a positive triangularity of +0.4. Concerning micro-instabilities, gyro-fluid and gyro-kinetic simulations indicate that TEM and ITG instabilities are at play in these plasmas. The good qualitative agreement between the observed experimental dependencies and the predictions of simulations suggests strongly that the TEM instabilities are involved in the transport of electron heat. The experimental study provides dependable scaling of the electron heat transport on plasma parameters that can now be used to test the prediction of transport simulations. New elements such as the saturation of electron heat transport at high values of the normalised temperature gradient and the decrease of electron heat transport towards negative triangularities have been demonstrated.

EPFL2006

Gyrokinetic particle-in-cell global simulations of ion-temperature-gradient and collisionless-trapped-electron-mode turbulence in tokamaks

Sébastien Jolliet

The goal of thermonuclear fusion research is to provide power plants, that will be able to produce one gigawatt of electricity. Among the different ways to achieve fusion, the tokamak, based on magnetic confinement, is the most promising one. A gas is heated up to hundreds of millions of degrees and becomes a plasma, which is maintained – or confined – in a toroidal vessel by helical magnetic field lines. Then, deuterium and tritium are injected and fuse to create an α particle and an energetic neutron. In order to have a favorable power balance, the power produced by fusion reactions must exceed the power needed to heat the plasma and the power losses. This can be cast in a very simple expression which stipulates that the product of the density, the temperature and the energy confinement time must exceed some given value. Unfortunately, present-days tokamaks are not able to reach this condition, mostly due to plasma turbulence. The latter phenomenon enhances the heat losses and degrades the energy confinement time, which cannot be predicted by analytical theories such as the so-called neoclassical theory in which the heat losses are caused by Coulomb collisions. Therefore, numerical simulations are being developed to model plasma turbulence, mainly caused by the Ion and Electron Temperature-Gradient and the Trapped-Electron-Mode instabilities. The plasma is described by a distribution function which evolves according to the Vlasov equation. The electromagnetic fields created by the particles are self-consistently obtained through Maxwell's equations. The resulting Vlasov-Maxwell system is greatly simplified by using the gyrokinetic theory, which consists, through an appropriate ordering, of eliminating the fast gyromotion (compared to the typical frequency of instabilities). Nevertheless, it is still extremely difficult to solve this system numerically due to the large range of time and spatial scales to be resolved. In this thesis, the Vlasov-Maxwell system is solved in the electrostatic and collisionless limit with the Particle-In-Cell (PIC) ORB5 code in global tokamak geometry. This Monte-Carlo approach suffers from statistical noise which unavoidably degrades the quality of the simulation. Consequently, the first part of this work has been devoted to the optimization of the code with a view to reduce the numerical noise. The code has been rewritten in a new coordinate system which takes advantage of the anisotropy of turbulence, which is mostly aligned with the magnetic field lines. The overall result of the optimization is that for a given accuracy, the CPU time has been decreased by a factor two thousand, the total memory has been decreased by a factor ten and the numerical noise has been reduced by a factor two hundred. In addition, the scaling of the code with respect to plasma size is presently optimal, suggesting that ORB5 could compute heat transport for future fusion devices such as ITER. The second part of this thesis presents the validation of the code with numerical convergence tests, linear (including dispersion relations) and nonlinear benchmarks. Furthermore, the code has been applied to important issues in gyrokinetic theory. It is shown for the first time that a 5D global delta-f PIC code can achieve a thermodynamic steady state on the condition that some dissipation is present. This is a fundamental result as the main criticism against delta-f PIC codes is their inability to deal with long time simulations. Next, the role of the parallel nonlinearity is studied and it is demonstrated in this work that this term has no real influence on turbulence, provided the numerical noise is sufficiently low. This result should put an end to the controversy that recently occurred, in which gyrokinetic simulations using different numerical approaches yielded contradictory results. Finally, thanks to the optimization of the code, the gyrokinetic model has been extended to include the kinetic response of trapped-electrons, in place to the usual adiabatic (Boltzmann) approximation. For the first time, global TEM nonlinear simulations are presented, and the role of the zonal flow on heat transport is analyzed. This study will help in acquiring some knowledge on the less-known TEM turbulence (as compared to ITG). In conclusion, this thesis is one of the main steps of the development of ORB5, which is now a state-of-the-art gyrokinetic code for collisionless ITG and TEM turbulence, and has brought several contributions to the understanding of these phenomena.

EPFL2009

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