Improving energy confinement in fusion plasmas by plasma shaping and current profile tailoring in the TCV tokamak

P. Angelino, A. Bottino¹, S. Brunner, Y. Camenen², A. Pochelon, O. Sauter, L. Villard
EPFL, CRPP, Association EURATOM - Confédération Suisse, Station 13, 1015 Lausanne

¹ present address: Max-Planck-Institut für Plasmaphysik, Postfach 1533, DE-85740 Garching bei München
² present address: International Institute for Fusion Science, University of Provence, Marseille, France

  Introduction The search for an abundant, sustainable, environmentally clean and climatically neutral energy source is more important than ever. Mankind is facing the challenge of a growing energy demand, particularly under the form of electricity, with a typical growth of a factor 4-5 expected from now to the end of the century [1, 2]. The fusion of light atomic elements like in the stars is an attractive option for several reasons, listed in box 1.

The most promising reaction is the fusion of deuterium and tritium (isotopes of hydrogen), which requires temperatures of the order of 10-20 keV (100-200 million °C) to bring the nuclei together. At these high temperatures, the gas is ionized, i.e. in a plasma state, macroscopically neutral but formed of electrically charged particles, which allows confining it with magnetic fields of a few Teslas. In magnetically confined fusion plasmas the density is very low and the pressure comparable to the pressure in a bicycle tire. For economic efficiency, it is important to have a confining "tire" that is not loosing too much pressure through heat and particle losses, keeping the plasma energy in "a well-insulating magnetic bottle". These losses have to be compensated by additional heating to maintain the required high temperatures. In standard plasmas, the losses are typically 50 times too high to be explained by inter-particle collisions alone and are due mainly to turbulence. But there are methods to reduce turbulence, as shown in the two next examples from the TCV tokamak experiment and numerical simulations in Lausanne.


Tokamak plasma shaping

The tokamak concept, invented in Russia in 1955, stands for "toroidalnaya kamera magnitnaya katushka", a toroidal chamber with magnetic coils. Tokamaks are used to create torus-like hot plasmas confined by magnetic fields. The plasma is embedded in a toroidal magnetic field (generated by a toroidal solenoid) and a poloidal field produced by a toroidal plasma current, itself induced by outside primary transformer coils. The superposition of these two fields forms a helical magnetic field, which confines the plasma. The magnetic helicity (or "rotational transform" of the field lines) is a key parameter for plasma confinement and was given the name of "safety factor" by early engineers.
Fusion energy basics
  • The most attractive reaction is the fusion of deuterium and tritium: D + T -> He + n + 17.6 MeV.
  • Deuterium is abundant, found in any water, thus well distributed geographically and practically inexhaustible. Doing the reaction with the deuterium found in one litre of water is equivalent to 300 litres of petrol.
  • Tritium is produced by breeding lithium, also abundant on Earth, using the fast neutrons (14 MeV) produced in the reaction above. Another comparison: the deuterium of your bath (45 l) and the lithium of the battery of your laptop is equivalent to 40 tonnes of petrol.
  • Interestingly, the reaction above is not a chain reaction. Thus, a “runaway” reaction and the resulting uncontrolled production of energy is not possible with fusion. Fusion reactions cannot be maintained uncontrollably: any disturbance or failure stops the reaction. This is why fusion is inherently safe.
  • Nuclear fusion reactors produce no high activity/long life radioactive waste. The "burnt" fuel is helium, a non-radioactive gas. Radioactive substances in the system are the fuel (tritium) and materials activated while the machine is running. The goal of the ongoing R&D programme is for fusion reactor material to be recyclable after 100 years.
  • Nuclear risks associated with fusion relate to the use of tritium, which is a radioactive form of hydrogen. However, the amount used is limited to few grams of tritium for the reaction and a few kilograms on site. During operation, the radiological impact of the use of tritium on the most exposed population is much smaller than that due to natural background radiation. For the reactor ITER, no accident scenario has been identified that would imply the need to take countermeasures to protect the surrounding population.
  • Risks of nuclear proliferation are extremely weak (no uranium, no plutonium, etc.). Fusion research was declassified in 1958 in the middle of the cold war.
  • Low CO2 emission.
Fig. 1). TCV facility. Left: view of the inside of the vacuum vessel covered by carbon tiles with ports for heating and diagnosing the plasma. Right: cut of the vacuum vessel with a schematics of the micro-wave heating beams at the 2nd and 3rd electron cyclotron harmonics, directed to a Dee-shaped plasma (i.e positive triangularity δ > 0)
Fig. 2). Effect of triangularity plasma shape on the poloidal cross-section of electrostatic potential cells from linear gyrokinetic calculations, for negative (left) and positive triangularity δ (right). This shows a reduction of the radial wavelength λ ~ 1/k at δ < 0, clearly identifiable in the midplane of left subplot, reducing transport [6]. The tokamak axis is on the left of the plasma cross-section, solid lines correspond to the plasma edge, dashed lines to the different magnetic surfaces in the plasma and n represents the toroidal mode number of the mode shown.
Fig. 3). Effect of current profile on the poloidal cross-section of the electrostatic potential for the n = 5 mode, in the case of
- a current density profile with a local off-axis maximum (indicated by the solid cyan line, left subplot) and
- a peaked, monotonic current density profile, from [9].
In the vicinity of the local off-axis maximum (left subplot), the electrostatic potential eddies are distorted and the radial extent of the mode is reduced, reducing transport.

These initial Russian tokamaks quickly showed excellent confinement properties but the results were not well known or accepted outside the USSR. In 1969, thanks to the development of high power lasers, the high electron temperature reached was confirmed by an international independent team using Thomson scattering. The news rapidly spread and led to the multiplication of tokamak experiments worldwide. In these initial tokamaks, the torus-like plasmas had a simple circular cross-section and the plasmas were heated Ohmically by the current flowing in the plasma.

The TCV device (Tokamak à Configuration Variable) at the Centre de Recherches en Physique des Plasmas (CRPP) at EPFL, which is part of the Association Suisse-EURATOM-Fusion, has been built to study the properties of shaped plasmas cross-sections, i.e. of different elongation, triangularity, and squareness [3]. TCV with its 16 independently supplied shaping coils, is a worldwide unique device for its flexibility in plasma shaping, see Fig.1. TCV for instance is the only tokamak that can generate negative triangularity plasma shapes (inverse Dee-shaped, symmetric to the plasma shape in Fig.1). In addition, the machine is equipped with an intense microwave power system [4], which allows heating the plasma through electron cyclotron resonance heating (ECRH). This heating method provides localized power deposition and is an important tool for controlling pressure and current profiles.

 

Improved heat confinement found at negative triangularity

Varying the plasma triangularity in TCV ECRH plasmas, it was found for the first time that the energy confinement time (the ratio of the energy in the plasma to the power used to heat it, a confinement figure of merit) nearly doubles from positive to negative triangularity δ [5]. ECRH, with its local power deposition property, enabled us to undertake local heat transport measurements over a large range of plasma parameters. This made it possible to separate the effects of plasma shape and collisionality on heat transport: transport decreases towards low triangularity and high collisionality [6]. The measurements show that heat transport is reduced by a factor up to two at mid-radius going from positive to negative triangularity. These experimental heat transport results are supported by linear and non-linear gyrokinetic simulations [7]. The micro-instabilities developing in these plasma conditions, the "trapped electron modes" (TEM), have essentially shorter radial wavelengths at negative triangularity, Fig.2a, compared to positive triangularity, Fig.2b. Recent non-linear calculations showed the relevance of such linear calculations in the specific case of TEM [8].

Let’s draw an analogy for the role of the radial wavelength on transport by comparing radial turbulent transport in tokamak plasma with the vertical transport in meteorology. In the summer, in the absence of wind, large vertical convection cells can develop – large cumulus - leading to vertical transport and finally storms. In the presence of wind, shearing the development of vertical cells, the vertical wavelength is reduced, leading to a reduction of the vertical transport: no storms are formed.

This analogy is suggestive and sheared flows are indeed known to reduce transport in tokamak plasmas. The full picture is however more complex and for instance the shearing of the magnetic field lines with radius, which by changing direction with radius, opposes like a braided net to radial convection, also has a beneficial effect on transport reduction.

Transport Barriers by current profile shaping

Internal transport barriers (ITB) are regions of reduced heat transport of energy and particles, which essentially allow for larger temperature and density gradients to develop while keeping the particles and the heat losses at a low level. The creation and control of ITBs are strongly related to the profile shape of the current density flowing along the torus, which determines the shearing of the magnetic field lines. In particular, ITBs tend to appear where the rotational transform profile has a local maximum (the cyan contour in Fig. 3a), At this radial location in the plasma, gyrokinetic simulations show that the radial size of the potential cells is largely reduced, Fig.3a, compared to the case of a monotonic rotational transform profile [9], Fig.3b. Appropriate current density profile tailoring can be used to control the position and strength of the ITB.

The TCV microwave power system can be employed for an accurate current density profile control through the generation of current from the waves, the so-called electron cyclotron current drive (ECCD). Thus, ITBs have been generated and studied in TCV in a variety of conditions [10], reaching electron temperatures up to 18 keV, typically one order of magnitude higher than the corresponding Ohmic condition. In TCV, ITBs can be generated in both electron temperature and density profiles.

ECCD was also used in TCV to produce for the first time in the world steady fully non-inductive discharges, in which all the current is driven by EC microwaves [11,12], an important step in the direction of the steady-state tokamak.

Conclusion

From the experience over the years in plasma shaping and current density profile tailoring, we have learned at TCV not only to develop interesting ways of improving confinement performance, but also – through comparison of modelling and experiment - found powerful tests of theoretical models of plasma stability and confinement [e.g. 13-15].

The comparison of plasma confinement, transport and, more recently, turbulence of TCV plasmas with results from gyrokinetic codes gives strong indications on the nature of turbulent transport. More crucially, this allows testing and improving the codes that will be applied to the computationally much more demanding cases of ITER, the International Thermonuclear Experimental Reactor in construction in the South of France [16] (see box 2).

There are various other issues that still need to be tackled in view of the development of a magnetic fusion reactor, such as the very large power loads on walls resulting from relaxation instabilities of edge plasma profiles. These aspects need skills from domains extending from plasma stability, transport, control, to materials. In this view, the flexibility in plasma shaping of TCV has allowed us to test for the first time a so-called "snow flake" divertor [17,18]. This new divertor concept has recently been proposed [19] with the aim of reducing heat loads on the first wall by spreading the power over a larger surface. This is an example of the various plasma geometries and concepts to test and develop, and which a machine like TCV can economically address owing to its relative small size and its flexibility in shaping and heating systems. Improved heating power systems, also delivering power to the ions instead of only to the electrons as presently, will open new fields of research relevant to reactor physics.

 

High power computation

The complexity of the systems under investigation implies that numerical computations are often the only way to make a quantifiable theoretical prediction. Stimulated by the fast growing performance of High Performance Computing (HPC) platforms, there has been in recent years a vigorous development in the development of codes devoted to fusion plasma physics. Among the still largely open and most challenging problems is the question of plasma turbulence. There has been dramatic improvement in first-principles based simulations of such phenomena. Thanks to the application of massively parallel algorithms that scale up to tens of thousands of processors computations of turbulence in the whole core of tokamaks such as TCV is now feasible. In order to be able to simulate the whole ITER plasma scalable computations up to hundreds of thousands of processors will be required. The CRPP is actively pursuing research in this field, notably in the frame of the HP2C initiative (Swiss Platform for High Performance High Productivity Computing).

 

References

[1] R. J. Goldston et al., Climate Change, Nuclear Proliferation and Fusion Energy, 23rd IAEA Fusion Energy Conf., Daejeon, Korea, Oct. 11-16, 2010, SEE/1-1a.
[2] L. Clarke et al., Int. climate policy architectures: Overview of the EMF 22 Int. Scenarios, Energy Economics 31 (2009) S64-S81.
[3] F. Hofmann, et al., Plasma Phys. Cont. Fusion 36 (1994) B277.
[4] T. P. Goodman and the TCV Team, Nucl. Fus. 48 (2008) 054011.
[5] A. Pochelon, T. P. Goodman, M. A. Henderson, et al., Nucl. Fusion 39 (1999) 1807.
[6] Y. Camenen et al., Nucl. Fusion 47 (2007) 510.
[7] A. Marinoni et al., Plasma Phys. Cont. Fusion 51 (2009) 055016.
[8] S. Jolliet, Thèse EPFL no. 4326 (2009).
[9] A. Bottino, O. Sauter at al., Plasma Phys. Cont. Fusion 48 (2006) 215–233.
[10] O. Sauter, S. Coda, T. P. Goodman, M. A. Henderson et al., Phys. Rev. Lett. 94 (2005) 105002.
[11] O. Sauter, et al., Phys. Rev. Lett. 84 (2000) 3325.
[12] S. Coda et al., Plasma Phys. Control. Fusion 42 (2000) B311.
[13] A. Fasoli for the TCV Team, Nucl. Fus. 49 (2009) 104005.
[14] A. Pochelon, Y. Camenen, A. Marinoni et al., Proc. of 22nd IAEA Fusion Energy Conference, Geneva, Switzerland, 12-18 Oct 2008, EX/P5-15.
[15] L. Villard et al., Plasma Phys. Cont. Fusion 52 (2010) 124038.
[16] ITER, International Thermonuclear Experimental Reactor
[17] F. Piras, et.al. Plasma Phys. Cont. Fusion 51 (2009) 055009.
[18] F. Piras et al., Phys. Rev. Lett. 105 (2010) 155003.
[19] D. D. Ryutov et al., Phys. Plasmas 15 (2008) 092501.

 

 

[Released: February 2012]