Discrepancy between modelled and measured radial electric fields in the scrape-off layer of divertor tokamaks: a challenge for 2D fluid codes?

Examination of radial electric field (E,.) profiles in the scrape-off layer (SOL) of ASDEX Upgrade (AUG) and JET revealed large discrepancies between 2D fluid edge modelling and experiment. Experimental profiles of plasma potential (V-p) in the outer (low field) side of the plasma, obtained with reciprocating Langmuir probes, decay radially with electron temperature, T-e, with the -eE(r)/del T-e ratio being > 1.5. In contrast, code simulated E-r are fairly low in most of the SOL (compared with -del T-e/e). Modelling with kinetic treatment of neutrals and drifts was performed using the SOLPS code for AUG cases and EDGE2D-Nimbus for JET cases. Mismatches between modelled and experimental E-r may be caused by the recently established tendency for the SOLPS code to underestimate T-e in the divertor of AUG. It was attributed to non-locality of parallel transport of supra-thermal, heat-carrying electrons originating upstream of the divertor, which are usually only weakly collisional and can penetrate, with few collisions, to the target. Ratios -eE(r)/del T-e obtained from the probe measurements in JET are of order 1.6, while in AUG, JT-60U and TCV they are of order 3. Such high values point to the possibility of fast electrons contributing, apart from target heat fluxes, also to the formation of the Debye sheath. The problem of the underestimation of E-r. in the codes must be closely related with the well-known problem of the underestimation of those parts of parallel ion flows in the SOL that are influenced by the toroidal field direction. It was demonstrated earlier that parallel ion flow at the outer midplane is dominated by the ion Pfirsch-Schluter flow, which in turn is partly driven by the radial electric field. The T-e and E-r discrepancies, as well as discrepancies between simulated and experimental parallel ion flows, raise a question of the validity of fluid codes for the plasma edge modelling and prompt the inclusion of kinetic effects into present-day 2D fluid codes which assume strong collisionality.

É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

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

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

Yann Camenen

EPFL2006

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

Andreas Pitzschke

EPFL2011

Extremely Shaped Plasmas to Improve the Tokamak Concept

Francesco Piras

EPFL2011