Diagnostics/ParticleFlux/RailProbe/Reports/Macha-2021-ECPP/main.tex

\documentclass{epsconf}
\usepackage{graphicx}
\usepackage{wrapfig}
\usepackage{color}
\title{Tokamak GOLEM for fusion education - chapter 12}

\author{\underline{P. Macha}$^{1,2}$, K. Hromasova$^{1,2}$, D. Kropackova$^{1}$, M. Lauerova$^{3}$, A. Socha$^{4}$, J. Malinak$^{1}$, \\ D. Cipciar$^{5}$, J. Cecrdle$^{1}$, V. Svoboda$^{1}$, J. Stockel$^{2}$, J. Adamek$^{2}$,  F. Papousek$^{1}´$ , L. Lobko$^{1}$ \\
$^{1}$ \textit{Faculty of Nuclear Sciences and Physical Engineering, CTU, Prague, Czech Republic} \\
$^{2}$ \textit{Institute of Plasma Physics, Czech Academy of Sciences, Prague, Czech Republic} \\
$^{3}$ \textit{Gymnazium Nad Aleji 1952, Prague, Czech Republic} \\
$^{4}$ \textit{Gymnazium a Stredni odborna skola Cihelni 410, Frydek-Mistek,} \textit{ Czech Republic} \\
$^{5}$ \textit{Faculty of Science, Masaryk University, Brno, Czech Republic}
}

\begin{document}

\maketitle

The GOLEM tokamak \cite{Svoboda, Svoboda2} is the oldest tokamak in the world. Currently, it serves mainly as an education device for students of tokamak physics. Remote control of the machine enables conducting experiments from all over the world. This contribution summarizes its main research topics of the last year.

% External Stabilization
\textbf{Due to stray fields from current drive windings, the GOLEM plasma is shifted vertically toward the upper wall of the chamber and horizontally toward the low field side.} Two external windings generating poloidal magnetic field were used to control the plasma position. A square waveform was selected as a predefined function and four KEPCOs connected in a parallel were used as power supply providing a current of 100 A. The resulting plasma position was determined from the signal of four Mirnov coils installed at the limiter.
%The influence of the horizontal magnetic field generated by the current in external coil can be observed as an effect on plasma column vertical displacement and the associated discharge duration.
This predefined stabilisation successfully increased the discharge duration. Experiments are ongoing to determine the optimal current in the external windings and the number of turns per coil.

%Long-range Langmuir probe correlations
\textbf{Two toroidally separated Langmuir probes (LPs) were used to study long-range correlations (LRC) in the edge plasma.} LRC were sought at the expected GAM frequency $f_{GAM} \approx 15$~kHz in both $V_\mathrm{fl}$ and $I_\mathrm{sat}$ regimes in the range $r = 65-75$ mm. A low-frequency LRC candidate for GAM was, indeed, found and it was investigated whether it has GAM properties, that is, whether the two LP signals have high coherence and near zero cross-phase. However, the explored frequency region did not exhibit near-zero cross-phase and coherence above $0.7$ at the same time at any radial position. This may suggest that GAM-like oscillations do not exist naturally in the explored plasma regime in the GOLEM tokamak. Future experiments will attempt to facilitate such a mode with edge plasma biasing.

% Fast ion temperature measurement in the GOLEM tokamak
\textbf{Ion temperature was measured with 5 $\mathrm{\mu s}$ temporal resolution using a swept ball-pen probe.} \cite{Adamek} The probe collector was swept from -30 V to +130 V at 100 kHz. The ion temperature $T_\mathrm{i}$ was obtained from the 4-parameter fit of the electron I-V characteristics branch

\begin{equation}
    I(V) = I_\mathrm{sat}^+ \left( \exp (\alpha_\mathrm{BPP}) \cdot \left[ 1+K \cdot \left( V - \Phi \right) \right] - \exp \left( \frac{\Phi - V}{T_\mathrm{i}} \right) \right),
\end{equation}
where $\alpha_\mathrm{BPP} = \ln \left( \frac{I^{-}_\mathrm{sat}}{I^{+}_\mathrm{sat}}  \right) = 0.25 \pm 0.09$ for GOLEM \cite{Diplomka}, $I^{-}_\mathrm{sat}$ and $I^{+}_\mathrm{sat}$ are the electron and ion saturation current, respectively, $\Phi$ is the plasma potential and $K$ is an empirical parameter facilitating linear increase of the electron current. The cut-off fitting technique (Fig. \ref{fig:dario} left) was applied to minimise the role of the cut-off voltage $V_{\mathrm{max}}$ on the fit results. Fig. \ref{fig:dario} middle shows the temporal evolution of $T_\mathrm{i}$ (measured with the described technique) and $T_\mathrm{e}$ (measured using the BPP+LP technique \cite{Adamek_bpp}) during a 1 ms time interval in main SOL. Fluctuations of the ion temperature ranging between 5 eV up to 40 eV reflect the turbulent behavior of the edge plasma and the resulting temperature histograms (Fig. \ref{fig:dario} right) are positively skewed.

\begin{figure}[h]
    \includegraphics[width=0.28\textwidth]{2150Example_full10.1507958_shot36709.png}
    \includegraphics[width=0.35\textwidth]{Ti_Te_2Discharges65mm_36709.png}
    \includegraphics[width=0.35\textwidth]{Histogram_CONF.png}
    \vspace{-0.5cm}
    \caption{Left) Cut-off technique. Mid) Evolution of $T_\mathrm{i}$ and $T_\mathrm{e}$. Right) Histograms of $T_\mathrm{i}$ and $T_\mathrm{e}$.}
    \label{fig:dario}
\end{figure}

\textbf{Lithization tests were conducted in a small vacuum cube.} A glow discharge ($U<1$ kV, $I=0.5$ A) in argon atmosphere between the chamber wall and a lithium electrode will be used in the future to create a homogeneous lithium film on the GOLEM wall. We have run tests of the lithization setup which revealed two critical problems. The first problem is the prompt oxidation of the metallic lithium during the installation of the electrode. The second problem is the risk of electrode melting. After about 15 minutes of the glow discharge, the electrode reached melting temperature. The solution is either to decrease the power in the circuit (preferably voltage) or to create a more homogeneous and unoxidized electrode to spread the power over a larger area. A positive result is that a spectral line of neutral Li has been observed in the glow discharge, meaning that lithium in, indeed, being sputtered from the electrode.

% Electron temperature measurements using rail probe on the tokamak GOLEM
% \textbf{In 2017, a new experimental type of flush-mounted probe called rail probe was used on Alcator C-Mode tokamak.} 


\textbf{A new rail probe (RP) was constructed and used to measure the electron temperature.} Rail probes are suitable for measurements in the tokamak divertor region as they can sustain high heat flux and reduce the sheath expansion effect. A probe head %(Fig. 2 top)
carrying a rail probe (length 40 mm, width 2 mm), a Langmuir probe (length 1.5 mm, diameter 1 mm) and a ball-pen probe \cite{Adamek_bpp} has been designed and installed on a new manipulator which can modify its inclination with
\begin{wrapfigure}{r}{0.4\textwidth}
    \centering
    %\includegraphics[width=0.25\textwidth]{FIG_NEW.png}
    \includegraphics[width=0.35\textwidth]{FIG2.png}
    \vspace{-0.5cm}
    \caption{%Top) The combined probe head. Bottom) 
    Comparison of $T_\mathrm{e}$ measured by LP, RP and BPP+LP.}
    \label{fig:RP}
\end{wrapfigure}
respect to the toroidal magnetic field within $\pm 10 ^\circ$ with a precision of $0.6 ^\circ$. First comparative measurements of the electron temperature $T_\mathrm{e}$ using a swept RP, a swept LP ($f = 5$ kHz) and the BPP+LP technique \cite{Adamek_bpp} while varying the probe head inclination.
The electron temperature was obtained from I-V characteristics using a four-parameter fitting formula with linear sheath expansion approximation, although the sheath expansion effect was observed smaller in the RP than in the LP. The temporal evolution of the electron temperatures of the RP, LP, and BPP is shown in Fig \ref{fig:RP} for two different probe head inclination angles. In the flat-top phase with the constant plasma current, good agreement of the three probe techniques is observed even for a very small inclination angle.

% Measurement of HXR radiation at the GOLEM tokamak generated by runaway electrons by using different scintillation detectors Lukáš Lobko 
\textbf{Four scintillation detectors were installed to measure runaway electron energy spectra using HXR spectrometry.} Two problems were encountered, shown in Fig. \ref{fig:HXR}. Electronics of standard photomultiplier tubes of scintillation detectors does not withstand too large HXR fluxes from REs generated in tokamaks. This problem is noticeable in Fig. \ref{fig:HXR} left, when NaI(Tl) detector signal (green color) dramatically drops at around 8 ms, while other detectors measure intense HXR signal. The second problem is the pile-up effect, illustrated in Fig. \ref{fig:HXR} right. Compare for each detector the areas where the signal has piled up and the areas where individual peaks are detected. These two problems are partially related. Our goal is to find optimal setup in which they will be negligible. This can be done by ensuring sufficient lead shielding and distance from the GOLEM tokamak.
\begin{figure}[h!]
    \centering
    \includegraphics[width=0.5\textwidth]{ffig1.png}
    \hfill
    \includegraphics[width=0.4\textwidth]{ffig2.png}
    \vspace{-0.5cm}
    \caption{Left) Comparison of HXR signals from 4 different scintillation detectors. Right) Comparison of piled-up signals and individual peaks.}
    \label{fig:HXR}
\end{figure}

%Measurements of edge turbulent structures using a Langmuir rake probe.
\begin{wrapfigure}{r}{0.35\textwidth}
    \centering
    %\includegraphics[width=0.25\textwidth]{FIG_NEW.png}
    \includegraphics[width=0.35\textwidth]{edited_picture3_small.PNG}
    \vspace{-1cm}
    \caption{Radial profile of ion saturated current skweness.}
    \label{fig:Isat}
\end{wrapfigure}
\textbf{The ion saturated current $I_\mathrm{sat}$, collected by a double rake probe, was used to demonstrate the presence of turbulent structures in the GOLEM edge plasma.} 
Its radial profile was measured in the range $r = 37-90$ mm, with $r = 85$ mm being the limiter radius. The $I_\mathrm{sat}$ fluctuations histograms were found skewed to the right, which indicates the presence of high positive fluctuations - blobs - in the edge plasma. $I_\mathrm{sat}$ skewness remained positive throughout the measured region, only hinting at a fall to zero or negative values at its inner edge, see Figure \ref{fig:Isat}. This suggests that the "blob birth zone" is located deeper in the chamber, in agreement with a basic LCFS reconstruction.
%This result was published in the form of a competitive high-school student research thesis (SOČ).

% Comparison of probe methods for Te measurement. 
% \textbf{An important part of thermonuclear fusion research is the study of edge plasma, a region available for measurements with electrostatic probes.} 
%
\begin{wrapfigure}{r}{0.35\textwidth}
    %\vspace{-1cm}
    \centering
    %\includegraphics[width=0.25\textwidth]{FIG_NEW.png}
    \includegraphics[width=0.35\textwidth]{ffig_small.png}
    \vspace{-1cm}
    \caption{Time evolution of $T_\mathrm{e}$.}
    \label{fig:Te}
\end{wrapfigure}

\textbf{A systematic comparison of two $T_e$ measurement techniques, using a swept Langmuir probe and the BPP+LP technique \cite{Adamek_bpp}, was carried out.} The edge electron temperature is commonly measured by a swept LP. This method has been verified by time and numerous experiments and models, but applying it is rather complicated and time-consuming. Conversely, the innovative combined method using a floating LP and BPP probe is simple to apply and offers higher time resolution, but its experimental and theoretical verification is somewhat lacking. The two methods were compared in a series of identical discharges. The four-parameter V-A characteristics fit was found more appropriate, as the ion current branch did not fully saturate. The two methods yielded identical results, see Figure \ref{fig:Te}. We therefore suggest that the combined method is suitable for the essential task of measuring the edge electron temperature, as it can make tokamak edge $T_e$ measurements readily available without the need for complicated and hand-checked data processing.

% \begin{figure}[h!]
%     \centering
%     \vspace{-0.5cm}
%     \includegraphics[width=0.275\textwidth]{edited_picture3_small.PNG}
%     \includegraphics[width=0.292\textwidth]{ffig_small.png}
%     \caption{Left) Radial profile of ion saturated current skweness. Right) Time evolution of $T_\mathrm{e}$.}
%     \label{fig:my_label}
% \end{figure}

\textbf{Acknowledgement: }This work was supported by the Grant Agency of the Czech Technical University in Prague, grant No. SGS19/180/OHK4/3T/14.
% [J. Adámek, Review of Scientific Instruments 87 (2016) 043510]

\begin{thebibliography}{99}
% \vspace{0.1cm}
\bibitem{Svoboda}
V. Svoboda et al, Fusion Engineering and Design \textbf{68} (2011) 1310-1314
% \vspace{-1.5em}
\bibitem{Svoboda2}
V. Svoboda et al, Journal of Fusion Energy \textbf{38} (2019) 253–261
% \vspace{-1.5em}
\bibitem{Adamek}
J. Adamek et al, Nuclear Fusion \textbf{61} (2021) 036023
% \vspace{-2.5em}
\bibitem{Diplomka}
D. Cipciar, \textit{Ion and electron temperature study in the edge plasma of the tokamak device}, MSc. thesis, Faculty of Nuclear Sciences and Physical Engineering, Czech Technical University, to be published in 2021
\bibitem{Adamek_bpp}
J. Adamek et al, Review of Scientific Instruments \textbf{87} (2016) 043510
\vspace{0cm}
\end{thebibliography}


\end{document}