Since 18 of December 2019 uses Nucleus credentials. Visit our help pages for information on how to Register and Sign-in using Nucleus.
May 10 – 15, 2021
Virtual Event
Europe/Vienna timezone
The Conference will be held virtually from 10-15 May 2021

Experimental validation of neutral beam current drive simulations in TJ-II plasmas

May 13, 2021, 8:30 AM
Virtual Event

Virtual Event

Regular Poster Magnetic Fusion Theory and Modelling P5 Posters 5


sadig mulas


The confinement of fast particles in fusion devices is shown to be strongly affected by Alfvén activity 1 and in order to model the associated wave particle interaction, a proper identification of the excited Alfvén Eigenmodes (AEs) is needed in the first place. The shear Alfvén spectrum is highly dependent on the iota profile of the magnetic configuration 2, which in turn depends on the net current present in the plasma. It has been proven that, in the TJ-II stellarator (R = 1.5 m, a = 0.2 m) 3 small changes on the iota profile, caused by the plasma current, can have a strong impact on the observed spectrum of shear Alfven waves [4, 5]. Ideally, a precise measurement of the iota profile, or failing this, an approximate knowledge of the different contributions to the plasma current is mandatory when it comes to obtaining reliable estimations of the shear Alfvén spectrum. While current drive in tokamaks has been extensively studied both experimentally, theoretically and numerically including comparisons with experiments for NBCD model validation, there have been only very few studies such as [6] of current drive in 3D non-axisymmetric configurations including stellarators so far.
Fig 1
The goal of the present work is the calculation of the neutral beam current drive (NBCD) for non-axisymmetric configurations and its validation in TJ-II NBI heated scenarios. Figure 1 shows the shear Alfven spectrum calculated with the STELLGAP code [7] for different plasma current profiles corresponding to different heating scenarios to illustrate the wide variety of conditions available in TJ-II for these studies.
In co-NBI scenarios (no current compensation), the dominant source of current is NBCD. This current is carried by the fast particles created during the NBI injection and its subsequent slowing down process. Theoretical models [8-10] that calculate the NBCD are available in the literature for different regimes and approximations and all of them rely on the assumption that the steady state fast ion distribution function is known. In order to obtain the fast ion distribution, numerical simulations must be performed. NBCD analysis in W7-AS [6] was carried out by means of a Fokker-Planck simulation code using first a simplified model to obtain the slowing down distribution function and then taking this distribution as input for theoretical NBCD models in the collisional [11] and collisionless [12] regimes. Preliminary estimations of NBCD (Fig. 2) were carried out for TJ-II plasmas [13] considering a more accurate calculation using Monte Carlo simulations. There, beam current drive was calculated using first the Monte Carlo code ASCOT [14, 15] to obtain the fast ion distribution function and NBCD was obtained using the model valid for the long mean free path regime [12].
Fig 2
The simulations performed with ASCOT did not include the charge exchange processes (CX), which have been shown to be very important in the range of energy of TJ-II plasmas; therefore, the NBCD results were corrected with simulations performed with the Monte Carlo code FAFNER [16], which is able to calculate the percentage of lost power due to CX. The calculations showed a reasonable agreement with the experimental data (Table 1).
In this work, neutral beam current drive is thoroughly investigated from numerical and experimental point of view, by carrying out simulations for different plasma parameters obtained in different heating scenarios (NBI with and without ECRH) and comparing with the large available data set. The comparison of the calculated NBCD with the one induced in the plasma when NBI is injected needs an estimation of the bootstrap current [17] contribution to the total measured plasma current (I_NBCD = I_P – I_BOOT). Where possible, a benchmark between FAFNER and ASCOT, both codes using now CX processes, will be presented in order to address the consistency of the theoretical-computational framework.
Table 1


1 W.W Heidbrink, Plasma Phys. Control Fus., 56 ,095030 (2014)
2 Ya. I. Kolesnichenko, Physics of Plasmas, 8, 491 (2001)
3 E. Ascasibar et al., Nuclear Fusion, 59, 112019 (2019)
[4] A. Cappa et al., P4.1040, 45th EPS, Prague, July 2018
[5] A. Rakha et al., Nuclear Fusion, 59, 056002 (2019)
[6] N. Marushchenko et al., D-5.003, 29th EPS, Montreux, June 2002
[7] D. Spong et al. Phys. of Plasmas 10, 3217 (2003)
[8] T. Ohkawa, Nuclear Fusion, 10, 185 (1970)
[9] D. F. H. Start et al., Plasma Physics, 22, 303-316 (1980)
[10] Y. J. Hu et al., Physics of Plasmas, 19, 034505 (2012)
[11] S.P. Hirshman, Physics of Fluids 21, p1295 (1978)
[12] N. Nakajima and M. Okamoto, NIFS-27, Japan, (1990)
[13] S. Mulas et al., P90, 22nd ISHW, Madison, September 2019
[14] J. Varje et al., Submitted to Comput. Phys. Commun. (2019)
[15] O. Asunta et al., Comput.Phys. Commun. 188, 33 (2015)
[16] G. Lister, Max-Planck-Institut für Plasmaphysik Technical Report IPP 4/222, 1985
[17] J.L. Velasco et al., Plasma Phys. Control. Fusion 53, 115014 (2011)

Affiliation CIEMAT
Country or International Organization Spain

Primary authors

Daniel López-Bruna (EsLNF) Ivan Calvo (Laboratorio Nacional de Fusion, CIEMAT) Dr Jose Luis Velasco Garasa (Laboratorio Nacional de Fusión, CIEMAT) Mervi Mantsinen (ICREA and Barcelona Supercomputing Center) Taina Kurki-Suonio (Aalto University) sadig mulas Álvaro Cappa (Laboratorio Nacional de Fusión CIEMAT) Joona Kontula

Presentation materials