Anisotropy and shaping effects on the stability boundaries of infernal ideal MHD modes in tokamak hybrid plasmas

Anisotropy and some limiting toroidal flow effects on the stability of nearly resonant ideal magnetohydrodynamic modes in hybrid shaped tokamak plasmas are investigated within the ideal MHD infernal mode framework. Such effects are found to alter the plasma magnetic well/hill, which can be interpreted as imparing the average curvature, and the strength of mode coupling. In line with previous results, it is found that better stability properties are achieved through deepening the magnetic well by special cases of uniform toroidal flow and parallel plasma anisotropy. Plasma shaping provides additional modifications to the magnetic well depth, whose global stabilising or destabilising effect depends on the mutual interplay of elongation, triangularity and toroidicity. Further stabilisation is achieved by weakening the mode drive in vertically elongated plasmas.

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

MHD stability limits in the TCV tokamak

Holger Reimerdes

Magnetohydrodynamic (MHD) instabilities can limit the performance and degrade the confinement of tokamak plasmas. The Tokamak à Configuration Variable (TCV), unique for its capability to produce a variety of poloidal plasma shapes, has been used to analyse various instabilities and compare their behaviour with theoretical predictions. These instabilities are perturbations of the magnetic field which usually extend to the plasma edge where they can be detected with magnetic pick-up coils as magnetic fluctuations. A spatially dense set of magnetic probes, installed inside the TCV vacuum vessel, allows for a fast observation of these fluctuations. The structure and temporal evolution of coherent modes is extracted using several numerical methods. In addition to the setup of the magnetic diagnostic and the implementation of analysis methods, the subject matter of this thesis focuses on four instabilities which impose local and global stability limits. Al1 of these instabilities are relevant for the operation of a fusion reactor and a profound understanding of their behaviour is required in order to optimise the performance of such a reactor. Sawteeth, which are central relaxation oscillations common to most standard tokamak scenarios, have a significant effect on central plasma parameters. In TCV, systematic scans of the plasma shape have revealed a strong dependence of their behaviour on elongation κ, and triangularity δ, with high κ and low δ leading to shorter sawteeth with smaller crashes. This shape dependence is increased by applying central electron cyclotron heating (ECH), which increases or decreases the sawtooth period, depending on the plasma shape. The response to additional heating power is determined by the role of ideal or resistive MHD in triggering the sawtooth crash. For plasma shapes where additional heating and consequently, a faster increase of the central pressure shortens the sawteeth, the low experimental limit of the pressure gradient within the q = 1 surface is consistent with ideal MHD predictions. The observed decrease of this limit with elongation is also in qualitative agreement with ideal MHD theory. Edge localised modes (ELMs), occurring in TCV Ohmic high-confinement mode discharges, were observed to be preceded by coherent magnetic oscillations. The precursors prior to small ELMs, believed to be of type III, and prior to larger ELMs, previously referred to as TCV large ELMs, show the same characteristics, which allows for an identification of both ELMs to be of type III according to the usual classification scheme. The detected poloidal and toroidal mode structures are consistent with a resonant flux surface close to the plasma edge. Unlike conventional MHD modes, these precursors start at a random toroidal location and then grow in amplitude and toroidal extent until they encompass the whole toroidal circumference. Thus, the asymmetry causing and maintaining the toroidal localisation of the ELM precursor, must be intrinsic to the plasma. Soft X-ray measurements show that the localised precursor always coincides with a central m = 1 mode, which can usually be associated with the sawtooth pre- or postcursor mode. A comparison of the phases indicates a correlation with the maximum of the central mode preceding the toroidal location of the ELM precursor and, therefore, a hitherto unobserved coupling between central modes and ELMs. Highly elongated plasmas promise several advantages, among them higher current and beta limits. During TCV experiments dedicated to an increasing of the plasma elongation, a new disruptive current limit, at values well below the conventional current limit corresponding to qa >2, was encountered for κ >2.3. This limit, which is preceded by a kink-type mode, is found to be consistent with ideal MHD stability calculations. The TCV observations, therefore, provide the first experimental confirmation of a deviation of the linear Troyon-scaling of the ideal beta limit with normalised current at high elongation, which was predicted over 10 years ago. In addition to the ideal beta limit, several other MHD events are observed in highly elongated plasmas. The axisymmetric mode causes vertical displacement events and thereby, imposes a lower current limit. This operational limit is in good agreement with theoretical predictions of the growth rate of the axisymmetric mode. Furthermore, minor disruption are occasionally observed. They are caused by bursts of tearing modes located close to the q = 1 surface, which are destabilised by flat central current profiles, typical for highly elongated plasmas. Neoclassical tearing modes (NTMs), which have been observed to limit the achievable beta in a number of tokamaks, arise from a helical perturbation of the bootstrap current caused by an existing seed island. Neoclassical m/n = 2/1 tearing modes have been identified in TCV discharges which are characterised by a low electron collisionality νe*, a medium ion collisionality ν*, a medium value of beta and strongly peaked pressure and current profiles. In contrast to other tokamak experiments, where sawteeth, fishbones or ELMs generate the seed island, the required seed island is provided by a conventional tearing mode. The island clearly shows a conventional and a neoclassical growth phase. The TCV results provide the first clear observations of such a trigger mechanism which could also explain the occurrence of "triggerless" NTMs observed in other experiments. The slowly growing seed island also allows a measurement of the critical seed island width ωcrit, which is observed to increase with increasing density. In conclusion, several local and a global stability limits are analysed. These instabilities can limit the pressure gradient and thereby, the performance of the plasma. The presented results reveal several previously unobserved features of commonly observed instabilities. Since the most of the new observations can be explained by theory, they improve the predictive capabilities with respect to new experiments. The experiments have also shown some new interactions among different instabilities, which add to the already crowded complexities of MHD phenomena in fusion plasmas.

EPFL2001

Bootstrap current compensation with electron-cyclotron waves in 3D reactor-size configurations

Sergi Ferrando I Margalet

In magnetic confinement devices, the inhomogeneity of the confining magnetic field along a magnetic field line generates the trapping of particles (with low ratio of parallel to perpendicular velocities) within local magnetic wells. One of the consequences of the trapped particles is the generation of a current, known as the bootstrap current (BC), whose direction depends on the nature of the magnetic trapping. The BC provides an extra contribution to the poloidal component of the confining magnetic field. The variation of the poloidal component produces the alteration of the winding of the magnetic field lines around the flux surfaces quantified by the rotational transform ι. When ι reaches low rational values, it can trigger the generation of ideal MHD instabilities. Therefore, the BC may be responsible for the destabilisation of the configuration. This thesis is divided into two parts. In the first part, we present a self-consistent method to calculate the BC and assess its effect on equilibrium and stability in general 3D configurations. This procedure is applied to two reactor-size prototypes (both with plasma volumes ∼ 1000m3): a quasi-axisymmetric (QAS) system and a quasi helically symmetric (QHS) system with magnetic structures that develop BC in opposite directions. The BC increases with the plasma pressure, therefore its relevance is enhanced when dealing with reactor-level scenarios. The behaviour of both prototypes at reactor level values of β ≡ (kinetic plasma pressure)/(magnetic pressure) is assessed, as well as its alteration of the equilibrium and stability. In the QAS prototype, BC-consistent equilibria have been computed up to β = 6.7% and the configuration is shown to be stable up to β = 6.4%. Convergence of self-consistent BC calculations for the QHS case is achieved only up to β = 3.5%, but the configuration is unstable for β ≥ 0.6%. The relevance of symmetry breaking modes of the Fourier expansion of the confining magnetic field on the generation of BC is studied for each prototype. This proves the close relationship between magnetic structure and BC. Having established the potentially dangerous implication of the BC, principally, in reactor prototypes, a method to compensate its harmful effects is proposed in the second part of the thesis. It consists of the modelling of the current driven by externally launched ECWs within the plasma to compensate the effects of the BC. This method is flexible enough to allow the identification of the appropriate scenarios in which to generate the required CD depending on the nature of the confining magnetic field and the specific plasma parameters of the configuration. Both the BC and the CD calculations are included in a self-consistent scheme which leads to the computation of a stable BC+CD-consistent MHD equilibrium. This procedure is applied in this thesis to simulate the required CD to stabilise the QAS and QHS prototypes introduced in the first part. The estimation of the input power required and the effect of the driven current on the final equilibrium of the system is performed for several relevant scenarios and wave polarisations providing various options of stabilising driven currents. Several scenarios have been devised for each prototype in order to drive current at the appropriate location and with the desired direction. Different polarisations and launching conditions have been employed to this purpose. In particular, a HFS launched X2 ECW with an input power of 1.5MW has been shown to drive sufficient current to maintain the rotational transform below the critical value 2/3 at β = 6.4% for the QAS reactor. Correspondingly, in the QHS reactor, an X3-mode ECW of 100KW was sufficient to drive the current required to push the rotational transform below unity near the magnetic axis at β = 3%. Thus, stabilisation of BC-driven instabilities with externally launched ECWs has been achieved for both contrasting configurations. The method proposed in this thesis allows also the utilisation of EBW in the generation of CD. The possible advantages of EBCD for compensation studies are described as well as their possible application to the two prototypes under consideration. The BC+CD procedure is particularly interesting to investigate new magnetic geometries as potential candidates for fusion reactors. With this numerical tool, it is possible to assess the implications of their consistent BC when operating at reactor level. It also allows to quantify how much power would be required to maintain the system MHD stable in these circumstances. Nevertheless, this method is flexible enough to be applicable to any configuration.

EPFL2006