Optimization of the snowflake diverted equilibria in the TCV tokamak

In support of the TCV experimental campaign aiming at studying H-mode plasmas with snowﬂake (SF) divertor, free boundary equilibrium and stability studies were performed with the SPIDER and KINX codes. Due to the high ﬂexibility of plasma shaping capabilities of TCV, SF divertor conditions can be reached for various plasma geometries. However, at high plasma current some conﬁgurations require poloidal ﬁeld (PF) coil currents close to the machine limit. This is particularly important when the equilibrium sensitivity to the edge pedestal proﬁles, which is higher than for standard X-point conﬁgurations, is taken into account. That is why the conﬁguration optimization should also include the proﬁle sensitivity study when planning the shot scenario.

Discrepancy between modelled and measured radial electric fields in the scrape-off layer of divertor tokamaks: a challenge for 2D fluid codes?

Jan Horacek, Richard Pitts, Marco Wischmeier

Examination of radial electric field (E,.) profiles in the scrape-off layer (SOL) of ASDEX Upgrade (AUG) and JET revealed large discrepancies between 2D fluid edge modelling and experiment. Experimental profiles of plasma potential (V-p) in the outer (low field) side of the plasma, obtained with reciprocating Langmuir probes, decay radially with electron temperature, T-e, with the -eE(r)/del T-e ratio being > 1.5. In contrast, code simulated E-r are fairly low in most of the SOL (compared with -del T-e/e). Modelling with kinetic treatment of neutrals and drifts was performed using the SOLPS code for AUG cases and EDGE2D-Nimbus for JET cases. Mismatches between modelled and experimental E-r may be caused by the recently established tendency for the SOLPS code to underestimate T-e in the divertor of AUG. It was attributed to non-locality of parallel transport of supra-thermal, heat-carrying electrons originating upstream of the divertor, which are usually only weakly collisional and can penetrate, with few collisions, to the target. Ratios -eE(r)/del T-e obtained from the probe measurements in JET are of order 1.6, while in AUG, JT-60U and TCV they are of order 3. Such high values point to the possibility of fast electrons contributing, apart from target heat fluxes, also to the formation of the Debye sheath. The problem of the underestimation of E-r. in the codes must be closely related with the well-known problem of the underestimation of those parts of parallel ion flows in the SOL that are influenced by the toroidal field direction. It was demonstrated earlier that parallel ion flow at the outer midplane is dominated by the ion Pfirsch-Schluter flow, which in turn is partly driven by the radial electric field. The T-e and E-r discrepancies, as well as discrepancies between simulated and experimental parallel ion flows, raise a question of the validity of fluid codes for the plasma edge modelling and prompt the inclusion of kinetic effects into present-day 2D fluid codes which assume strong collisionality.

2007

Snowflake Divertor Plasmas on TCV

Attila Bencze, Stefano Coda, Basil Duval, Federico Alberto Alfredo Felici, Ivo Furno, Jean-Marc Moret, Francesco Piras, Richard Pitts, Antoine Pochelon, Olivier Sauter, Costanza Zucca

Starting from a standard single null X-point configuration, a second order null divertor (snowflake (SF)) has been successfully created on the Tokamak à Configuration Variable (TCV) tokamak. The magnetic properties of this innovative configuration have been analysed and compared with a standard X-point configuration. For the SF divertor, the connection length and the flux expansion close to the separatrix exceed those of the standard X-point by more than a factor of 2. The magnetic shear in the plasma edge is also larger for the SF configuration.

PPCF2008

Pedestal Characteristics and MHD Stability of H-Mode Plasmas in TCV

Andreas Pitzschke

The tokamak à configuration variable (TCV) is unique in its ability to create a variety of plasma shapes and to heat the electron population in high density regimes using microwave power at the third harmonic of the electron cyclotron frequency. In the frame of this thesis, the impact of plasma shaping and heating on the properties of the edge transport barrier (ETB) in the high confinement mode (H-mode) was studied. This mode of operation is foreseen as one of the reference scenarios for ITER, the International Tokamak Experimental Reactor, which is being built to demonstrate the feasibility of thermonuclear fusion using magnetic confinement. A feature of H-mode regime operation are edge localized modes (ELMs), instabilities driven by the steep pressure gradients that form in the plasma edge region due to a transport barrier. During an ELM event, energy and particles are expelled from the plasma in a short burst. This will cause serious problems with respect to the heat load on plasma facing components in a tokamak of the size of ITER. Understanding of the phenomena associated with ELMs is thus required and dedicated investigations of their theory and experimental observations are carried out in many laboratories worldwide. This thesis presents several experimental and numerical investigations of tokamak behavior for configurations where the plasma edge plays an important role. From the experimental viewpoint, studies of transport barriers are challenging, as plasma parameters change strongly within a narrow spatial region. As part of the work presented here, the TCV Thomson scattering system was upgraded to meet the requirements for diagnosing electron temperature and density with high spatial resolution in the region of internal and external transport barriers. Simultaneously, the data analysis was significantly improved to cope with statistical uncertainties and alleviate eventual systematic errors. For measurements of the time evolution of density and temperature profile during the ELM cycle, the low repetition rate of the lasers used for Thomson scattering is a limiting. Although the system on TCV comprises 3 laser units that may be triggered in sequence with time separations down to 1 ms, time evolution over longer periods can only be reconstructed from repetitive events. In this context, an adjustment of the laser trigger to improve the synchronization with the ELM event is an advantage. A method was developed and implemented to generate a synchronizing trigger sequence, by a real-time monitoring of the D-alpha emission, which provides a marker for the ELM event. Recently, a "snowflake" (SF) divertor configuration, proposed as a possible solution to reduce the plasma-wall interaction by changing the divertor's poloidal magnetic field topology, was generated, for the first time, in TCV. A numerical code (KINX), based on a magnetohydrodynamic model (ideal MHD), was used to investigate the stability limits of this configuration under H-mode conditions and compare them with a similar standard single-null equilibrium. In a series of experiments, improved energy confinement was found and explained by improved stability of the edge region in the SF configuration. The influence of the pedestal structure in ELMy H-mode plasmas on the energy confinement and on ELM energy losses was investigated. The different ELM regimes found in TCV were analyzed, in particular the transition between type-III to type-I ELMs. The operational boundary of each ELM regime was characterized and verified by ideal MHD stability simulations for the ETB region. Recent studies on the scaling of the pedestal width with normalized poloidal pressure were confirmed. Using the capabilities of TCV, the influence of plasma shaping on pedestal parameters and MHD stability limits was investigated. In the past, models were developed to describe the onset of type-I ELMs, which are associated with modes in the ETB region arising from a coupling of pressure- and current-driven instabilities (coupled kink-ballooning modes). Experimental studies were performed to trace the temporal evolution of pedestal parameters characterizing the ETB during an ELM cycle. The results of these experiments were analyzed using information from MHD stability calculations. It is concluded that these models are capable of predicting limits as necessary conditions for ELM activity, but are not sufficient to fully explain ELM triggering.