@phdthesis{Kauffmann2012,
  author    = {Karla Kauffmann},
  title     = {Including Collisions in Gyrokinetic Tokamak and Stellarator Simulations},
  journal    = {Einbeziehung von St{\"o}{\"s}en in Gyrokinetische Simulationen von Tokamak und Stellarator Plasmen},
  url       = {https://nbn-resolving.org/urn:nbn:de:gbv:9-001221-0},
  year      = {2012},
  abstract  = {Particle and heat transport in fusion devices often exceed the neoclassical prediction. This anomalous transport is thought to be produced by turbulence caused by microinstabilities such as ion and electron-temperature-gradient (ITG/ETG) and trapped-electron-mode (TEM) instabilities, the latter ones known for being strongly influenced by collisions. Additionally, in stellarators, the neoclassical transport can be important in the core, and therefore investigation of the effects of collisions is an important field of study. Prior to this thesis, however, no gyrokinetic simulations retaining collisions had been performed in stellarator geometry. In this work, collisional effects were added to EUTERPE, a previously collisionless gyrokinetic code which utilizes the \δ f method. To simulate the collisions, a pitch-angle scattering operator was employed, and its implementation was carried out following the methods proposed in [Takizuka \& Abe 1977, Vernay Master's thesis 2008]. To test this implementation, the evolution of the distribution function in a homogeneous plasma was first simulated, where Legendre polynomials constitute eigenfunctions of the collision operator. Also, the solution of the Spitzer problem was reproduced for a cylinder and a tokamak. Both these tests showed that collisions were correctly implemented and that the code is suited for more complex simulations. As a next step, the code was used to calculate the neoclassical radial particle flux by neglecting any turbulent fluctuations in the distribution function and the electric field. Particle fluxes in the neoclassical analytical regimes were simulated for tokamak and stellarator (LHD) configurations. In addition to the comparison with analytical fluxes, a successful benchmark with the DKES code was presented for the tokamak case, which further validates the code for neoclassical simulations. In the final part of the work, the effects of collisions were investigated for slab and toroidal ITGs and TEMs in a tokamak configuration. The results show that collisions reduce the growth rate of slab ITGs in cylinder geometry, whereas they do not affect ITGs in a tokamak, which are mainly curvature-driven. However it is important to note that the pitch-angle scattering operator does not conserve momentum, which is most critical in the parallel direction. Therefore, the damping found in a cylinder could be the consequence of this missing feature and not a physical result [Dimits \& Cohen 1994]. Nonetheless, the results are useful to determine whether the instability is mainly being driven by a slab or toroidal ITG mode. EUTERPE also has the feature of including kinetic electrons, which made simulations of TEMs with collisions possible. The combination of collisions and kinetic electrons made the numerical calculations extremely time-consuming, since the time step had to be small enough to resolve the fast electron motion. In contrast to the ITG results, it was observed that collisions are extremely important for TEMs in a tokamak, and in some special cases, depending on whether they were mainly driven by density or temperature gradients, collisions could even suppress the mode (in agreement with [Angioni et al. 2005, Connor et al. 2006]). In the case of stellarators it was found that ITGs are highly dependent on the device configuration. For LHD it was shown that collisions slightly reduce the growth rate of the instability, but for Wendelstein 7-X they do not affect it and the growth rate showed a similar trend with collisionality to that of the tokamak case. Collisions also tend to make the ballooning structure of the modes less pronounced.},
  language  = {en}
}