Double rake probe spotlight

The following scripts introduces basic aspects of the double rake probe measurements.

The double rake probe diagnostic

This is the double rake probe:

Its wiki page is located here. The double rake probe has two rows of pins ("rakes"), out of which the first 6 in each row are connected to the data acquisition system. The pin signals are called DRP-R1 to DRP-R6 and DRP-L1 to DRP-L6.

Like any Langmuir probe, the rake probe pins can be electrically insulated (measuring the floating potential $V_f$), biased to a negative voltage -100 V (measuring the ion saturated current $I_{sat}$), their biasing voltage may be swept (measuring $V_f$, $I_{sat}$ and the electron temperature $T_e$ if the $I$-$V$ characteristic is fitted with an exponential function) or they can perform other measurements. During this campaign, all the pins measure the ion saturated current $I_{sat}$. This can be processed to gain information about the plasma turbulent transport.

Import basic libraries

To visualise and process the double rake probe data, first we import basic libraries, Numpy and Matplotlib. The line %matplotlib notebook enables the drawn figures to be interactive.

Load double rake probe data

The data directory of the double rake probe is http://golem.fjfi.cvut.cz/shots/{shot}/Diagnostics/DoubleRakeProbe/, where {shot} stands for the discharge number. Note that in different sessions, the data may be saved elsewhere and it might be needed to update the URL variable in the following text. The data directory may be located from the individual shot homepage, and tends to be the same within a single session.

In the following, we write a function to download the rake probe data.

Notice that in the last line, the data values are divided by 46.7. This is the constant which converts the raw DAS signal to the ion saturated current in A. Its physical meaning is the resistance of the measuring resistor in the "silver box".

Load the discharge beginning and end

Prior to plotting the double rake probe data, we'll write a little function which loads the time of the discharge beginning (plasma formation) and end (recombination). We'll use it to set the X axis limits later. Notice that t1 and t2 are in ms.

Plot all the signals from the double rake probe pins

Next, we load the double rake probe data for the current discharge and plot them.

Notice how the signal wanes on the pins with a higher number. This is because they are not as deep in the plasma column and they are surrounded by plasma of lower temperature and density, hence lower $I_{sat}$.

Load the rake probe position

Plot histogram of the $I_{sat}$ fluctuations

The basic method of investigating the presence of turbulent structures in the edge plasma is plotting the $I_{sat}$ histogram. A histogram is an approximation of the distribution function - the probability of a certain value occurring in the signal. It is plotted by the function pl.hist().

In the following, we plot the histogram of the $I_{sat}$ fluctuations in the test discharge shot. The fluctuations are sampled during a 1ms window from 8 ms to 9 ms (approx. middle of the discharge).

Notice that the histogram might be asymmetrical.

Calculate the distribution skewness

Skewness $S$ is the third statistical moment of a distribution function, the first two being the mean and the standard deviation. If the distribution is symmetric, its skewness is zero. For instance, the skewness of the normal (Gaussian) distribution is zero. The skewness of the $I_{sat}$ fluctuations in this case, however, isn't zero and the distribution is not normal. It has an abundance of large positive events (blobs) while lacking the same number of large negative events. This is evidence of the presence of blobs in the SOL.

In the following, we calculate the skewness using the function scipy.stats.skew.

Plot the skewness profile from all probes

In the following, we use all the probe pins and plot the profile of the $I_{sat}$ fluctuations skewness. The distance between the pins is 2.5 mm.