SOL RF physics modelling in Europe, in support of ICRF experiments

¹ CEA, IRFM, F-13108 Saint Paul Lez Durance, France.
² Institut Jean Lamour, UMR 7198, CNRS-University of Lorraine, F-54506 Vandoeuvre Cedex, France.
³ Department of Applied Physics, Ghent University, Belgium
⁴ Max-Planck-Institut fur Plasmaphysik, Garching, Germany.
⁵ LPP-ERM-KMS, TEC partner, Brussels, Belgium.
⁶ Laboratoire Jacques Louis Lions, UPMC-Paris VI, CNRS UMR 7598, Paris, France.

^* Corresponding author: laurent.colas@cea.fr

Published online: 23 October 2017

Abstract

A European project was undertaken to improve the available SOL ICRF physics simulation tools and confront them with measurements. This paper first reviews code upgrades within the project. Using the multi-physics finite element solver COMSOL, the SSWICH code couples RF full-wave propagation with DC plasma biasing over “antenna-scale” 2D (toroidal/radial) domains, via non-linear RF and DC sheath boundary conditions (SBCs) applied at shaped plasma-facing boundaries. For the different modules and associated SBCs, more elaborate basic research in RF-sheath physics, SOL turbulent transport and applied mathematics, generally over smaller spatial scales, guides code improvement. The available simulation tools were applied to interpret experimental observations on various tokamaks. We focus on robust qualitative results common to several devices: the spatial distribution of RF-induced DC bias; left-right asymmetries over strap power unbalance; parametric dependence and antenna electrical tuning; DC SOL biasing far from the antennas, and RF-induced density modifications. From these results we try to identify the relevant physical ingredients necessary to reproduce the measurements, e.g. accurate radiated field maps from 3D antenna codes, spatial proximity effects from wave evanescence in the near RF field, or DC current transport. Pending issues towards quantitative predictions are also outlined.