Particle transport in tokamak plasmas

The transport of particles in magnetically confined plasmas is of great importance for the development of fusion energy. It will determine techniques for fuelling, for controlling impurity concentrations and for the removal of the alpha particles produced by fusion reactions. The issues related to particle transport have received, until relatively recently, less attention than those regarding heat transport. Besides the greater experimental difficulty of measuring particle transport, this may explain why our understanding of this subject is still incomplete. The aim of this thesis is to document and tentatively interpret the experimental density profile behaviour in the TCV (Tokamak à Configuration Variable) and JET (Joint European Torus) tokamaks in the framework of different theoretical and semi-empirical models. The TCV tokamak is well suited to transport studies due to its extreme shaping capability, which allows the exploration of a wide range of different plasma conditions. This versatility is matched by the powerful and flexible electron cyclotron heating (ECH) system on TCV, which allows a control of the local power deposition profiles and current drive profiles. A study of particle transport on the JET tokamak has allowed us to compare the results of TCV to those of a much larger device and supplement the TCV study with the analysis in reactor relevant high confinement regime (H-mode). The experimental information was compiled into a database of density profiles in steady state, containing nearly 1000 samples for TCV and 600 samples for JET. The data analyzed covered a wide range of discharge conditions, including low confinement regime (L-mode) and H-mode discharges, ECH, ICH (ion cyclotron heating), LHCD (low hybrid), beam heated plasmas and include fully current drive discharges. The most relevant parameters which influence the density profiles were determined by regression. A detailed analysis of the particle sources showed that edge fuelling in TCV and JET cannot be responsible for density gradient in the plasma bulk, confirming the presence of particle convection or a 'pinch'. The existence of an anomalous pinch was unambiguously demonstrated both on JET and on TCV by the observation of peaked density profiles in stationary, fully relaxed, fully current driven discharges and hence in the absence of the neoclassical Ware pinch. An unexpected difference in the parameter dependencies was found in L-and H-modes. In TCV and JET, density gradient lengths (or profile peaking parameters) in L-mode were found to depend on magnetic shear with no dependence on collisionality. This lack of collisionality dependence in L- mode is inconsistent with current theoretical models. H-mode density profiles in JET, on the other hand, are clearly dependent on collisionality in agreement with theory and a prior observation on ASDEX-Upgrade, while exhibiting only a weak or no dependence on shear and temperature profiles. It was found that for TCV and JET, L mode density peaking can be interpreted as being due to turbulent equipartition, which assumes conservation of the magnetic moment and the longitudinal invariant during transport. The observation of a reduction of peaking in TCV with ECH supports drift wave turbulence theory, which predicts the appearance of outward particle convection, when trapped electron modes are destabilized. In JET H-mode, the weak secondary correlation of peaking with the electron-ion temperature ratio Te/Ti, may also be considered, at least quantitatively, as being supportive of drift wave turbulence theory.

Modelling of electron transport and of sawtooth activity in tokamaks

Transport phenomena in tokamak plasmas strongly limit the particle and energy confinement and represent a crucial obstacle to controlled thermonuclear fusion. Within the vast framework of transport studies, three topics have been tackled in the present thesis: first, the computation of neoclassical transport coefficients for general axisymmetric equilibria and arbitrary collisionality regime; second, the analysis of the electron temperature behaviour and transport modelling of plasma discharges in the Tokamak à Configuration Variable (TCV) ; third, the modelling and simulation of the sawtooth activity with different plasma heating conditions. The work dedicated to neoclassical theory has been undertaken in order to first analytically identify a set of equations suited for implementation in existing Fokker-Planck codes. Modifications of these codes enabled us to compute the neoclassical transport coefficients considering different realistic magnetic equilibrium configurations and covering a large range of variation of three key parameters: aspect ratio, collisionality, and effective charge number. A comparison of the numerical results with an analytical limit has permitted the identification of two expressions for the trapped particle fraction, capable of encapsulating the geometrical effects and thus enabling each transport coefficient to be fitted with a single analytical function. This has allowed us to provide simple analytical formulae for al1 the neoclassical transport coefficients valid for arbitrary aspect ratio and collisionality in general realistic geometry. This work is particularly useful for a correct evaluation of the neoclassical contribution in tokamak scenarios with large bootstrap current fraction, or improved confinement regimes with low anomalous transport and for the determination of the plasma current density profile, since the plasma conductivity is usually assumed neoclassical. These results have been included in the plasma transport code PRETOR. This code has been further extended and applied to the simulation of electron transport in TCV. In simulating the electron temperature profile of Ohmic sawtoothing plasmas, the proper description of the current density profile and the sawtooth activity play the dominant role and not the specific transport model, provided that a single parameter in the model is adjusted to match the global plasma performance. In TCV discharges with electron cyclotron heating (ECH), the behaviour of the electron temperature exhibits some characteristics which have been recently observed to be common to several tokamaks. In particular, with central heating the electron temperature profile is stiff outside the power deposition region, that is the gradient scale length is independent of the heating power and essentially constant along the minor radius. With off-axis heating, transport is strongly reduced in the central region of the plasma, whereas a steep increase of the heat conductivity is observed at the power deposition location. Although the semi-empirical Rebut-Lallia-Watkins (RLW) transport model does not involve a critical gradient scale length, as the experimental observations would suggest, rather a critical electron temperature gradient, we have shown that it allows simulations which reproduce the described experimental features with very good agreement. Due to the relatively low toroidal magnetic field of TCV, the experimental temperature gradient with ECH exceeds by far the threshold included in the model. It can thus be stated that the parametric dependence of the electron heat conductivity of this transport model is adequate to reproduce the electron transport for plasma parameters in the operation domain of TCV. PRETOR, interfaced with the experimental data and the code TORAY-GA for the computation of the ECH source, has hence been used as a reliable tool for transport analysis and planning of new experiments. This has contributed to the identification of an improved central electron confinement (ICEC) regime in TCV, characterized by a precise heating scheme (strong electron cyclotron current drive in the counter-direction in the plasma center, and localized off-axis heating), with a specific time sequence. Transport simulations and investigations of this regime, in particular dedicated to the reconstruction of the current density profile during the high performance phase, have motivated further experiments which have confirmed numerical predictions. As a consequence, the magnetic shear reversa1 has been identified as the crucial ingredient for ICEC. The powerful ECH system in TCV does not allow only strong global current profile modifications but also local tailoring which has significant effects on sawtooth activity. A model introduced for the prediction of the sawtooth period in the proposed International Thermonuclear Experimental Reactor (ITER) has been extended to be applicable to Ohmic and ECH discharges in TCV. The model has been found in agreement with the experimental observations and thereby we were able to identify the physical mechanisms which make ECH capable of controlling the sawtooth period. The parameter dependence of the relevant stability threshold has been found consistent with dedicated experiments demonstrating the effects of localized heating and current drive on the sawtooth period. The simulations have pointed out the effects of the magnetic shear and of the pressure gradient at the q = 1 surface. Moreover, the most efficient heating location to stabilize the sawtooth period has been identified as located outside the q = 1 surface before the sawtooth crash. The same period model has been used for the simulation of the sawtooth period in recent discharges with neutral beam injection (NBI) in the Joint European Torus (JET), performed to assess the role of beam ions on sawtooth stabilization. With the inclusion of an analytical expression for the fast ion contribution to the internal kink potential energy, validated by the hybrid kinetic/MHD code NOVA-K, the simulations have been found in remarkable good agreement with the experimental observations. This work has demontrated the role of beam ions in sawtooth stabilization and validated the stability threshold for the resistive internal kink which was predicted to be the sawtooth crash trigger relevant for ITER operation.

EPFL2001

É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

Analysis of plasma dynamic response to modulated electron cyclotron heating in TCV tokamak

Ilya Pavlov

The need of durable, economically acceptable and safe energy sources continues to stimulate studies in a field of thermonuclear fusion. The most successful solution for controlled magnetic fusion is the tokamak. The improvement of tokamak performance depends on the optimization of pressure and current density spatial distributions which can be modified by means of an auxiliary heating and a current drive. Electron cyclotron heating, in particular is a very important tool for the study and control of basic physical processes governing plasma confinement and stability, particularly because it allows the injection of highly localized intense power. ECH power deposition location plays a crucial role in sawtooth control and suppression, it is also important for tearing mode stabilization, and for implementation of closed loop systems for automatic control/suppression of magnetohydrodynamic activity. A part of the ECH power can be modulated (MECH), and used to identify where the ECH power has been deposited, and can also be useful in the experimental analysis of the electron transport in general. Nevertheless, despite the goal of MECH being a diagnostic and analysis tool, MECH can couple to plasma oscillations, such as sawteeth. MECH-sawtooth phase coupling adds significant complications in ECH deposition location and transport analysis, in some cases making the interpretations of results misleading. This is why it is important to get an insight into the phenomenon of MECH-sawtooth interaction, and to establish the boundaries where conventional types of modulation analysis can be successfully implemented. This thesis presents the analysis and interpretation of perturbative MECH experiments performed in the TCV tokamak with particular attention paid to the non-linear phase coupling of heat waves. TCV is equipped with a very flexible and high power ECH system. Two independent ECH systems permit MECH to be deposited at two different spatial locations, with two different types of modulating signals. This set-up was used to study simultaneous propagation of heat waves induced by MECH in non-sawtoothing plasmas, and in discharges with sawtooth activity. A new analysis method for the characterization of the plasma non-linear dynamic response to modulated heating was developed on the basis of Higher Order Spectral Analysis (HOSA) technique. This method applied to signals from different diagnostics, such as electron cyclotron emission and soft X-ray measurements, was extensively used to quantitatively characterize the effect of nonlinear phase coupling. In sawtooth free discharges a detailed analysis of the propagation of heat waves demonstrated that their phase coupling is solely related to properties of heat sources. It was demonstrated that if two heat waves are induced by non-coupled power sources (multi-beam MECH) then no phase coupling occurs. In the opposite case, when a source of perturbation (MECH) contains coupled harmonics, the corresponding heat waves demonstrate phase coupling. It was shown that these coupled heat waves loose their phase coherence while propagating in plasma. The dissipation of phase coupling is due to different phase velocities of heat waves and their diffusive damping. The new type of ECH power modulation accompanied with bicoherence analysis was proposed as a candidate for a reliable identification of EC power deposition location in a case of high frequency and low modulation depth MECH, including multi-beam heating. This type of MECH can be particularly important for real time control applications. In cases when MECH is applied to sawtoothing plasmas a direct experimental evidence of MECH-sawtooth non-linear phase coupling has been demonstrated using HOSA techniques, in particular bispectrum and bicoherence profiles. The detailed analysis presented here demonstrates a direct proof of periodic modification of sawtooth behavior by modulated ECH. It was shown that a simple diffusive model for the perturbed electron temperature with MECH source term can not be used for transport analysis in the presence of MECH-sawtooth coupling. Self-consistent time dependant transport simulations including a sawtooth relaxation model is capable of properly reproducing the effect of MECH-sawtooth phase coupling and can be used for transport analysis. On the basis of these simulations, an hypothesis for MECH-sawtooth coupling based on periodic modifications of magnetic shear was proposed. Additional information provided by HOS analysis enabled a correct interpretation of the plasma dynamic response to modulated heating in the presence of MECH-sawtooth coupling.

EPFL2008