Current work

Deputy Task Force Leader for the MST1 EUROfusion Workpackage

PhD Thesis

Tokamak electron heat transport by direct numerical simulation of small scale turbulence

In a fusion machine, understanding plasma turbulence, which causes a degradation of the measured energy confinement time, would constitute a major progress in this field. In tokamaks, the measured ion and electron thermal conductivities are of comparable magnitude. The possible sources of turbulence are the temperature and density gradients occurring in a fusion plasma. Whereas the heat losses in the ion channel are reasonably well understood, the origin of the electron losses is more uncertain. In addition to the the radial velocity associated to the fluctuations of the electric field, electrons are more affected than ions by the magnetic field fluctuations. In experiments, the confinement time can be conveniently expressed in terms of dimensionless parameters. Although still somewhat too imprecise, these scaling laws exhibit strong dependencies on the normalized pressure, β or the normalized Larmor radius, ρ*.

The present thesis assesses whether a tridimensional, electromagnetic, nonlinear fluid model of plasma turbulence driven by a specific instability can reproduce the dependence of the experimental electron heat losses on the dimensionless parameters β and ρ*. The investigated interchange instability is the Electron Temperature Gradient driven one (ETG). The model is built by using the set of Braginskii equations. The developed simulation code is global in the sense that a fixed heat flux is imposed at the inner boundary, leaving the gradients free to evolve.

From the nonlinear simulations, we have put in light three characteristics for the ETG turbulence: the turbulent transport is essentially electrostatic; the potential and pressure fluctuations form radially elongated cells called streamers; the transport level is very low compared to the experimental values.

The thermal transport dependence study has shown a very small role of the normalized pressure, which is in contradiction with the Ohkawa's formula. On the other hand, the crucial role of the electron normalized Larmor has been emphasized: the confinement time is inverse proportional to this parameter. Finally, the low dependence of turbulent transport with the magnetic shear and the inverse aspect ratio is also reported.

Although the transport level observed in the simulations is low compared to the experiments, we have tried a direct confrontation with Tore Supra results. This tokamak is well designed to study the electron heat transport. Keeping most of the parameters from a well referenced Tore Supra shot, the nonlinear simulation gives a threshold quite close to the experimental one. The observed turbulent conductivity is a factor fifty lower than the experimental one. An important parameter can not be matched: the normalized Larmor radius, ρ*.

This limitation has to be overcome in order to confirm this results. Finally, a rigorous confrontation between this result and gyrokinetic simulations has to conclude that the ETG instability cannot describe electron heat losses in tokamaks.

Keywords: thermonuclear fusion, tokamak, plasma, ETG turbulence, numerical simulations