Next Article in Journal
Symmetry and Quantum Features in Optical Vortices
Next Article in Special Issue
Thermonuclear Fusion Reactor Plasma-Facing Materials under Conditions of Ion Irradiation and Plasma Flux
Previous Article in Journal
Racemic and Meso Crystal Structures of an Axial-Chiral Spirobi-(dinaphthoazepin)ium Salt: Emergence of an S4-Symmetric Molecule
Previous Article in Special Issue
New Method for Calculation of Radiation Defect Dipole Tensor and Its Application to Di-Interstitials in Copper
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Evolution of Heavy Ion Beam Probing from the Origins to Study of Symmetric Structures in Fusion Plasmas

by
Alexander Melnikov
1,2
1
Kurchatov Center for Thermonuclear Energy and Plasma Technologies, National Research Center Kurchatov Institute, 123182 Moscow, Russia
2
Plasma Physics Department, National Research Nuclear University MEPhI (Moscow Engineering Physics Institute), 115409 Moscow, Russia
Symmetry 2021, 13(8), 1367; https://doi.org/10.3390/sym13081367
Submission received: 27 May 2021 / Revised: 24 June 2021 / Accepted: 9 July 2021 / Published: 27 July 2021
(This article belongs to the Special Issue Symmetry in Physics of Plasma Technologies)

Abstract

:
The overview discusses development of the unique fusion plasma diagnostics—Heavy Ion Beam Probing (HIBP) in application to toroidal magnetic plasma devices. The basis of the HIBP measurements of the plasma electric potential and processing of experimental data are considered. Diagnostic systems for probing plasma in tokamaks TM-4, TJ-1, TUMAN-3M and T-10, stellarators WEGA, TJ-II and Uragan-2M are presented. Promising results of the HIBP projects for various existing modern machines, such as TCV, TCABR, MAST, COMPASS, GLOBUS-M2, T-15 MD and W7-X and the international fusion tokamak reactor ITER are given. Results from two machines with similar size and plasma parameters, but with different types of the magnetic con-figuration: axisymmetric tokamak T-10 and helically symmetric stellarator TJ-II are compared. The results of studies of stationary potential profiles and oscillations in the form of quasimonochromatic and broadband fluctuations, turbulent particle flux, fluctuations of density and poloidal magnetic field are presented. The properties of symmetric structures—zonal flows and geodesic acoustic modes of plasma oscillations as well as Alfvén Eigenmodes excited by fast particles from neutral beam injection heating are described. General trends in the behavior of electric potential and turbulence in magnetized fusion plasmas are revealed.

1. Preface

In recent decades, diagnostics of high-temperature plasma has become an independent scientific discipline that has grown out of the physics of fusion plasma. Fusion plasma is a very unusual area of science. Unlike other areas of physics that study the properties of objects existing on Earth, plasma physics first creates an object of its research in terrestrial conditions—plasma with thermonuclear fusion temperatures higher than the temperature in the Sun core. Superconducting magnets, creating magnetic fields for plasma confinement, operate near the absolute minimum temperature, in conditions of space cold. The vacuum conditions in fusion machines for magnetic confinement of plasma are approaching to the conditions of cosmic space. The creation of such a complex and un-stable object as a fusion plasma in a terrestrial laboratory is in itself a formidable physical and technical issue. An equally difficult issue is plasma diagnostics—the development of tools and methods for studying plasma, reliable techniques and tools for measuring its properties and parameters, such as the temperature of its electrons and ions, density, stored energy, effective charge, and many others. The listed plasma parameters are not measured directly, which is obviously impossible in such a hot medium, but by their indirect signs, such as, for example, plasma radiation in different energy ranges. From the point of view of mathematics, the determination of the required plasma parameters from the quantities available for measurement is a classic example of inverse problems that require the development of special regularization methods for their solution [1,2]. In the last decades of the 20th century, one of the hottest areas of research has become the study of the electric field in the plasma of fusion machines with magnetic confinement. The reasons for this were the anomalous transport of particles and energy across the confining magnetic field, inherent in hot turbulent plasma; the identification of complex relationships between the electric field, plasma rotation, properties of turbulence and confinement. Over the past three decades, methods for diagnosing the potential of the electric fields in plasmas, as well as theoretical and numerical models of the transport processes associated with them, have been developed in the USSR and Russia. The present work is a brief overview of results of many years of research carried out by large teams under the leadership and with the active participation of the author at several fusion machines—tokamaks and stellarators.

2. Introduction

The study of electric fields is one of the urgent problems of modern physics of fusion plasma. These studies are related to both the fundamental problem—the emergence of an electric field in a quasineutral plasma [3], and an applied problem—the understanding of the electric field effects on the plasma transport processes in closed magnetic configurations [4]. It is known that the anomalous mechanisms linked with plasma turbulence are dominated in the particle and energy transport across confining electric field [5]. Currently, the hypothesis that the plasma turbulence and anomalous transport are stabilized by shear (radial inhomogeneity) of poloidal plasma rotation in crossed radial electric Er and toroidal magnetic Bt fields is generally accepted in magnetic confinement studies [6]. However, this hypothesis has both agreements and disagreements with experiment, so, the effect of Er to plasma confinement still remains unclear.
At present, it is believed that the problems associated with the electric field in tokamaks are becoming the principal ones in fusion research, and that for the success of the ITER project, it is necessary to change from studying the role of the electric field to controlling the anomalous transport with electric fields [7]. Measurement of the electric field in the plasma of modern fusion machines is a complex experimental issue [8]. At the plasma edge, the electric potential is measured with Langmuir probes. In the hot plasma core, the electric field, as a rule, is not measured directly, but it is determined from the measured plasma rotation velocity using spectroscopy or correlation reflectometry. The only direct method to find the electric potential in the plasma interior is the Heavy Ion Beam Probe (HIBP) diagnostics [9].
Diagnostic techniques for toroidal plasmas always have limitations concerning the access to the plasma. These constrains are dictated by the complexity of the toroidal devices and by the trend to reduce the ripple of the confining magnetic field by an increase in the number of toroidal field coils. This makes the diagnostic ports narrow in toroidal direction. Finally, most of the diagnostics observe just a part of the plasma vertical cross-section. They have to assume the symmetry of the measured parameter to complete the whole picture of its spatial distribution. The typical example of such an approach is the multi-chord plasma density measurement with a microwave interferometer as a function of the plasma radius that then assumes a poloidal symmetry of the density distribution.
On top of that, there are theoretical predictions of various asymmetries for plasma parameters. For example, there could be the so-called ballooning effect, which deforms the contour lines of plasma parameters at the Low-Field Side (LFS) of the torus vertical cross-section, breaking the in-out symmetry. Gyrokinetic simulations predict a strong increase of the broadband electrostatic turbulence in the LFS in respect to the High Field Side (HFS). Therefore, the issue of symmetry is an open one in many aspects of plasma studies. Resolving this issue is highly demanding for the diagnostic techniques.
This overview describes the HIBP diagnostics of toroidal plasma and its development over the past four decades for studying various phenomena in plasma including issues of their symmetry.

2.1. Heavy Ion Beam Probe—A Tool for Measuring Electrical Potential and Plasma Turbulence

The Heavy Ion Beam Probe (HIBP) is a unique method for plasma studies in fusion machines [10]. HIBP was first implemented in the late 1960s by Hickok and Jobs on the arc discharge and later on the ST tokamak [11]. For the first time in Russia (USSR), HIBP was implemented in the early 1980s by Krupnik and Nedzelskiy at the TM-4 tokamak [12,13]. Currently, HIBP is the only non-disturbing method for direct measurements of the electric potential in a hot plasma core [14]. In this paper, we consider the mathematical aspects of diagnostics, equipment and methods for measuring plasma parameters available for the modern HIBP [15].
To date, there are only two HIBP diagnostics operating in tokamaks in the world: Tuman-3M [16] (Ioffe Institute of Physics and Technology, St. Petersburg, Russia) and ISTTOK (Lisbon, Portugal) [17]. The T-10 tokamak [18] (Kurchatov Institute, Moscow, Russia) was terminated in 2018 due to creation of the next generation tokamak T-15MD. In addition, HIBP operates on the world’s largest stellarator LHD (National Institute for Fusion Science, Toki, Japan) [19,20,21], on the TJ-II stellarator (CIEMAT, Madrid, Spain) [22,23], as well as on the reversed field pinch MST [24] (Madison University, Madison, WI, USA) and on the tandem mirror GAMMA-10 [25] (Tsukuba University, Tsukuba, Japan).

2.2. Physical Principles of HIBP Measurements

To probe the plasma with the beam of heavy ions, the beam is injected across the confining magnetic field. When the beam particles fly through the plasma, some of them collide with plasma particles (mainly electrons) and lose one or more electrons. As a result, a fan of secondary ionized particles with higher charges is formed. The HIBP diagnostic scheme and an example of its realization in the T-10 tokamak are shown in Figure 1.
A small-aperture detector is usually placed outside the magnetic field that allows the secondary ion beam (part of the fan), formed at the secondary ionization point (sample volume or SV) on the probe beam trajectory, to enter the detector. The detected secondary ions carry information about the plasma parameters in the SV. The spatial resolution of the method λ is determined by the dimensions of the sample volume and mainly depends on the size of the aperture and the position of the detector. In real experiments on different machines λ = 0.2–2 cm.
HIBP allows us to simultaneously measure several plasma parameters: electric potential φ, electron density ne, and magnetic potential A (or poloidal magnetic field Bp or current density j). The position of the sample volume can be moved along the plasma cross section by varying the energy of the probing beam Eb or the injection angle of the beam into the plasma α. Thus, the set of points observed by the detector forms a two-coordinate (Eb, α) grid called the detector grid.
The trajectories of primary and secondary ions in a low-current tokamak lie in the meridional plane between the confining field coils or near it. For tokamaks with a high current density and for stellarators, the trajectories of the probing particles are shifted in the toroidal direction, and the detector grid becomes three-dimensional. This essential feature of the motion of probing particles should be taken into account during development projects of placing diagnostics on new machines installations, and in the analysis of measurement results.
Thus, a specific feature of HIBP diagnostics is that the values of plasma parameters are determined from measuring the characteristics of the secondary beam: intensity, energy, etc., while the position of the measurement point can only be found from calculations of the trajectories of probing particles. Therefore, before placing HIBP diagnostics on a new machine installation, not only the development of diagnostic equipment is required, but also preliminary calculations.
The position of the sample volume in the meridional plane depends on a large number of probing parameters, including: (i) geometric parameters: the coordinates of the injection xi, yi and detection points xD, yD, and the injection angle α (ii) the physical parameters: the toroidal magnetic field Bt, the energy Eb, the mass m, and the electric charge q of the probing particles.
Probing particles move in the magnetic field of the machine along a Larmor circle with a radius:
R L = c q e B t 2 m E b .
It is necessary to determine not the value of RL itself, but the values of the parameters on the right hand side of Equation (1), which are subject to independent restrictions. The equality RL(q, m, Eb, Bt) = const reduces the number of physical parameters from 4 to 3. So, for 2D optimization in the meridional plane, there are eight parameters: xI, yI, xD, yD, α, Bt, Eb, q. The need to pass the trajectories of particles through the ports of the vacuum chamber and the structural elements of the machine impose serious constraints on these parameters. These constraints are related to each other because the allowable tolerances for each parameter depend on the others. Such relations defy analytical description due to the complexity of the trajectory behavior. The configuration of the vacuum chamber in the meridional plane approximately determines their tolerances. The problem of optimizing experimental conditions includes the choice of optimal probing parameter values in the sense of the following goal functions:
  • Pass the beam through the existing vacuum vessel ports;
  • Find the detector line from the core to the edge of the plasma;
  • Find a detector grid covering the maximum part of the plasma cross section;
  • Optimize the range of beam energies.
In high-current tokamak discharges, as well as in stellarators, the magnetic field has a noticeable poloidal component Bp, which shifts the orbits of the probing particles out of the meridional plane and transforms them into spatial curves [26]. In this case, additional parameters appear that affect the sample volume. These are the toroidal coordinates of the injector and detector, ZI and ZD, and the injection angle between the initial velocity vector and the meridional plane β. This gives us a 12-dimensional optimization problem with implicit connections between the limits of parameters and with non-formalized goal functions. Figure 2 shows an example of solving this problem by ‘the shooting’ method in T-10. Examples of solving the optimization problem in other machines are given in Section 7 and Section 8. The first goal function was achieved in all considered machines, as the chosen probing schemes ensured the passage of trajectories through the vacuum vessel ports and allowed measurements. Maximization of the second and third goal functions provides the maximum possible radial and angular sizes of the observation area in plasma. The fourth function is not always so essential, because the energy range of modern accelerators is quite large. However, for large machines, such as the W-7X stellarator and the constructed ITER tokamak reactor, this task is of great importance. As will be shown in Section 7, the chosen injection schemes allow one to connect the core and the edge of the plasma with a detector line using the acceptable beam energy.

3. Mathematical Problems of Determining Plasma Parameters Using HIBP

3.1. Determination the Spatial Distribution of Electric Potential

The principle of HIBP measurements is based on the potentiality of the electric field, the conservation of the total energy of the charged particles in the plasma electric field (Figure 3). The probing beam, hereinafter referred to as the primary beam, enters the plasma with an initial energy Eb. At the ionization point, which is the plasma region of interest (SV), a particle of the primary beam loses an electron with potential energy e φ p l S V which is acquired by the secondary ionized probe ion. The total energy of secondary ions leaving the plasma is E d = E b + e φ p l S V . Therefore, the local potential in SV is equal to the energy difference:
φ p l S V = ( E d E b ) / e .
Thus, the problem of measuring the plasma electric potential does not pose serious mathematical difficulties. However, from the point of view of conducting physical measurements it is a very difficult task, as the dimensions of modern plasma machines (scale of 1–5 m) and the magnitude of the magnetic field (scale of 1–3 T) impose high requirements on the Larmor radius of probing particles, which must exceed the plasma size. According to (1), this requires the use of probing ions with the minimum charge number (q = 1), maximum possible mass (Cs+, Tl+, Au+) and energies (on a scale of several hundred keV—several MeV). On the other hand, the measured plasma potential has the scale of the electron temperature. In the plasma core, it ranges from several hundred eV to keV. At the plasma edge, it has a scale from several eVs to hundreds of eV. Thus, according to (2), the measurement of the potential claims the computation of a small difference of two large quantities. Such measurements require high accuracy of each quantity or highly stabilized high-voltage sources and the use of precision instruments to measure the energy of the secondary beam.

3.2. Determination of the Spatial Distribution of Plasma Density

The local values of the plasma density n(l) can be found from the relation:
i 2 = i 1 σ n ( l ) λ ( l ) ,
where i1 and i2 are the intensities of the primary and secondary beams in the SV; σ = σ υ e / υ b is the effective ionization cross section by electron impact, <συe> is the ionization rate averaged over the Maxwellian distribution of the electron velocity υe, itself a function of the electron temperature Te; υb is the initial velocity of the probing beam, υe > υb; l = l(x, y) is the coordinate on the detector line; and λ(l(E)) is the longitudinal dimension of SV, the length of the arc along the primary trajectory, from which the secondary ions reach the detector. The electron density n(l) is averaged over the length λ.
If the density is low and the dimensions of the machine are small, then we can assume that the intensities i1 and i2 are the intensities of the primary beam at the exit from the ion injector and the secondary beam at the detector. In this case, Relation (3) can be directly used to determine the local values of the plasma density. However, in modern fusion machines installations, both the primary and secondary beams are significantly attenuated, when passing through the plasma. Methods for determining plasma density using HIBP in the case of strong attenuation are described in [27,28], and they were used for interpretation of experimental data from the TM-4 tokamak.
The expression for the current of secondary ions to the detector It, taking into account the attenuation of along the beam trajectory, caused by elementary processes, including collisions of beam particles with plasma particles, can be written as follows:
I t ( l ) = 2 I b σ 12 ( l ) n ( l ) λ ( l ) exp ( I l n σ 12 ( s ) d s l D n σ 23 ( s ) d s ) ,
where Ib is the ion beam current leaving the injector, n is the electron density, σ12 and σ23 are the known effective cross sections for collisional ionization of the probing beam (e.g., Cs in the cases of TM-4 and TJ-II) by plasma electrons.
e + Cs + Cs 2 + + 2 e ,   e + Cs 2 + Cs 3 + + 2 e .
The unknown function in this nonlinear integro-functional equation is the plasma electron density n(l); functional dependence It(l) is determined experimentally; Ib is a known constant. The integration lines Il and lD, and the quantities λ(l) are found by calculating the trajectory and the detector line.
Equation (4) is an example of a nonlinear inverse problem with a nonclassical operator structure: the integral of the desired function is in the exponent. The limits of integration are variable; they depend on the coordinates of a point on the detector line l = l (x, y). This equation may be solved by a quasi-solution method: parametrization the desired function and finding the optimal values of the parameters by minimizing the discrepancy functional. The application of this method to the analysis of experimental data from TM-4 is shown in Figure 4. We see that the reconstructed density profile is consistent with the data of other diagnostics, so the method works well for TM-4 conditions. More recently, the same technique was applied for the TJ-II data. The results show that the reconstructed HIBP density profile lies within the error bars of the Thomson scattering data [29].

3.3. Determination of the Plasma Magnetic Potential (Field of Plasma Current)

The plasma current profile mainly determines the equilibrium of the plasma column, its stability, energy balance, and the characteristics of fusion plasma. At present, for measurements of the poloidal magnetic field, there is a very limited set of diagnostics [30], which do not have a high temporal and spatial resolution. Therefore, the problem of creating diagnostics for poloidal fields in toroidal plasma is still urgent. In this Subsection, we consider a method for determining the plasma current field from measurements of the shift of the probing beam caused by this field [31,32].
Let us consider the simplest model of a symmetric torus (without ripple), in which the toroidal magnetic moment is conserved for a probing charged particle. After integration the Lagrangian for the probing particle along the trajectories, we obtain the following expression for the desired function Aφ, the toroidal component of the magnetic potential:
( φ d φ 0 ) + e m c v [ q s 0 s i A φ r d s + ( q + k ) s i s d A φ r d s ] = = ( r 0 2 φ ˙ 0 + q e m c A φ 0 r 0 ) 1 v s 0 s d d s r 2 + k e m c v A φ ( s i ) s i s d d s r 2 r ( s i ) . ,
Here s0 is the injection point, si is the ionization point, and sd is the detection point, q and q + k are dimensionless charge numbers of primary and secondary beam particles. Equation (6) connects the unknown function Aφ with the angular displacement of the beam in the detector respect to the one of injector φdφ0, that is, with the measured quantity. Equation (6) is a nonlinear Volterra-type integral equation of the second kind. As the integration limits are variable, they depend on the coordinates of the ionization point in the plasma (points on the detector line). In contrast with the equation of the first kind, solving such equations is a well-posed problem. The function Aφ is unknown only inside the plasma column, for r < a. Outside the plasma, where r > a, Aφext is determined by its values inside the plasma as follows:
A φ e x t = 1 c j φ R d V = 1 4 π ( Δ A ) φ R d V = A φ e x t ( A φ ) ,
where the integration is carried out over the entire region of current flow, and R is the distance from the volume element dV to the position, where Aφext is computed. Note that the integrals in (6–7) are curvilinear integrals over 3D trajectories. The arc element ds along the particle trajectory depends on the unknown function Aφ: ds = ds(Aφ), therefore, Equation (7) is a nonlinear integral equation.
The solving of this equation is based on the iteration method. Figure 5 shows the results of the numerical solution of Equation (7) for the conditions of the TUMAN-3M tokamak.

4. Application of the HIBP to Fluctuation Measurements

The HIBP is a multi-purpose diagnostic tool capable to measure the local potential φ using the secondary beam energy, plasma density ne from the beam current It, and the poloidal magnetic field Bp from the toroidal angular displacement of beam φd or the linear beam shift ζ = φd R, where R is the distance from the detection point to the origin of the cylindrical coordinate system [33].
High temporal resolution (up to 0.5 μs) allows measurements of fluctuations of the parameters under study. Expressions for the links between the relative fluctuations of the measured quantities and the parameters under study are following:
For the variable component of the potential:
φ ˜ p l S V ( t ) = 2 U a n F δ i ( t ) ,
where δi is the normalized difference between the beam currents across the separated detector plates, Uan is the voltage applied to the analyzer, and F is its dynamic coefficient. Thus, there is a simple linear relationship between the desired value δ φ p l S V and the measured value δ i ( t ) .
If the spatial correlation length of density fluctuations is much less than the length of the beam trajectory, then the nonlocal terms in Equation (4) may be omitted due to integration of the oscillatory component δne over the beam trajectory. In this case:
δ n e ( ρ , t ) = I ˜ t ( ρ , t ) / I ¯ t ( ρ ) .
For large-scale density fluctuations δne, for which the radial correlation length is large and comparable to the plasma radius, the local data on oscillations in the sample volume can be contaminated by attenuation effects accumulating along the entire trajectory. For the experimental data, considered in this paper the large-scale modes are not typical, so the density oscillations are measured locally [34].
The oscillating component of the toroidal shift of the beam in the detector ζd, is described by the following expression [35]:
δ ζ d = δ ψ SV m L 2 d l R 2 1 m L 1 δ ψ d l R 2 2 m L 2 δ ψ d l R 2 + P ζ i n j m L 1 δ ψ d l R 2 + ζ i n j .
Here, ψ = Aφ/R, Pζ inj is the toroidal magnetic moment at the injection point, known constant. Note that the integrals in (10) are taken along the total paths L1 and L2 from the injector to the position of the sample volume (SV), and from SV to the detector. Equation (10) shows that the local term δψSV input the main contribution to the oscillations of beam shift. The contribution of integral terms is much smaller for oscillations with shorter correlation length (or radial wavelength) due to the path integration. It was shown in Reference [36] that the integral terms do not strongly affect the local terms from MHD oscillations, therefore HIBP provides practically local measurements of magnetic fluctuations. However, a nonlocal contribution to ζd from oscillations with a large correlation length, comparable to the length of the beam trajectory in plasma, is also possible.
In conclusion, we note that, despite the path integral effect, HIBP provides fairly local measurements of density and magnetic field fluctuations, taking into account the conditions specified above, as well as purely local measurements of electric potential fluctuations. All local measurements are averaged over the sample volume (SV), which size determines the spatial resolution of the HIBP diagnostic. The sample volume typically looks like the inclined elliptical disk with diameter equal to the diameter of the beam, and its thickness corresponds to the entrance slit of the analyzer. The typical radial dimension of SV ranges from 0.1 to 1 cm.

5. Mathematical Problems of Experimental Data Processing for Fluctuations of Plasma Parameters

5.1. Spectral Fourier Analysis of Oscillations

To analyze the oscillating component of an arbitrary signal x(t), we introduce the Fourier transform as Sx(f) and Sx*(f) as its complex conjugate. Thus the power spectral density (PSD) of oscillations is:
S x x ( f ) = S x ( f ) S x * ( f ) .
The cross-spectral oscillation power density (CSPD) for signals x(t) and y(t) is the following product of their Fourier transforms:
S x y ( f ) = S x ( f ) S y * ( f ) .
An example of calculating the PSD of oscillations of plasma parameters determined in a certain time domain is shown in Figure 6a. As a rule, when digitizing an experimental signal with a frequency of 1 μs, the PSD is calculated over a period of 0.25–1 ms, depending on the characteristic time of the change in the frequency of the processes under study.
The power spectra of characteristic peaks of monochromatic fluctuations corresponding to the Geodesic Acoustic Mode (GAM) [37,38] and magnetohydrodynamic (MHD) tearing mode with a poloidal mode number m = 2 are seen. In contrast to the solitary peak of the m = 2 mode, the frequency GAM peak is accompanied by a high-frequency satellite peak.
Figure 7 presents the physical picture of the zonal flows as symmetric torsional plasma oscillations [37]. While zonal flow by itself, also named residual zonal flow, typically has low-to-near zero frequency, GAM as a higher frequency branch of zonal flows has a typical frequency of around 20 kHz. Both GAM and tearing modes present examples of symmetric structure in the plasma toroidal magnetic configuration. GAM is poloidally and toroidally symmetric torsional plasma oscillation with n = m = 0, where n is the toroidal mode number [39]. MHD tearing mode is a chain of magnetic islands, that have helically symmetric structure, located on magnetic surface with rational value of safety factor q = m/n in the toroidal magnetic configuration [40].
This spectrogram, describing a process lasting 500 ms, consists of a sequence of elementary spectra, each of which is calculated in 1 ms. The potential spectrogram is dominated by the GAM peak, accompanied by a satellite peak [41]. Application of additional ECR heating in the time interval from 600 to 800 ms leads to an increase in the frequency of both the GAM and the satellite due to growth of the electron temperature [42].

5.2. Fourier Analysis of Coherency and Cross-Phase of Oscillations

Time evolution of the correlations between signals x(t) and y(t) are determined by spectrograms of their coherency Cohxy and cross-phase θxy, as following:
C o h x y ( f , t ) = | S x y | | S x x S y y | 1 / 2 ,   0   <   C o h x y   <   1 ,
θ x y ( f , t ) = tan 1 { Im ( S x y ) Re ( S x y ) } , π   <   θ x y   <   π
Coherence and cross-phase are modulus and argument of the complex function CPSD.
Figure 8 presents an example of calculating the coherence spectrogram. Values of statistically significant coherence indicate the Alfvén eigenmodes (AEs), which is another type of MHD mode. AEs are electromagnetic oscillations, propagating along the magnetic field line with Alfvén speed and possessing helically symmetric structure with finite m and n. In the vicinity of the magnetic flux surface with rational q value, equals to m/n for a specific mode, the dispersion relation for AE contains also a GAM component, causing GAM frequency to be a lower limit for AE frequency [43,44]. It presents a beautiful link between helically symmetric structure AE and toroidally symmetric structure GAM. Figure 7 shows that the frequency of AE decreases in time, caused by the density increase in accordance with Alfvén law f A E ~ B t n e 1 / 2 . The circles indicate different branches of the Alfvén eigenmodes with different poloidal mode numbers m = 4 and m = 3, calculated from the data of a set of magnetic probes.
Figure 9 shows an example of calculating the cross-phase spectrogram. We see that the magnitude of the cross-phase does not change, remaining equal to −π/2, although the mode frequency decreases twice with time, due to an increase in the density in this discharge. This constancy of the cross-phase values between the signals of the plasma density and the poloidal magnetic field indicates that the cross-phase is a stable individual characteristic for each Alfvén eigenmode.

5.3. Bispectral Fourier Analysis of Plasma Oscillations

Bispectral analysis is widely used to study nonlinear three-wave interaction processes. The quadratic coefficient of cross-bicoherence for three independent quantities x(t), y(t) and z(t) is defined as follows:
b x , y , z 2 ( f 1 , f 2 ) = | S x ( f 1 ) S y ( f 2 ) S z ( f 1 + f 2 ) * | 2 | S x ( f 1 ) S y ( f 2 ) | 2 | S z ( f 1 + f 2 ) | 2 ,   0   <   b 2 <   1 ,
where brackets denote the time averaging. Similar to the cross-phase, biphase is defined as an argument of the complex function:
θ x , y , z ( f 1 , f 2 ) = tan 1 Im S x ( f 1 ) S y ( f 2 ) S z ( f 1 + f 2 ) * Re S x ( f 1 ) S y ( f 2 ) S z ( f 1 + f 2 ) * , π   <   θ <   π .
Auto-bicoherence for random variable x(t) is defined as:
b x 2 ( f 1 , f 2 ) = | S x ( f 1 ) S x ( f 2 ) S x * ( f 1 + f 2 ) | 2 | S x ( f 1 ) S x ( f 2 ) | 2 | S x ( f 1 + f 2 ) | 2 ,   0   <   b x 2   <   1 .
Statistically significant bicoherence b points to three-wave interaction between oscillations at frequencies f1, f2 and f3 = f1 + f2. An example of calculating bicoherence and biphase is shown in Figure 10. The excess of the bicoherence coefficient over the background at the GAM frequency indicates a significant three-wave interaction between the GAM and the broadband turbulence of potential [46]. Note that, at the GAM frequencies the biphase has finite values, while outside them it has a stochastic character. Thus, analysis of HIBP data shows that GAM in T-10 linked with broadband turbulence [47,48], and is not excited by other mechanisms, such as, e.g., the effect of energetic particles.

6. Measurements by Several Spatial Channels

6.1. Measurement of Turbulent Particle Flux on the TJ-II Stellarator

A fundamentally new element in the HIBP operation is the use of a multichannel energy analyzer, which allows one to measure both the potential and density in several sample volumes. Figure 11 shows a cross-sectional view of a multichannel analyzer installed in T-10. On the TJ-II stellarator, where two diagnostic systems HIBP1 and HIBP2 are currently operating (see Section 7), on the second HIBP2 system, the same analyzer is installed, somewhat reduced in size. Let us consider the operation of a multichannel analyzer using the simplest example of a two-channel analyzer with two slots and two detectors. Such an analyzer is installed on the HIBP1 complex of the TJ-II stellarator. The detector line for the TJ-II two-channel analyzer is shown in Figure 12.
The probing parameters are adjusted in such a way that both sample volumes lie on the same magnetic surface, and they are shifted poloidally by ∆x ~1 cm, as shown in Figure 13. This method allows us to directly measure the poloidal component of the electric field Ep by the potential difference, Ep = (φ1φ2)/∆x. The finite distance between the investigated volumes limits the poloidal wave vector of oscillations kp = kθ < 3 cm−1. Direct measurement of Ep allows us to define the radial velocity of E × B drift Vr = Ep/Bt.
As oscillations of the total particle current or plasma density are measured simultaneously with the potential oscillations, we may directly calculate the radial turbulent particle flux:
Γ r ( t ) = n ˜ e V ˜ r = 1 / B t n ˜ e ( t ) E ˜ p ( t ) = Γ E × B ,
Such use of HIBP on stellarators was first applied on the TJ-II. To find Γ E × B ( t ) , the density fluctuations n ˜ e must be measured simultaneously with E ˜ p in the same spatial region. To analyze the frequency spectrum of the flux and its temporal dynamics, it is sufficient to measure the relative density fluctuations δ n e ( t ) = I ˜ t ( t ) / I ¯ t . In the low-density case, to estimate the absolute values of the flux Γ E × B ( t ) we may use the measured value I ˜ t ( t ) as n ˜ e .
In the high-density case, the beam attenuation should be accounted. The density should be normalized as n ˜ e = I ˜ t o t / I ¯ t o t n ¯ e , where the oscillating density component I ˜ t / I ¯ t is measured with HIBP, and the normalized factor n ¯ e is measured by another diagnostic, e.g., by the interferometry or by Thomson scattering.
Fourier transform for the time-dependent flux Γ E × B ( t ) , or the spectral function of flux is calculated as follows [49]:
Γ E × B ( f ) = 2 / B t Re ( S n E )
Here SnE is the Fourier cross-spectrum between density and poloidal electric field Ep oscillations.
Figure 14 shows an example of the calculated turbulent particle flux in discharge with NBI heating. The flux, linked with broadband turbulence, demonstrates the intermittent nature. It consists of a stochastic sequence of controversially directed bursts, which are mostly directed outwards. The spectrograms show Alfvén eigenmodes in the form of a family of quasi-monochromatic oscillations with a frequency decreasing in time. We see that AEs may contribute both to outward and inward flux, or also could produce no flux, depending on the phase relations between Ep and the density oscillations [50]. Usually, AEs contribute to the turbulent particle flux Γ E × B A E , making up an appreciable fraction of the total flux Γ E × B . Figure 14 shows that Γ E × B A E can be much larger than the broadband turbulent flux Γ E × B B B in the same frequency range [51]. It was shown that the value of the turbulent flux is comparable to the value of the total particle flux; therefore, ΓE×B plays an essential role in the particle balance. It was shown that the transition to the improved confinement mode is characterized by suppression of broadband fluctuations of the plasma potential and density, as well as by a decrease ΓE×B both at the plasma edge and in the hot core [52].

6.2. Measurement of Turbulent Particle Flux in the T-10 Tokamak

In contrast to the TJ-II stellarator, the geometry of the T-10 machine gives more uncertainties in determining the flux. The scheme of the experiment on measuring the flux with the 5-slit analyzer is shown in Figure 11. An example of calculating trajectories for the standard discharge: Bt = 2.4 T, and energy Eb = 220 keV is shown in Figure 15. We can see that the potential difference φkφj corresponds to the electric field E, directed along the segment kj. To measure the poloidal field, it is necessary to align the segment in the poloidal direction. This is achieved in the point of deepest penetration for the detector line of equal energy. [53]. Thus, we can measure the flux for each energy in one spatial point.
Figure 16 presents an example of the turbulent flux measurements in the standard T-10 discharge with ohmic and ECR heating. Only two frequency ranges with non-zero flux are seen in the spectral flux function: the quasi-coherent (QC) mode and the geodesic acoustic mode (GAM). In the remaining frequency ranges, the flux is stochastic. The QC mode dominates during the entire discharge, allowing particles to flow outward. The flux in the GAM frequency range Γ E × B G A M is by the factor of 5–6 less than the flux of the QC mode Γ E × B Q C [54].
Another example of the flux calculation for GAM is shown in Figure 17. Here we consider the ohmic T-10 discharge with slow density ramp-up. The density reaches the critical value at t = 600 ms, then the MHD instability develops, which leads to the disruption at t = 750 ms. Figure 17a shows that GAM dominates on the spectrogram of potential fluctuations, but its amplitude decreases with density increase. The spectrogram of fluctuations in the potential difference (Figure 17b) completely contrasts with Figure 17a, as in Figure 17b GAM does not dominate, but it is completely invisible due to the background noise. This observation agrees with the well-known fact that GAM on the potential has poloidal symmetry. Additionally, GAM was not observed on the spectrogram of flux [55].
An important conclusion follows from the above experiments: the radial turbulent particle flux at GAM frequencies is on the noise level or only slightly exceeds it due to the poloidal symmetry of potential oscillations at GAM frequencies.

6.3. Measurement of the GAM Poloidal Symmetry in the T-10 Tokamak

In the previous Subsection, we considered correlation measurements of potential, which are physically corrected as the potential measurements are completely local. In the low-density case, when the attenuation factors in the analysis of the variable component of the density may be neglected, the correlations of the density fluctuations may be also considered. The issue of the locality of density measurements is discussed in detail in the Reference [34]. Finally, as shown in Section 4, the density oscillations measurements are localized in the sample volume, so the density coherence is free from the non-local effects.
An example of a two-point measurement technique [56,57] is shown in Figure 18. The figure shows that the sample volumes are separated in the poloidal direction by 1.67–1.74 cm, but also have some radial shift of 0.81 cm.
The measurement results are shown in Figure 19. Time traces of the main discharge parameters are shown in Figure 19a. Figure 19b,c shows the low-frequency part of the power spectrograms of the potential oscillations measured using the center φ3 (Figure 19b) and edge φ5 (Figure 19c) slits. Figure 19d shows that the coherence coefficient between φ3 and φ5 in the GAM frequency domain reaches unity. Figure 19e shows that the phase shift between these poloidally shifted signals is zero, indicating that GAM is poloidal symmetric, i.e., m = 0.
Figure 20 shows the results of similar measurements carried out in the same mode on the radial interval of 18–26 cm. Figure shows that for closely located points 2 and 4, the cross-phase at the GAM frequency is close to zero over the entire measurement interval. The averaged error of cross-phase measurement, determined by the signal-to-noise ratio, is about 0.1 radians, increasing towards the edge to 0.2 radians. This gives the poloidal mode number m = 0 c with error ±1, and ±2 at the edge. Figure 21 shows the results of the cross-phase measurements for the most distant points 1 and 5 in the same shot, but in a wider radial range. We see that the cross-phase between the signals from the 1st and 5th points at the GAM frequency is also close to zero over the entire measurement interval. The error of phase measurement is a constant of about 0.1 radians on average, and the error in determining the poloidal mode number m = 0 decreased due to an increase in the distance between points S1 and S5 on average to ±0.5, except for the edge, where it remained at a level of ±2 [58]. Thus, phase measurements of the potential confirmed the theoretically predicted poloidal symmetry of potential oscillations for GAM.

7. HIBP Diagnostics and Their Outcome for Plasma Symmetry Study

7.1. TM-4 HIBP

The first HIBP complex in Russia (USSR) was on the TM-4 tokamak, created early in the 1980s by Krupnik, Bondarenko and Nedzelskiy [59]. The scheme of the complex is shown in Figure 22.
The energy of secondary particles was determined by a multigrid ion analyzer. On TM-4, the first measurements of the potential profiles were carried out at a fixed position of the sample volume (SV), depending on the discharge parameters, such as the average chord density, plasma current and toroidal magnetic field. These measurements were carried out in a series of reproduced shots, in which the SV scans plasma from the center to the edge. It was shown that the plasma potential at the core has a negative sign and a scale about several hundred volts. For the first time, the link between the potential and the average plasma density was established: a higher density corresponds to a larger absolute value of negative potential. The measurements of the plasma density profile with a HIBP are described in Section 3.2. The TM-4 tokamak was terminated in 1982 due to the construction of the T-15 tokamak. Creation and scientific exploitation of TM-4 HIBP complex made the physical and technological basis for further development of the HIBP diagnostics towards the tool to study symmetric structures in toroidal plasmas.

7.2. T-10 HIBP

The HIBP complex on the T-10 tokamak was created in the mid-1980s by Krupnik, Bondarenko, Khrebtov and Nedzelskiy [18,60]. In contrast with the TM-4 complex, it used a Cs ion beam accelerated up to an energy of 160 keV, and a ‘flat mirror’ analyzer [61] with a higher energy resolution than the grid analyzer of TM-4. As the complex was implemented into the existing machine, there were no convenient ports to conduct the beam through the plasma. For injection we use a standard wide vertical port, and to pull out the secondary beam, we have to cut the vacuum chamber and weld an additional port to it [62]. Such a significant reconstruction of the machine became possible due to the successful HIBP operation on the TM-4 and the obtaining of important physical results—first measurements of plasma potential profile. This important decision was adopted with support by Dnestrovskij, Razumova and Sokolov. The initial version of the HIBP complex is schematically shown in Figure 23.
Setting up the T-10 complex and conducting the ion beam into the analyzer encountered with an unexpected issue—a strong toroidal displacement of the beam, which hampered conducting the ion beam through the long narrow ports of the tokamak. To solve it, for the first time in the HIBP operation, the toroidal correction of secondary particles was applied and a secondary ion beamline was installed outside the vacuum chamber. The passage of the beam into the analyzer through narrow diagnostic ports of the tokamak became possible thanks to the assistance of the T-10 leading experimentalist Dr. Bobrovskiy. At the initial stage of the HIBP operation, the measurements were carried out in a series of repeated tokamak shots, in which the SV position moved from the plasma edge towards the center. The limitations on the beam energy and ion mass allowed us to carry out measurements at a reduced magnetic field Bt = 1.5 T in the range 0.6 < ρSV < 1 [63,64,65].
Later on, the T-10 HIBP complex has permanently evolved in the following directions [66]:
  • The Cs ions were replaced by heavier Tl ions. The energy range was consequently increased to 220, then to 280, 300, and finally to 330 kV. These changes allow to study discharges with microwave heating (ECRH) at 2.08 < Bt < 2.5 T. In discharges with Bt < 2.17 T, the radial region 0.2 < ρSV < 1 was investigated;
  • The ion beamlines were modernized in the following way:
    • Two serial pairs of steering plates were installed into the primary beamline. The first plates deflect the beam “counter-current” to the extreme toroidal position near the second plates. Thus, the beam approaches the second plates not axially, but maximally shifted. The second plates deflect the beam “co-current”. As a result, a zigzag trajectory is formed, and it became possible to use the full toroidal width of the port, not half of it, as in the axial injection of the beam, and maximize the toroidal correction angle β.
    • The secondary beamline was significantly expanded. The toroidal correction plates were placed inside the vacuum chamber to increase the angle of toroidal correction of the secondary trajectories β3. Correcting plates were equipped with a continuous baking system that maintains their operating temperature at the level of 220–250 °C. This allows us to completely avoid the deposition of hydrocarbon films on the surface of the plates, and to exclude the appearance of high-voltage breakdowns between the plates or from the plate to the grounded T-10 wall. Immediately before the working shot of the tokamak, the baking was turned off and the baking circuit was disconnected in order to exclude the possibility of the interaction of the closed loop with the current and the confining magnetic fields of the machine. The electromagnetic forces resulting from this interaction could deform the plates, making it impossible for the secondary beam to pass into the analyzer. The described upgrades made it possible to expand the operational limits of measurements up to the operational limits of the tokamak 140 kA < Ipl < 330 kA.
  • The technique of scanning the entry angle of the particle beam into the plasma was implemented. Additionally, the vertical correction plates for secondary trajectories were installed to optimize the particle entry angle into the analyzer. It allowed us to obtain fragments of the radial profile per discharge;
  • A control system for the beam was created, which allows us to select a complete set of control voltages in the primary and secondary ion beamlines in each subsequent tokamak shot based on the analysis of the beam position in the previous shot (injection and sweeping angles in primary and secondary beamlines α1, α2, α3, β1, β2, β3). Its implementation radically simplified and accelerated the process of selecting control voltages, to significantly increase the accuracy of beam positioning, thereby providing systematic measurements both with a fixed SV position or with radial scan.
    Besides these two traditional HIBP operating modes, the developed system made it possible to operate in new non-standard modes. The most popular among them:
    • Periodic variation of the SV between two spatial positions (“colon”),
    • Periodic change of the SV between several positions (“multiple points”),
    • Alternation of the scan and at point measurements during one shot (“scan + point”).
  • A feedback for toroidal displacement was created, which adjust the toroidal correction voltage in the secondary beamline (toroidal angle β3), depending on the measured toroidal displacement in the detector ζd. This allows us to automate the selection of the toroidal correction voltage from shot to shot, as well as during one shot, especially in shots with changing plasma current (ramp-up, ramp-down) as well as the start of the plasma discharge with sharp increase of the plasma current. The feedback system for toroidal displacement is described in detail in the Reference [67].
  • The initial single-channel energy analyzer was first replaced by a two-channel and then a five-channel one. As a result, we can carry out correlation measurements over several spatial channels. In particular, in the region of maximal beam penetration into the plasma, it became possible to measure the poloidal potential correlations or Ep, a turbulent particle flux ΓE×B, as well as the poloidal density correlations or rotation of turbulence. Closer to the plasma edge, it is possible to measure the velocity of radial correlations or the radial propagation of potential and density perturbations.
  • An emitter-extractor unit has been developed, which makes it possible to significantly increase the primary beam current from 2–20 μA to 100–130 μA. Its use has expanded the allowable density limit for HIBP towards both high densities (plasma core) and also ultra-low densities at the plasma edge and in the scrape-off layer (SOL) [68].
The final version of the complex and its photo are shown in Figure 1 and Figure 2. The T-10 tokamak was terminated in March 2018 due to the construction of the T-15MD tokamak. Creation and development of the T-10 HIBP complex provides comprehensive study of Geodesic Acoustic Mode (GAM)—the toroidally symmetric structures in tokamak plasmas, see Section 5.1 and Section 6.3.

7.3. TJ-I HIBP

The HIBP complex for the TJ-I tokamak in the CIEMAT Science Center in Madrid, Spain, was created early in the 1990s by Krupnik, Bondarenko, Khrebtov and Nedzelskiy [69]. It operated with a Cs+ ion beam, accelerated up to energy Eb = 60 keV. The scheme of the complex is shown in Figure 24. Continuing the T-10 experiment, the TJ-I used the correcting toroidal plates for the secondary beam. As well as in T-10, the analyzer of the ‘flat mirror’ type was used on the TJ-I. The TJ-I tokamak had rectangular toroidal field coils. The TJ-I was the first, where the radial scanning method was used. The results of its application for Eb = 20 keV are shown in Figure 25. We see that the plasma density has a bell-shaped profile with a pronounced maximum in the center, which shifts outward along the major radius during the discharge. The distribution of the plasma potential along the radius is inhomogeneous and also evolves with the density: the density decrease is accompanied by potential evolution towards the positive values.
The TJ-I tokamak was terminated in 1995 due to the construction of the TJ-II stellarator. Creation and scientific exploitation of TJ-I HIBP complex made the basis for the next-step TJ-II HIBP diagnostics, the most fruitful one in studies of the symmetric structures in toroidal plasmas.

7.4. TUMAN-3M HIBP

The HIBP complex on the small aspect ratio tokamak TUMAN-3M in the A.F. Ioffe Institute in St. Petersburg was created in the late 1990s by Krupnik, Bondarenko and Komarov [70]. It operates with a beam of K+ ions accelerated up to the energy Eb = 60 keV. The scheme of the HIBP complex is shown in Figure 26. Figure shows that the detector line connects the center and edge of the plasma. It is suitable for determination of radial profiles over a total radial interval.
Figure 27 shows the equatorial projection of the calculated trajectories of the probing particles. We see that the plasma current shifts the trajectories, as well as that variation of the injection toroidal angle β effectively changes the displacement of the trajectories. Choosing β, we can place the detector line close to the meridian plane of the torus, z = 0. On the TUMAN-3M, it is possible to measure the evolution of the potential and plasma density at fixed spatial points.

7.5. WEGA HIBP

Stellarator WEGA is a two-period torsatron with a low toroidal magnetic field and rotational transform ι/2π. The HIBP complex on the WEGA stellarator in the Max Planck Institute for Plasma Physics in Greifswald, Germany was created early in the 2000s by Krupnik, Bondarenko and Komarov [71,72,73]. It operates with a beam of Na+ ions accelerated up to the energy Eb = 50 keV. The scheme of the HIBP complex is shown in Figure 28, and the detector grid is shown in Figure 29.
Stellarator WEGA was terminated in 2006 due to the construction of the W7-X stellarator. The scientific exploitation of WEGA HIBP brings an additional experience of the HIBP operation in stellarators, so contributes to the next-step TJ-II HIBP diagnostics, the most fruitful one in studies of the symmetric structures in toroidal plasmas.

7.6. URAGAN-2M HIBP

The HIBP complex on URAGAN-2M stellarator at the Institute of Plasma Physics of the Kharkov Institute of Physics and Technology of the Academy of Sciences of Ukraine was created in the mid-2010th by Krupnik, Kozachek, Komarov, Zhezhera and Chmyga [74,75]. It operates with a Cs+ ion beam accelerated to energies up to Eb = 70 keV. The scheme of the HIBP complex is shown in Figure 30a, and the equatorial projection of trajectories is shown in Figure 30b. URAGAN-2M was recommissioned in 2006 after a long pause caused by the reconstruction. It is currently undergoing the stage of adjustment and reaching the design parameters.

7.7. TJ-II HIBP—The Most Advanced Diagnostics to Study Plasma Symmetric Structures

The TJ-II stellarator-heliac has been successfully operating from the early 2000s to the present at the National Fusion Laboratory of the CIEMAT Science Center in Madrid, Spain. From the start of its operation till the start of the world largest stellarator W7-X in 2016, the TJ-II was the world’s second largest stellarator after the LHD and the largest in Europe. Now it is the third stellarator after W7-X and LHD.
The initial version of the HIBP complex on the TJ-II was created in the early 2000s in collaboration with Krupnik, Nedzelskiy, Khrebtov and Bondarenko from the Kharkov Institute of Physics and Technology (Ukraine) and the National Research Center Kurchatov Institute (Russia) [76,77]. It operates with a Cs+ ion beam accelerated up to the energy Eb = 150 keV.
Photos of the HIBP complex on the TJ-II stellarator are shown in Figure 31. Figure 32 presents the scheme of the HIBP complex, and the magnetic configuration with detector line, along which the measurements were performed, is shown in frame.
The TJ-II stellarator turned out to be the unique machine, for which the HIBP complex was planned at the stage of machine design, rather than its operation, unlike other HIBPs. For this reason, HIBP design was optimized in advance aiming for the unique properties: the ability to obtain full profiles of plasma parameters by scanning voltage variation in one discharge.
An advanced HIBP in TJ-II have passed several steps to become a diagnostics capable to resolve the issue of the symmetry/asymmetry of the toroidal plasma parameters.
The first step was done in resolving the issue of LFS–HFS asymmetry. The detector line, started with radial interval −1 < ρ < 0 in the HFS and producing many important results in the characterization of the potential profile in biasing experiments [78,79], electron confinement study [80,81] and transition to the improved confinement regime [82,83], was then extended towards the LFS, 0 < ρ < 1, see the inset to Figure 32. This way the first experimental test for the potential profile symmetry was performed. The results are shown in Figure 33. This observation was the first experimental prove of the potential symmetry LFS–HFS in the toroidal devices. In this experiment, plasma density was varied by an order of magnitude, from 0.3 to 3 × 1019 m−3. It was shown that plasma potential evolves dramatically with the density raise, changing the shape and the sign [84]. In the whole observed domain of plasma parameters (density, temperature, collisionality, way of heating: ECRH, NBI) the evolving plasma potential is symmetric in terms of LFS–HFS. Later on, the symmetric structure of plasma potential profiles was proven for more extended domain of plasma parameters [85,86].
Operating with the full detector line −1 < ρ < 1, on the TJ-II, the first measurements in stellarators were carried out, in which the dependence of the potential on the line-averaged density for a fixed position of the sample volume (SV) was studied, as well as the potential profiles and their dependence on heating methods and plasma parameters were measured [87,88]. At this stage there was made a substantial contribution to the study of helically symmetric structure—Alfven Eigenmodes. Figure 34 shows an example of the Helicity Induced Alfven Eigenmode (HAE) detected and identified by HIBP combined with the MHD-code AE3D [51]. The next principally important step in the widening of the diagnostic capabilities of HIBP was the measurement of the two-dimensional (2D) spatial distribution of the plasma potential and turbulence. This 2D map was obtained by the movement of the detector line across the plasma with the Eb variation [89,90].
To construct a 2D map, the systematic measurements were performed in a series of the reproducible shots. The results are presented in Figure 35. This 2D map of plasma potential is the first experimental prove of the potential poloidal symmetry in the toroidal devices. Figure 35a shows that the contours of the plasma potential (equipotentials) are well in line with the magnetic flux surfaces [91]. In contrast, the 2D map of the plasma potential and density turbulence is not poloidally symmetric, it has higher level at the HFS, as presented in Figure 35b.
It is worthwhile to highlight that these observations are contrasting to the gyrokinetic expectations for toroidal devices, always predicting the contrary: higher level of turbulence at the LFS [92]. This experimental test should help the theorists to look deeper at the theoretical basis for the modelling.
The final step of the paramount importance is a duplication of the HIBP diagnostics in one machine and creation of the Dual HIBP [93,94]. This system is aimed for the study of the long-range toroidal correlations in plasma parameters, and, more specifically, zonal flow [95,96]. It consists of the former complex, now called HIBP-1 and the new HIBP-2, separated toroidally in a quarter of the torus as presented in Figure 36. The photo of the Dual HIBP is shown in Figure 37.
Long-range toroidal correlations were found with the Dual HIBP in low-frequency plasma potential oscillations. Their radial distribution is presented in Figure 38. Figure 38 shows that there is a statistically significant coherence between plasma potential oscillations, measured by HIBP1 and HIBP2 [97,98]. The only low-frequency (<30 kHz) oscillations are coherent and having a zero cross-phase. Therefore, the observed phenomenon presents toroidally symmetric (n = 0) low-frequency potential perturbations, which is, by definition, the zonal flow [99].
In addition to study zonal flows [100], the double HIBP was also focused on the long-range correlations driven by helically symmetric structures—Alfvén eigenmodes [101,102,103,104,105] and induced by pellet injection [106,107].
The development and scientific exploitation of TJ-II dual HIBP have shown an importance of the symmetry and symmetry braking studies at the flux surface in toroidal devices. The results obtained erects the interest to widen such studies in other tokamaks and stellarators. The next section describes the HIBP projects under consideration for various devices in operation or in construction.

8. Future Prospects to Study of Symmetric Structures in Toroidal Plasmas—Conceptual Design of the HIBP Diagnostics for Various Toroidal Devices

The technical demands to the future HIBP projects are basically formulated in Section 2, Section 3 and Section 4. An experience of the operation and physical results obtained during the HIBP operation, discussed in Section 5, Section 6 and Section 7, form the new and more challenging level for the new projects. It consists of the following items:
(a)
Simultaneous measurements of the all three signals for potential, density, magnetic oscillations for comprehensive analysis of the plasma phenomena, including turbulence;
(b)
Multichannel measurements for correlation studies, including plasma turbulence poloidal rotation and radial propagation, and turbulent particle flux;
(c)
Maximal extension of the detector grid for 2D mapping and study of the poloidal symmetry/symmetry braking;
(d)
Creation of the dual HIBP system for study of toroidal/helical symmetry.
The projects briefly presented below are oriented to address the critical issues of HIBP compatibility to the machines, and HIBP capability to reply to the new challenges (a)–(d).

8.1. HIBP Design for the TCABR Tokamak

The conceptual design of the HIBP complex for TCABR tokamak at the Plasma Physics Laboratory of the University of São Paulo in Brazil was developed in the mid-2000s by Kuznetsov and Krupnik on the initiative of Severo and Nasimento [108]. It is assumed that the diagnostics will operate with a Tl ion beam accelerated up to the energy Eb = 105 keV. The TCABR tokamak has rectangular toroidal field coils. The scheme of the HIBP complex is shown in Figure 39a, and the detector grid is shown in Figure 39b. At the TCABR tokamak, the polarization of the plasma edge (biasing) and heating by Alfvén waves [109] are studied. For both experimental programs, 2D data on the plasma potential and its fluctuations are required.

8.2. HIBP Design for the Globus-M2 Tokamak

The conceptual design of the HIBP complex for the Globus-M2 spherical tokamak with an extremely small aspect ratio (R/a = 0.36 m/0.24 m = 1.5), operating in the A.F. Ioffe Physical Institute, St. Petersburg, was based at the initial design studies for the earlier version Globus-M [70,110]. It is assumed that the diagnostics will operate with a Tl ion beam accelerated up to Eb ≤ 45 keV. Figure 40 shows the scheme of the HIBP complex. Figure 40a shows both the vertical and the equatorial projections of the calculated trajectories of the probing particles. Figure 40b shows the detector grid aiming for the study of the 2D structures at the outer half of plasma radius in the LFS.
At the Globus-M2 tokamak, studies of plasma confinement under conditions of an extremely small aspect ratio and plasma heating with NBI and lower hybrid waves are conducted with focus to Alfvén eigenmodes. Most of experimental programs require data on the plasma potential and turbulence.

8.3. HIBP Design for the COMPASS Tokamak

The small aspect ratio COMPASS tokamak, operating at the Institute of Plasma Physics of the Czech Academy of Sciences in Prague, is used to study the L-H transition in discharges with ohmic and NBI plasma heating. The studies are focused to Alfvén eigenmodes and GAMs [111,112]. These experimental programs require data on the plasma potential and its fluctuations.
Conceptual design of the HIBP for the COMPASS tokamak (R = 0.56 m, a = 0.23 m, Ip < 360 kA, Bt ≤ 2.1, PNBI (at 40 keV) = 2 × 0.3 MW) with a, was created on the initiative of Prof. J. Stökel. COMPASS has a single-null divertor plasma configuration with elongation k = 1.8, triangularity ∆ = 0.2, and horizontal plasma size of 16 cm. The detector line, calculated in the classical probing scheme for and rectangular toroidal field coils, does not pass through the plasma center. Therefore, we consider another possible combination of probing ports, with the injection of primary particles through the lower inclined port and registration of secondary particles through the upper inclined port. The calculation and design results are shown in Figure 41. We see that the detector line connects the center and the plasma edge. Thus, this probing scheme is preferable. However, for realization of this scheme, the direction of the magnetic field Bt should be reversed.
Figure 42 shows a detector grid corresponding to the scheme in Figure 41.

8.4. HIBP Design for the TCV Tokamak

The TCV (Tokamak á Configuracion Variable) tokamak at the Plasma Physics Laboratory of the Federal Polytechnic School in Lausanne, Switzerland, is used to study dependencies of plasma confinement depending on the shape of the magnetic configuration, the L-H transition under conditions of ohmic and NBI plasma heating. Particular attention is paid to the study of Alfvén eigenmodes and GAM. These experimental programs require data on the plasma potential and its fluctuations.
The conceptual design of the HIBP complex for the TCV was developed on the initiative of Tonetti and Weisen [113]. It is assumed that the diagnostics will operate with a Tl ion beam accelerated to Eb = 105 keV. The TCV has rectangular toroidal field coils. The scheme of the HIBP complex is shown in Figure 43a, and the detector grid is shown in Figure 43b.

8.5. HIBP Design for the MAST Tokamak

The MAST (Mega-Ampere Spherical Torus) tokamak at the Culham Science Center in Culham, UK, is used to study the plasma confinement in the configuration of a spherical tokamak, L-H transition under conditions of ohmic and NBI plasma heating, and divertor physics. Particular attention is paid to the study of the Alfvén eigenmodes and Edge Localized Modes (ELMs). These experimental programs require data on the plasma potential and its fluctuations.
The conceptual design of the HIBP complex for MAST spherical tokamak was developed in on the initiative of Winsor and Sharapov [114]. It is assumed that the diagnostics will work with the Tl ion beam accelerated to Eb = 105 keV. The MAST tokamak has parameters R = 0.85 m, a = 0.65 m, Bt = 0.6 T, Ipl = 1 MA, an extremely small aspect ratio and rectangular toroidal field coils. The vertical section of the tokamak vacuum chamber is also rectangular. The scheme of HIBP is shown in Figure 44a [115], and one of the probing options and the detector grid is shown in Figure 44b. Note that several trajectories (blue) have intersections that complicate interpretation of HIBP measurements, while other trajectories (red) have not intersections. The similar situation was in COMPASS (Figure 41).

8.6. HIBP Design for the T-15MD Tokamak

The T-15 MD tokamak in the National Research Center Kurchatov Institute, Moscow, Russia, is commissioned in 2021. T-15 MD intended to study the plasma confinement in conditions of a strong magnetic field (Bt ≤ 2 T) combined with a small aspect ratio (R/a = 1.48 m/0.67 m). It is planned to use ohmic, ECRH, ICRH and NBI plasma heating and current drive (CD). Particular attention will be paid to the study of the role of the electric field in L-H transitions, Alfvén eigenmodes, and GAMs. These experimental programs require 2D data on the plasma potential, density, magnetic field and its fluctuations.
Currently, conceptual design of the HIBP complex for T-15 MD is developed as a part of the diagnostic complex of the machine [116]. The HIBP will use the Tl ion beam with Eb = 400 keV. T-15 MD magnetic system was taken into account in the detailed trajectory calculations [117]. The scheme of the HIBP complex is shown in Figure 45a [118], and the detector grid is shown in Figure 45b.
In the future, it is planned to upgrade diagnostics and create a double HIBP, consisting of two identical complexes, shifted along the torus by about 90°. Thus, the dual HIBP on the T-15 MD tokamak will be basically similar to the dual HIBP on the TJ-II stellarator. The main task of the dual HIBP on T-15 MD is to study zonal flows and long-range correlations of potential and density at the frequencies of GAMs [119,120] and Alfvén eigenmodes [121] including the effects of the ECRH and CD [122]. On top of that, the special attention will be dedicated to the ICRH and CD [123,124].

8.7. HIBP Design for the W7-X Stellarator

The conceptual design of the HIBP complex for the world’s largest stellarator W7-X was developed in the early 2000s. This installation has the following design parameters: Bt = 3 T, R = 5.5 m, a = 0.53 m. W7-X was commissioned at the end of 2015 in the Max Planck Institute for Plasma Physics in Greifswald, Germany. W7-X intended to study the plasma confinement in the optimized 3D magnetic configuration. It uses microwave and NBI for plasma heating. Particular attention is focused on study the role of the electric field at L-H transitions and Alfvén eigenmodes. These experimental programs require data on the plasma potential and its fluctuations. A preliminary calculation of trajectories for the proposed W7-X plasma probing scheme is shown in Figure 46 [125,126].

8.8. HIBP Design for the International Experimental Tokamak Reactor ITER

The ITER, which is currently being assembled in the city of Cadarache in the south of France, aims to achieve the fusion power amplification factor Q = 5. The ITER diagnostics are subordinated to this engineering task. They are primarily focused to achieving the design parameters of the machine and provide the safety of its operation. So, HIBP is not included in the list of top-priority diagnostics for ITER, but it may be used in studying the physics of the L-H transition. Therefore, the HIBP allows one to study the periphery of the ITER plasma. Estimations have shown that HIBP will operate with a Tl ion beam accelerated to the energy of Eb = 3.5 MeV [127]. This is a standard energy range for large-scale machines in operation [19]. ITER has ovoid-shaped toroidal field coils. The scheme of the HIBP complex is shown in Figure 47, and the detector grid is shown in Figure 48a. In the standard ITER mode, the nominal plasma current will lead to a noticeable toroidal shift of the trajectories of the probing particles, shown in Figure 48b. We see that with the initial toroidal shift zI = 0 at the injection point I and the value of the toroidal angle β = 7°, the toroidal coordinate of the detection point D, zD = 10 cm, and easily fits into the dimensions of the horizontal port. Because in ITER the discharge will transit to the H-mode that is always accompanied by change in the electric field [52,128,129], the measurements of the electric potential and turbulence will be an important issue for the physics of ITER. Recently, it was first shown by HIBP in the T-10 tokamak that in the discharges with low collisionality—which mimic the ITER plasma condition—the core plasma potential has positive values and an enhanced level of turbulence, so the prediction for the ITER plasma was made [130]. Verification of this prediction could be a very important and challenging task for physics studies in ITER.

9. Summary

Over more than 30 years of development, diagnostics of fusion plasma using a beam of heavy ions has become a powerful and multifunctional tool for physical research at the magnetic confinement devices. Starting with a time-averaged potential signal measured in the single spatial point, it becomes the diagnostic for studying 2D spatial distributions of plasma potential and turbulence, and even 3D characteristics of turbulence. On modern operating medium-scale machines, this diagnostic makes a significant contribution to reveal the role of electric fields and turbulence in confinement and the transport of plasma energy and particles, physics of fast particles and Alfvén eigenmodes. These advances have allowed HIBP to study the spatial and frequency characteristics of toroidally symmetric structures: geodesic acoustic modes and zonal flows, and helically symmetric structures: Alfvén Eigenmodes. The future studies will be focused on the further development of operating HIBP capabilities and the installations to the other considered toroidal devices.

Funding

This research was funded by Russian Science Foundation, project number 19–12-00312 and supported in part by Competitiveness Program of MEPhI.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

I warmly acknowledge the long-term collaboration with my colleagues from HIBP group and T-10 tokamak in Kurchatov Institute, Moscow, Russia, HIBP group of Kharkov Institute of Physics and Technology, Kharkov, Ukraine, leaded by L.I. Lrupnik and A.S. Kozachek, and also colleagues from CIEMAT, Madrid, Spain, leaded by Carlos Hidalgo. Special words of thanks come to my teacher, Yu. N. Dnestrovskij and Sergey E. Lysenko for fruitful discussions and helpful comments.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Tikhonov, A.N.; Arsenin, V.V. Solutions of Ill Posed Problems; John, F., Translator, V.H., Eds.; Winston and Sons: Silver Spring, MD, USA; Wiley: New York, NY, USA, 1977. [Google Scholar]
  2. Dnestrovskij, Y.N.; Kostomarov, D.P.; Melnikov, A.V. Mathematical problems of active corpuscular diagnostics of plasmas. Proc. Moscov State Univ. Ser. 15 Appl. Math. Cybern. 1986, 3, 20–28. (In Russian) [Google Scholar]
  3. Melnikov, A.V. Applied and fundamental aspects of fusion science. Nat. Phys. 2016, 12, 386–390. [Google Scholar] [CrossRef]
  4. Van Oost, G.; Bulanin, V.V.; Donne, A.J.H.; Gusakov, E.; Kraemer-Flecken, A.; Krupnik, L.I.; Melnikov, A.; Nanobashvili, S.; Peleman, P.; A Razumova, K.; et al. Multi-machine studies of the role of turbulence and electric fields in the establishment of improved confinement in tokamak plasmas. Plasma Phys. Control. Fusion 2007, 49, A29–A44. [Google Scholar] [CrossRef]
  5. Kadomtsev, B.B. Collective Phenomena in Plasmas; Pergamon Press: Oxford, UK, 1982; ISBN 10-0080250182/13-978-0080250182. [Google Scholar]
  6. Burrell, K.H.; Carlstrom, T.N.; Doyle, E.J.; Finkenthal, D.; Gohil, P.; Groebner, R.J.; Hillis, D.L.; Kim, J.; Matsumoto, H.; Moyer, R.A.; et al. Physics of the L-mode to H-mode transition in tokamaks. Plasma Phys. Control. Fusion 1992, 34, 1859. [Google Scholar] [CrossRef]
  7. Tendler, M. Major achievements and challenges of fusion research. Phys. Scr. 2015, 90, 98002. [Google Scholar] [CrossRef]
  8. Melnikov, A.V. The study of electric potential of magnetically confined fusion plasmas. Comput. Nanotechnol. 2017, 2, 13–23. Available online: https://www.elibrary.ru/item.asp?id=29368049 (accessed on 1 January 2021). (In Russian).
  9. Melnikov, A.V. The evolution of the mathematical and technical aspects of heavy ion beam probe plasma diagnostics. Comput. Nanotechnol. 2018, 1, 31–57. Available online: https://www.elibrary.ru/item.asp?id=32775187n (accessed on 1 January 2021). (In Russian).
  10. Dnestrovskij, Y.; Melnikov, A.; Krupnik, L.; Nedzelskij, I. Development of heavy ion beam probe diagnostics. IEEE Trans. Plasma Sci. 1994, 22, 310–331. [Google Scholar] [CrossRef]
  11. Jobes, F.; Hickok, R. A direct measurement of plasma space potential. Nucl. Fusion 1970, 10, 195–197. [Google Scholar] [CrossRef]
  12. Bugarya, V.I.; Gorshkov, A.V.; Grashin, S.A.; Ivanov, I.V.; Krupin, V.A.; Krupnik, L.I.; Nedzel’skij, I.S.; Mel’nikov, A.V.; Razumova, K.A.; Sokolov, Y.A.; et al. Electric potential and toroidal and poloidal rotation velocities of a tokamak plasma. JETP Lett. 1983, 38, 404–408. [Google Scholar]
  13. Bugarya, V.; Gorshkov, A.; Grashin, S.; Ivanov, I.; Krupin, V.; Mel′nikov, A.; Razumova, K.; Sokolov, Y.; Trukhin, V.; Chankin, A.; et al. Measurements of plasma column rotation and potential in the TM-4 tokamak. Nucl. Fusion 1985, 25, 1707–1717. [Google Scholar] [CrossRef]
  14. Donné, A.J.H.; Melnikov, A.V.; Van Oost, G. Diagnostics for radial electric field measurements in hot magnetized plasmas. Czechoslov. J. Phys. 2002, 52, 1077–1096. [Google Scholar] [CrossRef]
  15. Melnikov, A.; Krupnik, L.; Eliseev, L.; Barcala, J.M.; Bravo, A.; Chmyga, A.; Deshko, G.; Drabinskij, M.; Hidalgo, C.; Khabanov, P.; et al. Heavy ion beam probing—diagnostics to study potential and turbulence in toroidal plasmas. Nucl. Fusion 2017, 57, 072004. [Google Scholar] [CrossRef]
  16. Askinazi, L.G.; Kornev, V.A.; Lebedev, S.V.; Tukachinsky, A.S.; Zhubr, N.A.; Dreval, N.B.; Krupnik, L.I. Heavy ion beam probe development for the plasma potential measurement on the TUMAN-3M tokamak. Rev. Sci. Instrum. 2004, 75, 3517–3519. [Google Scholar] [CrossRef]
  17. Cabral, J.A.C.; Malaquias, A.; Praxedes, A.; Van Toledo, W.; Varandas, C. The heavy ion beam diagnostic for the tokamak ISTTOK. IEEE Trans. Plasma Sci. 1994, 22, 350–358. [Google Scholar] [CrossRef]
  18. Bondarenko, I.S.; Gubarev, S.P.; Krupnik, L.I.; Mel’nikov, A.V.; Nedszel’skii, I.S.; Utkin, A.A.; Khrebtov, S.M. Plasma diagnostics using a heavy-ion beam in the T-10 tokamak. Sov. J. Plasma Phys. 1992, 18, 110. [Google Scholar]
  19. Fujisawa, A.; Iguchi, H.; Taniike, A.; Sasao, M.; Hamada, Y. A 6 MeV heavy ion beam probe for the Large Helical Device. IEEE Trans. Plasma Sci. 1994, 22, 395–402. [Google Scholar] [CrossRef] [Green Version]
  20. Ido, T.; Shimizu, A.; Nishiura, M.; Kato, S.; Nakano, H.; Ohshima, S.; Yokoyama, M.; Murakami, S.; Wakasa, A.; Nakamura, S.; et al. Development of 6-MeV Heavy Ion Beam Probe on LHD. Fusion Sci. Technol. 2010, 58, 436–444. [Google Scholar] [CrossRef]
  21. Melnikov, A.; Hidalgo, C.; Ido, T.A.; Shimizu, A.; Fujisawa, K.; Dyabilin, S.; Lysenko, S.E. Plasma Potential in Toroidal Devices: T-10, TJ-II, CHS and LHD. Plasma Fusion Res. 2012, 7, 2402114. [Google Scholar] [CrossRef] [Green Version]
  22. Bondarenko, I.S.; Chmuga, A.A.; Dreval, N.B.; Khrebtov, S.M.; Komarov, A.D.; Kozachok, A.S.; Krupnik, L.I.; Coelho, P.; Cunha, M.; Gonçalves, B.; et al. Installation of an advanced heavy ion beam diagnostic on the TJ-II stellarator. Rev. Sci. Instrum. 2001, 72, 583–585. [Google Scholar] [CrossRef]
  23. Melnikov, A.; Dyabilin, K.; Eliseev, L.; Lysenko, S.; Dnestrovskij, Y. Measurement and modelling of electric potential in the TJ-II stellarator. Probl. At. Sci. Technol. Ser. Thermonucl. Fusion 2011, 34, 54–73. (In Russian) [Google Scholar] [CrossRef]
  24. Lei, J.; Shah, U.; Demers, D.R.; Connor, K.A.; Schoch, P.M. Calibration and initial operation of the HIBP on the MST. Rev. Sci. Instrum. 2001, 72, 564–567. [Google Scholar] [CrossRef]
  25. Yoshikawa, M.; Sakamoto, M.; Miyata, Y.; Aoyama, M. Potential fluctuation study from the core plasma to end region in GAMMA 10. Nucl. Fusion 2013, 53, 073031. [Google Scholar] [CrossRef]
  26. Dreval, N.; Krupnik, L.; Chmyga, A.A.; Khrebtov, S.M.; Komarov, A.D.; Kozachek, A.S.; Hidalgo, C.; Melnikov, A.V.; Eliseev, L.G. Features of HIBP diagnostics application to stellarator-like devices. Probl. At. Sci. Technol. Ser. Plasma Phys. 2005, 2, 223–225. [Google Scholar]
  27. Dnestrovskii, Y.N.; Krupnik, L.I.; Mel’nikov, A.V.; Nedzel’skii, I.S. Density determination by probing with heavy-ion beams. Sov. J. Plasma Phys. 1986, 12, 130. [Google Scholar]
  28. Schwelberger, J.G.; Aceto, S.C.; Connor, K.A.; Zielinski, J.J.; Baylor, L.A.; England, A.C.; Isler, R.; Ma, C.H.; Murakami, M.; Uckan, T.; et al. Electron density profile measurement with a heavy ion beam probe. Rev. Sci. Instrum. 1992, 63, 4571–4573. [Google Scholar] [CrossRef]
  29. Khabanov, P.; Eliseev, L.; Melnikov, A.; Drabinskiy, M.; Hidalgo, C.; Kharchev, N.; Chmyga, A.; Kozachek, A.; Pastor, I.; De Pablos, J.; et al. Density profile reconstruction using HIBP in ECRH plasmas in the TJ-II stellarator. J. Instrum. 2019, 14, C09033. [Google Scholar] [CrossRef]
  30. Levinton, F.M. The motional Stark effect: Overview and future development (invited). Rev. Sci. Instrum. 1999, 70, 810–814. [Google Scholar] [CrossRef]
  31. Dnestrovskii, Y.N.; Mel’nikov, A.V. Determination of the current profile in a plasma by probing with heavy-ion beams. Sov. J. Plasma Phys. 1986, 12, 393. [Google Scholar]
  32. Weisen, H.; Melnikov, A.V.; Perfilov, S.; Lysenko, S. On the Possibility of Using a Heavy Ion Beam Probe for Local Poloidal Flux Measurements in a Tokamak. Fusion Sci. Technol. 2011, 59, 418–426. [Google Scholar] [CrossRef]
  33. Melnikov, A.V. Electric Potential in Toroidal Plasmas; Springer Nature Switzerland AG: Basel, Switzerland, 2019; p. 240. ISBN 978-3-030-03480-1. [Google Scholar]
  34. Khabanov, P.O. Computational analysis density fluctuations measurements with heavy ion beam probe. Comput. Nanotechnol. 2018, 1, 25–30. Available online: https://www.elibrary.ru/item.asp?id=32775186 (accessed on 1 January 2021). (In Russian).
  35. Melnikov, A.V.; Eliseev, L.G.; Jiménez-Gómez, R.; Ascasibar, E.; Hidalgo, C.; Chmyga, A.A.; Ido, T.; Khrebtov, S.M.; Könies, A.; Komarov, A.D.; et al. Study of Alfvén eigenmodes in the TJ-II stellarator. Plasma Fusion Res. 2010, 5, S2019. [Google Scholar] [CrossRef] [Green Version]
  36. Simcic, V.J.; Crowley, T.P.; Schoch, P.M. Internal magnetic and electrostatic fluctuation measurements of MHD modes in the TEXT tokamak. Phys. Fluids B 1993, 5, 1576. [Google Scholar] [CrossRef]
  37. McKee, G.R.; Fonck, R.J.; Jakubowski, M.; Burrell, K.H.; Hallatschek, K.; A Moyer, R.; Nevins, W.; Rudakov, D.L.; Xu, X. Observation and characterization of radially sheared zonal flows in DIII-D. Plasma Phys. Control. Fusion 2003, 45, A477–A485. [Google Scholar] [CrossRef] [Green Version]
  38. Zenin, V.N. Geodesic acoustic modes in the tokamaks. Comput. Nanotechnol. 2018, 1, 108–113. Available online: https://www.elibrary.ru/item.asp?id=32775193 (accessed on 1 January 2021). (In Russian).
  39. Zenin, V.N.; Eliseev, L.G.; Kozachek, A.S.; Krupnik, L.I.; Lysenko, S.E.; Melnikov, A.V.; Perfilov, S.V. Study of poloidal structure of geodesic acoustic modes in the T-10 tokamak with heavy ion beam probing. Probl. At. Sci. Technol. Ser. Plasma Phys. 2014, 6, 269–271. Available online: http://vant.kipt.kharkov.ua/TABFRAME.html (accessed on 1 January 2021).
  40. Eliseev, L.G.; Ivanov, N.V.; Kakurin, A.M.; Lysenko, S.E.; Maltsev, S.G.; Melnikov, A.V.; Perfilov, S.V.; Zenin, V.N. Study of the large-scale MHD mode and its effect on GAM in the T-10 tokamak. In Proceedings of the 42nd EPS Conference on Plasma Physics, Lisbon, Portugal, 22—26 June 2015; ECA. Volume 39E, p. P5.159. Available online: http://ocs.ciemat.es/EPS2015PAP/Pdf/P5.159.pdf (accessed on 1 January 2021).
  41. Melnikov, A.V.; Eliseev, L.G.; Lysenko, S.E.; Perfilov, S.V.; Shelukhin, D.A.; Vershkov, V.A.; Zenin, V.N.; Krupnik, L.I.; Kharchev, N.K. Correlation properties of Geodesic Acoustic Modes in the T-10 tokamak. J. Phys. Conf. Ser. 2015, 591, 012003. [Google Scholar] [CrossRef] [Green Version]
  42. Melnikov, A.V.; Eliseev, L.G.; Gudozhnik, A.V.; Lysenko, S.E.; Mavrin, V.A.; Perfilov, S.V.; Zimeleva, L.G.; Ufimtsev, M.V.; Krupnik, L.I.; Schoch, P.M. GAM evolution with ECRH Investigation of the plasma potential oscillations in the range of geodesic acoustic mode frequencies by heavy ion beam probing in tokamaks. Czech. J. Phys. 2005, 55, 349–360. [Google Scholar] [CrossRef] [Green Version]
  43. Breizman, B.N.; Pekker, M.S.; Sharapov, S.E. Plasma pressure effect on Alfvén cascade eigenmodes. Phys. Plasmas 2005, 12, 112506. [Google Scholar] [CrossRef] [Green Version]
  44. Eliseev, L.G.; Melnikov, A.V.; Ascasibar, E.; Cappa, A.; Drabinskiy, M.; Hidalgo, C.; Khabanov, P.O.; Kharchev, N.K.; Kozachek, A.S.; Liniers, M. Experimental observation of the geodesic acoustic frequency limit for the NBI-driven Alfvén eigenmodes in TJ-II. Phys. Plasmas 2021, 28. [Google Scholar] [CrossRef]
  45. Melnikov, A.V.; Eliseev, L.G.; Jiménez-Gómez, R.; Ascasibar, E.; Hidalgo, C.; Chmyga, A.A.; Komarov, A.D.; Kozachok, A.S.; Krasilnikov, I.A.; Khrebtov, S.M.; et al. Internal measurements of Alfvén eigenmodes with heavy ion beam probing in toroidal plasmas. Nucl. Fusion 2010, 50, 084023. [Google Scholar] [CrossRef]
  46. Melnikov, A.; Eliseev, L.; Lysenko, S.E.; Ufimtsev, M.; Zenin, V. Study of interactions between GAMs and broadband turbulence in the T-10 tokamak. Nucl. Fusion 2017, 57, 115001. [Google Scholar] [CrossRef] [Green Version]
  47. Gryaznevich, M.; Stockel, J.; Van Oost, G.; Del Bosco, E.; Svoboda, V.; Melnikov, A.; Kamendje, R.; Malaquias, A.; Mank, G.; Miklaszewski, R.; et al. Contribution of joint experiments on small tokamaks in the framework of IAEA coordinated research projects to mainstream fusion research. Plasma Sci. Technol. 2020, 22, 055102. [Google Scholar] [CrossRef]
  48. Riggs, G.A.; Nogami, S.; Koepke, M.; Melnikov, A.; Eliseev, L.; Lysenko, S.; Khabanov, P.; Drabinskij, M.; Kharchev, N.; Kozachek, A.; et al. Bispectral analysis of broadband turbulence and geodesic acoustic modes in the T-10 tokamak. J. Plasma Phys. 2021, 87. [Google Scholar] [CrossRef]
  49. Powers, E. Spectral techniques for experimental investigation of plasma diffusion due to polychromatic fluctuations. Nucl. Fusion 1974, 14, 749–752. [Google Scholar] [CrossRef]
  50. Melnikov, A.V.; Eliseev, L.G.; Alonso, A.; Ascasiba, E.; Cappa, A.; Chmyga, A.A.; Estrada, T.; Hidalgo, C.; Jiménez-Gómez, R.; Krasilnikov, I.; et al. Plasma potential profiles and oscillations in ECRH and NBI plasmas in the TJ-II stellarator. In Proceedings of the 37th EPS Conference on Plasma Physics, Dublin, Ireland, 21–25 June 2010; ECA. Volume 34A, p. Rep. P1.066. Available online: http://ocs.ciemat.es/EPS2010PAP/html/ (accessed on 1 January 2021).
  51. Melnikov, A.V.; Eliseev, L.G.; Ascasibar, E.; Chmyga, A.A.; Hidalgo, C.; Ido, T.; Jiménez-Gomez, R.; Komarov, A.D.; Kozachek, A.S.; Krupnik, L.I.; et al. Alfvén eigenmode properties and dynamics in the TJ-II stellarator. Nucl. Fusion 2012, 52, 123004. [Google Scholar] [CrossRef]
  52. Melnikov, A.; Eliseev, L.; Estrada, T.; Ascasibar, E.; Alonso, A.; Chmyga, A.; Hidalgo, C.; Komarov, A.; Kozachek, A.; Krupnik, L.; et al. The changes in plasma potential and turbulent particle flux in the core plasma measured by heavy ion beam probe during L-H transitions in the TJ-II stellarator. Nucl. Fusion 2013, 53, 092002. [Google Scholar] [CrossRef]
  53. Eliseev, L.G.; Lysenko, S.E.; Melnikov, A.V.; Krupnik, L.I.; Kozachek, A.S.; Zenin, V.N. Study of GAMs and related turbulent particle flux with HIBP in the T-10 tokamak. Probl. At. Sci. Technol. Ser. Plasma Phys. 2017, 107, 241–243. [Google Scholar]
  54. Eliseev, L.G.; Zenin, V.N.; Lysenko, S.E.; Melnikov, A.V. Measurement of geodesic acoustic modes and the turbulent particle flux in the T-10 tokamak plasmas. J. Phys. Conf. Ser. 2017, 907, 12002. [Google Scholar] [CrossRef]
  55. Eliseev, L.G.; Melnikov, A.V.; Lysenko, S.E.; Khabanov, P.O.; Zenin, V.N.; Drabinskij, M.A.; Kharchev, N.K.; Kozachek, A.S.; Krupnik, L.I.; Team, H. Evaluation of Turbulent Particle Flux by Heavy Ion Beam Probe in the T-10 Tokamak. Plasma Fusion Res. 2018, 13, 3402106. [Google Scholar] [CrossRef] [Green Version]
  56. Eliseev, L.; Melnikov, A.; Perfilov, S.; Hidalgo, C. Two point correlation technique for the measurements of poloidal plasma rotation by heavy ion beam probe. Plasma Fusion Res. 2012, 7, 2402064. [Google Scholar] [CrossRef] [Green Version]
  57. Drabinskiy, M.A.; Eliseev, L.G.; Khabanov, P.O.; Melnikov, A.V.; Khartchev, N.K.; Sergeev, N.S.; Grashin, S.A. Radial structure of quasi-coherent mode in ohmic plasma of the T-10 tokamak. J. Phys. Conf. Ser. 2019, 1383, 012004. Available online: https://iopscience.iop.org/article/10.1088/1742-6596/1383/1/012004 (accessed on 1 January 2021). [CrossRef]
  58. Melnikov, A.; Eliseev, L.; Perfilov, S.; Lysenko, S.; Shurygin, R.; Zenin, V.; Grashin, S.; Krupnik, L.; Kozachek, A.; Solomatin, R.; et al. The features of the global GAM in OH and ECRH plasmas in the T-10 tokamak. Nucl. Fusion 2015, 55, 063001. [Google Scholar] [CrossRef]
  59. Bondarenko, I.S.; Krupnik, L.I.; Nedzel’skii, I.S.; Mel’nikov, A.V. Plasma probing with a heavy ion beam in the TM-4 tokamak. Sov. J. Tech. Phys. 1986, 31, 1390. [Google Scholar]
  60. Melnikov, A.V.; Bondarenko, I.S.; Efremov, S.L.; Kharchev, N.K.; Khrebotv, S.M.; Krupnik, L.I.; Nedzelskij, I.S.; Zimeleva, L.G.; Trofimenko, Y.V. HIBP diagnostics on T-10. Rev. Sci. Instrum. 1995, 66, 317. [Google Scholar] [CrossRef]
  61. Melnikov, A.V.; Kharchev, N.K.; Eliseev, L.G.; Bondarenko, I.S.; Krupnik, L.I.; Khrebtov, S.M.; Nedzelskij, I.S.; Trofimenko, Y.V. Calibration of the heavy ion beam probe parallel plate analyzer using the gas target and reference beam. Rev. Sci. Instrum. 1997, 68, 308–311. [Google Scholar] [CrossRef]
  62. Krupnik, L.; Melnikov, A.; Nedselskij, I. Development of beam probe diagnostics and recent measurements on the TJ-1 and T-10 tokamaks. Fusion Eng. Des. 1997, 34–35, 639–644. [Google Scholar] [CrossRef]
  63. Melnikov, A.V.; Krupnik, L.I.; Nedzelskij, I.S.; Kharchev, N.K. Plasma potential measurements in the T-10 tokamak. In Proceedings of the 18th EPS Conference on Control. Fusion and Plasma Physics, Berlin, Germany, 3–7 June 1991; Volume 15C, p. 221. [Google Scholar]
  64. Melnikov, A.V.; Zimeleva, L.G.; Eliseev, L.G.; Bondarenko, I.S.; Krupnik, L.I.; Khrebtov, S.M.; Nedzelskij, I.S.; Kharchev, N.K. HIBP measurements of the plasma electric potential on T-10. In Proceedings of the 23rd EPS Conference on Controlled Fusion and Plasma Physics, Kiev, Ukraine, 24–28 June 1996; Volume I, p. 434. [Google Scholar]
  65. Melnikov, A.V.; Tarasyan, K.N.; Vershkov, V.A.; Dreval, V.A.; Sushkov, A.V.; Lysenko, S.E.; Sannikov, V.V.; Gorshkov, A.V.; Krupnik, L.I.; Nedzelskij, I.S. Space and time evolution of plasma potential in T-10 under variation of main gas influx. IEEE Trans. Plasma Sci. 1994, 22, 363–368. [Google Scholar] [CrossRef]
  66. Melnikov, A.; Drabinskiy, M.; Eliseev, L.; Khabanov, P.; Kharchev, N.; Krupnik, L.; De Pablos, J.; Kozachek, A.; Lysenko, S.; Molinero, A.; et al. Heavy ion beam probe design and operation on the T-10 tokamak. Fusion Eng. Des. 2019, 146, 850–853. [Google Scholar] [CrossRef]
  67. Kharchev, N.K. Correction of toroidal displacement of diagnostic heavy ion beam. Comput. Nanotechnol. 2018, 1, 58–61. Available online: https://www.elibrary.ru/item.asp?id=32775188 (accessed on 1 January 2021). (In Russian).
  68. Drabinskiy, M.A. The study of the electric potential of magnetically confined fusion plasmas. Comput. Nanotechnol. 2018, 1, 62–70. Available online: https://www.elibrary.ru/item.asp?id=32775189 (accessed on 1 January 2021). (In Russian).
  69. Bondarenko, I.S.; Khrebtov, S.M.; Krupnik, L.I.; Nedzelskij, I.S.; Gordeev, O.A.; Kharchev, N.K.; Melnikov, A.V.; Tarasyan, K.N.; Zimeleva, L.G.; Hidalgo, C.; et al. Heavy ion beam probe diagnostics on TJ-1 tokamak and the measurements of the plasma potential and density profiles. Rev. Sci. Instruments 1997, 68, 312–315. [Google Scholar] [CrossRef]
  70. Melnikov, A.V.; Zimeleva, L.G.; Krupnik, L.I.; Nedzelskij, I.S.; Trofimenko, Y.V.; Minaev, V. Heavy ion beam probe systems for tight aspect ratio tokamaks. Rev. Sci. Instrum. 1997, 68, 316–319. [Google Scholar] [CrossRef]
  71. Krupnik, L.I.; Deshko, G.N.; Zhezhera, A.I.; Chmyga, A.A.; Komarov, A.D.; Kozachek, A.S.; Melnikov, A.V.; Perfilov, S.V.; Otte, M.; Shubert, M. The heavy ion beam probing development for WEGA stellarator. Fusion Sci. Technol. 2006, 50, 276–280. [Google Scholar] [CrossRef] [Green Version]
  72. Podoba, Y.; Otte, M.; Wagner, F.; Bondarenko, I.; Chmyga, A.; Deshko, G.; Komarov, A.; Kozachek, A.; Krupnik, L.; Khrebtov, S.; et al. First HIBP results on the WEGA stellarator. AIP Conf. Proc. 2008, 993, 235–238. [Google Scholar] [CrossRef] [Green Version]
  73. Bondarenko, I.S.; Chmyga, A.A.; Deshko, G.N.; Komarov, A.D.; Kozachok, A.S.; Krupnik, L.I.; Khrebtov, S.M.; Zhezhera, A.I.; Podoba, Y.Y.; Otte, M.; et al. HIBP results on the WEGA stellarator. Probl. At. Sci. Technol. Ser. Plasma Phys. 2009, 15, 28–30. [Google Scholar]
  74. Krupnik, L.I.; Mel’nikov, A.V.; Nedzel’skii, I.S.; Samokhvalov, N.V.; Schwelberger, J.; Connor, K.A. Plasma diagnostics using a heavy-ion beam in the U-2M stellarator. Plasma Phys. Rep. 1994, 20, 173–175. [Google Scholar]
  75. Krupnik, L. 50 years of hot plasma diagnostic with heavy ion beam probing (HIBP) at the Kharkov institute of physics and technology. Probl. At. Sci. Technol. Ser. Plasma Phys. 2021, 154–162. [Google Scholar] [CrossRef]
  76. Bondarenko, I.S.; Chmyga, A.A.; Dreval, N.B.; Khrebtov, S.M.; Komarov, A.D.; Kozachok, A.S.; Krupnik, L.I.; Melnikov, A.V.; Yudina, O.A.; Coelho, P.; et al. Installation of the advanced heavy ion beam probing diagnostic on the TJ-II stellarator. Czechoslov. J. Phys. 2000, 50, 1397–1412. [Google Scholar] [CrossRef]
  77. Krupnik, L.; Bondarenko, I.; Chmyga, A.; Dreval, M.; Khrebtov, S.; Komarov, A.; Kozachok, A.; Hidalgo, C.; García-Cortés, I.; Rodriguez-Rodrigo, L.; et al. The first operation of the advanced heavy ion beam probing diagnostic on the TJ-II flexible heliac. Fusion Eng. Des. 2001, 56–57, 935–939. [Google Scholar] [CrossRef]
  78. Melnikov, A.V.; Hidalgo, C.; Chmyga, A.A.; Dreval, N.B.; Eliseev, L.; Khrebtov, S.M.; Komarov, A.D.; Kozachok, A.S.; Krupnik, L.I.; Pastor, I.; et al. Plasma potential measurements by the heavy ion beam probe diagnostic in fusion plasmas: Biasing experiments in the TJ-II stellarator and T-10 tokamak. Fusion Sci. Technol. 2004, 46, 299–311. [Google Scholar] [CrossRef]
  79. Hidalgo, C.; Pedrosa, M.A.; Dreval, N.; McCarthy, K.J.; Eliseev, L.; Ochando, M.A.; Estrada, T.; Pastor, I.; Ascasíbar, E.; Calderón, E.; et al. Improved confinement regimes induced by limiter biasing in the TJ-II stellarator. Plasma Phys. Control. Fusion 2003, 46, 287–297. [Google Scholar] [CrossRef]
  80. Estrada, T.; Krupnik, L.; Dreval, N.; Melnikov, A.; Khrebtov, S.M.; Hidalgo, C.; van Milligen, B.; Castejón, F.; Ascasíbar, E.; Eliseev, L.; et al. Electron internal transport barrier formation and dynamics in the plasma core of the TJ-II stellarator. Plasma Phys. Control. Fusion 2003, 46, 277–286. [Google Scholar] [CrossRef] [Green Version]
  81. Estrada, T.; Alonso, A.; Chmyga, A.A.; Dreval, N.; Eliseev, L.; Hidalgo, C.; Komarov, A.D.; Kozachok, A.S.; Krupnik, L.; Melnikov, A.V.; et al. Electron internal transport barriers, rationals and quasi-coherent oscillations in the stellarator TJ-II. Plasma Phys. Control. Fusion 2005, 47, L57–L63. [Google Scholar] [CrossRef] [Green Version]
  82. Estrada, T.; Happel, T.; Eliseev, L.; López-Bruna, D.; Ascasíbar, E.; Blanco, E.; Cupido, L.; Fontdecaba, J.M.; Hidalgo, C.; Jiménez-Gómez, R.; et al. Sheared flows and transition to improved confinement regime in the TJ-II stellarator. Plasma Phys. Control. Fusion 2009, 51. [Google Scholar] [CrossRef]
  83. Van Milligen, B.P.; Pedrosa, M.A.; Hidalgo, C.; Carreras, B.A.; Estrada, T.; Alonso, J.A.; de Pablos, J.L.; Melnikov, A.; Krupnik, L.; Eliseev, L.G.; et al. The dynamics of the formation of the edge particle transport barrier at TJ-II. Nucl. Fusion 2011, 51, 113002. [Google Scholar] [CrossRef]
  84. Melnikov, A.V.; Alonso, A.; Ascasíbar, E.; Balbin, R.; Chmyga, A.A.; Dnestrovskij, Y.N.; Eliseev, L.; Estrada, T.; Fontdecaba, J.M.; Fuentes, C.; et al. Plasma potential evolution study by HIBP diagnostic during NBI experiments in the TJ-II stellarator. Fusion Sci. Technol. 2007, 51, 31–37. [Google Scholar] [CrossRef]
  85. Melnikov, A.; Hidalgo, C.; Eliseev, L.; Ascasibar, E.; Chmyga, A.; Dyabilin, K.; Krasilnikov, I.; Krupin, V.; Krupnik, L.; Khrebtov, S.; et al. Plasma potential and turbulence dynamics in toroidal devices (survey of T-10 and TJ-II experiments). Nucl. Fusion 2011, 51. [Google Scholar] [CrossRef]
  86. Krupnik, L.I.; Melnikov, A.V.; Hidalgo, C.; Eliseev, L.G.; Chmyga, A.A.; Komarov, A.D.; Kozachok, A.S.; Perfilov, S.V.; Lysenko, L.G.; Zhezhera, A.I. Recent measurements of the electric potential profile and fluctuations in ECRH and NBI in TJ-II stellarator. Probl. At. Sci. Technol. Ser. Plasma Phys. 2009, 1, 31–33. [Google Scholar]
  87. Krupnik, L.; Alonso, A.; Ascasibar, E.; Estrada, T.; Hidalgo, C.; Van Milligen, B.; Ochando, M.; Pedrosa, M.; De Pablos, J.L.; Tribaldos, V.; et al. Radial electric fields and confinement in the TJ-II stellarator. Czechoslov. J. Phys. 2005, 55, 317–339. [Google Scholar] [CrossRef]
  88. Melnikov, A.V.; Eliseev, L.G.; Grashin, S.A.; Gudozhnik, A.V.; Lysenko, S.E.; Mavrin, V.A.; Perfilov, S.V.; Vershkov, V.A.; Krupnik, L.I.; Chmyga, A.A.; et al. Study of the plasma potential evolution during ECRH in the T-10 tokamak and TJ-II stellarator. Czechoslov. J. Phys. 2005, 55, 1569–1578. [Google Scholar] [CrossRef]
  89. Sharma, R.; Khabanov, P.O.; Melnikov, A.V.; Kharchev, N.K.; Sánchez, E.; Chmyga, A.A.; Deshko, G.N.; Eliseev, L.G.; Hidalgo, C.; Khrebtov, S.M.; et al. Poloidal 2D scans to investigate potential and density profiles in the TJ-II stellarator using heavy ion beam probe. In Proceedings of the 45th EPS Conference on Plasma Physics, Prague, Czech Republic, 2–6 July 2018; ECA. Volume 42A, p. P5.1061. Available online: http://ocs.ciemat.es/EPS2018PAP/pdf/P5.1061.pdf (accessed on 1 January 2021).
  90. Sharma, R.; Khabanov, P.O.; Melnikov, A.V.; Hidalgo, C.; Cappa, A.; Chmyga, A.; Eliseev, L.G.; Estrada, T.; Kharchev, N.K.; Kozachek, A.S.; et al. Measurements of 2D poloidal plasma profiles and fluctuations in ECRH plasmas using the heavy ion beam probe system in the TJ-II stellarator. Phys. Plasmas 2020, 27, 062502. [Google Scholar] [CrossRef]
  91. Eliseev, L.G.; Melnikov, A.V.; Drabinskiy, M.A.; Khabanov, P.O.; Kharchev, N.K.; Lysenko, S.E.; Hidalgo, C.; Barcala, J.M.; Cappa, A.; Liniers, M.; et al. 2D distributions of potential and density mean-values and oscillations in the ECRH and NBI plasmas at the TJ-II stellarator. In Proceedings of the 28th IAEA Fusion Energy Conference (FEC 2020), 10–15 May 2021; EX/P6-17. Available online: http://iaea.org/events/fec-2020 (accessed on 1 January 2021).
  92. Garbet, X.; Idomura, Y.; Villard, L.; Watanabe, T.H. Gyrokinetic simulations of turbulent transport. Nucl. Fusion 2010, 50, 043002. [Google Scholar] [CrossRef]
  93. Melnikov, A.; Krupnik, L.; Barcala, J.M.; Bravo, A.; Chmyga, A.; Deshko, G.; Drabinskij, M.; Eliseev, L.; Hidalgo, C.; Khabanov, P.O.; et al. Heavy ion beam probing—a tool to study geodesic acoustic modes and Alfven eigenmodes in the T-10 tokamak and TJ-II stellarator. Probl. Atom. Sci. Technol. Ser. Plasma Physics 2017, 107, 237–240. [Google Scholar]
  94. Chmyga, O.O.; Ascasibar, E.; Barcala, J.; Drabinskiy, M.A.; Eliseev, L.G.; Hidalgo, C.; Khabanov, P.O.; Khrebtov, S.M.; Komarov, O.D.; Kozachok, O.S.; et al. A dual heavy ion beam probe diagnostic on the TJ-II stellarator. Probl. Atom. Sci. Technol. Ser. Plasma Phys. 2019, 1, 248–251. [Google Scholar]
  95. Fujisawa, A.; Itoh, K.; Iguchi, H.; Matsuoka, K.; Okamura, S.; Shimizu, A.; Minami, T.; Yoshimura, Y.; Nagaoka, K.; Takahashi, C.; et al. Identification of zonal flows in a toroidal plasma. Phys. Rev. Lett. 2004, 93, 165002. [Google Scholar] [CrossRef]
  96. Fujisawa, A.; Ido, T.; Shimizu, A.; Okamura, S.; Matsuoka, K.; Iguchi, H.; Hamada, Y.; Nakano, H.; Ohshima, S.; Itoh, K.; et al. Experimental progress on zonal flow physics in toroidal plasmas. Nucl. Fusion 2007, 47, S718–S726. [Google Scholar] [CrossRef]
  97. Melnikov, A.V.; Krupnik, L.I.; Ascasíbar, E.; Cappa, A.; Chmyga, A.A.; Deshko, G.N.; Drabinskij, M.A.; Eliseev, L.G.; Hidalgo, C.; Khabanov, P.O.; et al. ECRH effect on the electric potential and turbulence in the TJ-II stellarator and T-10 tokamak plasmas. Plasma Phys. Control. Fusion 2018, 60, 084008. [Google Scholar] [CrossRef]
  98. Zhezhera, A.I.; Chmyga, A.; Deshko, G.N.; Eliseev, L.G.; Hidalgo, C.; Khrebtov, S.M.; Komarov, A.D.; Kozachek, A.S.; Krupnik, L.I.; Melnikov, A.V.; et al. New capabilities of plasma potential and density measurements using a dual heavy ion beam probing (HIBP) diagnostic in the TJ-II stellarator. Probl. Atom. Sci. Techn. Ser. Plasma Phys. 2017, 107, 261–264. [Google Scholar]
  99. Fujusawa, A. A review of zonal flow experiments. Nucl. Fusion 2009, 49, 013001. [Google Scholar] [CrossRef]
  100. Hidalgo, C.; Chmyga, A.; Chercoles, J.; Deshko, G.N.; Eliseev, L.G.; Khabanov, P.O.; Khrebtov, S.M.; Komarov, A.G.; Kozachek, A.S.; Krupnik, L.I.; et al. On the influence of ECRH on neoclassical and anomalous mechanisms using a dual heavy ion beam probe diagnostic in the TJ-II stellarator. In Proceedings of the 26th IAEA Fusion Energy Conference, Kyoto, Japan, 17–22 October 2016; Rep. EXC/P7-44. Available online: https://nucleus.iaea.org/sites/fusionportal/Shared%20Documents/FEC%202016/fec2016-preprints/preprint0209.pdf (accessed on 1 January 2021).
  101. Nagaoka, K.; Ido, T.; Ascasíbar, E.; Estrada, T.; Yamamoto, S.; Melnikov, A.V.; Cappa, A.; Hidalgo, C.; Pedrosa, M.A.; van Milligen, B.P.; et al. Mitigation of NBI driven Alfven eigenmodes by electron cyclotron heating in the TJ-II stellarator. Nucl. Fusion 2013, 53, 072004. [Google Scholar] [CrossRef]
  102. Jiménez-Gómez, R.; Könies, A.; Ascasíbar, E.; Castejón, F.; Estrada, T.; Eliseev, L.G.; Melnikov, A.V.; Jiménez, J.A.; Petty, D.G.; Jiménez-Rey, D.; et al. Alfvén eigenmodes measured in the TJ-II stellarator. Nucl. Fusion 2011, 51, 033001. [Google Scholar] [CrossRef]
  103. Melnikov, A.; Eliseev, L.; Castejón, F.; Hidalgo, C.; Khabanov, P.; Kozachek, A.; Krupnik, L.; Liniers, M.; Lysenko, S.; De Pablos, J.L.; et al. Study of NBI-driven chirping mode properties and radial location by the heavy ion beam probe in the TJ-II stellarator. Nucl. Fusion 2016, 56, 112019. [Google Scholar] [CrossRef]
  104. Garcia-Munoz, M.; Sharapov, S.E.; Van Zeeland, M.A.; Ascasibar, E.; Cappa, A.; Chen, L.; Ferreira, J.; Galdon-Quiroga, J.; Geiger, B.; Gonzalez-Martin, J.; et al. Active control of Alfvén eigenmodes in magnetically confined toroidal plasmas. Plasma Phys. Control. Fusion 2019, 61, 054007. [Google Scholar] [CrossRef]
  105. Yamamoto, S.; Nagasaki, K.; Nagaoka, K.; Varela, J.; Cappa, A.; Ascasibar, E.; Castejón, F.; Fontdecaba, J.M.; Regaña, J.M.G.; González-Jerez, Á.; et al. Effect of ECH/ECCD on energetic-particle-driven MHD modes in helical plasmas. Nucl. Fusion 2020, 60, 066018. [Google Scholar] [CrossRef]
  106. Alonso, J.A.; Sánchez, E.; Calvo, I.; Velasco, J.L.; McCarthy, K.J.; Chmyga, A.; Eliseev, L.G.; Estrada, T.; Kleiber, R.; Krupnik, L.I.; et al. Observation of oscillatory radial electric field relaxation in a helical plasma. Phys. Rev. Lett. 2017, 118, 185002. [Google Scholar] [CrossRef] [Green Version]
  107. McCarthy, K.J.; Panadero, N.; Combs, S.K.; Tamura, N.; Ascasíbar, E.; Calvo, M.; Chmyga, A.; Estrada, T.; Fontdecaba, J.M.; García, R.; et al. The impact of fast electrons on pellet injection in the stellarator TJ-II. Plasma Phys. Control. Fusion 2018, 61, 014013. [Google Scholar] [CrossRef]
  108. Chmyga, A.A.; Deshko, G.N.; Dreval, N.B.; Melnikov, A.V.; Perfilov, S.V.; Galvao, R.M.; Kuznetsov, Y.K. Heavy ion beam probe design study for TCABR. Probl. At. Sci. Technol. Ser. Plasma Phys. 2003, 1, 160–162. [Google Scholar]
  109. Kuznetsov, Y.; Nascimento, I.; Silva, C.G.; Figueiredo, H.; Guimarães-Filho, Z.; Caldas, I.; Galvão, R.; Severo, J.; Toufen, D.; Ruchko, L.; et al. Long-distance correlations in TCABR biasing experiments. Nucl. Fusion 2012, 52. [Google Scholar] [CrossRef]
  110. Melnikov, A.V.; Perfilov, S.V. Heavy Ion Beam Probe (HIBP) diagnostics design study for GLOBUS-M tokamak. In Proceedings of the 26th EPS Conference on Plasma Physics and Controlled Fusion, Maastricht, The Netherlands, 14–18 June 1999; ECA, 23J. p. P4.116. [Google Scholar]
  111. Melnikov, A.V.; Markovic, T.; Eliseev, L.; Adamek, J.; Aftanas, M.; Bilkova, P.; Boehm, P.; Gryaznevich, M.; Imrisek, M.; Lysenko, S.E.; et al. Quasicoherent modes on the COMPASS tokamak. Plasma Phys. Control. Fusion 2015, 57. [Google Scholar] [CrossRef]
  112. Seidl, J.; Krbec, J.; Hron, M.; Adamek, J.; Hidalgo, C.; Markovic, T.; Melnikov, A.; Stockel, J.; Weinzettl, V.; Aftanas, M.; et al. Electromagnetic characteristics of geodesic acoustic mode in the COMPASS tokamak. Nucl. Fusion 2017, 57, 126048. [Google Scholar] [CrossRef] [Green Version]
  113. Siegrist, M.; Tonetti, G.; Weisen, H.; Melnikov, A.; Perfilov, S.V.; Krupnik, L.; Karpushov, I. The Heavy Ion Beam Diagnostic Project for the TCV Tokamak. In Proceedings of the 21st IEEE/NPS Symposium on Fusion Engineering SOFE 05, Knoxville, TN, USA, 26–29 September 2005; pp. 1–4. [Google Scholar] [CrossRef]
  114. Melnikov, A.V.; Perfilov, S.V.; Sharapov, S.E. Heavy ion beam probe for studying potential and turbulence on MAST. In Proceedings of the 37th EPS Conference on Plasma Physics, Dublin, Ireland, 21–25 June 2010; ECA. Volume 34A, p. P5.120. Available online: http://ocs.ciemat.es/EPS2010PAP/html/ (accessed on 1 January 2021).
  115. Melnikov, A.V.; Perfilov, S.V. Effect of the sample volume splitting in heavy ion-beam probing (HIBP) for spherical tokamaks (by the example of HIBP for mega ampere spherical tokamak). Rev. Sci. Instrum. 1999, 70, 1402–1406. [Google Scholar] [CrossRef]
  116. Melnikov, A.V.; Sushkov, A.V.; Belov, A.M.; Dnestrovskij, Y.N.; Eliseev, L.G.; Ivanov, D.P.; Kirneva, N.A.; Korobov, K.V.; Krupin, V.A.; Lysenko, S.E.; et al. Physical program and conceptual design of the diagnostics of the T-15 upgrade tokamak. Fusion Eng. Design 2015, 96, 306–310. [Google Scholar] [CrossRef]
  117. Ilin, A.M.; Khabanov, P.O.; Melnikov, A.V. Probing ion trajectory simulations for the HIBP diagnostics at the T-15MD tokamak. J. Phys. Conf. Ser. 2019, 1383, 012006. [Google Scholar] [CrossRef]
  118. Drabinskiy, M.; Melnikov, A.; Khabanov, P.; Eliseev, L.; Kharchev, N.; Ilin, A.; Sarancha, G.; Vadimov, N. Conceptual design of the heavy ion beam probe diagnostic for the T-15MD tokamak. J. Instrum. 2019, 14, C11027. [Google Scholar] [CrossRef]
  119. Melnikov, A.; Eliseev, L.; Lysenko, S.; Perfilov, S.; Shelukhin, D.; Vershkov, V.; Zenin, V.; Krupnik, L.; Kozachek, A.; Kharchev, N.; et al. Long-distance correlations of geodesic acoustic modes in T-10. Probl. At. Sci. Technol. Ser. Thermonucl. Fusion 2015, 38. [Google Scholar] [CrossRef]
  120. Melnikov, A.V.; Vershkov, V.A.; Eliseev, L.; Grashin, S.A.; Gudozhnik, A.V.; Krupnik, L.I.; Lysenko, S.E.; Mavrin, V.A.; Perfilov, S.V.; Shelukhin, D.A.; et al. Investigation of geodesic acoustic mode oscillations in the T-10 tokamak. Plasma Phys. Control. Fusion 2006, 48, S87–S110. [Google Scholar] [CrossRef]
  121. Melnikov, A.V.; Eliseev, L.G.; Ascasíbar, E.; Cappa, A.; Castejón, F.; Hidalgo, C.; Ido, T.; Jiménez, J.A.; Kozachek, A.S.; Krupnik, L.I.; et al. Effect of magnetic configuration on nonlinear evolution of NBI-driven Alfvén modes in TJ-II. Nucl. Fusion 2016, 56, 076001. [Google Scholar] [CrossRef]
  122. Cappa, Á.; Varela, J.; Lopez-Bruna, D.L.; Ascasíbar, E.; Liniers, M.; Eliseev, L.; Fontdecaba, J.; García-Regaña, J.; González-Jerez, A.; Kharchev, N.; et al. Stability analysis of TJ-II stellarator NBI driven Alfvén eigenmodes in ECRH and ECCD experiments. Nucl. Fusion 2021, 61, 066019. [Google Scholar] [CrossRef]
  123. Ongena, J.; Messiaen, A.; Melnikov, A.; Ragona, R.; Kazakov, Y.O.; Van Eester, D.; Dnestrovskii, Y.N.; Khvostenko, P.; Roy, I.; Romannikov, A. Conceptual study of an ICRH system for T-15MD using traveling wave antenna (TWA) sections. Fusion Eng. Des. 2019, 146, 787–791. [Google Scholar] [CrossRef]
  124. Melnikov, A.V.; Ongena, J.; Messiaen, A.M.; Ragona, R.; Sushkov, A.V.; Kazakov, Y.O.; Van Eester, D.; Dnestrovskii, Y.N.; Khvostenko, P.P.; Roy, I.N. Conceptual study of an ICRH traveling wave antenna (TWA) for T-15MD at 60 MHz. AIP Conf. Proc. 2020, 2254, 070007. [Google Scholar] [CrossRef]
  125. Perfilov, S.; Melnikov, A.; Krupnik, L.; Hartfuss, H.J. Applicability of Heavy Ion Beam Probing for Stellarator W7-X. AIP Conf. Proc. 2006, 812, 199–202. [Google Scholar] [CrossRef] [Green Version]
  126. Perfilov, S.; Melnikov, A.; Krupnik, L.; Hartfuss, H.J. Applicability of Heavy Ion Beam Probing for Stellarator W7-X. Fusion Sci. Technol. 2007, 51, 38–45. [Google Scholar] [CrossRef] [Green Version]
  127. Melnikov, A.V.; Eliseev, L.G. Optimized heavy ion beam probing for International Thermonuclear Experimental Reactor. Rev. Sci. Instrum. 1999, 70, 951–954. [Google Scholar] [CrossRef]
  128. Melnikov, A.V.; Eliseev, L.G.; Razumova, K.A.; Krupnik, L.I.; Chistyakov, V.V.; Myalton, T.B.; Dnestrovskij, Y.N.; HIBP Team. On the sensitivity of the plasma electric potential to the edge transport barrier formation. Jpn. J. Plasma Fusion Res. Series 2000, 4, 203–208. [Google Scholar]
  129. Alikaev, V.V.; Borshchegovskij, A.A.; Vershkov, V.A.; Volkov, V.V.; Gorshkov, A.V.; Gott, Y.V.; Grashin, S.A.; Dremin, M.M.; Eliseev, L.G.; Esipchuk, Y.V.; et al. Investigation of the H-mode during ECRH in the T-10 tokamak. Plasma Phys. Rep. 2000, 26, 917. [Google Scholar] [CrossRef]
  130. Melnikov, A.V.; Eliseev, L.G.; Grashin, S.A.; Drabinskiy, M.A.; Khabanov, P.O.; Kharchev, N.K.; Krupin, V.A.; Lysenko, S.E.; Ryzhakov, D.A.; Soloviev, N.A.; et al. Evolution of the electric potential and turbulence in OH and ECRH low-density plasmas at the T-10 tokamak. In Proceedings of the 28th IAEA Fusion Energy Conference (FEC 2020), 10–15 May 2021; Conference was held online. EX/5-6 (Nucl. Fusion under review). Available online: https://conferences.iaea.org/event/214/contributions/17104/ (accessed on 1 January 2021).
Figure 1. (a) HIBP experimental set-up at T-10 [15]: (1) accelerator, (2) primary beamline, (3) secondary beamline, (4) analyzer, (5) trajectory of primary Tl+ ions (red curve), (6) trajectory of secondary Tl++ ions (yellow curve), (7) ionization point or sample volume (SV). (b) Photo of the HIBP diagnostic hardware. Reproduced courtesy of IAEA. Copyright 2017 IAEA.
Figure 1. (a) HIBP experimental set-up at T-10 [15]: (1) accelerator, (2) primary beamline, (3) secondary beamline, (4) analyzer, (5) trajectory of primary Tl+ ions (red curve), (6) trajectory of secondary Tl++ ions (yellow curve), (7) ionization point or sample volume (SV). (b) Photo of the HIBP diagnostic hardware. Reproduced courtesy of IAEA. Copyright 2017 IAEA.
Symmetry 13 01367 g001
Figure 2. Detector grid for T-10 with the field Bt = 1.55 T: lines of equal beam energy (Eb) are marked in green, lines of equal injection angle α, labeled with scanning voltage Uscan are blue. The asterisks denote the nodes of the detector grid available for observation with all slits of multi-slit detector.
Figure 2. Detector grid for T-10 with the field Bt = 1.55 T: lines of equal beam energy (Eb) are marked in green, lines of equal injection angle α, labeled with scanning voltage Uscan are blue. The asterisks denote the nodes of the detector grid available for observation with all slits of multi-slit detector.
Symmetry 13 01367 g002
Figure 3. Principle of plasma potential measurement.
Figure 3. Principle of plasma potential measurement.
Symmetry 13 01367 g003
Figure 4. Plasma density profiles in TM-4: (a) Ohmic heated (OH) plasma, n ¯ e = 1.5 × 1019 m−3; dash line and × is a solution of Equation (3), without attenuation; solid line and ○ is solution of Equation (4) with a quasi-solution method. (b) OH ( n ¯ e = 1.3 × 1019 m−3) and ECR heated plasmas; ○ is HIBP, + is interferometer; □ are probes and Δ is Thomson scattering data.
Figure 4. Plasma density profiles in TM-4: (a) Ohmic heated (OH) plasma, n ¯ e = 1.5 × 1019 m−3; dash line and × is a solution of Equation (3), without attenuation; solid line and ○ is solution of Equation (4) with a quasi-solution method. (b) OH ( n ¯ e = 1.3 × 1019 m−3) and ECR heated plasmas; ○ is HIBP, + is interferometer; □ are probes and Δ is Thomson scattering data.
Symmetry 13 01367 g004
Figure 5. (a) Toroidal shift of trajectories as function of the path traversed by the probing particle in TUMAN-3M (Bt = 0.5 T; Ipl = 100 kA; R = 54 cm; a = 24 cm), with parabolic profile of current; sα1 and sα2 are distances from the injector to the plasma boundary; sd is the distance to the detector point; curves with different colors shows trajectories, born in different sample volumes, marked with black rectangles. (b) Shift of beam particles in the detector as a function of energy at a total plasma current Ipl = 100 kA with different profiles of current density, marked with numbers. The width of the beam track for each current profile corresponds to error in definition of the beam position relatively the meridional plane of 1 cm.
Figure 5. (a) Toroidal shift of trajectories as function of the path traversed by the probing particle in TUMAN-3M (Bt = 0.5 T; Ipl = 100 kA; R = 54 cm; a = 24 cm), with parabolic profile of current; sα1 and sα2 are distances from the injector to the plasma boundary; sd is the distance to the detector point; curves with different colors shows trajectories, born in different sample volumes, marked with black rectangles. (b) Shift of beam particles in the detector as a function of energy at a total plasma current Ipl = 100 kA with different profiles of current density, marked with numbers. The width of the beam track for each current profile corresponds to error in definition of the beam position relatively the meridional plane of 1 cm.
Symmetry 13 01367 g005
Figure 6. (a) PSD of the potential and plasma density fluctuations measured with a heavy ion beam, and the PSD of magnetic oscillations measured with a magnetic probe in the typical T-10 ohmic discharge [15]. (b) Spectrogram of plasma potential oscillations measured in the T-10 discharge with additional ECR heating from 610 to 800 ms. Reproduced courtesy of IAEA. Copyright 2017 IAEA.
Figure 6. (a) PSD of the potential and plasma density fluctuations measured with a heavy ion beam, and the PSD of magnetic oscillations measured with a magnetic probe in the typical T-10 ohmic discharge [15]. (b) Spectrogram of plasma potential oscillations measured in the T-10 discharge with additional ECR heating from 610 to 800 ms. Reproduced courtesy of IAEA. Copyright 2017 IAEA.
Symmetry 13 01367 g006
Figure 7. Conceptual scheme of a zonal flow, a radially localized, poloidally and azimuthally symmetric (m = 0, n = 0) electrostatic potential structure, along with its resulting radial electric field Er, and turbulence stabilizing poloidal E × B flow.
Figure 7. Conceptual scheme of a zonal flow, a radially localized, poloidally and azimuthally symmetric (m = 0, n = 0) electrostatic potential structure, along with its resulting radial electric field Er, and turbulence stabilizing poloidal E × B flow.
Symmetry 13 01367 g007
Figure 8. Spectrogram of the coherency between the total HIBP current (plasma density) oscillations, measured in the plasma core at r/a = 0.16 and a MP signal [45]. Reproduced courtesy of IAEA. Copyright 2010 IAEA.
Figure 8. Spectrogram of the coherency between the total HIBP current (plasma density) oscillations, measured in the plasma core at r/a = 0.16 and a MP signal [45]. Reproduced courtesy of IAEA. Copyright 2010 IAEA.
Symmetry 13 01367 g008
Figure 9. Spectrogram of the cross-phase between signals It (~ne) and poloidal magnetic field Bp (~ζd) measured in the same SV in discharge with AEs. Insets show the histograms of the cross-phase computed over the areas marked by rectangles; the insets show the histograms for the cross-phases [45]. Reproduced courtesy of IAEA. Copyright 2010 IAEA.
Figure 9. Spectrogram of the cross-phase between signals It (~ne) and poloidal magnetic field Bp (~ζd) measured in the same SV in discharge with AEs. Insets show the histograms of the cross-phase computed over the areas marked by rectangles; the insets show the histograms for the cross-phases [45]. Reproduced courtesy of IAEA. Copyright 2010 IAEA.
Symmetry 13 01367 g009
Figure 10. Auto-bicoherency (a) and biphase (b) for the plasma potential in the ohmic discharge of the T-10 tokamak [46]. Reproduced courtesy of IAEA. Copyright 2017 IAEA.
Figure 10. Auto-bicoherency (a) and biphase (b) for the plasma potential in the ohmic discharge of the T-10 tokamak [46]. Reproduced courtesy of IAEA. Copyright 2017 IAEA.
Symmetry 13 01367 g010
Figure 11. Multi-slit energy analyzer for direct measurements of particle flux: 5-S—five entrance slits; D—detector; G—grid; GP—ground plate; HVP—high-voltage plate; W—adjustment window.
Figure 11. Multi-slit energy analyzer for direct measurements of particle flux: 5-S—five entrance slits; D—detector; G—grid; GP—ground plate; HVP—high-voltage plate; W—adjustment window.
Symmetry 13 01367 g011
Figure 12. Detector line for the two-slit energy analyzer [15]. The upper slit is marked in blue, the lower one—in red. Reproduced courtesy of IAEA. Copyright 2017 IAEA.
Figure 12. Detector line for the two-slit energy analyzer [15]. The upper slit is marked in blue, the lower one—in red. Reproduced courtesy of IAEA. Copyright 2017 IAEA.
Symmetry 13 01367 g012
Figure 13. Scheme of the poloidal electric field measurements in the plasma using the double-slit analyzer. The potential φ1 and φ2 are measured by the first and second slit correspondingly; ∆x is distance between the points of measurements; Ep is poloidal electric field; Bt is toroidal magnetic field, Vr is the velocity of radial drift in crossed Ep × Bt fields.
Figure 13. Scheme of the poloidal electric field measurements in the plasma using the double-slit analyzer. The potential φ1 and φ2 are measured by the first and second slit correspondingly; ∆x is distance between the points of measurements; Ep is poloidal electric field; Bt is toroidal magnetic field, Vr is the velocity of radial drift in crossed Ep × Bt fields.
Symmetry 13 01367 g013
Figure 14. The time evolution of the NBI-sustained discharge with AEs in TJ-II. (a) Frequency resolved particle flux ΓE×B (in a. u.) measured by HIBP at ρ = −0.5, kθ < 2 cm−1; red color means positive outward flux, and blue means negative inward flux. (b) Cross-phase θEp-ne between Ep and ne oscillations. Only the points with high Cohne- Ep > 0.3 are shown. The color bar is in radians. Three chosen branches of the AE are marked by color ovals. (c) The histograms of the cross-phase for each branch is marked with corresponding colors, indicating the flux direction; left box (red): θEp-ne = −3/4π, corresponds to the outward flux; central box (green): θEp-ne = −π/2, corresponds to zero flux; right box (blue): θEp-ne ~0 corresponds to the inward flux. (d) PSD of the turbulent particle flux, taken at some typical time instant, averaged over 1 ms. Three frequency peaks related to the AE branches identified above are marked with corresponding colors. Adapted from [51] Copyright 2012 IAEA.
Figure 14. The time evolution of the NBI-sustained discharge with AEs in TJ-II. (a) Frequency resolved particle flux ΓE×B (in a. u.) measured by HIBP at ρ = −0.5, kθ < 2 cm−1; red color means positive outward flux, and blue means negative inward flux. (b) Cross-phase θEp-ne between Ep and ne oscillations. Only the points with high Cohne- Ep > 0.3 are shown. The color bar is in radians. Three chosen branches of the AE are marked by color ovals. (c) The histograms of the cross-phase for each branch is marked with corresponding colors, indicating the flux direction; left box (red): θEp-ne = −3/4π, corresponds to the outward flux; central box (green): θEp-ne = −π/2, corresponds to zero flux; right box (blue): θEp-ne ~0 corresponds to the inward flux. (d) PSD of the turbulent particle flux, taken at some typical time instant, averaged over 1 ms. Three frequency peaks related to the AE branches identified above are marked with corresponding colors. Adapted from [51] Copyright 2012 IAEA.
Symmetry 13 01367 g014
Figure 15. Scheme for measuring the electric field in plasma using the 5-slit analyzer. The potential φk is measured in the point Sk by slit with number k = 1–5; Ep is the poloidal component of the electric field required to calculate Vr; the radial Ep × Bt drift velocity in the crossed fields, and Er is its radial component.
Figure 15. Scheme for measuring the electric field in plasma using the 5-slit analyzer. The potential φk is measured in the point Sk by slit with number k = 1–5; Ep is the poloidal component of the electric field required to calculate Vr; the radial Ep × Bt drift velocity in the crossed fields, and Er is its radial component.
Symmetry 13 01367 g015
Figure 16. Time evolution of the T-10 discharge with additional ECR heating. (a) Flux spectral function ΓE×B (rel. un.) measured at ρ = 0.78, kθ < 2; the red color corresponds to the outward flux, positive, blue—to the inward flux, negative. The time-frequency regions of the quasi-coherent mode (QC) and geodesic acoustic mode (GAM) are marked; (b) the total particle fluxes averaged over GAM ( Γ E × B G A M ) and QC ( Γ E × B Q C ) frequency ranges, and the line-average density n ¯ e [54]. Reproduced with permission IoP.
Figure 16. Time evolution of the T-10 discharge with additional ECR heating. (a) Flux spectral function ΓE×B (rel. un.) measured at ρ = 0.78, kθ < 2; the red color corresponds to the outward flux, positive, blue—to the inward flux, negative. The time-frequency regions of the quasi-coherent mode (QC) and geodesic acoustic mode (GAM) are marked; (b) the total particle fluxes averaged over GAM ( Γ E × B G A M ) and QC ( Γ E × B Q C ) frequency ranges, and the line-average density n ¯ e [54]. Reproduced with permission IoP.
Symmetry 13 01367 g016
Figure 17. Time evolution of Ohmic T-10 discharge with slow density ramp-up. (a) Power spectral density (PSD) of potential oscillations measured by the central slit φ3, SV is located at r/a = 0.57; (b) PSD for difference of potentials between the central and edge slits φ1φ3 = Ep Δx; (c) the frequency resolved spectral function of flux ΓE×B. The frequency domain of GAM is marked by rectangle [54]. Reproduced with permission IoP. Copyright IoP.
Figure 17. Time evolution of Ohmic T-10 discharge with slow density ramp-up. (a) Power spectral density (PSD) of potential oscillations measured by the central slit φ3, SV is located at r/a = 0.57; (b) PSD for difference of potentials between the central and edge slits φ1φ3 = Ep Δx; (c) the frequency resolved spectral function of flux ΓE×B. The frequency domain of GAM is marked by rectangle [54]. Reproduced with permission IoP. Copyright IoP.
Symmetry 13 01367 g017
Figure 18. Two-point technique of correlation measurements of the plasma electric potential using the 5-slit analyzer.
Figure 18. Two-point technique of correlation measurements of the plasma electric potential using the 5-slit analyzer.
Symmetry 13 01367 g018
Figure 19. Correlational measurements of electric potential with 5-slits analyzer. Traces of plasma current Ip and line-averaged density n ¯ e in discharge with Bt = 2.3 T (a); power spectral densities of potential oscillations measured in two most distant SVs by upper slit φ1 (b); and by lower slit φ5 (c); their coherency (d); and cross-phase (e) at r = 21 cm. GAM dominates in all PSDs. Adapted from [58] Copyright IAEA 2015.
Figure 19. Correlational measurements of electric potential with 5-slits analyzer. Traces of plasma current Ip and line-averaged density n ¯ e in discharge with Bt = 2.3 T (a); power spectral densities of potential oscillations measured in two most distant SVs by upper slit φ1 (b); and by lower slit φ5 (c); their coherency (d); and cross-phase (e) at r = 21 cm. GAM dominates in all PSDs. Adapted from [58] Copyright IAEA 2015.
Symmetry 13 01367 g019
Figure 20. The phase shift between the potential oscillations and the poloidal mode number m measured by intermediate slits of the analyzer in Ohmic discharge with Bt = 2.3 T, ne = (0.9–2.4) × 1019 m−3, Ip = 220 kA.
Figure 20. The phase shift between the potential oscillations and the poloidal mode number m measured by intermediate slits of the analyzer in Ohmic discharge with Bt = 2.3 T, ne = (0.9–2.4) × 1019 m−3, Ip = 220 kA.
Symmetry 13 01367 g020
Figure 21. The phase shift between the potential oscillations measured by most distant slits 1 and 5 of the analyzer for the plasma scenario shown in Figure 19. Adapted from [58] Copyright IAEA 2015.
Figure 21. The phase shift between the potential oscillations measured by most distant slits 1 and 5 of the analyzer for the plasma scenario shown in Figure 19. Adapted from [58] Copyright IAEA 2015.
Symmetry 13 01367 g021
Figure 22. Scheme of the TM-4 HIBP complex: 1—source power unit; 2—high-voltage source for accelerator; 3—ion source; 4—accelerating tube; 5—high-voltage divider; 6—ion beamline; 7—vacuum valves; 8—primary beam pick-up; 9—plasma; 10—tokamak chamber; 11—grid analyzer for secondary ions; 12—secondary-electron multiplier. Adapted from [13] Copyright IAEA 1985.
Figure 22. Scheme of the TM-4 HIBP complex: 1—source power unit; 2—high-voltage source for accelerator; 3—ion source; 4—accelerating tube; 5—high-voltage divider; 6—ion beamline; 7—vacuum valves; 8—primary beam pick-up; 9—plasma; 10—tokamak chamber; 11—grid analyzer for secondary ions; 12—secondary-electron multiplier. Adapted from [13] Copyright IAEA 1985.
Symmetry 13 01367 g022
Figure 23. Initial scheme of the heavy ion beam probe on the T-10 tokamak.
Figure 23. Initial scheme of the heavy ion beam probe on the T-10 tokamak.
Symmetry 13 01367 g023
Figure 24. Scheme of the heavy ion beam probe on the TJ-I tokamak. (a) Machine and diagnostics; (b) detector grid with lines of equal injection angle from −12° to +12° (dashed), and the equal beam energy Eb from 6 to 60 keV (solid lines).
Figure 24. Scheme of the heavy ion beam probe on the TJ-I tokamak. (a) Machine and diagnostics; (b) detector grid with lines of equal injection angle from −12° to +12° (dashed), and the equal beam energy Eb from 6 to 60 keV (solid lines).
Symmetry 13 01367 g024
Figure 25. The profiles of the plasma potential (a) and density (b) for three spatial scans in TJ-I.
Figure 25. The profiles of the plasma potential (a) and density (b) for three spatial scans in TJ-I.
Symmetry 13 01367 g025
Figure 26. Scheme of the HIBP on the TUMAN-3M tokamak, the major axis of the torus is on the left. Thin curves mark the trajectories of probing particles, I—injector, D—detector. Asterisks indicate points of the detector line available for observation. The calculation was carried out for a beam of Cs+ ions with the beam energies Eb = 4–13 keV, the point with Eb = 10 keV is located in the plasma center.
Figure 26. Scheme of the HIBP on the TUMAN-3M tokamak, the major axis of the torus is on the left. Thin curves mark the trajectories of probing particles, I—injector, D—detector. Asterisks indicate points of the detector line available for observation. The calculation was carried out for a beam of Cs+ ions with the beam energies Eb = 4–13 keV, the point with Eb = 10 keV is located in the plasma center.
Symmetry 13 01367 g026
Figure 27. Numerical study of the sensitivity of the toroidal shift of the trajectories to the plasma current and the initial toroidal angle β of the TUMAN-3M. Projections of the trajectories onto the equatorial plane of the torus (x, z) for different values of the plasma current I and the initial toroidal angle β (a–d); XI and XD injector and detector coordinates; asterisks mark the points of the detector line available for observation for energies Eb = 4, 5, 6, 8, and 10 keV.
Figure 27. Numerical study of the sensitivity of the toroidal shift of the trajectories to the plasma current and the initial toroidal angle β of the TUMAN-3M. Projections of the trajectories onto the equatorial plane of the torus (x, z) for different values of the plasma current I and the initial toroidal angle β (a–d); XI and XD injector and detector coordinates; asterisks mark the points of the detector line available for observation for energies Eb = 4, 5, 6, 8, and 10 keV.
Symmetry 13 01367 g027
Figure 28. Scheme of the HIBP on the WEGA stellarator. Bt = 0.5 T, ι/2π = 0.18.
Figure 28. Scheme of the HIBP on the WEGA stellarator. Bt = 0.5 T, ι/2π = 0.18.
Symmetry 13 01367 g028
Figure 29. Detector grid for Na+ ions in the WEGA stellarator. The beam energy Eb and the injection angle (expressed in scanning voltage Uscan) are indicated.
Figure 29. Detector grid for Na+ ions in the WEGA stellarator. The beam energy Eb and the injection angle (expressed in scanning voltage Uscan) are indicated.
Symmetry 13 01367 g029
Figure 30. Scheme of the HIBP on the URAGAN-2M stellarator. The optimal detector line is marked by blue asterisks, the green solid line is the trajectory; I is injector, D is detector, and P is plasma; (a) vertical plane, (b) equatorial plane.
Figure 30. Scheme of the HIBP on the URAGAN-2M stellarator. The optimal detector line is marked by blue asterisks, the green solid line is the trajectory; I is injector, D is detector, and P is plasma; (a) vertical plane, (b) equatorial plane.
Symmetry 13 01367 g030
Figure 31. Photo of the HIBP in TJ-II: (a) ion energy analyzer; (b) ion beam injector.
Figure 31. Photo of the HIBP in TJ-II: (a) ion energy analyzer; (b) ion beam injector.
Symmetry 13 01367 g031
Figure 32. Scheme of the HIBP in TJ-II. Inset shows plasma cross-section and detector line for Eb = 127 keV, scanning voltage limits are indicated: −6 kV at HFS = −1) and +4 kV at LFS (ρ = +1).
Figure 32. Scheme of the HIBP in TJ-II. Inset shows plasma cross-section and detector line for Eb = 127 keV, scanning voltage limits are indicated: −6 kV at HFS = −1) and +4 kV at LFS (ρ = +1).
Symmetry 13 01367 g032
Figure 33. Plasma potential profile evolution in TJ-II. The first proof of the LFS–HFS symmetry [85]. Reproduced courtesy of IAEA. Copyright 2011 IAEA.
Figure 33. Plasma potential profile evolution in TJ-II. The first proof of the LFS–HFS symmetry [85]. Reproduced courtesy of IAEA. Copyright 2011 IAEA.
Symmetry 13 01367 g033
Figure 34. 3D structure of HAE mode m/n = 4/7; (a) top view; (b) poloidal cut [51]. Reproduced courtesy of IAEA. Copyright 2012 IAEA.
Figure 34. 3D structure of HAE mode m/n = 4/7; (a) top view; (b) poloidal cut [51]. Reproduced courtesy of IAEA. Copyright 2012 IAEA.
Symmetry 13 01367 g034
Figure 35. The spatial distribution of the plasma potential (a) and its oscillations (b) over the vertical cross-section in TJ-II. (a) Potential contours cover about one half of the plasma poloidal cross-section: the maximum of 2D potential distribution (~1000 V) is located at the magnetic axis, equipotentials are consistent with vacuum magnetic flux surfaces, and 2D potential distribution is symmetric (LFS–HFS and up–down). (b) The maxima of 2D RMS distribution are located at the magnetic axis (~30 V) and at the edge (~35–40 V), equipotentials are basically consistent with vacuum magnetic flux surfaces, 2D potential distribution is not fully symmetric: at the mid-radius ~15 V at LFS vs. ~20 V at HFS, i.e., 1.3 times higher. [91].
Figure 35. The spatial distribution of the plasma potential (a) and its oscillations (b) over the vertical cross-section in TJ-II. (a) Potential contours cover about one half of the plasma poloidal cross-section: the maximum of 2D potential distribution (~1000 V) is located at the magnetic axis, equipotentials are consistent with vacuum magnetic flux surfaces, and 2D potential distribution is symmetric (LFS–HFS and up–down). (b) The maxima of 2D RMS distribution are located at the magnetic axis (~30 V) and at the edge (~35–40 V), equipotentials are basically consistent with vacuum magnetic flux surfaces, 2D potential distribution is not fully symmetric: at the mid-radius ~15 V at LFS vs. ~20 V at HFS, i.e., 1.3 times higher. [91].
Symmetry 13 01367 g035
Figure 36. Dual HIBP on TJ-II—the tool to study plasma symmetric structures. (a) Top view of the TJ-II diagnostics layout. (b) Side view of HIBP-2.
Figure 36. Dual HIBP on TJ-II—the tool to study plasma symmetric structures. (a) Top view of the TJ-II diagnostics layout. (b) Side view of HIBP-2.
Symmetry 13 01367 g036
Figure 37. Photo of the dual HIBP on TJ-II.
Figure 37. Photo of the dual HIBP on TJ-II.
Symmetry 13 01367 g037
Figure 38. Long-range toroidal correlations in plasma potential oscillations in TJ-II [97]. Plasma potential long-range-correlations radial distribution in the ECRH discharge #39894, PECRH = 2 × 240 kW. (a) Potential coherency spectrogram coh(φ1(ρSV1), φ2(ρSV2)) vs. radius ρSV1 by HIBP-1; ρSV1 scans from −1 at HFS to +1 at LFS, while HIBP-2 measures at fixed position ρSV2 = −0.63; (b) spectrogram of the cross-phase between φ1(ρSV1) and φ2(ρSV2), calculated for the coherency, exceeding noise level. (c) Radial distribution of potential coherency averaged by 0 < f < 20 kHz. Reproduced courtesy of IoP. Copyright 2018 IoP.
Figure 38. Long-range toroidal correlations in plasma potential oscillations in TJ-II [97]. Plasma potential long-range-correlations radial distribution in the ECRH discharge #39894, PECRH = 2 × 240 kW. (a) Potential coherency spectrogram coh(φ1(ρSV1), φ2(ρSV2)) vs. radius ρSV1 by HIBP-1; ρSV1 scans from −1 at HFS to +1 at LFS, while HIBP-2 measures at fixed position ρSV2 = −0.63; (b) spectrogram of the cross-phase between φ1(ρSV1) and φ2(ρSV2), calculated for the coherency, exceeding noise level. (c) Radial distribution of potential coherency averaged by 0 < f < 20 kHz. Reproduced courtesy of IoP. Copyright 2018 IoP.
Symmetry 13 01367 g038
Figure 39. (a) Scheme of the HIBP at TCABR tokamak; (b) Detector grid of the TCABR tokamak for Tl ions with beam energy Eb = 55–106 keV and the injection angle 65–76°. The grid calculated with a step of 5 keV and 1°.
Figure 39. (a) Scheme of the HIBP at TCABR tokamak; (b) Detector grid of the TCABR tokamak for Tl ions with beam energy Eb = 55–106 keV and the injection angle 65–76°. The grid calculated with a step of 5 keV and 1°.
Symmetry 13 01367 g039
Figure 40. (a) Scheme of the HIBP at Globus-M2 tokamak. HIBP primary (black) and secondary (red) trajectories for a detector line of the probing scheme using 78° input port, Eb = 40 keV, with the secondary beamline. Primary beamline angles are α = 45°, β = −5°. Secondary beamline angles are α = 0°, β = 20°. Top: side view, bottom: bottom view. Green star denotes the detection point, blue diamond denotes plasma center. (b) Detector grid for the probing scheme shown in (a).
Figure 40. (a) Scheme of the HIBP at Globus-M2 tokamak. HIBP primary (black) and secondary (red) trajectories for a detector line of the probing scheme using 78° input port, Eb = 40 keV, with the secondary beamline. Primary beamline angles are α = 45°, β = −5°. Secondary beamline angles are α = 0°, β = 20°. Top: side view, bottom: bottom view. Green star denotes the detection point, blue diamond denotes plasma center. (b) Detector grid for the probing scheme shown in (a).
Symmetry 13 01367 g040
Figure 41. Design of the heavy ion beam probing scheme for the COMPASS tokamak.
Figure 41. Design of the heavy ion beam probing scheme for the COMPASS tokamak.
Symmetry 13 01367 g041
Figure 42. Detector grid for the Tl ions at Bt = 2.1 T and beam energies Eb = 130–250 keV. The grid is calculated with a step of 10 keV. The line Eb = 210 keV allows us to connect the center and edge of the plasma, 0 < r/a < 1.
Figure 42. Detector grid for the Tl ions at Bt = 2.1 T and beam energies Eb = 130–250 keV. The grid is calculated with a step of 10 keV. The line Eb = 210 keV allows us to connect the center and edge of the plasma, 0 < r/a < 1.
Symmetry 13 01367 g042
Figure 43. (a) Schematic of the HIBP complex on the TCV tokamak: Bt = 1.07 T, Ipl = 135 kA; (b) detector grid for Tl ions, Eb = 120–250 keV and the injection angle is 65–80°. The grid is drawn with a step of 10 keV and 1°.
Figure 43. (a) Schematic of the HIBP complex on the TCV tokamak: Bt = 1.07 T, Ipl = 135 kA; (b) detector grid for Tl ions, Eb = 120–250 keV and the injection angle is 65–80°. The grid is drawn with a step of 10 keV and 1°.
Symmetry 13 01367 g043
Figure 44. (a) Schematic of the HIBP complex on the MAST spherical tokamak: Bt = 0.6 T, Ipl = 1 MA; I—injector; A, B and C in inset are various detector positions. Inset demonstrates intersection of secondary trajectories; (b) detector grid for Cs ions. Areas without intersection the secondary trajectories are marked in red, with the crossing—in blue. Lines of equal angle are solid, lines of equal energy are dashed.
Figure 44. (a) Schematic of the HIBP complex on the MAST spherical tokamak: Bt = 0.6 T, Ipl = 1 MA; I—injector; A, B and C in inset are various detector positions. Inset demonstrates intersection of secondary trajectories; (b) detector grid for Cs ions. Areas without intersection the secondary trajectories are marked in red, with the crossing—in blue. Lines of equal angle are solid, lines of equal energy are dashed.
Symmetry 13 01367 g044
Figure 45. (a) Schematic of the HIBP trajectories on the T-15 MD tokamak: Bt = 2 T, Ipl = 2 MA; (b) part of the poloidal section of T-15 MD chamber and the detector grid. The calculations were performed for Tl ion beam in reduced magnetic field Bt = 1 T. Lines of equal energy Eb = 80–220 keV are green.
Figure 45. (a) Schematic of the HIBP trajectories on the T-15 MD tokamak: Bt = 2 T, Ipl = 2 MA; (b) part of the poloidal section of T-15 MD chamber and the detector grid. The calculations were performed for Tl ion beam in reduced magnetic field Bt = 1 T. Lines of equal energy Eb = 80–220 keV are green.
Symmetry 13 01367 g045
Figure 46. Scheme of HIP for the W7-X stellarator. (1)—primary trajectories injected through the port A11 at the bottom of the figure; (2)—fan of secondary trajectories directed to the analyzer through the port N11; (3)—total fan of secondary trajectories.
Figure 46. Scheme of HIP for the W7-X stellarator. (1)—primary trajectories injected through the port A11 at the bottom of the figure; (2)—fan of secondary trajectories directed to the analyzer through the port N11; (3)—total fan of secondary trajectories.
Symmetry 13 01367 g046
Figure 47. Schematic of the HIBP complex on the fusion reactor ITER: Bt = 6 T, Ipl = 15 MA. The inset shows the detector grid and fragments of ion beam lines. The asterisks mark the detector line, which allows one to study the radial profiles of the plasma potential; I—injection point, D—detection point, dimensions in cm.
Figure 47. Schematic of the HIBP complex on the fusion reactor ITER: Bt = 6 T, Ipl = 15 MA. The inset shows the detector grid and fragments of ion beam lines. The asterisks mark the detector line, which allows one to study the radial profiles of the plasma potential; I—injection point, D—detection point, dimensions in cm.
Symmetry 13 01367 g047
Figure 48. (a) The detector grid for ITER. Calculated for Eb = 2.6–3.4 MeV with steps of 0.1 MeV and angle 0.5°. The line for α = 13° allows one to probe a fragment of the potential profile at the plasma periphery, 0.76 < r/a < 0.98, by changing the energy from shot to shot. The line Eb = 3.3 MeV allows one to obtain a fragment of the potential profile 0.76 < r/a < 1 by changing the injection angle β. (b) Toroidal shift in the equatorial projection of the detector line with α = 14°.
Figure 48. (a) The detector grid for ITER. Calculated for Eb = 2.6–3.4 MeV with steps of 0.1 MeV and angle 0.5°. The line for α = 13° allows one to probe a fragment of the potential profile at the plasma periphery, 0.76 < r/a < 0.98, by changing the energy from shot to shot. The line Eb = 3.3 MeV allows one to obtain a fragment of the potential profile 0.76 < r/a < 1 by changing the injection angle β. (b) Toroidal shift in the equatorial projection of the detector line with α = 14°.
Symmetry 13 01367 g048
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Melnikov, A. Evolution of Heavy Ion Beam Probing from the Origins to Study of Symmetric Structures in Fusion Plasmas. Symmetry 2021, 13, 1367. https://doi.org/10.3390/sym13081367

AMA Style

Melnikov A. Evolution of Heavy Ion Beam Probing from the Origins to Study of Symmetric Structures in Fusion Plasmas. Symmetry. 2021; 13(8):1367. https://doi.org/10.3390/sym13081367

Chicago/Turabian Style

Melnikov, Alexander. 2021. "Evolution of Heavy Ion Beam Probing from the Origins to Study of Symmetric Structures in Fusion Plasmas" Symmetry 13, no. 8: 1367. https://doi.org/10.3390/sym13081367

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop