\documentclass[12pt,a4paper]{article}
\usepackage[utf8x]{inputenc}
\usepackage[USenglish]{babel}
\usepackage[T1]{fontenc}
%\usepackage{inputenc}
\usepackage{graphicx}
\usepackage{lmodern} %Type1-font for non-english texts and characters
\usepackage{ucs}
\usepackage{amsmath}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{color}
\usepackage{url}
\hoffset=-1cm
\voffset=-2cm
\oddsidemargin=0cm
\evensidemargin=0cm
\marginparwidth=0cm
\marginparsep=0cm
\textwidth=18cm
\textheight=24cm
\graphicspath{{figs/}}
\begin{document}
%\pagestyle{empty}
\begin{center}
\vspace{3 cm}
\Large
{\bf Instructions for student measurements on the GOLEM tokamak}\\
~\\*[0.5cm] \normalsize
Dr. Gergo Pokol$^{(1)}$, Laszlo Horvath$^{(1)}$, Csaba Buday$^{(1)}$, Daniel Imre Refy$^{(1)}$
\\
\vspace{1 cm}
{
\it \small $^{(1)}$ Department of Nuclear Techniques, Budapest University of Technology and Economics, Association EURATOM, M\H uegyetem rkp. 3-9., H-1111 Budapest, Hungary\\
}
\vspace{2 cm}
Last update: \today
\end{center}
\vspace{4 cm}
\begin{abstract}The instructions present a 6 hour long measurement program making use of the remote measurement potential of the GOLEM tokamak located at the Faculty of Physical and Nuclear Engineering of the Czech Technical University (CTU). The purpose of the measurement is to demonstrate the very basics of tokamak operation, and to get the students acquainted with basic properties and operational limits.
\end{abstract}
\newpage
\section{Introduction}\label{sec:introduction}
Participants of the present student measurement are strongly encouraged to read chapters 1 and 3 of John Wesson's The science of JET \cite{wesson99science} to be downloaded from
\url{http://www.iop.org/Jet/fulltext/JETR99013.pdf}.
This reading gives an introduction to the basic concept of thermonuclear fusion and magnetic confinement.
\begin{figure}[ht]
\centerline{\resizebox{100mm}{!}{\rotatebox{0}{\includegraphics{tokamak.pdf}}}}
\caption{Structure of a tokamak: (1) vacuum chamber, (2a) poloidal field coil / vertical field coil, (2b) toroidal field coil, (2c) transformer coil, (3) plasma, (4) plasma current, (5) magnetic field line, (7) radial direction ($r$).}
\label{fig:tokamak}
\end{figure}
Basic build up of a tokamak can be seen on Figure \ref{fig:tokamak}. The confining magnetic field structure is the result of the superposition of the toroidal magnetic field ($B_t$) generated by external coils and the poloidal magnetic field ($B_p$) generated by a strong toroidal plasma current ($I_{pl}$) induced by the transformer coil.
In the resulting magnetic geometry, field lines are winding helically around a torus surface, which is called the magnetic surface. The tokamak magnetic field consists of such nested magnetic surfaces. The helical structure at each magnetic surface is described by the safety factor ($q$). It gives the number of toroidal turns necessary for the magnetic field line at the given magnetic surface to reach its original position poloidally. On large aspect ratio circular tokamaks (like GOLEM), where the major radius ($R$) is much larger than the minor radius ($r_0$), it can be approximated by:
\begin{equation}
q(r,t)= \dfrac{r}{R} \dfrac{B_t(t)}{B_p(r,t)},
\label{eq:q_def}
\end{equation}
where $R$ is the major radius of the magnetic axis. An illustration for the meaning of the safety factor can be seen on Figure \ref{fig:q}.
\begin{figure}[ht]
\centerline{
\resizebox{80mm}{!}{\rotatebox{0}{\includegraphics{q_small.pdf}}}
\resizebox{80mm}{!}{\rotatebox{0}{\includegraphics{q_high.pdf}}}
}
\caption{Magnetic field lines in a tokamak for different safety factors.}
\label{fig:q}
\end{figure}
An important concept regarding the energy balance of the tokamak fusion reactor is the energy confinement time ($\tau_E$). It is the characteristic time of energy loss:
\begin{equation}
P_{loss}= \dfrac{W_{pl}}{\tau_E},
\label{eq:tau}
\end{equation}
where $P_{loss}$ is the power lost and $W_{pl}$ is the total plasma energy. The energy confinement time is a global parameter of the confined plasma, and reaching higher values is of central interest of tokamak research.
The Lawson criterion is a simple threshold for self-sustained thermonuclear fusion plasma burn at optimum temperature, and it also includes the energy confinement time along with plasma density ($n$):
\begin{equation}
n \tau_E > 10^{20} \;\rm s m^{-3}.
\label{eq:Lawson}
\end{equation}
More general information on fusion power production and the tokamak concept can be found at the following sites:
\url{http://www.magfuzio.hu}, \url{http://www.iter.org}, \url{http://www.jet.efda.org/}
\subsection{GOLEM}\label{sec:golem}
The GOLEM tokamak is a tokamak with full remote control capability and educational purpose. It is a small sized tokamak device equipped with basic controls and diagnostics having dimensions:
\begin{itemize}
\item Major radius at the magnetic axis: $R_0=0.4\;\rm m$.
\item Minor radius: $r_0=0.1\;\rm m$.
\item Radial position of the limiter: $a=0.085\;\rm m$.
\end{itemize}
The device was originally called TM1. Designed and constructed in Kurchatov Institute of Nuclear Research (Soviet Union), it was one of the first operational tokamaks in the world. The original concept of the device did not include poloidal field coils of stabilization however, it was believed that having a multiple layer metallic chamber enclosing layers of vacuum would help to achieve better stability of plasma column. This resulted in the realization of a liner inside the vacuum vessel that also acts as a plasma limiter. The capacitor battery for toroidal field coils and transformer filled several rooms.
\begin{figure}[ht]
\centerline{\resizebox{100mm}{!}{\rotatebox{0}{\includegraphics{golem_photo.pdf}}}}
\caption{Photo of the GOLEM tokamak}
\label{fig:golem_photo}
\end{figure}
Some time later, there was a microwave heating system integrated and the device was renamed to TM1-MH. The microwave heating, in addition to ohmic heating had to heat the plasma further. After the device was moved to Institute of Plasma Physics, Czech Academy of Sciences (IPP CAS) in September 1977, thanks to cooperation between the Kurchatov Institute and IPP, some changes in the engineering took place. The microwave heating system was left in Russia, alongside with the most of the oil capacitors of the toroidal field generation, since there was not enough room in the new tokamak hall. A few years later, the device went under major reconstruction. The vacuum vessel was replaced for a new one, the layer of vacuum between the liner and coating were fully removed and a feedback stabilization system was integrated instead. The power supply was substituted by a stronger one, and the ignition was replaced by a glow discharge. Between the years 1977 and 2007 there were several small changes over the device, such as the use of new diagnostics sensors.
In the end of 2007 the device was transferred to the Faculty of Physical and Nuclear Engineering of the Czech Technical University (CTU). Work on the reoperation started on the 14th of July 2008 with limited capabilities, and improvements are still underway. Further upgrade of GOLEM is envisaged in a near future - an increase of $B_t$, $I_p$ and the discharge duration. Dynamic plasma position stabilization is under present consideration and investigation. Basic diagnostics will be enriched with the plasma density measurement (microwave interferometer), $H_{\alpha}$ and X-ray radiation measurement will be installed in a near future. Investigation of plasma edge physics with the help of the various probe measurements is planned, as the previous version of the GOLEM tokamak, the CASTOR had a very good inspiring tradition in this field of interest.
More information can be reached via the home site of the GOLEM tokamak:
\begin{center}
\url{http://golem.fjfi.cvut.cz}
\end{center}
Detailed experimental arrangement can be seen at the following link (\emph{This should be studied in detail before the measurement!}):
\begin{center}
\url{http://golem.fjfi.cvut.cz/?p=tokamak}
\end{center}
\textbf{The parameters to be set remotely:}
\begin{itemize}
\item Toroidal magnetic field ($B_t$) through the voltage of the toroidal field capacitor bank ($U_B=U_{C\_Bt}$), range: $400-1400 \;\rm V$.
\item Toroidal electric field ($E_t$) through the capacitor bank for the current drive ($U_E=U_{CD}$), range: $100-600 \;\rm V$.
\item The time delay between the triggers of the toroidal magnetic field and the current drive ($T_{CD}=\tau_{OH}$), range: $0-20000 \;\rm \mu s$.
\item Hydrogen gas pressure ($p_{H2}$), range: $0-100 \;\rm mPa$.
\item Preionization ON/OFF
\end{itemize}
Figre \ref{fig:delays} shows the effect of time delay parameter.
\begin{figure}[ht]
\centerline{\resizebox{120mm}{!}{\rotatebox{0}{\includegraphics{kresba.pdf}}}}
\caption{Time delay parameters.}
\label{fig:delays}
\end{figure}
\textbf{The diagnostics used during the session to be accessed online:}
\begin{itemize}
\item Time resolved measurement of loop voltage ($U_l$).
\item Time resolved measurement of total toroidal current by Rogowski coil ($I_t$).
\item Time resolved toroidal magnetic field by coil measurement ($B_t$).
\item Time resolved measurement of plasma radiation by photodiode.
\item Vacuum chamber pressure ($p_{ch}$).
\item The temperature of the vacuum chamber ($T_{ch}$).
\end{itemize}
\section{Measurement procedure}\label{sec:procedure}
This section summarizes the technical procedures necessary for the remote control of the GOLEM tokamak.
\subsection{Communication with local support}
Primary real time communication to the local support (Dr. Vojtech Svoboda) is through Skype instant messaging. For this purpose a user (name: nti.hallgato, passwd: vendeg0 or name: nti.hallgato2, passwd: vendeg1) has been created. Before starting the experiment an instant messaging conference is to be started including users nti-hallgato, gergo\_pokol (Dr. Gergo Pokol, Hungarian supervisor) and tokamak.golem (Dr. Vojtech Svoboda, Czech support). After the session the conversation log should be saved and appended to the measurement log.
\subsection{Remote control}
Measurements are to be set up and shots initiated using the web interface of GOLEM tokamak, which can be seen on figure \ref{fig:remote_control}. The exact url address of it is provided by Dr. Vojtech Svoboda just at the beginning of the session.
\begin{figure}[ht]
\centerline{\resizebox{100mm}{!}{\rotatebox{0}{\includegraphics{remote_control.pdf}}}}
\caption{Remote control interface of GOLEM tokamak.}
\label{fig:remote_control}
\end{figure}
\subsection{Remote data access}
All the recorded data and the settings for each shot are available at the GOLEM website. The root directory for the files is:
\begin{center}
\url{http://golem.fjfi.cvut.cz/operation/shots/<shotnumber>/}
\end{center}
Basic data of the present shot series are collected at a page to be reached at:
\begin{center}
\url{http://golem.fjfi.cvut.cz/operation/currentsession/}
\end{center}
during the session and later archived at:
\begin{center}
\url{http://golem.fjfi.cvut.cz/operation/tasks/TrainingCourses/HUNTRAIC/}
\end{center}
In order to facilitate the procedure of data analysis, a MATLAB package is available for basic data processing. (This package is also compatible with the OCTAVE freeware software.) The task is to build a proper work flow using these building blocks. It should be noted that these routines do not cover the whole procedure, some additional programs are supposed to be written by the students. The routines are listed in the Table below:
\begin{tabular}{l|l|l}
File name&Input parameters&Description\\
\hline
GOLEM\_get\_data.m & shot\_nr&\begin{tabular}[x]{@{}c@{}}Loads raw data from database\\into the MATLAB workspace\end{tabular} \\
GOLEM\_plot\_rawdata.m & shot\_nr&\begin{tabular}[x]{@{}l@{}}Makes plots of the time varying\\raw data\end{tabular} \\
GOLEM\_offset\_correction.m & \begin{tabular}[x]{@{}l@{}}raw\_signal, time\_vector, \\t1,t2\end{tabular} &\begin{tabular}[x]{@{}l@{}}Makes offset correction\\for raw data\end{tabular} \\
GOLEM\_cut\_data.m &\begin{tabular}[x]{@{}l@{}}raw\_signal, time\_vector, \\t1,t2\end{tabular}&Crop the given signal \\
GOLEM\_integrate.m &time\_vec, signal&\begin{tabular}[x]{@{}l@{}}Integrates the given signal\end{tabular} \\
GOLEM\_chamber\_current.m &\begin{tabular}[x]{@{}l@{}} time\_vec, \\$I_t$, $U_l$, $R_{ch}$, $L_{ch}$\end{tabular}&\begin{tabular}[x]{@{}l@{}}Calculates chamber current\\integrating equation \eqref{eq:vac} \end{tabular} \\
GOLEM\_diff.m &$x$, $y$& Calculates $dx/dy$ \\
\end{tabular}
\paragraph{GOLEM\_get\_data.m}
The return value of GOLEM\_get\_data.m contains then \textbf{rawdata} structure with the following elements:
\begin{itemize}
\item \textbf{nr:} shotnumber
\item \textbf{timedata:} structure, contains vectors of time signals
\begin{itemize}
\item \textbf{t:} time axis vector in [s]
\item \textbf{U\_l:} loop voltage measurement raw signal vector in [V]
\item \textbf{dB\_t:} toroidal filed coil raw signal vector in [V]
\item \textbf{dI\_t:} Rogowski coil raw signal vector in [V]
\item \textbf{Photo:} photodiode raw signal vector in [V]
\end{itemize}
\item \textbf{N:} number of data points
\item \textbf{samplerate:} samplerate of the measurements in [Hz]
\item \textbf{pressure:} pressure of vacuum chamber in [mPa]
\item \textbf{T\_ch:} temperature of the chamber in [K]
\item \textbf{trigger:} time delay between starting diagnostics and toroidal magnetic field drive in [s]
\item \textbf{time\_delay:} time delay between toroidal field and inductive current drive in [s]
\item \textbf{Bt\_calibration:} calibration factor of toroidal magnetic field diagnostic in [T/Vs]
\item \textbf{Rogowski\_calibration:} calibration factor of plasma current diagnostics in [A/Vs]
\item \textbf{U\_loop\_calibration:} calibration factor of loop voltage diagnostic [V/V]
\end{itemize}
Elements of structures can be referenced as e.g. \textbf{rawdata.timedata.U\_l}. \emph{Measured signals are saved in the timedata structure, but these are raw signals needing further processing to produce the physical quantities measured!} Signal processing steps are described in the next section.
\section{Measurement tasks, method of evaluation}\label{sec:tasks}
This section starts with a short description on how to reconstruct the measured plasma parameters from the raw signals returned by GOLEM\_get\_data.m. Measurement tasks are detailed in the later subsections.
The sampling rate of the time resolved measurements (samplerate), time delay between starting diagnostics and toroidal magnetic field drive (trigger) and time delay between toroidal field and inductive current drive (time\_delay) are returned by GOLEM\_get\_data.m, and these are to be used whenever needed instead of the examples provided in this description.
The simplest signal to be reconstructed is the loop voltage ($U_l$): The measurement loop of the loop voltage is connected to a voltage divider, therefore the signal must be multiplied by a calibration factor (U\_loop\_calibration) as plotted in Figure \ref{fig:bd_U_l}.
\begin{figure}[ht]
\centerline{\resizebox{100mm}{!}{\rotatebox{0}{\includegraphics[trim=0cm 8.5cm 9cm 4cm, clip=true]{bd_U_l.pdf}}}}
\caption{Block diagram showing the steps of data processing for loop voltage measurement.}
\label{fig:bd_U_l}
\end{figure}
The toroidal magnetic field ($B_t$) and the total current ($I_{tot}$) raw signals must be integrated before multiplying by calibration factors (Bt\_calibration and Rogowski\_calibration). The reason for this is that the voltage measured is induced in these diagnostic loops and coils by the changing of the toroidal and poloidal magnetic field respectively.
Integrated magnetic measurements are very sensitive to the DC bias of the measurement circuit, which needs to be corrected for. If the sampling rate is 1 MHz, and the shot starts at 5 ms, we have 5000 samples from the background noise. It is better to exclude a few samples around the swithing time point. This is important, because these samples measure the bias, and we can correct the integrated values with this factor.
Figure \ref{fig:bd_dB_t} shows the block diagram for the necessary steps of processing for the toroidal magnetic field signal. Routines for all the steps are ready, they should just be parametrized and linked.
\begin{figure}[ht]
\centerline{\resizebox{160mm}{!}{\rotatebox{0}{\includegraphics[trim=1cm 8.5cm 0cm 4cm, clip=true]{bd_dB_t.pdf}}}}
\caption{Block diagram showing the steps of data processing for toroidal magnetic field measurement.}
\label{fig:bd_dB_t}
\end{figure}
The block diagram for the total current measured by the Rogowski coils is only slightly more complicated: Switching the toroidal magnetic field on causes an offset in the toroidal current measurement, which has to be corrected by subtracting the average value measured in the $\tau_{OH}$ long interval before switching on the toroidal electrical field from the integrated current value.
\begin{figure}[ht]
\centerline{\resizebox{160mm}{!}{\rotatebox{0}{\includegraphics[trim=1cm 5cm 0cm 4cm, clip=true]{bd_dI_t.pdf}}}}
\caption{Block diagram showing the steps of data processing for total plasma current measurement.}
\label{fig:bd_dI_t}
\end{figure}
\subsection{Determination of vacuum chamber parameters}\label{sec:vacuum}
In GOLEM, part of the toroidal current always flows in the vacuum vessel, which has to be taken into account during the interpretation of experimental results. In a vacuum shot, when no plasma is formed, it is possible to determine the resistance of the vacuum vessel: all the current measured by the Rogowski-coil flows in the vessel. This is an important parameter for further evaluations.
Let us denote the loop voltage with $U_{l}$, the resistance of the chamber by $R_{ch}$, the total current (which is the chamber current ($I_{ch}$) in this case) with $I_{tot}$ and the inductance of the chamber by $L_{ch}$.
The circuit equation is then
\begin{equation}
U_{l}(t)= R_{ch} \cdot I_{tot}(t) + L_{ch} \dfrac{d I_{tot}}{d t}.
\label{eq:vac}
\end{equation}
Using the loop voltage measurement and the Rogowski-coil, we have both $U_{l}$, $I_{tot}$ and $dI_{tot}/dt$ measured, so $R_{ch}$ and $L_{ch}$ can be determined.
A simple method is the following: Just after switching on the toroidal electric field, the toroidal current is still close to zero ($I_{tot}\approx 0$), so $U_{l}\approx L_{ch} dI_{tot}/dt$, so $L_{ch}$ can be determined. On the other hand, at the flat top of the current curve ($dI_{tot}/dt \approx 0$) equation \eqref{eq:vac} simplifies to $U_{l}\approx R_{ch} \cdot I_{tot}$, so $R_{ch}$ can be estimated.
A more sophisticated method is a 2D least squares linear fit making use of all data points $(U_{l},I_{tot},dI_{tot}/dt)$. Since we have only two independent parameters $R_{ch}$ and $L_{ch}$, the fitted plane has to pass through the origin. If we divide equation \eqref{eq:vac} by $I_{tot}$, we can simplify the task to a 1D least squares linear fit, which can be easily implemented in MATLAB (OCTAVE), using \textit{polyfit} function.
Values of $R_{ch}$ and $L_{ch}$ should be calculated for about \textbf{5} discharges having different parameters, and the results should be compiled to a single best estimate for both parameters. Estimation should be performed by both methods described above, and the results of the method giving the more precise estimates should be used in the further steps.
\subsection{Plasma breakdown}\label{sec:breakdown}
After measuring the vacuum chamber properties, we can make the next step towards creating a tokamak plasma: we can let $H_2$ gas into the chamber before initiation of the toroidal electric field. The $p_{H2}$ value, which can be set as a discharge parameter, is a control parameter for the inlet valve. The actual value of the pre-discharge gas pressure is measured by a vacuummeter $p_{ch}$.
As we will see, letting $H_2$ gas into the chamber is not always sufficient to produce a plasma. The toroidal electric field must also reach a critical value for mass ionization, in other words plasma breakdown.
The task is to plot the $p_{ch}$ against the maximum of the loop voltage spikes in the beginning of the discharge for several discharges, and indicate the plasma breakdown by the shape of the symbols. Shots should be concentrated around the critical line separating breakdown and non-breakdown shots. Detailed scan should be performed for a given magnetic field and the effect of the magnetic field should be studied with a few discharges. During this exercise the pre-ionization should be turned on to produce more reproducible results, but the effect of turning it off could also be studied. About a total of \textbf{30} discharges are available for this exercise.
\subsection{Estimation of main plasma parameters}\label{sec:parameters}
If plasma breakdown occurs, plasma parameters can be determined - with different accuracy - from the measured parameters. The aim of this task is to investigate the effect of different parameters on the performance of the discharge, and reach discharges with the highest central temperature, plasma energy or energy confinement time. This task should result in about \textbf{25} discharges.
\subsubsection{Plasma current}
A simple electrical model for the inductive current drive is a time-varying voltage source ($U_l(t)$) connected to the plasma and the vacuum chamber in parallel can be seen on Figure \ref{fig:circuit}. Both the vacuum chamber and the plasma are modeled by LR circuits. The main difference is, that while the internal inductance and resistance of the chamber are constant, and thus they can be measured separately, the parameters of the plasma differ in each discharge.
\begin{figure}[ht]
\centerline{\resizebox{120mm}{!}{\rotatebox{0}{\includegraphics{circuit.pdf}}}}
\caption{Model of the inductive current drive circuit}
\label{fig:circuit}
\end{figure}
The basic circuit equations are:
\begin{align}
\label{eq:chamberLR}
U_l(t) &= R_{ch} \cdot I_{ch}(t) + L_{ch} \dfrac{d I_{ch}(t)}{d t}\\
\label{eq:plasmaLR}
U_l(t) &= R_{pl}(t) \cdot I_{pl}(t) + L_{pl} \dfrac{d I_{pl}(t)}{d t}\\
I_{tot}(t)&= I_{pl}(t) + I_{ch}(t)
\end{align}
The chamber parameters have already been determined according to Section \ref{sec:vacuum}. Integration of the \eqref{eq:chamberLR} circuit equation using the initial condition $I_{tot}(t=0)=I_{ch}(t=0)$ is implemented in the routine GOLEM\_chamber\_current.m to arrive to $I_{ch}(t)$. This can then be used to determine the plasma current, as $I_{pl}(t)=I_{tot}(t)- I_{ch}(t)$ as shown in Figure \ref{fig:block_diagrams_current}. Plasma resistivity can be determined in turn from equation \eqref{eq:plasmaLR}.
\begin{figure}[ht]
\centerline{\resizebox{160mm}{!}{\rotatebox{0}{\includegraphics[trim=1cm 5cm 0cm 4cm, clip=true]{block_diagrams_current.pdf}}}}
\caption{Block diagram showing the steps of data processing for the plasma current measurement.}
\label{fig:block_diagrams_current}
\end{figure}
Having calculated the plasma current, a threshold can be defined significantly exceeding the calculation accuracy to safely determine the beginning and end of the plasma discharge. Using this threshold in the "find" function, one can cut the time signals to the extent of the discharge for further processing. Time duration of the dsischarge is also an imprtant parameter.
\emph{This task needs some programming that should be done parallel to the task described in Section \ref{sec:breakdown}!}
It can be attempted to investigate the effect of the $L_{pl}\approx 0 \;\rm H$ approximation by a more careful integration of choosing $L_{pl}\approx L_{ch}$ in the time region with plasma. If significant differences are found, this latter approximation has to be implemented for all further data processing.
A suitable threshold in plasma current can be used to determine the discharge duration and cut out the interval of the measured signals relevant for plasma diagnostics.
Plasma current has to be calculated for all discharges with plasma and the maximum value and the discharge duration have to be included in the shot summary table.
\subsubsection{Plasma heating power}
In the GOLEM tokamak the only heating mechanism of the plasma is ohmic heating resulting from current flowing in a conductor with finite resistivity. The ohmic heating power can be calculated as:
\begin{eqnarray}
P_{OH}(t) = R_{pl}(t) \cdot I_{pl}^2(t)
\end{eqnarray}
Ohmic heating power has to be calculated for all discharges with plasma and the maximum value has to be included in the shot summary table.
\subsubsection{Central electron temperature}
Specific resistivity of a fully ionized plasma only depends on its electron temperature ($T_e$) and effective charge number ($Z_{eff}$). This dependence is quantified by the Spitzer formula \cite{nrl_formulary}. It has to be noted that the ion temperature can be very much different from electron temperature. The effective charge number is determined by the amount, composition and state of impurities in the $H_2$ plasma, and we can take value $Z_{eff}\approx 2.5$ for GOLEM plasmas.
Center of the plasma has higher temperature, and lower resistivity with higher current density, which makes the estimation of the electron temperature ambiguous from an integrated value of resistivity ($R_{pl}(t)$). However, if we use an equilibrium temperature profile \eqref{eq:temp_prof} (Figure \ref{fig:temp_profile}), measured in more detailed measurements \cite{brotankova}, we can estimate one parameter of the profile, which is in this case the central electron temperature ($T_{e0}(t)$):
\begin{equation}
T_{e}(r,t)=T_{e0}(t)\left( 1-\dfrac{r^2}{a^2} \right)^2
\label{eq:temp_prof}
\end{equation}
\begin{figure}[ht]
\centerline{\resizebox{100mm}{!}{\rotatebox{0}{\includegraphics{temp_profile.pdf}}}}
\caption{Equilibrium temperature profile used in the estimation of central plasma temperature}
\label{fig:temp_profile}
\end{figure}
The central electron temperature ($T_{e0}$) is then calculated from equation (3.20) of \cite{brotankova}, which itself is based on Spitzer's resistivity formula:
\begin{equation}
T_{e0}(t)=\left(\dfrac{R_0}{a^2}\dfrac{8 Z_{eff.}}{1544}\frac{1}{R_{pl}(t)}\right)^{2/3},
\end{equation}
where $R_{pl}(t)$ is in Ohms, distances are in meters and we get $T_{e0}(t)$ in electronvolts.
It has to be noted that plasma in the GOLEM tokamak is only fully ionized in the central region, $Z_{eff}$ can be estimated with large uncertainty and even the $a$ plasma small radius might change in an unmonitored way due to the lack of plasma stabilization. All these factors make the estimation of the central electron temperature quite uncertain.
Nevertheless, central electron temperature has to be calculated for all discharges with plasma and the maximum value has to be included in the shot summary table.
\subsubsection{Electron density}
In its current state, the GOLEM tokamak does not have any density measurements. However, as electron density is needed for further calculations, we estimate its order of magnitude from the state law of ideal gases.
For the average density it is assumed, that it is constant during the discharge, apart from the dissociation of the hydrogen gas. There is a 30 second delay between the gas filling and the actual shot, which is enough for the gas to reach thermal equilibrium with the chamber wall. Chamber temperature is monitored with respect to the room temperature, and the difference is normally be zero, but should be checked. (If chamber temperature is not measured, room temperature can be used instead.) The ideal gas law is used to give an order of magnitude estimate of the electron density (in particle/$m^3$):
\begin{equation}
n_{avr} = \dfrac{2p_{ch}}{k_B T_{ch}}.
\end{equation}
We have to note that this is a very rough estimate basically for two reasons:
\begin{enumerate}
\item Plasma in the GOLEM tokamak is not fully ionized, which makes us overestimate the electron density.
\item Due to the plasma-wall interaction, adsorbed gases are released from the surface of plasma facing components during the discharge. These atoms enter the plasma and can be ionized, thus making us underestimate the electron density.
\end{enumerate}
The order of magnitude estimate of the average electron density has to be calculated for all discharges with plasma and included in the shot summary table.
\subsubsection{Plasma energy}
The total energy content can be simply calculated from the temperature, density and volume ($V$), based on the ideal gas law, taking into account the assumed \eqref{eq:temp_prof} temperature profile:
\begin{equation}
W_{pl}(t)=V \dfrac{ n_{avr} k_B T_{e0}(t)}{3}.
\end{equation}
The information that the magnetic field reduces the degrees of freedom of the particles to two has been used to derive this formula.
Uncertainty of this formula is dominated by the uncertainty of our density estimate, which makes it good only for an order of magnitude estimate. Qualitative time trace reflects that of the electron temperature and thus is more reliable.
Nevertheless, plasma energy has to be calculated for all discharges with plasma and the maximum value has to be included in the shot summary table.
\subsubsection{Energy confinement time}
Having an estimate for the plasma energy, the energy confinement time can be estimated. The loss power can be estimated from the energy balance:
\begin{equation}
P_{loss}(t)=P_{OH}(t)-\dfrac{d W_{pl}}{d t}
\end{equation}
We then have to just substitute it into the definition (\ref{eq:tau}) of the energy confinement time:
\begin{equation}
\tau_E(t)=\frac{W_{pl}(t)}{P_{loss}(t)}.
\end{equation}
Given the uncertainty of the input parameters, maximum value for the energy confinement time should be taken with care. Nevertheless, its maximum should be included in the shot summary table.
\subsection{$q=2$ disruptions}\label{sec:disruptions}
When the plasma current grows so strong that the edge safety factor, defined by \eqref{eq:q_def}, reaches the value of 2, a plasma instability resonant to the $q=2$ rational surface destabilizes, and a discharge terminating disruption occurs. This limit of operation is to be attempted to be reached in this task using about \textbf{5} dedicated shots.
We can calculate the poloidal field at the edge (for large aspect ratio circular tokamaks) using Amp\`{e}re's law, as the enclosed current is the total plasma current:
\begin{equation}
B_p(a,t)= \dfrac{\mu_0}{2\pi} \dfrac{I_{pl}(t)}{a},
\end{equation}
where $a$ is the plasma minor radius. Substituting this expression into formula \eqref{eq:q_def}, the safety factor at the edge can be estimated as:
\begin{equation}
\label{eq:q}
q(a,t)= \dfrac{a^2}{R_0} \dfrac{2 B_t(t) \pi}{\mu_0 I_{pl}(t)}.
\end{equation}
Discharges aiming to reach a low $q(a,t)$ need as large plasma current as possible. As we have very limited control over the evolution of plasma current in GOLEM, we can also set the $\tau_{OH}$ time delay to set up a discharge at the declining phase of the toroidal magnetic field, which will constantly decrease the edge safety factor.
In order to monitor the success of our efforts, the evolution of the discharges should be plotted on the Hugill diagram. The Hugill diagram positions a discharge on the plane of two parameters:
\begin{itemize}
\item Inverse edge safety factor: $\dfrac{1}{q(a,t)}$
\item Murakami parameter (normalized density): $\dfrac{n_{avg} R_0}{B_t(t)}$
\end{itemize}
The Hugill diagram serves as an operation envelope for tokamaks. If either the Murakami parameter is too high or the inverse edge safety factor reaches the value of 0.5, the plasma disrupts.
First, the temporal evolution of the dedicated shots aiming $q=2$ disruptions should be plotted on the Hugill diagram. Afterwards, all previous shots could be plotted to check that none reach the region $1/q>0.5$.
\section{Requirements for the measurement logbook}\label{sec:requirements}
There should be a single measurement log written in English language. The measurement log should include:
\begin{itemize}
\item Exact method of executing the measurement.
\item Exact method of derivation of final results.
\item Results of the measurement tasks.
\item Time traces of shot parameters for some characteristic shots.
\item Shot summary table indicating all calculated parameters of all shots executed. (Unsuccessful shots also need to be recorded.)
\item Skype log as appendix
\end{itemize}
All comments and proposals regarding the measurements are welcome.
It would be nice, if the participants of the session would write some thankful words to the \url{http://golem.fjfi.cvut.cz/hodnoceni.html} homepage, and/or send a postcard to the following address: Tokamak GOLEM, Brehova st. 7, Prague 1, Czech Republic. These activities should also be indicated in the logbook.
\bibliographystyle{unsrt}
\bibliography{refs1}
\end{document}
