Studies of High Voltage Breakdown Phenomena On ICRF (Ion Cyclotron Range of Frequencies) Antennas
Studies of High Voltage Breakdown Phenomena On ICRF (Ion Cyclotron Range of Frequencies) Antennas
Studies of High Voltage Breakdown Phenomena On ICRF (Ion Cyclotron Range of Frequencies) Antennas
Volodymyr Bobkov
Coupling of ICRF (Ion Cyclotron Range of Frequencies) power to the plasma is one of the
standard methods to heat plasmas in toroidal devices with magnetic confinement. However
voltage limits on the ICRF antenna used to launch the waves sometimes lead to a limitation of
the power. These limits are related to a variety of high voltage breakdown phenomena in the
presence of plasma that depend, in particular, on spatial charge effects and particle fluxes to
the electrodes.
An ICRF probe has been developed to study the high voltage phenomena. The open end of
a coaxial line models the high voltage region of the antenna. The voltage limits were studied in
well defined conditions in a test facility without magnetic field and in the real conditions of the
peripheral plasma of the ASDEX Upgrade divertor tokamak.
The ICRF probe was installed in the test facility and conditioned in vacuum by high power
pulses to reliable operation with 60 kV, 200 ms or 80 kV, 20 ms pulses. During the conditioning,
vacuum arcs occur mainly at the probe head. The arcs appear often when dark field emission
currents are measured. The presence of a plasma density of 10 15 m−3 (delivered by a high
aperture ion source) does not affect the voltage stand-off of the probe unless the pressure of
working gas is increased beyond a critical level: a semi-self-sustained glow discharge is ignited
at a pressure of 0.15 Pa for He and 0.03 Pa for air. These pressures are about one order of
magnitude lower than the pressures required for ignition of a self-sustained glow discharge at 80
kV. Cathode spots on the surface of the inner conductor are formed in the semi-self-sustained
discharge and often lead to the formation of the arc discharge.
When the ICRF probe is installed in ASDEX Upgrade and is well conditioned (to the
maximal voltages achieved in the test facility), high voltage breakdown on the probe often
correlates with activity of edge localized modes (ELMs). The breakdown characteristics are
similar to that of the cathode spots formation in the semi-self-sustained discharge glow discharge.
The maximal RF voltage on the ICRF probe increases from shot to shot, i.e. an additional
conditioning effect is observed during plasma operation. The voltage limit of the probe can be
increased by application of a positive DC bias to the inner conductor while at the same time
the rectified current associated with the collection of ions across magnetic field is suppressed.
It was found that the appearance of ELMs and other intermittent events in the scrape-off-layer
(SOL) plasma in the region of the probe head lead to a local dissipation of a high fraction of
RF power.
The role of ELMs as RF breakdown trigger is confirmed by observations during operation
of the full-size AUG ICRF antenna. A reliable arc detection system is required for the ICRF
antennas (not every breakdown triggered by ELMs is easy to detect), otherwise the overall
performance of the antennas degrades due to appearance of quasi-stationary arc discharges.
The antennas operates more reliably when the antenna conductors are conditioned with plasma.
Measures to improve the antenna voltage stand-off in the presence of plasma are suggested:
an optically closed Faraday screen; glow discharge conditioning; a form of antenna conductors
to minimize ion collection across the magnetic field and minimize asymmetry of electrodes along
the field; neutral density reduction inside the antenna. Further work should be focused on the
choice of the antenna materials, parasitic absorption of the RF power and the antenna-plasma
interaction for different DC boundary conditions of the antenna circuit.
Contents
1 Introduction 1
1.1 Fusion research . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2 Fusion of energetic particles . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.3 Tokamak concept and plasma heating . . . . . . . . . . . . . . . . . . . . . 5
1.4 Heating of plasma with ICRF . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.5 Power limitations of ICRF antennas . . . . . . . . . . . . . . . . . . . . . . 11
1.6 Outline of the thesis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2 Phenomenology of RF breakdown 14
2.1 Main parameters . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.2 Power transmission . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.3 Breakdown development on the ICRF antenna . . . . . . . . . . . . . . . . 18
2.4 Gas discharge phenomenology . . . . . . . . . . . . . . . . . . . . . . . . . 20
2.4.1 DC discharges . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.4.2 RF discharges . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
2.4.3 RF discharges responsible for voltage limitation . . . . . . . . . . . 26
2.5 RF vacuum arc ignition . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
2.5.1 Field emission . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
2.5.2 Conditioning by high voltage . . . . . . . . . . . . . . . . . . . . . 28
2.5.3 Spark stage of RF vacuum breakdown . . . . . . . . . . . . . . . . 29
2.6 Charge particles in electrode gap in vacuum . . . . . . . . . . . . . . . . . 30
2.6.1 Particle motion in vacuum at high RF voltage . . . . . . . . . . . . 31
2.6.2 Particle flux focusing on the microscale . . . . . . . . . . . . . . . . 33
2.6.3 Thermal desorption and skin-effect . . . . . . . . . . . . . . . . . . 35
2.6.4 Particle stimulated desorption . . . . . . . . . . . . . . . . . . . . . 37
2.6.5 Secondary emission processes . . . . . . . . . . . . . . . . . . . . . 37
2.6.6 Multipactor in vacuum . . . . . . . . . . . . . . . . . . . . . . . . . 38
2.6.7 Mean free-pass and cross-sections of ionization processes . . . . . . 38
2.7 Self-sustained RF glow discharge . . . . . . . . . . . . . . . . . . . . . . . 39
2.7.1 Role of inductively coupled discharge . . . . . . . . . . . . . . . . . 39
2.7.2 Capacitively coupled discharge . . . . . . . . . . . . . . . . . . . . . 39
2.7.3 Multipactor plasma discharge (multipactor affected by gas) . . . . . 41
2.7.4 Pressure hysteresis for RF discharge existence . . . . . . . . . . . . 42
2.7.5 RF gas discharge conditioning . . . . . . . . . . . . . . . . . . . . . 43
2.8 Summary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
iii
3 Plasma in the electrode gap 45
3.1 Approach to a DC sheath . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
3.2 A RF sheath: frequency ranges . . . . . . . . . . . . . . . . . . . . . . . . 48
3.2.1 Comparing ω0 with ωpe . . . . . . . . . . . . . . . . . . . . . . . . . 48
3.2.2 Comparing ω0 with ωpi . . . . . . . . . . . . . . . . . . . . . . . . . 50
3.3 Plasma screening properties . . . . . . . . . . . . . . . . . . . . . . . . . . 51
3.3.1 Electrical field for the thin sheath (s < d) . . . . . . . . . . . . . . 52
3.3.2 Basic dynamics of the thin high-voltage RF sheath (s < d) . . . . . 55
3.3.3 Basic dynamics of the thick high-voltage RF sheath (s > d) . . . . 62
3.3.4 Surface electrical field and a transition to s ∼ d . . . . . . . . . . . 65
3.3.5 Role of ponderomotive force for density reduction . . . . . . . . . . 66
3.4 Influence of a magnetic field . . . . . . . . . . . . . . . . . . . . . . . . . . 67
3.4.1 Confinement of particles in the electrode gap . . . . . . . . . . . . . 67
3.4.2 Charging of the plasma in the magnetic field . . . . . . . . . . . . . 68
3.4.3 Multipactor conditions . . . . . . . . . . . . . . . . . . . . . . . . . 69
3.4.4 Effect on the effective interelectrode distance . . . . . . . . . . . . . 70
3.5 Summary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
4 Experimental approach 72
4.1 Concept of the experiment . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
4.1.1 RF and DC power generators . . . . . . . . . . . . . . . . . . . . . 75
4.2 Experimental device . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
4.3 Setup of the experiment . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
4.3.1 Setup in the test facility . . . . . . . . . . . . . . . . . . . . . . . . 82
4.3.2 Setup in ASDEX Upgrade . . . . . . . . . . . . . . . . . . . . . . . 83
4.4 Measurements . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
4.4.1 RF measurements . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
4.4.2 DC measurements . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
A Abbreviations 132
Introduction
The reduction or complete substitution of energy produced from fossil fuels needs reli-
able alternatives. Nuclear fusion, in particular, thermonuclear fusion, is one of the most
promising candidates despite the fact that further research and development is needed.
1
2 CHAPTER 1. INTRODUCTION
The first reaction product of (1.1) He4 (α-particle) is completely inert. Furthermore
the energy of the α-particles can be used to sustain fusion reaction (see next section). The
second product of (1.1) are fast neutrons which are collected by a blanket around the
plasma [2, 3]. Neutrons heat the blanket and the heat can be converted into electricity.
However the first wall of the fusion device can be activated by energetic neutrons. To
minimize the negative influence of neutrons, low-activation materials should be used in
the blanket. There exists a positive aspect of the neutron production: neutrons from
the D-T reaction can be used for the production of T in a tritium breeding blanket by
reacting with Li (lithium):
D and Li can be isolated from seawater in significant amounts. Therefore there would
be no problem to get easily accessible and cheap resources. Also for environment there
is no CO2 production and no danger of an uncontrolled reactor. The products of fusion
reactions are not radioactive. The radioactive waste can come from the activation of the
walls of a fusion reactor but is much smaller than that for a fission reactor.
1. energy of the reacting particles, should be high to overcome the Couloumb barrier
in thermonuclear fusion, particles get energy by increasing of temperature of plasma
(at the required temperatures matter is in plasma state);
3. energy confinement time τE (introduced by Lawson [1] for a pulsed fusion device),
a characteristic timescale of the loss of energy carried by the plasma particles after
external energy sources are switched off.
Pf us
Qf us = (1.2)
Pext
1.2. FUSION OF ENERGETIC PARTICLES 3
The condition:
Qf us = 1 (1.3)
is defined as the so-called ”break-even” – when the power released during the fusion
reactions in plasma becomes equal to the power delivered to plasma by external heating
systems.
To get an idea about values of Qf us required for total production of power we introduce
an efficiency of energy conversion from fusion energy to electric energy η thermal and an
efficiency of conversion from electric energy to energy delivered by heating systems η ext .
The total efficiency of conversion of the fusion power to electricity η total is given by:
1
η total =η thermal − (1.4)
Qf us · η ext
The efficiency η thermal has typically values of 0.3 to 0.4. Value of η ext depends on the type
of the fusion device and heating systems used. For equality of the energy produced in the
fusion device that is converted into electricity and the energy consumed by the heating
systems (η total = 0), Qf us ≈ 10 is required for systems with magnetic confinement while
Qf us ≈ 100 is required for systems with inertial confinement and laser heating.
Only 20% of the fusion energy released from a single fusion reaction goes to the
energy of α-particles α while 80% of the energy is carried by neutrons (e.g. for Qf us = 1
α-particle power Pα,f us constitutes 20% of the heating power Pext ). The concept of a
steady-state fusion device deals with the so-called ignition or burn condition. This means
that the fusion reaction is self-sustained, i.e. the required temperature of fuel is sustained
solely by α-particles heating. The power balance in the steady-state fusion device can be
written as:
Pα,f us + Pext = Ploss (1.5)
where Ploss is the power of losses from the plasma characterized by the empirical energy
confinement time τE . For the burn condition Pext = 0 (or Qf us = ∞) and:
Pα,f us > Ploss (1.6)
For a 50% mixture of deuterium and tritium one gets the condition:
1 2
n hσf us vi α V > 3ni Ti V /τE (1.7)
4 i
where ni is the density of the deuterium/tritium mixture, Ti – temperature, V – volume
where the reaction takes place, α – energy of α-particle produced by the fuion reaction.
The expression is transformed to:
12 Ti
ni τ E > (1.8)
hσf us vi α
Within 10% accuracy one can represent the reaction rate in the temperature range 10-20
keV as: hσf us vi ≈ 1.1 · 10−24 · Ti2 m3 s−1 (Ti in keV). Using α = 3.5 MeV one gets an
approximate condition for ignition:
ni Ti τE > 3 · 1021 m−3 keV s (1.9)
4 CHAPTER 1. INTRODUCTION
Getting condition (1.9) we assumed that temperature and density are the same for the
whole volume where the fusion reaction takes place.
Even if ignition is not achieved, the conditions for total power production may be
achieved when Qf us is sufficiently high. The total energy conversion efficiency (1.4) should
be as close to the efficiency of fusion energy conversion η thermal as possible. There are
different experimental approaches which are distinguished by types of particle confinement
and heating systems.
2. Laser fusion and ion fusion. In these systems [5, 6] the conditions for the thermonu-
clear burn are achieved by bombarding a target (a D − T pellet) by a laser beam
or beam of ions. Fueling is realized via repetitive usage of many targets, i.e. this
type of systems are pulsed. The target is a capsule with frozen or liquid deuterium-
tritium. The target should be symmetrically compressed by implosion process to
the conditions of high energy and temperature and very high density (≈1000 times
the density of liquid D −T ). The confinement times are of the order of nanoseconds.
This leads to the high power requirements of the ”driver” (beam) systems which
should deliver more than 1 MJ of energy during the confinement time. Because
of low efficiency of driver systems, in particular laser heating systems: η ext < 0.1,
the so-called target gain (that is actually an analog of Qf us in terms of energy)
is required to be higher than 100 to achieve positive energy gain. Use of particle
accelerators to heat the target can improve the efficiency: η ext = 0.2 − 0.35. There
are two approaches to optimize target gain: a direct drive - where a target is di-
rectly bombarded by laser beam or ion beam and heated; an indirect drive - when
the energy from beams is first absorbed by a high-Z enclosure and up to 80% of
it is converted to the energy of X-rays which drive the implosion. During the last
decade a large progress has been made in understanding the underlying physics of
inertial confinement and targets with high gain have been designed. Experiments
are planned, in particular at the National Ignition Facility, U.S.A., to reach ignition
and fusion power Pf us >100 MW during a quasi-steady-state operation (repetition
frequency of ≈ 10 Hz).
1.3. TOKAMAK CONCEPT AND PLASMA HEATING 5
OH-Transformer
Vertical field
coils
Particle motion
The basic principles of particle movement are shown in Fig. 1.2. The charged particles in
a magnetic field move along the magnetic field lines and have a circular orbit in the plane
perpendicular to the field lines. The radius of the orbit is called gyroradius.
The toroidal magnetic field created by external (poloidal) coils is non-uniform and is
decreasing from the center of the torus to the outer wall (low-field side). Therefore if only
the toroidal field is present, a vertical drift of particles appears having opposite directions
for different charge polarity. The separated charges produce a vertical electric field which
results in an additional drift of the plasma radially towards the low field side. The scheme
of this charge separation and the plasma drift is illustrated in Fig 1.2a.
To avoid the radial plasma drift, a poloidal component of the magnetic field is created
by a toroidal plasma current (Fig. 1.2b). The magnetic lines connect top and the bottom
sides of the torus. Since the plasma has high electrical conductivity along the magnetic
field lines, the regions of different excess charges are short circuited, thus the vertical
electric field is reduced and the radial drift is decreased.
A particle moving in a non-uniform magnetic field changes the distribution of its
energy between the velocity perpendicular and velocity parallel to the magnetic field
depending on the magnitude of the field. The parallel velocity decreases when particle
goes to the high field side. Therefore particles having high initial parallel velocity may
pass the high field side and circulate continually around the torus. These are called
passing. The particle with initially low parallel velocity are reflected to the low field side
and are trapped. The trajectories of guiding centers (centers of gyromotion) of a passing
and a trapped particles are shown in Fig. 1.2b. The diagrams for the movement projected
into the poloidal plane are shown in the right part of the figure. The trapped particles
1.3. TOKAMAK CONCEPT AND PLASMA HEATING 7
a) b)
Trapped particle
B ter
cen
ing
guid orbit Passing particle
→ ter flux
B cen
ing surface
→ guid orbit
B
+ + +
→ → →
E E B Magnetic
+
- - - flux surfaces
Radial drift
of plasma
→
B + Ip - toroidal
plasma current
Figure 1.2: Main particle motion in a tokamak: a) only toroidal magnetic field, b) in the
presence of toroidal and poloidal magnetic fields.
move along the so-called banana orbit. The presented picture of particle movement is
realized when collisions are rather rare, so passing particles make more than one toroidal
turn and trapped particles execute more than one banana orbit before they collide with
other species.
Heating systems
For a fusion reactor, in particular based on the tokamak principle, the heating systems [2,
3] are important for:
In most present experiments the achieved fusion power is negligible compared to the
heating power. Therefore it is important to simulate the conditions at which fusion
reactors will operate. Use of the heating systems in the fusion experiments allows to:
- simulate α-particle heating to get a valuable experience for operation of the future
fusion reactors;
- test systems for stabilization of MHD instabilities, current drive and avoidance of
impurity accumulation in the central plasma.
The plasma heating in magnetically confined plasmas is realized by using heating systems
which may be classified as follows:
1. ohmic heating
joule heating by the induced toroidal current in tokamaks; usually not sufficient to
reach ignition conditions;
3. heating by launching waves in different frequency ranges into plasma (no fuelling of
the plasma), the most important ones are:
In fusion reactors working in the conditions when the power of the external heating
amounts to a small fraction of produced fusion power, it will unfortunately not be possible
to provide a central plasma fuelling. Therefore for a correct simulation of the reactor
operation a heating system without particle input should be used. α-particles in the
reactor heat mainly electrons. Since α-particle heating power density is proportional
to the fusion reaction rate hσf us vi, and therefore to Ti2 , the heating profiles are peaked.
Thus for the correct simulation of α-particle heating a central electron heating is required,
e.g. electron-cyclotron or ICRF heating. In addition to electron heating, ICRF heating
produces energetic ions which firstly, can help to reach ignition parameters, and secondly
allow to model existence of fast α-particles and α-particle heating simultaneously.
Let us describe basic principles and problems of ICRF heating.
where Zi is ion charge number, e is the elementary charge, mi is ion mass. The frequency
range used for operation in the magnetic fields typical for the fusion devices with magnetic
confinement is typicaally from 20 (fundamental frequency of H at 1.3 Tesla) to 120 MHz
(second harmonic of H at 4 T). The vacuum wavelength of RF waves at these frequencies
is λ = 15 m and λ = 2.5 m respectively.
Principle of heating
The heating principle relies on the excitation of the waves at the plasma edge that propa-
gate into the plasma [3, 15]. The power of the waves is absorbed inside the plasma. Two
types of waves can be used for the heating [15]:
Heating by fast waves is used often in the so-called minority heating scenario. Heating
of ions requires left-hand polarization of the electric field. When one has a plasma of one
sort of ions, the fast wave at the fundamental frequency is right hand polarized and the
wave is poorly absorbed in the ion cyclotron resonance layer. To overcome the problem of
polarization, it is possible to fill in a small fraction of another sort of ions (minority ions)
and use fundamental frequency for this sort of ions. The polarization is mainly determined
by majority ions while at the resonance layer for minority ions the wave contains a right
hand polarized component and the minority ions absorb energy efficiently.
There are also other heating scenarios involving, e.g. a heating by fast waves on the
second harmonic or a conversion of the fast wave into ion-Bernstein wave (an electrostatic
wave that is efficiently damped on electrons) in the region of ion-ion hybrid resonance
when the concentration of the second sort of ions is relatively high. The minority heating
(e.g. by minority hydrogen in deuterium plasmas) has become a routinely used ICRF
heating scenario due to its relative simplicity and reliability.
Typically the fast wave is launched by an antenna that situated at the low field side of
a tokamak (see Fig. 1.3 showing ICRF antenna in ASDEX Upgrade). The electromagnetic
fields corresponding to the fast wave field are excited by one or many poloidal straps. Use
of a Faraday screen with rods oriented along the external magnetic field reduces electric
field of the waves along the magnetic field and therefore reduces the excitation of slow
waves.
10 CHAPTER 1. INTRODUCTION
power input
edge plasma
0.28
a
scrape-off layer
plasm
core
plasma R
edge
−0.84
divertor tiles
antenna straps
Faraday screen
(inside the antenna box)
−1.40
0.80 1.16 1.52 1.88 2.24 2.60
Radial position R, [m]
The geometry of the ASDEX Upgrade experiment and the ICRF antennas is shown in
Fig. 1.3. The antenna includes two current-carrying straps (the right part of the figure
shows that each strap consists of two loops) made of stainless steel and separated by a
septum. The straps are mounted in an antenna box with a back plate playing a role of
the RF return conductor (behind the straps in Fig. 1.3). Limiters (plasma facing graphite
tiles, see the antenna limiter contour in the figure) protect the stainless steel and copper
components from the plasma. A Faraday screen is constructed from molybdenum rods
typically with TiC coating on the rods.
The strap is grounded and has a length from λ/8 to λ/4 depending on the frequency
of operation. This means that a spatial distribution of RF voltage exists along the strap.
For this reason a range of the voltages is present on the antenna strap, starting from
zero at the position of the ground to a maximum at the side where it is connected to the
transmission line.
The plasma geometry in Fig. 1.3 corresponds to AUG discharges of a special plasma
shape, the so-called edge optimized configuration (EOC) characterized by relatively small
distance (from 2.5 to 4 cm) between the separatrix and the antenna limiter along the
antenna contour. This plasma shape fits the antenna contour relatively well. Many
plasma discharges of ASDEX Upgrade which were performed during this work used the
EOC plasma shape.
1.5. POWER LIMITATIONS OF ICRF ANTENNAS 11
R0 2
P = V (1.11)
2Z02
where Z0 is the characteristic impedance of a transmission line, R - the coupling
resistance, V - the peak value of the maximum voltage on the transmission line. One can
see that for a given R0 , an increase of the coupled power should be made by an increase of
voltage, furthermore the dependence on the voltage is quadratic. For example, for ASDEX
Upgrade antennas the resistance of 2 Ω and delivery of 2 MW power to plasma correspond
to a voltage on the transmission line of 25 kV. Since for typical frequencies used for ICRF
heating the antenna strap length is somewhat larger than λ/4, the corresponding voltage
on the antenna straps is a few kV less than 25 kV. The voltage on the ASDEX Upgrade
antenna is presently limited for reliable operation to about 20-25 kV. Higher voltages
should allow higher power and thus higher power densities.
Therefore voltage stand off, i.e. the maximum voltage sustained without a breakdown,
is an important issue to increase the power launched on the antenna and/or to decrease
an area occupied by the ICRF launchers in the fusion devices with magnetic confinement.
In experiments it is often observed that the maximal voltage on the antenna is affected
by the presence of plasma. The plasma conditions in the region of the ICRF antenna in
the fusion device are very specific due to complex stationary and transient processes which
take place at the plasma boundary.
As the confinement of plasma is not ideal, particles and heat are transported across the
magnetic field. The total transport in the devices with magnetic confinement has three
contributions [3]:
(a) low collisionality banana regime, the particles move as described in 1.3, the
trapped particles are displaced by the width of the banana orbit due to a colli-
sion and contribute to transport dominantly, diffusion grows with collisionality;
(b) highly collisional regime, transport is defined by resistive plasma diffusion in
torus, the collisions are so frequent that the trapping of particles is negligible,
diffusion grows with collisionality;
(c) plateau regime, a transition region from low to high collisionality regimes, is
characterized by an approximately constant diffusion coefficient;
2. turbulent, is associated with the appearance of instabilities, not yet well understood;
3. ripple losses, the losses caused by particles toroidally trapped in a magnetic field
ripple between toroidal field coils.
The neo-classical transport is higher than transport described by the classical diffusion
- a diffusion across the magnetic field in a cylindrical plasma where diffusion coefficient
is defined by collision time and gyroradius (the latter is the appropriate scale of particle
displacement after a collision in the classical case). However the measured diffusion coeffi-
cients are usually even higher than neo-classical, mainly due to the anomalous transport.
In the tokamak different instabilities may be observed which lead to the anomalous trans-
port and are described by magneto-hydrodynamic and kinetic theory. The instabilities
lead to the changes of effective transport coefficients in the core of plasma and at the
plasma edge leading to losses of particles and heat from the region of confined plasma.
Thus the edge instabilities affect the plasma density profile. Since the plasma density
profile is an important parameter for coupling of waves by the ICRF antenna, the insta-
bilities affect the efficiency of launching the waves into the plasma. The instabilities may
produce plasma formations (blobs) [16, 17] which can penetrate fast into the tokamak
scrape-off-layer and to the antenna. One of the experimentally observed instabilities is
the so-called edge localized mode (ELM, see section 5.3 of this thesis for a more detailed
description). When this instability takes place, the plasma boundary is distorted in a
way that is analogous to the distortion during formation of plasma protuberances on the
sun. Strong turbulent particle fluxes, electric and magnetic fields interact and dense, hot
plasma formations can be ejected from the main plasma boundary.
Therefore, due to uniqueness and complexity of the processes occurring at the plasma
edge of tokamak, it is impossible to reproduce the real plasma conditions in a small
laboratory experiment. For this reason the real plasma conditions of the plasma boundary
of the fusion device with magnetic confinement (in our case - tokamak ASDEX Upgrade)
are essential for the experimental studies of high voltage breakdown phenomena on the
ICRF antennas.
1.6. OUTLINE OF THE THESIS 13
Phenomenology of RF breakdown on
the ICRF antenna
Arcing in a fusion device with magnetic confinement is usually associated with excitation
of arcs on a plasma-wall boundary [18, 19, 20]. These arcs (called unipolar arcs) lead to
an increase of erosion rate of the first wall and a release of impurities to the edge plasma.
We are interested in arcing that occurs on the ICRF antennas when an external RF
power is applied. The arcing can be considered as bipolar arcs. However arcing is a final
stage of a developed high voltage breakdown. Let us start with a review of the possible
processes that may be responsible for the voltage limits of the ICRF antenna.
There are several types of breakdown which can be viewed as candidates leading to
the voltage limitation on the ICRF antenna:
- high-voltage breakdown in vacuum or vacuum arc;
- Townsend breakdown of neutral gas leading to formation of a glow discharge [21];
- multipactor discharge which can be subdivided in two types: pure electron resonance
discharge [21, 22] and multipactor discharge with gas ionization [22].
Multipactor is possible only if RF voltage is applied to the electrodes. The basic mech-
anisms leading to the development of a vacuum arc and a glow discharge allow them to
appear both with RF and DC voltages applied.
In real conditions of interaction of a RF antenna with the peripheral plasma of a fusion
device with magnetic confinement two possible approaches exist:
1. First-order approach. Effects are defined by local physical conditions. Plasma edge
parameters are considered as defined by external plasma inflow and its interaction
with electrode surfaces, RF voltage and RF current on ICRF antenna. The inter-
action does not the affect the parameters of the central plasma.
2. Second-order approach. The ICRF power affects globally the central plasma. The
plasma edge parameters are affected by the changes in the central plasma. This
results in an additional non-linear dependence on the ICRF power.
In the framework of this thesis the second-order approach is not considered.
14
2.1. MAIN PARAMETERS 15
Local plasma parameters. Presence of a plasma near the antenna is essential. The
plasma can screen external electrical field by forming sheaths of space charge near
the electrodes. The space charge effect may lead to an increased local electric field on
the electrode surface compared to the electric field in vacuum. The plasma potential
should be taken into account for the description of the electric fields. The electron
temperature is an important parameter defining ion flows from quasineutral plasma.
The parameters are connected self-consistently.
RF voltage/current on the antenna. The ICRF antenna operates with high currents
and high voltages. High RF voltages are mainly discussed in the work.
DC boundary condition for the antenna. The fact whether the system is capable
to conduct DC current or not defines the boundary condition which affects plasma
parameters. The reason for this is sheath rectification and self-biasing effects which
appear when RF voltage is applied to an electrode in plasma.
Magnetic field. The presence of a magnetic field affects almost every parameter in the
list: local plasma parameters, effective geometry of the system of electrodes, sheath
effects, multipactor conditions.
Electrode material and surface properties. Surface processes play a decisive role in
the ignition of a vacuum arc and the startup of glow and arc discharges.
VRF
Z0 Z0 R0
antenna
RF plasma
generator transmission matching resonant
line circuit line
DC-break
2
Vmax
Pmax = (2.1)
2Z0
where Vmax is voltage stand-off of the transmission line. The typical values of the charac-
teristic impedance of the transmission lines vary from 20 Ω to 75 Ω and differ from typical
plasma R values. For this reason a wave carrying energy from a RF generator to the load
(the forward wave) is partly reflected. A matching circuit represents an inhomogeneity
which produces an additional reflection wave towards the RF generator with the same
amplitude but out of phase with the wave reflected from the load. The two reflected
waves interfere and in the ideal case only the forward wave exists between the matching
circuit and the RF generator. Due to the presence of both forward and reflected waves a
standing wave is formed (see Fig. 2.1) between the antenna and the matching circuit. A
2.2. POWER TRANSMISSION 17
standing wave is characterized by voltage standing wave ratio higher than unity:
Vmax
k= >1 (2.2)
Vmin
.
It is convenient to use the relations for the generalized voltage reflection coefficient ρ.
According to the definition this complex coefficient is the ratio between the phasor value
of voltage of the forward wave and the phasor value of voltage of the reflected wave. It
can be written as:
Z − Z0
ρ= (2.3)
Z + Z0
where Z = V /I is a complex impedance along a line (V and I the phasor values of voltage
and current) and Z0 the characteristic impedance of the line. The value of the impedance
Z at any point of the line is often referred as input impedance.
If the RF system is matched perfectly the input impedance at the end of the matching
circuit connected to the transmission line equals exactly to characteristic impedance Z 0
of the transmission line. This case is often referred as ”traveling wave mode” (k = 1)
for the transmission line. In reality the effective impedance R0 of antenna-plasma system
varies in time and no perfect matching is instantaneously obtainable. For this reason a
standing wave is often formed between the matching circuit and the RF generator.
As a result of the mismatch higher voltages along the lines exist compared to the
”traveling wave mode”.
Quality of resonance is a measure which helps to understand how the resonant (un-
matched) part of RF system reacts on the load changes. It is defined as:
Wstored
Q= one cycle
(2.4)
Wdissipation
The Q-value determines how much energy is dissipated during one RF oscillation com-
pared to the amount of reactive energy stored in the system. It is important to distinguish
between Qno load without load connected and Qload for the resonant line with load. Applied
to the ICRF system, the former notion is used for operation in vacuum with Wdissipation
defined by inherent RF losses while the latter corresponds to operation with high tem-
perature plasma in the toroidal chamber. The Q-value without load connected can be
estimated from:
βp
Qno load = (2.5)
2αat
where βp – phase factor and αat – attenuation factor of a wave. For a coaxial transmission
line with the radii a and b of the inner and the outer electrodes respectively one has an
expression for the attenuation factor defined by electrical resistance of the conductors:
Rsq 1 1
αat = + (2.6)
4πZ0 a b
where Rsq = 1/σδ – resistance of a square meter of the surface defined by the electrical
conductivity σ and the skin depth δ.
18 CHAPTER 2. PHENOMENOLOGY OF RF BREAKDOWN
- antenna electrodes
this work is mainly concentrated on studying the antenna voltage stand-off since the
influence of plasma on the voltage limit is often essential.
With a higher level of the reflected power the generator reacts on the external load by
altering its output RF power from the requested values and often by starting to produce
modulation of the power on low frequencies. Reflected power at the generator is decreased
by:
2.3. BREAKDOWN DEVELOPMENT ON THE ICRF ANTENNA 19
matching
circuit
RF Pfeed
ρb
generator
x
x=0 x = xb
feeding breakdown
point location
- any ”advanced” ICRF system is designed in a way to reduce reflected power at the
generator.
This helps to obtain the value of the input impedance on the generator side close to 50 Ω.
For the ASDEX Upgrade ICRF system, a 3 dB hybrid power splitter on the antenna
side [28, 29, 30] is used (see section 5.6). The reflected power at the generator side in
this system is mostly dependent on the asymmetry of the reflected power from the two
antennas. If the reflection from the two antennas is approximately the same, the reflected
power at the generator is not influenced by the absolute value of the reflected power from
a single antenna. Online matching circuits, such as Fast Ferrite Tuners (FFT) [31, 32]
can be used to match a change in the load impedance during RF pulse as well. Thus the
generator can sometimes continue to deliver power after the breakdown on the antenna
has evolved into arcing.
Now we assume that the generator proceeds to deliver the same power after the break-
down. For simplicity the circuit shown on the Fig. 2.2 is considered. The type of the
breakdown defines the initial value of reflection coefficient ρb at the load. According to
transmission line theory an amplitude voltage-current (V-I) characteristic of an ”external
power source” at the breakdown location xb is set by the formulas:
v
u
u 2 · P f eed (1 − |ρ f eed | 2 ) exp(−γx ) − ρ
b f eed exp(γx b )
Ib = |Ix=xb | = t ·
(2.8)
1+ρf eed
Re Z0 · 1−ρf eed 1 − ρ f eed
1 + ρ
Ub = |Ux=xb | = Ix=xb · Z0 ·
b (2.9)
1 − ρb
ρf eed = 0.5 · ρb · exp(γxb ) + ρmatch (2.10)
20 CHAPTER 2. PHENOMENOLOGY OF RF BREAKDOWN
where Pf eed , ρf eed – forward power and reflection coefficient at the feeding point, ρmatch
– reflection coefficient of the matching circuit at the feeding point, γ = iβ + α. The case
ρf eed = 0 represents a line perfectly matched for the load characterized by ρb .
The relations (2.8)–(2.10) for Ib and Ub should be satisfied as as well as a relation
between Ib and Ub set by local physical conditions. External circuit and physical conditions
at the breakdown location set the actual scenario for the development of a discharge. The
amplitude V-I characteristic of the external circuit can limit the current at the load as the
system is not specially matched for the loads typical for gas discharges. It may influence
the discharge load impedance since the load (defined by gas discharge properties) depends
on the power applied non-linearly. However equations (2.8)–(2.10) do not describe the
reaction of the external circuit on transient (very fast) changes of the load. In this case
the stored (circulating) RF power is available to be dissipated in a discharge. Therefore
the stored power characterized by the Q-factor may influence the discharge evolution
(see below in section 2.4.2). Thus the breakdown can lead either to a steady-state gas
discharge (the existence of which is defined by external circuit with equations (2.8)–(2.10)
or to a pulsed discharge (a discharge with fast termination the external circuit being
characterized by the stored reactive energy).
The load impedance of gas discharge is one of the most interesting issues discussed in
gas discharge physics. Taking in account that load impedance is a ratio between current
and voltage on electrodes it is important to discuss the phenomena which may appear in
an electrode gap for different neutral gas pressures.
V
Rbal Imax
N
Vc
F
M
B C
Vt
D E G
O H P
A
K
I max I
Figure 2.3: Stationary voltage-current characteristics of the electrode gap for DC voltage.
Black curve - vacuum, p < 0.1 Pa; gray curve - low-pressure, p ∼ 0.1 ÷ 100 Pa [33, 21].
Dashed gray curve - load curve defined by the external DC circuit.
2.4.1 DC discharges
A qualitative picture summarizing the discharge phenomena with DC voltage between two
parallel plates can be presented with the help of DC V-I characteristic of the electrode
system [21, 33, 34]. In Fig. 2.3, the stationary characteristics for DC voltage are shown.
Each point on the solid curves is a stationary point, i.e. the change along the curves
happens with a characteristic time much longer than time of the dissipation of energy
stored in reactive elements of the power circuit.
The black curve in Fig. 2.3 corresponds to the vacuum case or to pressures p < 0.1
Pa. The gray curve represents the characteristics of self-sustained discharges at the low-
pressure case which is defined:
glow discharge.
An operational point on the stationary V-I curve is chosen by the so-called load curve
(dashed gray line in Fig. 2.3):
- the whole glow discharges collapses into a current channel (the so-called contraction)
leading to high current densities and an arc.
2.4. GAS DISCHARGE PHENOMENOLOGY 23
~
V
~
Vc N
F
M
~ C E
Vt D
B
L G
H
A O P
K ~
I
Figure 2.4: Stationary voltage-current characteristics of the electrode gap for RF voltage.
Black curve - vacuum, P < 0.1 Pa; gray curve - low-pressure, P & 0.1 Pa. Constructed
from [22, 35, 36].
The arc discharge (segment GH) represents a respawning group of the cathode spots
having certain properties depending mainly on the cathode material and the maximal
current the circuit can deliver.
2.4.2 RF discharges
In Fig. 2.4 qualitative stationary amplitude V-I characteristics for the two-electrode sys-
tem with RF voltages applied are shown. The black curve corresponds to the vacuum
case P < 0.1 Pa. The characteristic shown by the gray curve is constructed based on the
RF discharge phenomena at low-pressures from the noted range of interest (the condi-
tion (2.11)), e.g. P ∼ 10 Pa for He. The experimental results from [35, 36] are reflected
in the curve qualitatively. At the low pressures the glow discharge can be created at a
high voltage which is of interest for us since we are looking for the voltage limits in the
low-pressure range.
The real discharge in the electrode gap is dependent on the system details like geom-
etry, materials of the electrodes, gas and may not include every segment from Fig. 2.4
or may include some features (see below). However we try to summarize the effects ob-
servable in the double electrode RF system at low pressures. As we already know, a real
stationary amplitude V-I characteristic is defined by a load curve derived from equations
(2.8)–(2.10). In comparison to the DC case the load curve is not linear and depends on
the load impedance. Moreover the curve contains more parameters of the system.
Let us consider the vacuum case (black curve, Fig. 2.4). With increasing the RF
voltage on the electrodes a small RF current is registered proportional to the voltage
because of a non-zero capacitance between the electrodes. The possibility to conduct
24 CHAPTER 2. PHENOMENOLOGY OF RF BREAKDOWN
capacitive current through the electrode gap is one of the principal distinction of an RF
discharge to the DC discharges to be discussed later.
At a RF voltage typically less than 1 kV, multipactor ([22], section 2.6.6) appearance
is probable. The effect is observed as an increase of RF current in a very narrow voltage
range (L) and may be a purely electron discharge or a gas amplified discharge. At high
voltages (segment MN) field emission currents become detectable and at some critical
voltage Ṽc a vacuum arc similar to the DC vacuum arc is ignited (segment NO). Segment
OP represents RF-sustained arc discharge on the electrodes.
The gray curve in Fig. 2.4 shows the behavior of the electrode system at the pressure
when gas ionization starts playing a role.
When the RF voltage is increased, the RF current increases proportionally untill the
voltage satisfies resonant conditions for multipactor which are influenced by gas ([22],
section 2.7.3). Ionization of neutral particles by electron impact assists multipactor and
makes an additional source of electrons. Gas in the electrode gap is partly ionized and
the discharge is not purely electron. Gas desorption from the electrodes under condi-
tions of electron bombardment can lead to an increase of the current (and voltage, see
section 2.7.3) illustrated by the dashed gray curve in Fig. 2.4. At higher pressure than
showed in Fig. 2.4, the voltage required to have conditions for the multipactor effect is
of the order of the voltage Ṽt to start an RF glow discharge. Therefore ignition of the
RF glow discharge (CD) (paragraph 2.7) at higher pressures may be accompanied by
multipactor.
At the pressures P < 100 Pa the glow discharge may exist in two forms: α- (DE)
and γ- (EF) discharges. In the α-discharge the main fraction of the RF current on the
electrodes has a displacement nature due to formation of sheaths near the electrodes. The
α-discharge is not dependent on the material of the electrodes and is sustained only by
the ionization of gas in the gap. The γ-regime of the discharge depends on processes of
recycling of particles on the electrode surfaces and therefore on the electrode material. The
transition from α to γ discharge for the range of low pressure is usually monotonical and
rather smooth for the light gasses (e.g. He [35, 36, 37]). Experiments with heavy inert
gases show that the transition form α to γ regimes may be accompanied by a hysteresis
on the I-V characteristic (e.g. for Ar [38]). A more detailed discussion about the features
of α and γ discharges will follow in section 2.7.
~ a) ~ b)
V V
Q1 < Q2 Q1 < Q2
1
free
e
rg
n-
ha
2
breakdow
disc
γ-
e
α - dis
charg 1
2
actor arc discharge arc discharge
multip
~ ~
I I
Figure 2.5: Qualitative stationary load curves for two different systems with Q > 1
matched for high impedance operation: system 1 has a lower Q than system 2: a) vacuum
operation, b) low-pressure, P & 0.1 Pa.
accurately [45]. However for the initial phase of cathode spot development we consider
pure field emission, i.e. equation (2.13) is applicable.
Taking into account the sinusoidal alternating voltage on the electrode a good ap-
proximation to equation (2.13) is made [42, 43, 44] for ERF as peak electric field on the
electrode surface:
5/2 5/2
M · βRF ERF −6.83 · 109 φ3/2
IRF = ARF · · exp (2.14)
φ7/4 βRF ERF
where M = 5.27 · 10−12 (V/m)−1/2 · (eV)3/4 is a constant derived from the integration over
the RF period. Experimentally measured dependence of field emission current is usually
fitted to equation (2.14) to derive ARF and βRF [44]. As it can be seen from (2.14), the
emitted electron current is exponentially dependent on (βRF ERF ) and φ. These quantities
are very sensitive to the surface conditions on the electrodes:
1. Surface material. In the literature different values can be found for the work func-
tions [46]. One of the reasons lies in the different measurement techniques (by field
emission, photoemission or contact potential difference). Another reason is in dif-
ferent finishing of the surface of the samples before the measurements are made
(electrolytical or mechanical polishing, baking at different temperatures, plasma
cleaning) which affect the surface purity and homogeneity (impurity inclusions, for-
mation of oxide or nitride films) as well as surface microroughness which influences
electrical field amplification. High φ is required to shift the effect of field emission
to the higher electric fields. For the ICRF antenna the material should be chosen
which do not form oxides, nitrides and carbons which have lower work function
than pure bulk material. However such metals as stainless steel, Cu, Ag, Au, N i
are often used as surface electrode materials (for pure samples φ ≈ 4 − 5 eV). A
majority of these materials have tendency to form oxide films on the surface and a
reduction of φ by 1-2 eV may be the result. For instance, W and T a as well their
oxides and carbons have the high work function of φ > 4.5 eV .
2. Surface microroughness (small-scale geometry). Small microprotrusions (sometimes
called whiskers) with characteristic sizes from 1 nm to 1 mm lead to electrical field
amplification factor of β = 10−500. These regions on the surface are called emission
sites and the main fraction of the field emission current is carried by them. Very
high local current densities (∼ 106 − 1014 A/m2 ) may be present on these sites.
Surface microroughness is characterized by an effective height of microprotrusions
Rz and written in µm. Rz is affected by the methods of electrode finishing.
3. Surface homogeneity. The boundaries between regions with different material den-
sities bring an undesirable increased diffusion of impurities from bulk towards the
surface along these boundaries and its accumulation on the surface. The growth of
microprotrusions on the surface is often observed on grain boundaries.
4. Surface purity. Contamination of the surface by impurities leads to high rates of
gas desorption during high power operation and also to appearance of the emit-
ting sites [51]. Chemical and physical bonded contaminants, e.g. oxide layers, dust
28 CHAPTER 2. PHENOMENOLOGY OF RF BREAKDOWN
the breakdown initiation. After a certain time of such operation, the surface roughness
parameter decreases and field emission currents tend to be reduced too.
Similarly an already developed cathode spot may disappear without ignition of other
cathode spots if the conditions on the surface are not favorable. In other words the
arcing can be terminated after some time even if the power applied to the electrode was
initially sufficient for the arc sustainment. For example, the long-sustained arcing cleans
the surface and the arc may vanish if it was sustained partly by the fact that impurities
were present on the surface. However, even the cleaning effect of the long-term arcing
does not compensate for the surface erosion produced by cathode spots. The erosion is
definitely unfavorable for the high voltage operation and leads to the developement of
a microstructure of the surface and at least to transient sparking, when high voltage is
applied again.
We will discuss other general aspects and methods of the ICRF conditioning in the
subsections 2.6.2, 2.6.3, 2.7 of this chapter.
- vaporization of one of the emitting sites on the cathode by Joule heating during
field emission, temperature of the sites is influenced by Nottingham [55, 45] (heat-
ing/cooling) effect;
metal
where ωpi – metal ion plasma frequency. However the value of the plasma densities of
20 −3
about 10 m may be still present near the electrode after the half of RF period when the
plasma is assumed to consist of single-charged ions of iron. On the other hand the inverse
polarity of the voltage favors further ionization of the metal vapors and heating of the
surface by electron beam, created by the field emission effect on the other electrode. This
may lead to the further amplification of the plasma density and makes the formation of an
arc on the cathode probable for the RF case. Due to the change of the voltage polarity in
a time comparable to that of the spark stage, the development of the RF spark can take
longer as the combined field-thermal emission (the cathode process) should take place to
sustain the arc.
The transition to a (quasi)-steady-state vacuum arc discharge happens when two main
conditions are realized (both for DC and RF cases):
1. external circuit delivers a current exceeding a critical value at a voltage high enough
to sustain cathode spots;
- heating of emitting cites: Joule heating by field emission currents and heating
by ion impact;
- gas desorption due to the heating of the emitting sites;
- ion stimulated desorption;
- emission of secondary electrons and ions by impact of ions and fast neutrals
born in the electrode gap.
2. Volume processes:
3. Anode processes:
- gas desorption due to local heating by the electron beam and electron stimu-
lated desorption (ESD) of gas from the anode;
2.6. CHARGE PARTICLES IN ELECTRODE GAP IN VACUUM 31
where e – elementary charge, V0 – amplitude value of RF voltage across the gap, ω0 – gen-
erator frequency, me – electron mass, d – gap distance. The corresponding displacement
of electrons is:
eV0
de = 2 (2.17)
ω 0 me d
At relatively low RF voltages (below 5 kV) the displacement of electrons during half of
RF period can be equal to the interelectrode distance d. If these conditions are fulfilled
in a system with electrodes having high secondary electron emission coefficient (see below
section 2.6.5), the number of electrons can be multiplied. The effects in connection to the
electron multiplication are described below in sections 2.6.6, 2.7.3, 3.4.3.
Let us consider the case of high RF voltages. Electrons leave the gap in a time much
smaller than RF half-period when tpass π/ω0 or:
de
1, (2.18)
d
i.e. when the electron displacement during RF period is much larger than the gap distance.
For ω0 < 120 MHz, d = 1 cm criterion (2.18) is fulfilled for voltages V0 > 5 kV. Assuming
that the most part of electrons is emitted when field emission currents are at maximum
(high-voltage phase of RF oscillation) the noted condition exists always and is valid both
for field emission electrons and for the secondary emitted electrons. Therefore in vacuum
for the voltages of interest (V0 > 10 kV) electrons leave the electrode gap fast. The
maximum velocity an electron can obtain is:
r
max V0
ve = (2.19)
2me
32 CHAPTER 2. PHENOMENOLOGY OF RF BREAKDOWN
+
H
20 MHz
O+
Ar+
2eV0 V0 , [kV]
di = 2
= 0.485 , [m] (2.20)
ω 0 mi d µi (f0 , [MHz])2
where mi – ion mass, µi = mi /mp – ion mass expressed in proton mass, f0 – linear
frequency. The simple estimation (2.20) allows to determine the reaction of different ions
on high voltage in Fig. 2.6. The figure shows that light ions (H + and D + ) may pass the
whole electrode gap, i.e. di ∼ d. Heavy ions (µ > 10) can be in principle assumed as
motionless in vacuum during RF cycle. However displacement during half of RF cycle
may be larger than di . For ions having an initial velocity in direction to the electrode with
lower voltage
pcomparable or higher than the velocity an ion can achieve in the voltage
max
gap vi = (2V0 )/(mi ), the effective displacement di is larger. Having initial velocities
of the order of vimax in phase with external electric field means that ion is either back-
scattered at the time close t = nπ/ω, n = ±1, ±2, ±3, ... or it is born by charge exchange
with fast neutrals.
Thus we can conclude that light ions of hydrogen isotopes for the parameter range
considered can follow RF oscillations of electrical field in vacuum. Typical energies of
these ions vary from few hundred eV to eV0 . For heavy ions which can appear in the
gap by gas desorption from the electrodes (e.g. after exposure in atmosphere) one can
assume that di d though a finite ion current onto the cathode exists always. It can be
estimated that the maximal energies of ions collected by cathode are of the order of few
keV for 10 kV< V0 <80 kV.
The electron current onto an electrode is much larger during RF half-cycle than the
ion current for the case when plasma density in the gap and ion space charge are small.
2.6. CHARGE PARTICLES IN ELECTRODE GAP IN VACUUM 33
a) b)
qV
qEn0 = d
d h
h
qE t=β t
qV qV
qE n=β d
vn=√ 2eV
m
d
qV
qEn=-β d
vn = √ 2eV
m
dt
Change of the distribution of electrical field by ion space charge in the gap can be neglected
and the distribution can be assumed as in vacuum as long as:
r
20 V
d (2.21)
eni
where ni – ion density in the gap, V – voltage on the electrodes. Assuming ni =const the
ion density represents the density of plasma inside the gap. The screening properties of
plasma are discussed in details below in chapter 3.
surface, where h is the microprotrusion height, and ending on the surface. Furthermore
we assume that E has a direction corresponding to acceleration of ions to the surface at
a distance larger than h from the surface. At a distance smaller than h from the surface
the field is assumed retarding. The arrows on the Fig. 2.7b present the forces qE acting on
the particle, where q – particle charge. Since h d, the particle velocity
p at the boundary
of appearance of E and on the surface is considered to be v0 = (2V )/(m). Fields E
and E are presented as functions of the main vacuum field En0 = V /d by using field
dt
h
βt
=1
1.2 β
Tangential displacement
βt
=0.5
1.0 β
βt
=0.1
0.8 β
βt
=0.01
0.6 β
0.4
0.2
e Et 2 v2
dt = tn = βt 0 t2n , (2.22)
m d
where tn – a time of flight of the particle to the surface in the fields E and E . The
time is written:
s
1 d 2β h βh
tn = 3F1 (z) − F1 (z) 1 − 2 −2 , z= (2.23)
β v0 F1 (z) d d
where
1 √
F1 (z) = 1 + 1 + 2z (2.24)
2
After substitutions we get an expression for dt :
F12 (z) βt
dt = h (2.25)
2z β
which does not depend on the particle energy and mass. The estimation of dt is presented
in Fig. 2.7b versus z.
Values of βt can be comparable with the values of the field amplification factor β
(section 2.5) which in turn can reach 100 or even higher. For the 1 cm gap this would
mean h = 100 µm that is already a realistic scale for a microprotrusion or a microparticle.
The closeness of z = V h/d to 1 means that the potential drop over the length h in the
region of protrusion is close to V .
We can see from Fig. 2.8 that the focusing of the particle flux is an effect acting on
the surface area close to (π h2 ), i.e. within the microprotrusion. The charged particle
flux has a maximal value in the center of the spikes on the electrode surface which from
2.6. CHARGE PARTICLES IN ELECTRODE GAP IN VACUUM 35
one hand makes the spikes even more dangerous sites for initiation of breakdown, from
the other - there exists another mechanism of the spikes removal in addition to the field
emission effect, i.e. an additional option for conditioning.
Often an ion flux is considered to describe the cathode processes. The total effect of
the ion focusing depends on the ion flux flowing towards the electrode surface. The effect
can be enchanced when a high ion density is present near electrodes and one needs to
account for space charge. Calculations in [56] which considered ion space charge showed
that ion density in the regions of a microprotrusion is about one order of magnitude
higher compared to the density in the regions of flat surface. The increase of ion density
is very crucial for the start-up of breakdown processes: the electric field is enchanced in
the region of microprotrusions and an increased heat flux is delivered to the tips. Based
on considering of heating of the surface tips, an approach for description of processes on
the electrode surface in the presence of high electric field was applied in [57]. It was shown
that evaporation of the microprotrusions and consequent breakdown is possible even only
by heating of the tips by ions, not accounting for field emission currents.
On the other hand, the ion bombardment of the surface tips may lead to an enchance-
ment of field emission by excitation of metal electrons which finally results in a higher
probability of the electron tunneling through the potential barier on the metal-vacuum
boundary.
Thus we can draw the conclusion that the presence of the ion flux (especially a dense
flux, e.g. of plasma ions) to the surface of the electrodes increases the probability of a
breakdown in vicinity of the surface microprotrusions.
a) b)
100 100
Temperature, [grad C]
80
Temperature, [grad C]
t = 10 s 80
= 3 mm
60 t=8s 60
t=6s
40 40 = 5 mm
t=4s
20 20
t=2s = 10 mm
= 5 mm, copper
0 0
0 1 2 3 4 5 0 0.5
distance from the center of the strap, [mm] distance from the center of the strap at t= 5 s
Figure 2.9: Temperature profiles for the pulse of 1 kA RF current: a) different pulse lengths
for stainless steel electrode; b) different strap thickness for 5 s RF pulse, black curves -
stainless steel electrode, gray curve - stainless steel electrode with current carrying copper
layer on the surface.
on the solution of a heat diffusion equation from Appendix B, one can get the surface
temperatures of an antenna strap for a pulse corresponding to the peak antenna current
of 1 kA. This corresponds approximately to the antenna operation at its voltage limits.
The temperature rise is shown in Fig. 2.9 for a strap width of 20 cm. Fig. 2.9a and
black lines in Fig. 2.9b correspond to an antenna strap completely made of stainless steel.
The gray curve correspond to a strap with a current carrying copper layer (coating of Cu)
on the stainless steel strap.
Poor heat conduction of stainless steel needs to be considered and leads to a high initial
temperature of the strap: usually 40-80 o C. The initial temperature should be added to
the temperature in the figure. Fig. 2.9a shows temperature for the pulses of 1 kA current
with different lengths when half of the strap thickness equals to 5 mm (strap thickness
for AUG ICRF antenna is about 10 mm). It is seen that the temperature profile keeps
the form after 1 s pulse and increases in the absolute value.
Fig. 2.9b presents the dependence of the temperature on the strap thickness. Tem-
perature on the surface of the strap is found to be very sensitive to the strap thickness,
especially at l < 5 mm.
When a RF-current-carrying material is changed from stainless steel to copper (e.g.
by deposition of a copper coating of several skin depths on the stainless steel electrode),
the heat flux onto the surface decreases from 136 kW/m2 to 23 kW/m2 . The change of
the surface temperature decreases proportionally (see solution (B.6)).
The calculated surface temperature of the antenna may already reach the maximal
values achieved when torus is baked before operation. Furthermore the rise happens dur-
ing an operational pulse of the ICRF antenna and under the conditions of plasma presence
during ICRF operation. The charged particles of the plasma contribute significantly to
the gas desorption.
2.6. CHARGE PARTICLES IN ELECTRODE GAP IN VACUUM 37
takes place at the energies from 50 to 100 keV. The maximum yield depends very much
on the type of incident ions and the surface material but is usually below 10.
neutral gas consists not only of hydrogen isotopes, but of impurity gases with a higher
mass. In other words impurity atoms (C, Ar, ...) and molecules (N2 , Cα Hβ , Cα Dβ ...)
have more degrees of freedom resulting in a high probability of the step-wise ionization
due to the interaction of excited atoms or molecules in metastable states with other heavy
species.
A comprehensive description of elementary collisional processes and ways of mea-
surements of probabilities of the processes can be found in [64, 65, 66]. The data on
electron-neutral collisions at low energies, ionizing electron-neutral collisions and ionizing
proton-neutral collisions are presented in Appendix C.
high voltage
1000 RF discharge
(hysteresis)
RF amplitude [V]
multipactor RF discharge
electron multipactor plasma with Paschen law
for ignition
RF
dis
cha
rge
100
ex
ist
s
if ignited
-5 -4 -3 -2 -1 0 1 2 3
10 10 10 10 10 10 10 10 10
Gas pressure [Pa]
Figure 2.10: Pressure dependence of voltages of existing forms of RF gas discharge in the
region of multipactor [22, 63].
may be present only for a short period before the transition to the γ-discharge. In the γ
discharge the mechanism of electron birth by secondary emission plays a significant role.
At the pressures p ∼10-100 Pa (we consider only pressure from the range p ≈ 0.1−100 Pa)
the γ discharge is similar to the DC glow with a cathode potential drop. This means that
the cathode sheath is self-adjusted in order to have optimal conditions for the generation of
charged particles in it by ionization and ion-electron secondary emission. The conditions
correspond to the operation in the minimum of Paschen curve for the sheath. The averaged
plasma density profile normally has two maxima in vicinity of each of the electrodes.
A RF discharge at the left branch of the Paschen curve and at substantially high volt-
age exists in the γ-form and is characterized by other similarities with DC glow discharges.
Electrons pass the electrode gap in a time shorter than characteristic time of the change
of external electric field, i.e. the voltage is quasi-constant. However neither in DC nor
in RF case a creation of the effective cathode sheath is possible in the left branch of the
Paschen curve. The mean free pass is comparable with the interelectrode distance and
the sheath thickness required to provide optimal conditions for ionization is larger than
the distance. To ignite such a discharge both secondary emission and volume ionization
are important and the Paschen curve at high voltages has a threshold dependence on
pressure. The threshold pressure depends very much on the gas [69] for the high voltages.
Most of the experimental data exists for voltages below 1 keV. From the experiments
in [70, 71] it is clear that, at the left branch of the Paschen curve, effects connected with
secondary emission processes can take place. At voltages of few hundred volts there often
exists a region where the ignition voltage has a weak dependence on the pressure (ignition
voltage is almost constant for a range of pressures). The dependence does not follow the
classical Paschen curve behavior at the left branch. By insulating the electrodes with
glass plates [70] it was possible to determine that the reason for the behavior is really the
secondary emission. The γ processes participate in the ignition and sustainment of the
discharge by multiplication of the electrons in the gap. According to section 2.6.6, for
the RF voltage the movement of the electrons may be in a resonance with the external
electric field. The effect may be coupled to collisions of electrons with neutrals.
values above unity at high energies and partly compensates for the low ion density, and
therefore plays a role as well as electron-electron secondary emission (but in the energy
range of low emission yield).
2.8 Summary
The power transmission from RF generator to plasma by ICRF transmission lines and
antenna can be limited by the parasitic RF breakdown in the RF system and, in particular,
on the antenna. There exists a number of effects which may assist or lead to the high
voltage ICRF breakdown that should be taken into account simultaneously. RF vacuum
arc, RF glow discharge - self-sustained and semi-self-sustained are the candidates that
may lead to arcing on the antenna.
High voltage operation makes restrictions on the surface state of the electrodes. Ini-
tiation of breakdown happens typically at the emission centers where the work function
is low (or where a non-metallic impurity inclusion exists which makes the electron tun-
neling easier through interfaces) and/or microprotrusions are present which increase field
emission currents from the surface.
We also have concluded that such ions as H + , H2+ , He+ can not be considered com-
pletely motionless during the RF period of the high electrical fields.
When external plasma appears in the electrode gap, a few additional mechanisms can
be triggered simultaneously. We have reviewed the most important ones:
- increase of charged particle fluxes to the surface, local focusing of the fluxes to
microprotrusions;
The processes listed make one think about an improvement of vacuum conditions in
the region of the ICRF antenna that carries high RF currents and is heated.
One of the main processes, the amplification of electrical field on the surface by space
charge effects, is important but not listed above. It is discussed separately in the next
chapter.
Chapter 3
When plasma is present in the electrode gap, the electrical fields on the electrodes may
become significantly higher than the vacuum fields. Ignition of arcs on the surface of the
electrodes is directly dependent on the surface electric field. This makes it important to
clarify the characteristic behavior of the fields during a RF period and the influence of
the RF field on the densities of ions/electrons in the sheath and on the density of the
quasineutral plasma.
In section 2.6 some attempts were made to understand the influence of the presence
of charged particles on the formation of an arc when the influence of space charge is
negligible (i.e. condition (2.21) is fulfilled). Now we initially assume that space charge
can be present in the electrode gap and the condition (2.21) is not fulfilled. However we
consider the neutral pressure to be low enough for all types of particle-particle collisions
to be negligible in the near-electrode regions (sheaths) of the gap where spatial charge is
present. The potential distribution and the flows of charged particles to the electrodes
depend strongly on the density of a plasma which is sustained by a plasma source. The
latter is considered external, though the following description is also suited for the plasma
sustained by volume ionization in the region of quasineutral plasma of the electrode gap.
At first we discuss sheath theory in the case of no magnetic field. That is practically
the same as in the case when magnetic field lines are not nearly parallel to the electrode
surfaces since plasma diffusion parallel to the electrodes is not considered in the 1-D
models being used.
45
46 CHAPTER 3. PLASMA IN THE ELECTRODE GAP
electrons are much faster. The electron mean energy often defines the plasma properties.
Maxwellization of the electron energy distribution is associated with electron-electron col-
lisions. We assume that the external plasma source contains initially thermalized electrons
and ions. Therefore the electron-electron collision time τee can be evaluated by applying
the concept of electron temperature Te (energy distribution function is initially close to
the maxwellian):
11 (Te , [eV])3/2
τee ≈ 3.45 · 10 , [s]; ln Λ = 24 − 0.5 ln ne − ln Te , (3.1)
ln Λ · ne , [m−3 ]
where Te &10 eV, ne – electron density, ln Λ - Coloumb logarithm. For Te = 15 eV,
ne =1016 m−3 one estimates τee ≈ 62 µs that is much longer than the RF period (e.g. 20
ns for 50 MHz). Thus non-maxwellian electrons can be thermalized only if a configuration
of electric field exists in the gap which confines electrons for several RF periods.
If electron collisions with the heavy particles are elastic, the energy transfer from
electrons to the heavy particles is very small, though a one-dimensional component of the
velocity is broadened (but does not thermalize).
When Te & Ti (Ti – ion temperature), electrons are lost on the walls more frequently.
The loss of electrons charges the plasma to some positive potential. The plasma tends
to be quasineutral and a voltage drop is formed on a thin sheath of positive (ion) spatial
charge between the plasma and the electrode forming a ”bell-shaped” configuration of
potentials that partially confines electrons inside the gap. The sheath has a thickness of:
r s
0 kTe 3 Te , [eV]
λDe = = 7.43 · 10 , [m] (3.2)
ne e 2 ne , , [m−3 ]
Here λDe - is the so-called electron Debye length. For typical electron temperatures of
Te = 15 eV and plasma densities ne > 1014 m−3 , the Debye length λDe <0.29 cm and
is thus significantly smaller than the gap distance of few cm. When an electrode is
electrically insulated from any surface with a defined potential, the electrode is floating.
The potential drop between the electrode and plasma (assuming singly charged ions with
adiabatic index of 3) is [73]:
kTe m e Ti
Vsh ≈ − ln 2π 3 +1 (3.3)
2e m i Te
For cold (Ti Te ) deuterium ions one gets Vsh ≈ 3.2 kTe /e. For a system with insulated
electrode, no steady-state electric current flows through the gap. This means that the
potential drop (3.3) appears in order to reflect electrons from the electrodes and make
the electron current to the electrodes equal to the ion current.
If the circuit is closed, electron currents may flow to each of the electrodes (see
Fig. 3.1a). The voltage drop on the sheaths is reduced compared with expression (3.3),
but is still of the order of Te /e. The net electron current is zero if both electrodes of the
gap are completely identical (Fig. 3.1a). If one of the electrode has a larger area, as in
Fig. 3.1b, it collects more electrons and the plasma induces a small net current between
two electrodes. A positive net current flows externally from the small electrode to the
large one.
3.1. APPROACH TO A DC SHEATH 47
ϕ a) ϕ b)
~λ De ~λ De ~λ De ~λ De
kT kT
∼ ee ∼ ee
0 -
-
- 0 - -
- - -
- - -
I th = 0 I th -
- - I th
Figure 3.1: a) Completely identical electrodes, no net electron current; b) electrode areas
are different, net electron current exists.
initially high velocities at the boundary between the quasineutral plasma and the sheath:
r
kTe
vi > v B = (3.7)
mi
where vB is the Bohm velocity. Expression (3.7) is called the Bohm criterion1 . Acceler-
ation of ions to the required velocity occurs in the so-called presheath. The presheath is
quasineutral but a finite electric field exists inside it which ”drives” ions to the velocity
vB . The dimensional scale of the presheath is much larger than λDe and is of the order of
d [74].
Thus the criterion (3.7) should be fulfilled if a sheath of ion spatial charge is formed.
However the formation of the sheaths should be discussed separately when a RF voltage is
applied to the electrodes, though approximately the same approach can be used further.
1
Expression (3.7) is written for the assumption Te Ti (see, e.g. [73]). We assume here and below
that the condition Te Ti is fulfilled, though we keep in mind that a higher ion temperature Ti can
result in a higher ion velocity (and current) at the boundary of the quasineutral plasma and therefore in
a higher ion density in the sheath.
3.2. A RF SHEATH: FREQUENCY RANGES 49
ω0 < ωpe
The plasma is capable to screen the electric field. If Esh > Ep , plasma electrons re-
spond to external RF electric field of arbitrary magnitude, and therefore the exter-
nal electrical field does not penetrate into plasma deeper than the screening sheath
thickness. Also electron density and fluxes are quasi-stationary. Appropriately this
means densities above 5 · 1012 − 1014 m−3 for f0 =20-120 MHz. As we will see in
chapter 5 the condition for the density is fulfilled for the location of ICRF antennas
on ASDEX Upgrade and for the experiments conducted in the frameworks of this
thesis.
Densities above 8 · 1012 m−3 correspond to the case λDe < d, where d ≈ 1 cm –
interelectrode distance, Te is taken 15 eV. However the sheath thickness s may be
either smaller than d or larger than d depending on instantaneous voltage applied to
the electrodes. The discussion on the comparison of s with d follows in section 3.3.
ω0 > ωpe
The screening capabilities of the plasma are reduced. If Esh < Ep , the electrical field
starts penetrating into the plasma since electrons are too slow compared with the
time of change of the electric field. For the frequency range of interest 20-120 MHz
the condition means that electron density should be lower than 5 · 1012 − 1014 m−3
respectively.
Densities below 8 · 1012 m−3 correspond appropriately to λDe > d and s > d for
the usual parameters. The charged particles are moving in the gap according to
the vacuum electrical field distribution and spatial charge effects are neglected. The
latter statement is provided by s > d and is applicable for the higher plasma densities
unless voltage on the electrodes is high enough (see also condition (2.21)).
50 CHAPTER 3. PLASMA IN THE ELECTRODE GAP
2
di ωpi
≈ 2 (3.10)
s ω0
Therefore the ratio ω0 /ωpi defines the displacement of an ion compared with the sheath
thickness.
The condition ωpi = ω0 for f0 = 20 − 120 MHz means in terms of ion density in the
sheath:
ni , [m−3 ] = (9 · 1015 − 3.3 · 1017 ) · µi (3.11)
where µi - ion to proton mass ratio. Higher densities correspond to the condition ω0 < ωpi ,
lower to ω0 > ωpi .
ω0 < ωpi
Ions acquire energy on the RF timescale for arbitrary RF electric fields. Positive space
charge prevails in the near-electrode region. In this case the sheath evolves as a quasi-
stationary (DC) sheath and Langmuir probe theory is applicable. This condition means
also smallness of the capacitive RF current through the sheath compared to the current
carried by charged particles.
ω0 > ωpi
In this range ions can be assumed almost motionless on the RF scale since a displacement
of a single ion during RF cycle is lower than the sheath thickness and therefore much lower
than the gap interelectrode distance. In this case ions respond mainly to time-averaged
electric fields. If ω0 ωpi , many theoretical models use only the time-averaged electric
fields to describe the ion density.
The condition ω0 ωpi is fulfilled for RF glow discharges for heavy (µi > 10) gases.
For hydrogen and its isotopes the condition may not be so strict as for the heavy gases, e.g.
Ar used in the etching/coating reactors based on the capacitive RF glow discharges [68].
Moreover in the limiter shadow of ASDEX Upgrade, where ICRF antenna is situated, the
range of the densities is broad (see section 5.3.1) and the condition may be fulfilled only
for certain plasmas or for a limited time interval.
3.3. PLASMA SCREENING PROPERTIES 51
f pi = 120 MHz
100 100
kV kV
10
kV
s>d
10
s, [m]
}
1k kV
}s
V
0.1
kV 1k
V ~d
Sheath thickness
0.1
kV
s<d
f pi = 20 MHz
f pi = 50 MHz
Ion density, [m-3 ]
Figure 3.2: Sheath thickness versus ion density. Solid vertical lines - boundaries ω 0 = ωpi .
Solid lines - sheath thickness for different voltage over the sheath for n i =const and (3.12).
Dashed lines - sheath thickness according to (3.18), Te =15 eV, when the ion density is
reduced.
- the high density case, when the maximal sheath thickness is much smaller than the
gap distance s < d, the thin sheath;
- the low plasma case, when the sheath is larger than s > d for the most part of RF
cycle, the thick sheath.
The estimation of the ion density for each of these two cases can be understood from
Fig. 3.2 where the sheath thickness (3.12) is presented versus ni for ni =const by solid
lines. One can see that the case s > d is only possible when ω0 ωpi . The case s < d
takes place both for ω0 > ωpi and for ω0 < ωpi .
For the case s < d the vacuum electrical field is strongly changed by the field spatial
charge. This leads to an increase of the electrical field on the electrode surface. One can
estimate this increase by using simple sheath models.
The absolute value of the electrical field on the electrode surface with the potential V on
it is determined for the high-voltage case by:
r r
dΦ 2enp kTe 2e(Vp − V ) 1/4
= |Eel | = · · (3.16)
dx Φ=V 0 e kTe
The last expression 2 shows that, as a result of the reduced ion density across the sheath,
the dependence of the electric field on the electrode surface grows as a root of the fourth
order with the voltage, i.e. very slowly compared to the vacuum case. The sheath thick-
ness in this case is equal to:
√ 2e(V − V ) 3/4
2 p
s= · λDe (3.18)
3 kTe
that is larger than for the model with the constant ion density. The sheath thickness
according to (3.18) is presented in Fig. 3.2 by the dashed lines for different voltages
(Vp − V ).
Similarly (for example, in [34]) it is possible to derive a more general expression - the
so-called Child-Langmuir law for the spatial charge limited instantaneous current ji from
the plasma. The current ji is not necessarily equal to j0 = eni · Vi = ens vB (where ns is
the ion density on the boundary of the sheath and the quasineutral plasma) like it was for
the previous case. The relation connects the ion current which is emitted by the plasma
with the voltage (Vp − V ) on the sheath:
r
2e (Vp − V )3/2
ji = K 1 · 0 · · (3.19)
mi s2
In the ”worst case” when γie is about 10 and one considers hydrogen ions, the resulting change in the
electric field will be below 13%. We can neglect the effect as the uncertainty connected with measurements
of the plasma density np in the scrape-off-layer is higher (see section 5.3).
54 CHAPTER 3. PLASMA IN THE ELECTRODE GAP
case of high plasma density (the plasma completely screens the external electric field)
depends rather weakly on the applied voltage (∼ V 1/4 ). The space charge limited current
for the high voltage sheath is also weakly dependent on the voltage (it is close to the
Bohm current j0 ). However the energy of the ions which carry the current increases
accordingly to the voltage between the plasma and the electrodes. The strong dependence
of the current and the electric field on the plasma density plays an important role in the
conditions of the scrape-off layer plasma where the plasma density can be often changed
by orders of magnitudes.
It is important to consider asymmetry of the electrodes of the gap. The electrodes can
construct either a symmetrical system (electrodes are identical and have the same area)
or an asymmetrical system (one of the electrodes is larger). For the last case one has
a ”large” electrode and a ”small” electrodes characterized by the large and small areas
correspondingly.
Furthermore we will rely on the following assumptions in the further discussion:
(a) quasineutrality of the plasma in all regions of the gap but sheaths;
(b) the sheath is collisionless;
(c) the sheath is a high voltage sheath and is ”electron free”: eV0 kTe leads to
exp(e(Φ − Φp )/kTe ) 1 and the sheath boundary is sharp;
(d) presence of the sheaths with the voltage drop on them ∼ kTe /e means ”elec-
trical contact” between plasma and the electrode;
(e) a change of the ion density in the sheath takes place in a time much longer
than a change in electron density: electric field in the sheath changes instantly
(without a delay) according to the changes of the external field;
(a) quasineutrality of the plasma in terms of the time-averaged densities n̄e = n̄i
in the whole gap;
(b) the gap is collisionless;
(c) voltage distribution over the electrode gap is close to the distribution in vac-
uum.
3.3. PLASMA SCREENING PROPERTIES 55
ϕ
ωt = π/2
V0
V=V0 sin(ω 0 t)
ωt = 3π/2
V=0
0
-V0
x
0 s d-s d
Figure 3.3: Spatial variation of potential at maximal positive and negative voltage (solid
lines) and time-averaged potential (dashed line) for symmetrical electrodes in a high den-
sity case.
capacitance (the resistance of the capacitor is much lower than plasma resistance) is used
in the circuit of RF electrode either in symmetrical or asymmetrical systems. In other
words the case with the blocking capacitor means that the RF electrode is DC-floating.
For the symmetrical system there exists no difference between the case of grounded or
DC-floating RF electrode. For the asymmetrical case (see connections in Fig. 3.4a) the
RF electrode can acquire a potential Vdc .
From the condition of zero net current we can derive the time averaged plasma poten-
tial V̄p [68] assuming the potential of the plasma is:
where ∆Vp is obviously a RF component of the plasma potential. The ion current density
on the sheath boundary
√ 3 is the Bohm current ji = enp vB /K2 , where K2 is a constant
(K2 = np /ns ≈ ē) , relating to a reduction of the ion density in the presheath [80].
In a collisionless sheath near a grounded electrode the Bohm current ji is constant. The
electron current near the ground electrode is determined with the Langmuir formula:
s
1 2eV̄p eV
p
je = · enp · · exp − (3.21)
4 me kTe
R 2π R 2π
Equality of the time-averaged ion ( 0 ji d(ω0 t)) and electron ( 0 je d(ω0 t)) currents leads
to an expression for the plasma potential:
kTe e∆Vp
V̄p = |Vf | + ln I0 (3.22)
e kTe
where |Vf | = −kTe /2e ln(ēmi /2πme ) is the potential of the insulated probe (floating
potential) and I0 ((e∆Vp )/kTe ) is the modified Bessel function of the first order.
For the high voltages ((e∆Vp )/kTe ) 1) expression (3.22) is transformed to:
kTe 2πe∆Vp
V̄p = ∆Vp + |Vf | − ln (3.23)
2e kTe
then by exchanging in (3.21) Vp → Vp − V and V̄p → |V̄p − Vdc |, using (3.23), one can
achieve the relation:
2πe(V − ∆V )
V0 + Vdc kTe 2πe∆Vp 0 p
V̄p = + |Vf | − ln + ln (3.25)
2 4e kTe kTe
Expression (3.25) (see also Fig. 3.3) confirms that for the symmetrical system the time-
averaged plasma potential is close to V0 /2. It was noted in the beginning of this chapter
3
We use ē for the Euler number 2.7182... in order to differentiate it with the elementary charge e.
3.3. PLASMA SCREENING PROPERTIES 57
that for two equivalent electrodes with closed external circuit the plasma potential is
positive due to the fact that a certain electron charge leaves the gap forever. Here we
see that when a RF voltage is applied to the electrode, the plasma potential is charged
more positively since the charge initially going to the electrodes increases. Application
of the RF field leads to a higher plasma potential than (3.3). In this context we can
speak about the so-called sheath rectification effect [77]: a time-averaged sheath voltage
near every electrode is much higher than the voltage in (3.3). The high plasma potential
leads to high ion DC currents towards the both electrodes. The sheath rectification is
accompanied by a self-biasing potential Vdc in the asymmetrical system with the opened
external DC circuit.
System with the RF driven electrode small and the external circuit opened
In Fig. 3.4a the small electrode is RF driven. Also for this system which includes a blocking
capacitor, the ion current is compensated by the electron current for each electrode. This
requires that the electron density should be in contact with the electrodes for some time
during RF period.
At ω0 t = π/2 the plasma is charged to Vdc + V0 when it reaches the small electrode.
At ω0 t = 3π/2 the plasma reaches the large electrode and its potential becomes close to
0. Thus the time-averaged plasma potential should be indeed close to (Vdc + V0 )/2 (see
expression (3.25)).
Self-bias of the RF electrode Vdc is negative if the RF driven electrode is smaller than
the ground electrode. This can be shown from a simple electrical model which represents
the sheath near left electrode as a capacitor C1 of fixed value connected in parallel with
a resistance R1 of fixed value. The capacitor describes the RF displacement current
flowing through the sheath while the resistance is a model for the response of the ion
conductive current in the sheath to the RF voltage. For the right electrode we use C2
and R2 respectively. The capacitors are connected in series with the blocking capacitor
C so that it does not limit RF current: C C1 and C C2 . From the equality of RF
current flowing through the left and the right electrodes:
1 1
ω0 C 1 + V1 = ω 0 C 2 + V2 (3.26)
R1 R2
where the RF components of the voltage drops near the electrode are respectively: V 1 = ∆Vp
and V2 = V0 − ∆Vp , one gets:
C2 + 1/(ω0 R2 )
∆Vp = · V0 (3.27)
C1 + C2 + (R1 + R2 )/(ω0 R1 R2 )
Equation (3.27)4 for the RF voltage drop V2 at the grounded electrode can be also
used for the voltage drop V1 at the RF driven electrode by exchanging C2 to C1 and
R2 to R1 . The equation describes the voltage drops on the hybrid capacitive-resistive
4
Equation (3.27) can be applied to any system as it is derived from equation (3.26) for the RF currents
which is a general relation and is independent of the DC boundary condition and the electrode area.
58 CHAPTER 3. PLASMA IN THE ELECTRODE GAP
ϕ
a) ϕ b)
C1 C2 C1 C2
R1 R2 R1 R2
V0 ωt = π/2
V0
Vp
Vp
ωt = π/2 C ωt = 3π/2
0
0 ωt = 3π/2
-V0
~ -V0 ~
Vdc< 0
x x
0 d 0 d
Figure 3.4: Spatial variation of potential at maximal positive and negative voltage (solid
lines) and time-averaged potential (dashed line) for the system with the small RF electrode:
a) without blocking capacitor; b) with blocking capacitor.
voltage divider that is the linear model for the sheath. It is important to note that in the
limit ω0 ωpi the impedances of capacitances are very high so that capacitive current is
much lower than resistive: ω0 C1,2 1/R1,2 and the voltage divider becomes completely
resistive. For ω0 ωpi respectively ω0 C1,2 1/R1,2 and the voltage divider is completely
capacitive.
For the system from Fig. 3.4a the minimal plasma potential is 0, the maximal is 2V̄p ,
therefore V̄p = ∆Vp . This system has also V̄p ≈ (Vdc + V0 )/2. Hence one gets an expression
for the self-bias Vdc of the right electrode:
C2 − C1 + (R1 − R2 )/(ω0 R1 R2 )
Vdc = · V0 (3.28)
C1 + C2 + (R1 + R2 )/(ω0 R1 R2 )
It is reasonable to consider that C1 > C2 since capacitance of a parallel plate capacitor
is proportional to the surface area A and the left electrode area A1 is larger than the area
of the right one A2 . Also for the resistances: R1 and R2 are inversely proportional to the
surface of a conductor. The role of the conductor here is played by the plasma and one
has R1 < R2 . Thus from equation (3.28) one concludes that Vdc is negative. One also sees
that both types of the dividers: capacitive and resistive describe the rectification effect
(existence of negative Vdc ), hence for ω0 ∼ ωpi the rectification also occurs.
One can write for V1 and V2 and the system considered in Fig. 3.4a:
V0 + Vdc
V1 = ∆Vp = V̄p =
2
V0 − Vdc
V2 = V0 − ∆Vp = V0 − V̄p =
2
As Vdc < 0 and |Vdc | < V0 (see equation (3.28)), one gets the requirement that the RF
component of voltage on the sheath near the large electrode V1 < V0 /2 whereas near the
small electrode V2 > V0 /2. The time averaged plasma potential of the system V̄p < V0 /2.
3.3. PLASMA SCREENING PROPERTIES 59
C1 a) C2 C1 b) C2
ϕ R1 R2 ϕ R1 R2
ωt = π/2 ωt = π/2
V0 Vdc > 0
V0
Vp Vp
C
0 0
ωt = 3π/2 ωt = 3π/2
~ -V0
~ -V0
x x
0 d 0 d
Figure 3.5: Spatial variation of potential at maximal positive and negative voltage (solid
lines) and time-averaged potential (dashed line) for the system with the large RF electrode:
a) without blocking capacitor; b) with blocking capacitor.
System with the RF driven electrode small and the external circuit closed
The ICRF antennas are usually DC grounded. The asymmetrical case for the small RF
electrode is shown in Fig. 3.4b. An asymmetry in the electrode areas leads to the battery
effect in the presence of plasma. Though the electrodes have the same DC potential, a
DC current may be induced in the circuit. The mechanism of self-biasing does not work
for the electrode but it applies to the plasma potential.
The plasma (electron density) should reach at least one of the electrodes, otherwise
a continuous loss of the ions from the plasma would make the system non-stationary.
Since electrodes can not acquire DC potential and the ion and electron currents do not
compensate each other (DC current is present), the plasma can not reach both electrodes.
The alternating component of sheath voltage is lower at the large electrode and the plasma
does not touch the large electrode. An ion spatial charge is always present near the large
electrode. This means that the plasma potential Vp is always higher than 0. If one
accounts for the electron temperature this means that the minimal value of the plasma
potential is higher than (3.3). The time averaged plasma potential V̄p > V0 /2. The effect
can be considered as sheath rectification accompanied by self-biasing of plasma.
The asymmetrical systems shown in Fig. 3.5 can be treated in the same way. For these
systems the RF electrode is the large electrode and the small electrode is grounded. For
the case with the blocking capacitor (Fig. 3.5a) we achieve Vdc > 0. For the case of the
closed external circuit (Fig. 3.5b) the plasma touches only the grounded electrode. For
both cases the plasma potential Vp is 100% modulated and V̄p > V0 /2.
60 CHAPTER 3. PLASMA IN THE ELECTRODE GAP
Blocking capacitor and V0 < VDmax < 2V0 V0 > Vpmax > 0 V0 /2 > ∆Vp > 0
∆Vp
small RF driven electrode V0 /2 < V̄D < V0 V0 /2 > V̄p > 0 V̄p
≈ 100%
(Fig. 3.4a)
Blocking capacitor and V0 > VDmax > 0 V0 < Vpmax < 2V0 V0 /2 < ∆Vp < V0
∆Vp
large RF driven electrode V0 /2 > V̄D > 0 V0 /2 < V̄p < V0 V̄p
≈ 100%
(Fig. 3.5a)
No blocking capacitor and V0 < VDmax < 2V0 Vpmax ≈ V0 V0 /2 > ∆Vp > 0
small RF driven electrode V0 /2 < V̄D < V0 V0 /2 < V¯p < V0 100% > ∆VV̄p
p
> 0%
(Fig. 3.4b)
No blocking capacitor and VDmax ≈ V0 V0 < Vpmax < 2V0 V0 /2 < ∆Vp < V0
large RF driven electrode V0 /2 < V̄D < V0 V0 /2 < V¯p < V0 ∆Vp
V̄p
≈ 100%
(Fig. 3.5b)
Table 3.1: Maximal instantaneous and averaged voltages on the sheaths and plasma poten-
tials for different DC boundary conditions and asymmetry. Left limit corresponds to the
small asymmetry A1 /A2 ≈ 1, right limit corresponds to the large asymmetry A1 /A2 1.
Summary table
We can see that the high density plasma screens the externally applied electric field by
forming a sheath near the electrodes as well as adds the plasma potential Vp to the voltage
drops across the sheaths. The sheath voltages, and in particular their rectified compo-
nents, depend strongly on the asymmetry of the system. The table 3.3.2 summarizes the
rectified component and maximal values of the voltage drop across the sheath near the RF
and ground electrodes. The left limit of the values in the table reflects the symmetrical
case (Fig. 3.3) while the right one is the limit of the highest asymmetry or A1 /A2 1.
The maximal values VDmax = (Vp − V )max and Vpmax define the maximal values of
electrical field on the surfaces of RF driven and grounded electrodes respectively. The
maximal field ultimately influences the probability of arc ignition.
The rectified components V̄D and V̄p define the stationary ion fluxes on the surface of
the electrodes. These fluxes can be essential for the sputtering on the antenna surface.
For the case with a magnetic field, the rectified plasma potentials are also important for
the electrical field distribution in the area of the antenna (see below section 3.4.2).
The modulation of the plasma potential are also described in the table 3.3.2.
The case with the blocking capacitor and small RF electrode for highly asymmetrical
3.3. PLASMA SCREENING PROPERTIES 61
system is characterized by an almost double value of the maximal voltage drop on the RF
electrode compared to the symmetric case and by a high DC self-bias of the RF electrode.
The latter allows for surface processing of the RF electrode surface by bombardment
of ions having energies close to eV0 . The small DC potential of the plasma (non-zero
when accounting for electron temperature) results in low ion fluxes towards the grounded
surfaces, that have large area and could be sources of undesirable impurities. These
properties are used when building the systems for coating deposition and etching on
the base of RF glow discharge or RF biasing of the electrode (substrate) in externally
sustained plasma. The systems are widely used in science and industry. There exists
practically no difference between the DC floating and dielectric substrate for heavy gases
(ω0 >> ωpi ). This enlarges the variety of the materials which can be used for substrate
or coatings. In these systems the power is mainly consumed to sustain a discharge and to
accelerate ions towards the RF electrode. Only a small fraction of the power goes for the
acceleration of ions towards the grounded electrodes. The other 3 cases in the table are
characterized by the formation of high DC ion fluxes towards the grounded electrodes.
The most interesting cases for ASDEX Upgrade antennas are the cases with the DC
grounded electrodes. In the case of the small RF driven electrode, the RF voltage drop
between the RF driven electrode and the plasma has the maximal value between V0 and
2V0 . The same situation is observed between the ground electrode and the plasma when
the ground electrode is smaller than the RF driven one. The sign of the maximal RF
voltage drop is the same for both cases and corresponds to the electric field direction from
plasma to the electrodes, i.e. electron field emission and ion bombardment take place at
that voltage. Due to the almost doubled potential drop for a strongly asymmetric system
compared to a symmetrical system, fluxes of ions with energies twice higher than eV0
may be collected by the electrodes when ω0 . ωpi . Consequently the system of electrodes
should be as close to the symmetrical as possible to reduce the ion energies.
As we will see in 3.4, a magnetic field influences the effective geometry of the electrodes
and requires a different approach for the electrode surfaces oriented differently with respect
to the magnetic field.
One should also note that the modification [76] of the Child-Langmuir law for the
time averaged ion current stated in section 3.3.1 is not applicable for all the systems. The
modification used the model which assumes that a boundary of electron density reaches
the electrode during RF period. The latter condition is not correct, in particular, for the
case of the asymmetric systems without a blocking capacitor (see also Fig. 3.4b, 3.5b)
where the electrons do not reach the large electrodes and the oscillations of the boundary
of the electron density in the sheaths near the large electrodes are smaller, when the
asymmetry is higher. Therefore the Child-Langmuir law modified in [76] should be used
with care for the RF case.
62 CHAPTER 3. PLASMA IN THE ELECTRODE GAP
ϕ
V0 ωt = π/2
V=V0 sin(ω 0 t)
V=0 0
ωt = 3π/2
-V0
x
0 d
Figure 3.6: Spatial variation of potential at positive and negative voltage (solid lines) and
time-averaged potential (dashed line) for symmetrical electrodes in low plasma density
case.
To derive the time-averaged density one uses (3.36) for the first half-period and (3.35)
for the second half-period:
Z Z
1 π 0→π 1 π 0→π
n̄e = n (x) d(ω0 t) + n (d − x) d(ω0 t) (3.37)
π 0 e π 0 e
The integrated expression is:
√ √
ξ0 d √
n̄e = C0 p · x+ d−x (3.38)
(eV0 )/me
where C0 is a constant:
√ Z
2 π
d(ω0 t) Γ2 ( 41 )
C0 = p = √ (3.39)
π 0 sin(ω0 t) π π
64 CHAPTER 3. PLASMA IN THE ELECTRODE GAP
From expression (3.38) one can see that the averaged electron density has a distribution
with a maximum in the middle of the gap x = d/2:
√ ξ0 d
n̄max
e = 2 C0 p (3.40)
(eV0 )/me
and it is inversely proportional to the square root of the amplitude value of the voltage
V0 . The averaged electron density equals the ion density according to (3.29). Therefore
one can speak of a reduction of plasma density by application of a high voltage in the
case when the density is small enough, so the spatial charge effects do not change the
distribution of the electrical field during most part of the RF cycle.
The described reduction of the plasma density may very well affect the neutral pressure
balance in the region of the ICRF antenna. The neutral gas can be ionized by electron
impact and be transferred to the electrodes thus reducing the neutral pressure near the
antenna.
Rectified current
In the symmetric case the RF current is sinusoidal. If the areas of the electrodes are
different, the currents are asymmetric. With the approximation (2.21) one can neglect
the non-linear sheath motion [72, 76]. The current is transmitted mostly by electrons, as
the current density at the two electrodes can be considered to be the same: the density
is the same and defined by the plasma injection parameter ξ0 and the electrode voltage.
The current is higher for the half-period when the large electrode collects electrons. In
a time-averaged picture a direct current flows externally from the large electrode to the
small one if the external DC-circuit is closed. In the first approximation, assuming that all
electrons injected with the external plasma immediately reach the electrodes, the current
can be estimated as:
IDC = ξ0 d (A1 − A2 ) (3.41)
where A1 and A2 – the surface area of the large and small electrodes respectively. This
expression presents an actual saturation DC current which can be reached at high voltages
for the case of the thick sheath.
3.3. PLASMA SCREENING PROPERTIES 65
17 -3
n p= 1⋅10 m
1
uum
vac
0.1
1 10 100
Voltage amplitude, [kV peak]
Figure 3.7: Electrical field on the electrode surface versus voltage between electrode and
plasma: solid black line - reference field in vacuum for the 1 cm gap, dashed black line
- surface electrical field for ni =const (3.13), gray lines - surface electrical field (3.16),
Te =15 eV, from the model which accounts for a reduction of ni . The arrows show that for
the thick sheath the electrical field approaches the vacuum field as the density in the gap
is reduced.
1. the capabilities of plasma to localize electric fields are small: ωpe < ω0 or s & d at
ωpe > ω0 for a given RF voltage and interelectrode distance 7 ;
and
5
The ponderomotive force applies similarly for electrons in RF fields of an electromagnetic wave. In
our case we do not consider RF magnetic field and the Lorentz force associated with this field.
6
Here we do not consider the cases when waves with strong RF electric fields can exist in the plasma.
7
Use Fig. 3.2 for values of the density.
3.4. INFLUENCE OF A MAGNETIC FIELD 67
- a change of the effective collection areas of electrodes with respect to charged particle
movement which is related to the confinement of charged particles as discussed in
section 3.4.1;
- a large difference between the movement of electrons and ions in the RF electric fields
with frequencies close or above the ion cyclotron frequency, the drift approximation
can be applied only to electrons, see section 3.4.1;
- a change of plasma potential leads to consequences for the electric field distribution
since the plasma is confined: a build-up of a complex 3D distribution of electric field
is briefly discussed in section 3.4.2;
- at least two spatial dimensions should be used to treat electron trajectories in mag-
netic field for multipactor conditions: a discussion in section 3.4.3 is devoted to
the features of the electron parallel and perpendicular movement in the external
magnetic field and the influence on multipactor conditions.
where ωce , ωci are the electron and ion cyclotron frequencies respectively.
In addition to the gyromotion, the particles (or the so-called guiding centers) drift in
the direction perpendicular to magnetic field with velocities:
E × B ∇pi,e × B
v i,e,⊥ = + (3.44)
B2 qi,e · nB 2
where E - external electric field, qi = e, qe = −e. The first term is the so called E-cross-B
drift which has the same direction for electrons and ions. The second term is a diamagnetic
drift with the opposite sign for electron and ions and gives an averaged current connected
with gyromotion and the drift of a single particle in the non-homogeneous magnetic field:
2
mi,e v⊥
v d⊥ = · (B × ∇B) (3.45)
2qi,e B 2
68 CHAPTER 3. PLASMA IN THE ELECTRODE GAP
The first order correction to the velocity in (3.45) can be derived by applying the operator
(∂/∂t + ∇ · V ) to (3.45) as a zero-order expression and filling in the first term of it. As
a result we obtain the polarization velocity:
mi,e ∂E
v i,e,pol = · (3.46)
qi,e B 2 ∂t
The particles are magnetized, hence the gyroradii are much smaller than the gap
distance and electrode dimensions. This results in a change of the effective geometry of
the electrodes, i.e. the surfaces should be accounted that are intersected by magnetic field
lines. This makes the radius of curvature of the electrodes important. There exists no big
difference in terms of particle fluxes compared to the case without magnetic field for the
surfaces to which the magnetic field line have rather large angle of incidence α (α > 5◦ ).
Equations (3.44) and (3.45) are valid for the approximation ω0 ωc . The magnetic
moment:
1 2
2
mi,e v⊥
µ= (3.47)
B
is constant in this frequency range. The energy is distributed between v⊥ and vk depending
on the initial conditions and B.
The heating of ions is possible by the RF electric fields at ion cyclotron range of
frequencies. This means ω0 ∼ ωci and µ 6= const which results in non-validity of the drift
approximation for ions, though for electrons one still can use it.
The concept of the ponderomotive force described in section 3.3.5 can be applied
for electron motion along the magnetic fields lines. For ions, a perpendicular non-linear
motion can be also present.
E × B drift multipactor?
The real geometry of the antenna and its surroundings allows for the existence of RF
electric fields in various directions with respect to B. The E × B velocity changes the
sign according to the polarity of the electric field and in principle can play the role of
the velocity parameter in the resonant conditions. This is only possible if the electrons
are not lost during RF cycle along the magnetic field lines, i.e. the electric field along B
is close to zero, because the ratio between the velocity along the magnetic field line and
E × B is very high:
v Bk ωce mi
= ≈ (3.48)
v ×
ω0 me
For the same reason v × requires relatively small gaps and high RF voltages to fit
the resonant conditions. Compared to the multipactor along magnetic field lines which
require V0 ∼ 200 V for d ∼ 1 cm, fulfilling the resonance condition with v × would
require V = 800 kV for d ∼ 1 cm. Such voltages are irrelevant and therefore the resonance
distance for the maximal voltages on the ICRF antennas is about 1/30 cm, much smaller
than the typical antenna dimensions. Furthermore electron energies associated with the
E × B movement are low with respect to the required (∼ 200 eV for stainless steel) to
get the secondary emission coefficient above unity. This makes the E × B-multipacting
improbable.
1. Electrons move along magnetic field lines and are confined, i.e. the perpendicular
diffusion of electrons is limited.
2. Magnetic field lines intersect the wall wich has a complicated 3D geometry. Con-
sequently, it is likely that a system of two electrodes with a interelectrode distance
close to the resonant can be formed.
3. Gyromotion of electrons changes the incident angle modifying the secondary emis-
sion yield which depends on v⊥ , vk (the yield is typically increased since gyromotion
increases vk ).
70 CHAPTER 3. PLASMA IN THE ELECTRODE GAP
The conditions listed make multipator possible at very unexpected regions of the vacuum
vessel in the presence of magnetic field [85]. The multipactor effect is especially dangerous
when a high gas desorption (electron stimulated desorption and desorption due to surface
heating by RF current) leads to high neutral densities which can be ionized. We have
discussed in chapter 2 (section 2.7.3) that the so-called ”multipactor plasma” can be
sustained at low pressures (p & 3 · 10−2 Pa) [63].
3.5 Summary
The electrostatic effects were considered for the description of the plasma - RF antenna
interaction. For a plasma density corresponding to ωpe > ω0 , the external electrical field
can be screened by the ion space charge in the case without magnetic field or in the case
with large angles of incidence of magnetic field to the surface of the electrodes.
If the sheath associated with the spatial charge near the electrodes is thinner than
the electrode gap, the surface electrical field is amplified. The field depends on the the
voltage (V − Vp ) between the quasineutral plasma and the electrode as (V − Vp )1/4 and
√
on the plasma density np as np . This makes variations of the plasma density essential
for the surface electric field. The ion current from the plasma to the electrode is space-
charge limited and energies of the ions which contribute to the current are defined by
the DC voltage drop between the plasma and the electrodes in the case ω0 ωpi or by
the RF potential drop if ω0 . ωpi . The RF and DC voltages between the plasma and
the electrode in a system with different electrode areas can reach 2V0 , where V0 is the
amplitude of the RF voltage between the electrodes. For a system of electrodes with
equivalent areas, the corresponding DC potential is about only V0 /2 while the maximal
RF potential drop plasma-electrode (V − Vp ) ∼ V0 . This makes the consideration of the
antenna electrode asymmetry important in order to reduce the voltage drops near the
antenna electrodes and to reduce the energy of ions which bombard the surface of the
electrodes. The presence of the external magnetic field should also be taken into account
as it influences the effective area which collects charged particles.
If the sheath is thick compared with the interelectrode distance, the applications of
a high RF voltage to the electrodes reduces the plasma density in the gap. For the
3.5. SUMMARY 71
intermediate case, when sheath thickness is of the order of the gap distance, the plasma
density is also reduced. In addition the neutral pressure balance is affected, since the
neutral gas can be ionized and transferred to the electrodes.
The magnetic field has additional consequences for the RF breakdown: multipactor
and RF glow gas discharges by simplifying the ignition of both. Geometrical consider-
ations are important to prevent the discharge ignition: avoiding the distances between
two conducting surfaces along the magnetic field lines which fit the resonant conditions
for the multipactor (in practice this can be only partly realized by using simple geome-
try with well-defined interelectrode distances); minimizing the effective electrode distance
d∗ = V gas /S, where V gas – geometrical volume of the neutral gas, S – surface limited by
the intersection of the electrodes by magnetic field lines.
Chapter 4
Experimental approach
The ICRF antennas in a fusion device with magnetic confinement are situated in the
scrape-off layer region. For ASDEX Upgrade the antennas are placed at the low field
side. Many mechanisms of charged particle transport from the core plasma to the outer
wall (see also chapter 1) exist in a tokamak. The particle transport affects properties of
the SOL plasma. The presence of the intermittent plasma fluxes towards the wall due
to the appearance of instabilities at the plasma edge, for example edge-localized modes
(ELMs), is an inherent property of the fusion experiments. Similar plasma conditions can
not be reproduced in a test facility.
The ICRF antennas operate under these conditions, but can not provide an explicit
knowledge on the high voltage breakdown for several reasons:
1. The antenna has a complex 3D-geometry with various alignments of the electrodes
with respect to the magnetic field. This makes it difficult to analyze experimental
data on the local phenomena and predict the critical aspects of the antenna in terms
of the RF breakdown.
2. The whole range of the RF voltages from zero to the maximal one (the latter be-
ing defined by the RF power coupled to the antenna) is present on the antenna
components which face the SOL plasma.
3. The antenna is a tool for the plasma heating and is usually not diagnosed specifically
to get a detailed knowledge on the RF breakdown.
4. The plasma heating by the ICRF antennas affects the tokamak central plasmas as
well as the edge plasmas. This makes the problem strongly non-linear and very
difficult to treat (as we already noted in the beginning of chapter 2 we apply a ”first
order” approach in the frameworks of this thesis trying to avoid the global plasma
changes).
72
4.1. CONCEPT OF THE EXPERIMENT 73
a) b)
high RF voltage
high RF voltage
RF current
flowing
Figure 4.1: Experimental modeling of the ICRF antenna high voltage region: a) inductive
loop of the ICRF antenna, b) model: coaxial open end of a RF probe.
and should not lead to the global plasma changes as occur during operation of the ICRF
antennas.
This chapter is devoted to the description of an experimental device - an active ICRF
probe which has been developed as a part of this thesis in order to gain data on the local
effects taking place during the interaction of a high RF voltage with the SOL plasma.
RF generator
~
measurements feeding line
vacuum of RF voltage
probe head matching element
feedthrough and current
(open end) (short circuit)
Plasma
V - +
coaxial resonant line + -
DC circuit
Standing wave pattern:
V
Figure 4.2: Concept of the ICRF probe experiment. Two standing wave patterns cor-
respond to two resonant frequencies for the resonator, the solid line corresponds to the
frequency used for the matching.
The matching of the resonant system uses the following scheme: by adjusting the
frequency, the input impedance at the point of the connection of the feeding line is fitted
to an impedance equal to Z0 + iY , where Z0 is the characteristic impedance of the feeding
line, Y is the imaginary part of the impedance which is compensated by the stub tuner
reactive impedance Ystub = −Y at the same location. Therefore the input impedance
of the whole resonator equals to Z0 which corresponds to the matching. The matching
position can be changed only between RF pulses, therefore for a single RF pulse one is
limited to one matching position.
In Fig. 4.2 two standing waveforms are shown. Both correspond to the matching of
the resonator. Two frequencies can be used for the ≈ 5λ/4 resonator with a difference of
about 2 MHz. The highest frequency of 51.52 MHz (solid gray line) is chosen to have a
lower RF voltage at the position of vacuum feedthrough (see also the next section for the
description of the setup).
Because of the sheath effects described in chapter 3, it is important:
Therefore the inner conductor facing the plasma is DC-disconnected from the rest of the
resonant line by an inner DC break (a capacitor for the inner conductor). The stub tuner
with the short circuit allows a cable connection from the inner conductor at the plasma
side to an external DC circuit.
4.1. CONCEPT OF THE EXPERIMENT 75
Measurements of the RF voltage and current are available . λ/2 from the open end.
The directional couplers in the feeding line measure the forward power from the generator
and the reflected power. A more detailed description is found in section 4.4.1.
Thus the concept for the RF probe working in the ion cyclotron range of frequencies
allows for:
At first we describe the power generators used in the experiments. The experimental
device of the ICRF probe is subsequently discussed in more detail.
DC power supply
Four DC power supplies KEPCO BOP 100-4M [87] were used for experiments where
the inner conductor was set at a certain DC bias. The supplies can work in pair in
master/slave modes. One power supply generates 100 V DC voltage with a maximal
current of 4 A.
76 CHAPTER 4. EXPERIMENTAL APPROACH
coaxial output 2 MW
output lines
variable short
input matching
grid circuit
~
20 mW >
100 W
>
4 kW
>
100 kW
oscillator
amplifier amplifier amplifier
∅ 89 mm
∅ 55 mm
Probe head
front view:
∅ 4 mm threads for heating wire
mounting of tips
6 mm
1) Probe head
The drawing of the probe head is shown in Fig. 4.4. The probe head consists of a coaxial
inner and outer conductor. The outer conductor is a 3 mm thick stainless steel tube
with outer diameter of 89 mm. The inner conductor has a diameter of 55 mm. The
corresponding characteristic impedance is 25 Ohm.
The inner conductor consists of an exchangeable head close to the plasma, the in-
termediate piece (stainless steel) with a heating wire for outgassing, and the long inner
conductor (stainless steel tube) of the resonant part of the line.
The head was manufactured from stainless steel and was mechanically polished to the
microroughness parameter Rz < 0.5 (see [86] for a detailed description). The head design
allows for the installation of tips of different shape with a height of 4 mm for a local
increase of electric field on the electrode surface. The tips were prepared in the same way
as the head [86] and can be mounted with the help of a thread in 4 mm diameter holes in
the head. The four holes are distributed evenly around the probe head (i.e. 90o between
each hole). The holes can be adjusted with respect to a direction of an external magnetic
78 CHAPTER 4. EXPERIMENTAL APPROACH
pressure gauge
vacuum
vacuum
feedthrough: feedthrough
massive
cylinder for the cables
with 3 holes quartz rods
for the cables
space for
a vacuum
feedthrough
ceramic disc for optical fibers
(not used in the
frameworks of
metallic this thesis)
potential rings
field. The tips were not installed for the experiments in a magnetic field performed in the
frameworks of this thesis.
ode
electr
ded inner
DC-groun
ode
electr
inner tching
ontro
lled
eak to ma
DC-c DC br circuit
ak
wer le
RF po
outer RF
sma 4
to pla conductor
screening
massive shield capacitors
insulator of connection to the inner RF at the cable
3 the DC break inner conductor conductor output
1 screening (DC ground)
DC cable 2 metallic
capacitor
connection screening shield
to the inner capacitor insulators
conductor of the cables
Figure 4.6: Screening of the cable placed in the RF inner conductor to close the circuit
for the RF power leak to the inner space of the DC break of the inner conductor.
3) DC cabling in RF resonator
The voltage strength of the DC-break and the cables inside the inner conductor should be
high (at least > 5 kV) to withstand the DC voltages. In addition there exists the problem
that RF power can be coupled to the inner conductor of the main coaxial line through
the inner DC break (Fig. 4.6).
The estimation of the RF power can be done as follows. We know that electrical field
in a capacitor is proportional to a voltage on the capacitor electrodes. In our case the
voltage across the capacitor is a voltage drop on the DC-break defined by the current
flowing on the outer surface of the inner conductor and RF impedance of the capacitor
1/ω0 CDC−break , where CDC−break =5 nF is the DC-break capacitance. The capacitance
is relatively low for the frequency used and the RF voltage drop across the DC-break
can be rather high. Furthermore due to geometrical restrictions and mechanical design
the DC-break is located in the vicinity of the maximum of the RF current. The current
maximum is estimated to be about 2.4 kA for 60 kV voltage maximum in the main probe
line. The corresponding voltage drop on the DC-break is ≈ 1.5 kV. The voltage across
the DC-break can be considered as the local voltage on the parasitic transmission line
formed by the inner surface of the inner conductor and an electrical cable inside it. If
the power reflected back from the parasitic line is assumed to be small and characteristic
impedance of the cables is ≈ 50 Ω, the resulting coupled power is somewhat higher than
20 kW (this is the upper limit estimation).
Even at a coupled power of a few kW, existence of a resonant condition inside the
uncertain configuration of the electrodes in the inner conductor can lead to very high
voltages and damage of the cables. Therefore the proper RF screening of the cable is very
80 CHAPTER 4. EXPERIMENTAL APPROACH
important.
The scheme of the screening is illustrated in Fig. 4.6. Several screening techniques
were used (see the corresponding numbers in the figure):
1. the cable is screened by a metallic shield with a massive electrical connection to the
DC grounded part of the inner electrode;
3. a RF filtering high-voltage capacitor is mounted between the shield and the DC-
controlled inner electrode at the boundary between shielded and unshielded cable;
4. RF filtering high-voltage capacitors are mounted at the output of the cable from
the matching resonator line (a scheme of the short circuit is shown in Fig. 4.6).
In reality two DC cables are used in the inner space of the inner conductor which allow
simultaneously to control the DC boundary condition and to heat the inner conductor
head. One cable is connected directly to the DC-controlled side of the inner conductor,
the other is connected to the same electrode through the heating spiral shown in Fig. 4.4.
With this scheme the inner conductor can be heated by an AC 50 Hz current from an
insulation transformer. At the same time one can control DC conditions or perform DC
measurements with the help of the cable directly connected to the inner conductor. For
the scheme of two cables the shield (number (1) in Fig. 4.6) was used for both cables
(twisted pair) and two additional screening capacitors were mounted: one connecting the
both wires of the cables at the boundary between shielded and unshielded cables and one
connecting the second wire to the DC-controlled inner electrode.
When the RF power leak was not suppressed by the screening, arcing between the
DC cables and the shield occurred as well as arcing between the shield and the inner side
of the inner conductor was observed: the cables and the inner electrode were damaged.
The arcing inside the DC cables occurred in the regions of the cable bending. After the
described screening techniques were used, arcs did not appear and the cables were not
damaged during the long term experiments. Also RF power was screened well at the
input/output of the resonator matching circuit (RF amplitude was much less than 1 V
on the end of DC cables disconnected from the external power supplies) and therefore RF
power did not penetrate into the external DC power supplies which can be also damaged
if the power level is high.
4) Matching circuit
A 9-inch line segment with a massive variable short is used as a matching element (see
Fig. 4.7). The 9-inch line was used instead of a 6-inch line in order to increase the reliability
of operation of the pressurized lines under high voltages. The matching concept allows
for an output of the cables from the inner space of the inner conductor. The cables, after
being RF filtered at this output, are connected to the DC measurement (control) lines
(see section 4.4.2) and to the AC insulation transformer for heating respectively.
4.3. SETUP OF THE EXPERIMENT 81
matching line
connection to
the probe line
feeding line
power input
connection to the
flexible coax. cable
DC break
(inner+outer) RF probes
5) Feeding line
The feeding line is a line with 50 Ω characteristic impedance. The line is equipped with
a DC-break for inner and outer coaxial conductors to avoid a ground loop between the
grounds of RF generator and toroidal chamber. The RF pick-up probes mounted on the
line allow for measurements of forward and reflected power. A 3-inch flexible high power
cable of about 6 m long was used to allow some radial movement of the probe by midplane
manipulator motor drives. A 6-inch transmission line to RF generator with iron inner
conductors was used to maximize attenuation between the probe and the generator. The
attenuation is used to minimize influence of the fast changes of the load of the probe on
the generator.
CCD camera
metallic grid
The RF pick-up probes were calibrated using special lines. In the case of the RF probes
mounted . λ/2 from the open end the same parts of the main coaxial lines were mounted
separately and terminated by 50 Ω resistance as a load. The transmission measurements
were conducted by a Hewlett-Packard 8753D [88] network analyzer in the frequency range
30-60 Mhz. The corrections were made for the line with the standing wave ratio k=2
(characteristic impedance of the probe lines used for calibration is 25 Ω). The coupling
factors for the pick-up RF probes for the forward and the reflected power in the feeding
lines were measured in a similar setup without the corrections for k since the characteristic
impedance of the feeding lines equals 50 Ω and k = 1.
The probe can be mounted in the test facility field (see Fig. 4.8). A 6-inch tube is used
as a vacuum vessel. The outer conductor of the probe was supported by a support made
of an insulator. The parasitic coaxial line between the outer conductor of the probe and
the vacuum vessel (the tube) was screened from RF power by a shielding ring (shield in
the figure). A vessel with multiple ports was installed for pumping, for the installation
of a plasma source and for optical observation of the ICRF probe head. The port with a
glass window in front of the probe head was screened by a metallic grid.
A Hall-type ion source was used as a plasma source. Ions in the source are accelerated
by the anode layer of the axisymmetric discharge in crossed electric and magnetic fields
[89]. The discharge exists in the circular electrode gap with 10 cm diameter and 4 mm
thickness. The ion source produces a tubular beam with a total current of 5-300 mA,
average current density 0.3-20 mA/cm2 , ion density 1014 -1016 m−3 and average ion energy
200-2000 eV. The tubular ion beam is injected perpendicular to the open end of coaxial
4.3. SETUP OF THE EXPERIMENT 83
line. The mode of operation allows to have the ion energy distribution function close to
a maxwellian with the maximum at about 30% of the accelerating voltage.
The plasma density was measured near the ICRF probe head by a Langmuir probe
collecting beam ions. The density was estimated from the ion saturation current using
the ion velocity corresponding to the maximum energy of the initially known ion energy
distribution function.
In addition a CCD camera was installed with a time resolution of 20 ms for observation
of the ICRF probe head (see Fig. 4.8).
Data acquisition
A Nicolet 460 oscilloscope system [90] measured the amplitudes of the RF signals with a
time resolution of several tens of milliseconds using well-calibrated linear RF detectors.
Fast Acqiris DC 265 8-bit digitizers [91] with a bandwidth of 150 MHz and a sampling
rate of 500 Msamples/s allow to save the RF signals digitally. The RF amplitudes and
the phase were obtained by fitting the waveform with sinusoidal functions of time.
When a RF breakdown leads to ignition of an arc, the arcing generates a noise with a
broad spectrum of frequencies. A high fraction of the energy is carried by low frequency
modes which often lead to saturation of the saved signals. To prevent the saturation,
the RF signals passed through analog high-pass filters of the first order with a resonance
frequency of 1 MHz.
The data acquisition system was saving the data of the RF detected signals and DC
voltage/current and was integrated into the shotfile system of ASDEX Upgrade, except
for the Acqiris system which was equipped with the high-pass filters and saved the last 2
to 4 ms of the RF voltage, current and power before a breakdown. The breakdown was
detected by setting the critical reflected power level as the trigger for the data acquisition.
4.4 Measurements
To get an information on the processes occurring at the probe head a few options for
the measurements are available. Measurements by RF probes and the DC circuit are the
most important to characterize RF breakdown phenomena.
4.4.1 RF measurements
Because of geometrical constraints set by the construction of the midplane manipulator
the RF voltage and current are measured only a distance of . λ/2 from the probe head.
The location of the measurements is few tens of cm closer to the open end than the voltage
maximum corresponding to ≈ λ/2. This allows to measure RF current of a relatively high
amplitude: one has no problem of a small signal.
However the values of the voltage and the current at the probe head should be recon-
structed from these measurements. The reconstruction is done by calculating the standing
wave pattern with a constant load impedance at the head. The load impedance for the
ICRF probe in vacuum is modeled as: a resistive part Rr and a stray capacitance Cstray
connected in parallel. If Rr Z0 , where Z0 = 25 Ω is the characteristic impedance, the
voltage standing wave pattern has the same form for any value of the input power and
4.4. MEASUREMENTS 85
where Vrfmeas is the measured value of the RF voltage and A = 1.359 is a reconstruction
coefficient. From the distant measurements of RF voltage and current it is possible to
deduce the information on the load impedance of the open end. However when using the
amplitudes of the RF signals, one needs to deal with a few problems. Values of RF voltage
depend strongly on the matching position, hence one needs to account for the mismatch.
Very fast processes resulting in fast load changes have a small influence on the RF voltage
because a high reactive power is stored in the high-Q RF circuit.
The RF voltage on the open end is controlled by a waveform of the input power. Thus
the resonator can be represented as a wide-band amplifier. For the wide-band amplifier
transient characteristics are often described in terms of a rise-time trise , the time it takes
the envelope of modulation to increase from 10% to 90% of its final value [92]. The rise
time can be estimated from:
0.7
trise ≈ (4.2)
B3dB
where B3dB is bandwidth of the amplifier [92] which can be expressed in terms of RF
frequency f0 and Q-value:
f0
B3dB = (4.3)
Q
Assuming a constant Q-factor, for the parameters noted above trise ≈3 µs. This time
resolution is not enough to resolve a vacuum arc.
The output power of the RF generator may be strongly affected if the reflected power
on the generator side exceeds the critical level (2.7). For this reason the RF voltage often
becomes an unreliable physical magnitude. A good option for a reliable measurement of
the load impedance is to use the phase value between voltage and current at the point of
measurement. For an electrotechnical model with a constant stray capacitance one can
obtain a relation between the phase and the load resistance.
If we assume the stray capacitance to be Cstray = 5 pF, the voltage and the phase
depend on the open end resistance according to Fig. 4.10. One can see that the phase
has a low sensitivity to a change of the load resistance and the evolution of the load can
not be detected by the phase measurements if the load changes too little. However we are
interested in the effects connected to the change of the load resistance by orders of magni-
tude (transition to RF glow discharge, RF sustained arcing) and the phase measurements
are applicable.
The method of time-resolved voltage and current measurements has been used in more
simple systems to measure the power coupled to a discharge, see e.g. [93]. In our case it
is convenient to use Q measurements to estimate the power coupled to the load.
To measure the power coupled to the plasma Pcoup by the ICRF probe for load changes
happening slower than the time (4.2), one can apply measurements of Q-factor and use
the earlier measured Q without a load. This requires the assumption that the power going
to the intrinsic losses in the conductors is the same for the same voltage for the cases with
86 CHAPTER 4. EXPERIMENTAL APPROACH
Figure 4.10: Sensitivity at the point of measurements for 300 kW net power. Voltage
curves: 1 - matched for 1 MΩ, 2 - matched for 5 kΩ, 3 -matched for 1 kΩ.
where Q is the Q-factor measured with the load, Qno load is measured without the load
(in vacuum). This simple equation may be applied to the ICRF antenna as well. One can
see here that the ideal case Pcoup = Pnet (where Pnet is the net RF power) corresponds to
the case Qno load → ∞, since Q ≥ 1. Equation (4.4) may be written in terms of energy
conservation:
Q
Pnet = Pcoup + Ploss , Ploss = Pnet (4.5)
Qno load
meaning that the net RF power goes partly to the load and partly to the losses.
4.4.2 DC measurements
The DC external circuit may be changed as shown in Fig. 4.2: open DC circuit (inner
conductor is DC floating), a low resistance can be connected (closed DC circuit) or a DC
bias can be applied. The first configuration allows to measure the rectified voltage, the
second - to measure precisely the rectified current, the third - to affect the current by
connecting an external power supply with a certain polarity.
Most of the experiments in the test facility in the presence of plasma were performed
with the inner conductor DC-floating. This configuration in the test facility allows one
to have the measurement of the rectified current from the charging characteristics of the
DC-break capacitor if the plasma conditions can be assumed constant.
Rectification effects were discussed in chapter 3. The rectified self-bias voltage can
reach the value of the order of RF amplitude in the DC-floating configuration for the inner
4.4. MEASUREMENTS 87
conductor when the plasma density is high. For the ASDEX Upgrade the high self-biasing
voltage on the inner conductor is possible during the high density transient phases (e.g.
ELMs). Technically it is not allowed to have voltages above 5 kV on the inner conductor.
For this reason, the majority of experiments with the ICRF probe in ASDEX Upgrade
were carried out with the resistance Rdc = 10 Ω (see Fig. 4.2) connected setting the DC
boundary condition to the DC short circuit. This configuration is very similar to that of
the ASDEX Upgrade ICRF antennas.
The DC biasing experiments were conducted by connecting the DC power supplies
(described in section 4.1.1) in series with a balast resistance Rbal = 50 − 100 Ω . The
resistance Rbal was used to limit the DC current when an arc appears in the probe and
transforms a high resistance load to a short circuit.
Chapter 5
In the following chapter the experimental results both from the test facility and from
ASDEX Upgrade are presented. The studies were performed in three steps:
2. experiments with a Langmuir and the ICRF probes on the midplane manipulator
of ASDEX Upgrade;
3. studies of the high voltage operation of the AUG ICRF antenna and comparison
with the ICRF probe.
88
5.1. TEST FACILITY RESULTS 89
top). However the same values of voltages as for the probe without the tips were achieved
by intense high-voltage conditioning in vacuum.
Since the energy is mostly dissipated in the system as intrinsic RF losses, the match-
ing position changes during the pulse due to heating of the resonator components (sec-
tion 2.6.3). RF pulses of 200 ms duration at a voltage of 50 kV are already critical for the
system matching. Using the formula (B.6) in Appendix B, one can estimate the temper-
ature increase of the conductors. The estimation results in about 50 o C for a current of
2 kA with a duration of 200 ms on a conductor with 3 mm thickness. The temperature
coefficient of the resistance of the stainless steel is about 5 · 10−3 (o C)−1 . The surface
resistance Rsq at 50 o C is 25% higher than the resistance at the room temperature. Ac-
cording to the theory briefly reviewed in section 2.2 the change of the phase along the
transmission line due to the losses is assumed to be zero. If one assumes that the observed
mismatch happens only because of the changed resistance of the stainless steel compo-
nents and makes an estimation in terms of the above mentioned theory, the corresponding
change in the absolute value of the input reflection coefficient |ρf eed | is 8 · 10−4 . This is
too small to explain the experimentally observed increase of the reflected power of about
10-20% of the forward power at the feeding point that is proportional to |ρf eed |2 . Thus the
change of the reflected power is connected either to a change of the phase characteristics
of the stainless steel lines between the open end and the feeding point or to a change of
characteristics of the matching circuit. The change of the geometry (e.g. lengths) due to
the heating of the electrodes can affect both. A further estimation of these effects is not
considered in the frameworks of this thesis.
To minimize the reflected power level, the ICRF probe was operated in the test facility
mostly with short pulses, while for the long pulses in the test facility and in ASDEX
Upgrade the matching was chosen to minimize the reflected power level at the time when
a RF breakdown was anticipated.
The other aspect of high voltage in the resonator with stainless steel electrodes, is a
high gas desorption which leads to even higher importance of the pulse length during the
conditioning and operation.
2 µs
800
[kW]
400
forward power
0
[kW] 400 reflected power
200
0
120
[kV]
80
voltage at the head
40
0
0
phase
[degree]
-40
-80
follows the initial phase change and is finished when the arc is developed. The duration
of the spark stage of a few hundred nanoseconds agrees with the expectations from the
discussion in section 2.5.3.
In addition to the phase jump of about 900 , the increase of the reflected power and an
increased light emission observed by the CCD camera confirm the vacuum arc ignition at
the probe head. It is also observed that the RF voltage degrades very fast during the arc
development. For a constant Q-factor the time required for the voltage to change should
be about 3 µs (see section 4.4.1). The observation of a faster change of voltage contradicts
to the assumption Q=const and therefore indicates a fast decrease of Q. This means that
the vacuum arc after the spark stage dissipates the energy very efficiently.
As was discussed in section 2.4.2, the evolution of the arc discharge is dependent on
the Q-factor. For Fig. 5.1 the switch-off of the RF power by the generator corresponds
to the time -0.5 µs. The increase of the reflected power and the decrease of the voltage
before the switch-off is a demonstration of the self-adjusting mechanism of RF arc. It
should be also noted that the RF voltage shown in Fig. 5.1 after the arc ignition may be
overestimated since it includes the contribution from the arc produced noises in the broad
frequency range. The plotted voltage is the amplitude of the sinusoidal fit to the data
collected with the help of the fast data acquisition system (described in section 4.3.1).
For the whole series of experiments with vacuum arc (which also included variation
of the critical power level for the RF generator switch-off), the damage of the electrodes
was observed only as single points. Thus the vacuum arc for the ICRF probe is self-
suppressed by a mismatch which follows the transformation of high resistance load to the
short circuit.
5.1. TEST FACILITY RESULTS 91
a) high arcing
probability
b) high arcing
probability
Rectified current, [µA] 400
10
200
100 1
0 0.1
20 40 60 80 100 20 40 60 80 100
Voltage, [kV peak] Voltage, [kV peak]
Figure 5.2: Rectified field emission current:a) linear graph, b) logarithmic graph. Dia-
monds - measured current, solid curve - fit to the formula (2.14) by varying A RF and
βRF .
circuit).
The dependence on the voltage shows clear exponential behavior peculiar to the field
emission (see Fig. 5.2a). The measured current is actually a time-averaged difference
between the current emitted by the inner electrode and the current emitted by the outer.
For the field amplification of 7.5 for the inner electrode one can neglect the contribution
to the current by the outer electrode because of the exponential part dependence of the
current on the electrical field in equation (2.14).
A standard fit [44] by the two parameters: the emitting area ARF and the local
amplification factor βRF (see section 2.5.1) can be successfully applied to the experimental
data when the amplification β0 · βRF = 7.5 · βRF is accounted for the inner electrode. A
Poisson-weighted least-squares fit results in ARF = 1.99 · 10−16 m2 , βRF =334 and the
corresponding solid curves in Fig. 5.2a,b.
For the voltages above 80 kV many RF pulses longer than 10 ms are terminated by
arcing at the probe head. Before the data presented in Fig. 5.2 was taken the probe
and the tips have been conditioned to have the arcing probability low enough for making
about 10 consequent arc-free RF short (3 ms) pulses possible.
In section 2.6.2 of chapter 2 the focusing of the ions on the microprotrusions was
discussed. The focused ion flux may in principle influence the emission characteristics of
the protrusions and the rectified current. To look for this effect an externally sustained
ion density (plasma) presence is required.
400
200
100 ion
m iss
field e
+
0 a
s m
a n
pl io
c e iss ted
u r m
e t r ac
-100 so d bs
n el su
io fi
sat. current
-200
0 20 40 60 80 100
Voltage, [kV peak]
Figure 5.3: Rectified field emission current and plasma current. Squares - measured recti-
fied current, black solid curve - fit to the formula (2.14). Stars - measured rectified current
in plasma of the ion source, black dotted curve - polynomial fit to stars, gray solid curve
- fit to the formula (2.14) accounting the saturation current. Black dashed line - differ-
ence between the fits for the currents by a subtraction of the solid black and the solid gray
curves.
of the positive component with plasma and the positive current without plasma. The
increase is observed also with other plasma injection parameters. We can note that the
positive rectified current is changed when an external plasma is present. However it is
rather difficult to determine the exact reasons for the current increase since a few effects
act simultaneously. At least three processes can contribute to the rectification current:
1. Ion bombardment and modification of the surface.
Ions modify the surface structure (also due to focusing to the tops of microprotru-
sions (section 2.6.2)). Eventually this results in the decrease of the field amplifi-
cation factor βRF . However during the bombardment, an increased field emission
yield may occur (charged particles can pre-heat the emission centers and excite
metal electrons).
2. Secondary ion-electron emission.
At high ion energies the emission has a yield comparable to or even higher than unity
(see section 2.6.5). The total yield of the emission is proportional to the integral
ion flow to the surface and is not directly connected to its density (or surface area).
However in the gap between the coaxial electrodes, the electric field is higher near
the inner electrode. If the injection of ions into the gap is assumed homogeneous,
the ion flux is higher to the inner electrode due to the distribution of electric field.
This may lead to the positive rectified current. Also an electric field distribution
in the 3D geometry of the ”parasitic electrodes” - vacuum chamber and the front
surface of the inner conductor affects the rectified current by different total yields
of the secondary emission.
3. Volume ionization.
If only ionization by electron impact is considered, it is stronger in the RF period
corresponding to the field emission from the inner electrode since electrons emitted
by electric field contribute to the ionization in addition to externally injected elec-
trons. In principle this can result in a contribution to the positive rectified current.
On other hand an increase of the time-averaged plasma density results in the nega-
tive rectified current. Ion impact ionization should contribute more to the negative
rectified current by producing charged particles in the volume.
An increase of the positive current is observed when one switches from the case of
no external plasma to the case with plasma. Though the reasons of the increase can not
be distinguished exactly (we can only say a contribution of ionization is less probable
than a contribution from the processes occurring on the electrodes since also no light
emission - no excited atoms were observed from the volume), this observation confirms
that the processes of generation of charged particles in a high voltage electrode gap are
significantly enchanced when an external plasma is supplied.
The existence of a plasma density of the order of 1015 m−3 at low neutral pressure
pHe 6 0.1 Pa does not influence the high voltage RF breakdown. Arcing is initiated with
no detectable difference compared to the vacuum arc. Also the Q-factor does not change
(the change is within the error bars, the error being about 10% for the Q-value). In
Fig. 5.4 the measurements of voltage depending on the net RF power are presented by
5.1. TEST FACILITY RESULTS 95
100
30 Ω
=2 =2 Ω
Q R0
00 =4
=1 R0
Q 0
=5
Q 2.5
=1 5
Q .2
=6
Q
1
1 10 100 1000
Net power, [kW]
Figure 5.4: Operation with and without ion source. Black circles - measurements in
vacuum, black dashed line - fit for Q=230, gray circles - measurements at low neutral
pressure with the presence of the external plasma, black triangles - measurements in the
semi-self sustained RF glow discharge in He, black lines correspond to Q = 12.5 (ICRF
antenna operation with R0 =2 Ω) and to Q = 6.25 (R0 =4 Ω).
black circles for the vacuum case and by gray circles for the case with the external plasma
presence at the pressures of helium below 0.15 Pa. For comparison with the voltages and
powers of the real ICRF antenna operation, the maximal voltages on the antennas are
shown for two values of the effective resistance R0 : R0 = 2 Ω and R0 = 4 Ω.
Self-sustained discharge
Ignition of a self-sustained RF glow discharge follows the Paschen curve. A stable thresh-
old He pressure for ignition of a glow discharge can be obtained only by intensive con-
ditioning in the He glow discharge at low RF voltages. The pressure required to ignite
the glow discharge in He before the conditioning at 80 kV is pHe = 0.5 Pa. After the
conditioning the threshold pressures for 80 kV are pHe = 2 Pa for He and pair = 0.18 Pa
for air.
The resonator of the ICRF probe does not allow for matching the high voltage self-
sustained glow discharge. To reach the high voltages required for the glow ignition one
96 CHAPTER 5. RESULTS AND DISCUSSION
1 ms
80
40
0
15
reflected power
[kW]
10
5
100 µs
0
22
30 voltage at the head
[kV]
20 18
10 14
0 10
-82 phase
[degree]
-82
-84 -84
-86 -86
needs the vacuum matching position, at which the ignition of the discharge results in
almost total reflection of the power to the RF generator. The ICRF probe does not allow
for fast change of matching position (see section 4.1). For this reason the transition of RF
glow discharge to arc discharge can not be studied. However when an external plasma
source is switched on, the semi-self-sustained discharge (see section 2.7.4) can be ignited
which can be matched rather well.
At neutral pressures of pHe = 0.15 Pa and pair = 0.03 Pa, and corresponding plasma
densities nHe
e = 9 · 1014 m−3 and nair
e = 2 · 1015 m−3 measured by the Langmuir probe,
the semi-self-sustained RF glow discharge is ignited for the whole range of operational
voltages: from 2 to 30 kV. The maximal achievable voltage is restricted by a mismatch
of the matching circuit and its non-linear dependence on the power (see Fig. 5.4). No
resolvable dependence of the threshold pressure on the voltage applied is observed in the
above mentioned range of voltages. As expected from chapters 2 and 3, the ignition of the
discharge is accompanied by a plasma density increase and leads to an increased negative
charge collected by the inner conductor.
The Q factor of the semi-self sustained discharge is substantially lower than that
in vacuum. Voltage measurements for the He discharge are presented in Fig. 5.4 by
triangles. It is also observed that voltage dependence versus power does not fit the
5.1. TEST FACILITY RESULTS 97
quadratic law (1.11)1 . The quality factor of the resonance at a net power below 50 kW
is about 100 and corresponds to the power dissipation of about 56 % of the net power at
the open end according to equation (4.4). At a power of 300 kW Q ≈ 52 corresponds to
77% power dissipation level. The non-linear dependence on the power is a typical feature
of a glow discharge (see I-V characteristics in section 2.4.2).
The semi-self-sustained discharge has another interesting feature. Cathode spots ap-
pear on the inner electrode at RF voltage higher 5 kV. Furthermore for the voltage lower
than 20 kV (peak) the cathode spots do not evolve into an arc discharge. In Fig. 5.5 mea-
surements for this effect are shown. The measurements include forward power, reflected
power, voltage at the probe head and the phase. The spots in the RF semi-self-sustained
discharge have a characteristic time of formation from 1 to 3 µs. The phase change during
a spot ignition is about 2◦ . This indicates a local character of the change of the increased
discharge current. The system recovers from the first spike (corresponding to the cathode
spot) to the initial value after 30 µs (see e.g. phase measurements in the figure). The last
spike is probably an ignition of several spots which lead to a more pronounced reaction.
However no full arc is developed for the pulse shown in Fig. 5.5. The damage of the
electrodes after the appearance of the spots on the inner electrode is arc tracks (spread
erosion) and can be observed by the CCD camera during the discharge as well as visually
after the discharge.
800 µs
multipactor plasma
200 forward power
[kW]
100
0
20
reflected power Phase pattern of
[kW]
10
the terminating breakdown:
0
3 µs
40
RF voltage
[kV]
-80
[degree]
20
-90
0 phase
0 -100
DC self-bias -1×10
-6 -7
-5×10 0 5×10
-7 -6
1×10
-6
1.5×10
-6
2×10
-400
[V]
Figure 5.6: Experiment with a high voltage multipactor plasma (gray region).
through the inner conductor. This shows that the charged particles are either confined or
generated by the multipactor effect. After the DC voltage reaches about 1 kV a breakdown
occurs. This breakdown is accompanied by a noise pattern on the phase signal indicating
that the breakdown does not occur in the electrode gap, but occurs on the insulation
between the DC cable and the inner conductor at the grounded side.
The experiment illustrates that the multipactor plasma is already a high voltage phe-
nomena when measured at increased neutral pressures with the plasma presence. It also
shows that complications may arise with the parasitic RF power coupled to the DC circuit
when an ICRF system with a varying DC boundary condition is developed.
field emission currents (an increase) is observed in the presence of plasma. In this case the
field emission currents compete with the currents carried to the electrodes by externally
injected charged particles. At neutral pressures higher than the above mentioned values
the semi-self-sustained discharge appears. The latter is a form of RF glow discharge and
is ignited at pressures about one order of magnitude lower than a pressure required for
ignition of a self-sustained glow discharge (discharge without external plasma injection).
The semi-self-sustained discharge can be transformed into an arc discharge by formation
of cathode spots on the surface of the inner conductor if the RF voltage exceeds 20 kV.
At lower voltages the appearance of single cathode spots does not necessarily lead to the
glow-to-arc transition.
plasma structures which may appear as a result of the development of a wide class of drift
instabilities occurring in the SOL plasma [17]. The radial component of the electric field
E is responsible for the poloidal drift of the turbulent plasma structures while the exis-
tence of poloidal component of the electric field moves plasma formations radially. The
existence of a strong poloidal electrical field can be explained by a ∇B polarization of a
plasma formation (also called ”blob” or ”filament”) [102, 103] which is initially detached
from the edge of the main plasma after development of an instability. The presence of
the high electrical poloidal fields in the SOL was treated in [16] where it was shown that
the fields appear in the SOL from different electrostatic plasma potentials of the flux
tubes, the latter being intersected by divertor and limiter surfaces. The different plasma
potentials may exist due to the existence of the space charge sheaths which screen these
potentials from the zero potential of the wall (see chapter 3). This effect is often referred
as a ”finite sheath resistivity”.
100 CHAPTER 5. RESULTS AND DISCUSSION
16 Thomson
Langmuir probes at midplane 10 scattering
manipulator (sector 8): # 15685
movable probes
14
10 0 50 100
R - Rsep, [mm]
limiter position ICRF probe
position
sector 8
Figure 5.7: Thomson scattering and Langmuir probe systems for measuring density and
temperature profile of the edge plasma on AUG. Profile of electron density with respect to
separatrix position.
The turbulence in the SOL results in the convective E × B transport. The turbulent
structures are sometimes also considered as a reason to increase the spatial scale of trans-
port because the perturbed magnetic and electric fields can lead to an enhanced transport
of charged particles (not trapped inside the turbulent structures) across these structures.
An existence of coherent turbulent structures in the SOL and a plasma interchange be-
tween them may further increase the radial transport [101].
Edge localized modes (ELMs) [104, 105, 106] are observed at the plasma boundary
in H-modes. ELMs are initiated by magnetohydrodynamic instabilities and lead to the
very steep increase of the plasma density in the scrape-off-layer. ELMs are registered
practically by every diagnostic system for the edge plasma, for example, Mirnov magnetic
coils, Dα divertor radiation, Langmuir probes, etc.
Therefore the plasma density in the limiter shadow shows very non-stationary behav-
ior. This affects the measured profiles of the edge plasma density and electron temperature
which can be measured with the help of the diagnostics available on AUG.
Nd:YAG-laser beam on the electrons of plasma. Measurements by the system show usually
small gradients of plasma density behind the limiter but the plasma density values are
very scattered because of the small sensitivity of the diagnostic to the typical densities
present behind the limiter.
The second system is a set of Langmuir probes that can be mounted on the midplane
manipulator [108]. Since the ICRF probe and Langmuir probes can not be mounted si-
multaneously on the manipulator, the same configuration of the AUG plasma and heating
power was used to get information on plasma parameters for ICRF probe operation.
Though the measuring systems are situated in different toroidal sectors of ASDEX
Upgrade (Thomson scattering - sector 3, midplane manipulator - sector 8, overall number
of sectors - 16), we assume the axial symmetry of the tokamak with equal conditions and
edge plasma parameters in the different tokamak sectors.
The profile of the electron density measured by the diagnostics is presented in Fig. 5.7
for the discharge with Edge Optimized Configuration (EOC) (see the description in chap-
ter 1), 5 MW neutral injection heating power, 1.1 MW ICRH power, and a line-averaged
plasma density of 9·1019 m−3 . The measurements presented in Fig. 5.7 were made between
ELMs. The measurements of electron temperature in the limiter shadow result in a value
of ≈15-20 eV almost independent of the position in the region 30 mm < R − Rsep < 60 mm.
One can see that in the limiter shadow (R − Rsep > 22 mm) the density profile is
rather flat to about R − Rsep ≈ 50 mm. The plasma density here is modulated by the
intermittent transport by about 50 % (see also [98]). The plasma density for the positions
R − Rsep > 50 mm becomes even more scattered since it is sustained here mostly by
the intermittent transport. Approximately at the position corresponding to the most
frequently used position of the ICRF probe (Fig. 5.7), the averaged plasma density for
this discharge starts to decay further. This is connected to the fact that other components
of the wall start acting as limiters. At the probe position (taking into account about 5 mm
uncertainty of the inner conductor collecting area) the density varies from 2 · 1016 m−3 to
4 · 1017 m−3 .
If one compares the measured densities with the diagram 3.2 from chapter 3 at the
frequency of about 50 MHz, one finds out that the densities correspond mostly to the
region when the sheath thickness is comparable or higher than the 14 mm interelectrode
distance of the probe. This situation is changed when an ELM comes to the outer wall
of the tokamak. During this time the plasma density in the limiter shadow increases by a
factor of 10-100. The sheath thickness during ELMs becomes small compared to the gap.
Therefore the plasma density during ELMs is increased to the values that correspond to
the regime of the thin sheath (s < d, from section 3.3.2) while for the densities between
ELMs, the regime of the thick sheath is mostly realized (s > d, from section 3.3.3).
As an illustration to the intermittent character of the density changes in the limiter
shadow, the measurements by the Langmuir probe with a radial position almost constant
in time are presented in Fig. 5.8. These were performed by H. W. Müller et al. [99] on
the AUG midplane manipulator. The Langmuir probe was placed 3.6 - 4 cm behind the
limiter in a plasma configuration with a distance of about 5 cm between the separatrix
and the limiter.
102 CHAPTER 5. RESULTS AND DISCUSSION
# 15884
2.0 400
Mod_even, [a.u.]
0
40
ϕ floating, [V]
20
0.1
Ion saturation
current, [A]
0.01
0.001
A1
A1
A2
A2
B
no magnetic field with magnetic field
Figure 5.9: Geometrical asymmetry of the probe without and with magnetic field.
5.3. EXPERIMENTS IN ASDEX UPGRADE 103
a) b)
R DC
positive positive
current current
- if the DC circuit is closed and the rectified current is measured, the current directions
are defined as in section 5.1.3 (Fig. 5.10a).
The DC I-V characteristic for the H-mode discharge of the same configuration as for
the discharges used for the density profile measurements (Fig. 5.7) is presented in Fig. 5.11.
104 CHAPTER 5. RESULTS AND DISCUSSION
0.02
Current, [A]
0
-0.02
Figure 5.11: DC voltage-current characteristic for the approximately coaxial probe in mag-
netic field. Gray - experimental data, black - the averaged data.
The characteristic was measured between the ELMs. From the characteristic one can see
that the saturation currents do not differ much. Therefore the characteristic is close to
the characteristic of a double Langmuir probe and the probe asymmetry is rather small
(A1 and A2 are close to each other). The value of negative (ion) saturation current 14
mA is less than the value of the positive saturation current of 20 mA. This indicates that
the inner conductor area A1 is still smaller than the outer conductor area A2 . From the
overview of the asymmetries in chapter 3 one can anticipate that when an RF voltage is
applied to the probe, the DC self-bias voltage will be negative in the case of the opened
external DC circuit. In the case of the closed DC circuit, a negative DC current will flow
from the large electrode to the small one according to the definition in Fig. 5.10a.
The measured characteristic gives an electron temperature of Te ≈ 15 eV in agreement
with conventional Langmuir probe measurements (section 5.3.1). Taking a density of
2 · 1016 m−3 one gets an ion saturation current density of about 17 A/m2 . For the 14 mA
ion saturation current measured by the ICRF probe this results in the effective collecting
area of the inner conductor Aef f ≈4 cm2 . For the cylinder of 1 cm height and 5.5 cm
diameter one has about 35 cm2 . The small measured Aef f is obviously connected with
the limiter configuration of the RF probe and the fact that charged particles are collected
by the probe (a large Langmuir probe) and in this configuration the plasma density is
decreased. The outer conductor plays the role of a limiter since the inner conductor is
placed 6 mm deep inside the outer (see Fig. 4.4).
5.3. EXPERIMENTS IN ASDEX UPGRADE 105
# 15739
2.0 800
-1.0 0
1.0
0.8
0.6
0.4
0.2
0
3.46 3.48 3.50 3.52
time, [sec]
Figure 5.12: Upper signals: Mirnov coils, Dα radiation level from the inner divertor, the
rectified current of the ICRF probe.
broad frequency range. In Appendix D, frequency spectra for the signals of the rectified
current and phase are presented for the type I ELM precursors and type I ELMs. A time
resolved picture of ELMs is described below in section 5.3.7.
a) b)
B B
e e
α
E α
E
i
i
6 mm 6 mm
Figure 5.13: Collisionless movement of charged particles to the front surface of the ICRF
probe. a) case of DC voltage on the probe, b) case of high voltage RF voltage on the probe
at ω0 & ωci .
The ratio is much less than unity for a realistic range of parameters: E < 10 MV/m,
B=1-4 Tesla. However the ion gyroradius at Ti,⊥ <1 kV is smaller than 3 mm, so the ions
do not pass the 6 mm distance to the inner conductor (if the angles α discussed previously
are sufficiently low), and the outer conductor still works as a limiter.
Ions can be collected by the front surface of the inner conductor when higher RF
voltages are applied to the inner electrode (Fig. 5.13b, see also [73]). The perpendicular
ion energy is changed by the application of high RF voltages, because the electrical field is
changing during the period of gyromotion. In this case we can not use the concept of ion
gyroradius nor the concept of the perpendicular temperature Ti,⊥ (see also the discussion
in section 3.4).
If the RF voltage is high enough, one can use an estimation of the ion displacement as a
displacement without magnetic field according to equation (2.20). This kind of estimation
can be applied if the displacement is much larger than ion gyroradius. In this case one
neglects the change of the velocity direction of an ion due to the Lorentz force defined
by the presence of the external magnetic field . The approximation is applicable for a
voltage V0 > 10 kV. Taking the displacement of the electron due to the polarization drift
108 CHAPTER 5. RESULTS AND DISCUSSION
from 5.2, the ratio of electron and ion displacements during half of RF period is:
de,⊥ me ω02
= · 2 1 (5.6)
di,⊥ mi ωci
The ratio is also smaller than unity as the frequency of the RF probe is close to ω0 ≈ 4ωci .
Therefore for both cases we have the perpendicular ion current much larger than the
perpendicular electron current. Furthermore, for high RF voltage the ion displacement
can be of the order of 1 cm (see also Fig. 2.6). This allows ions far from the inner
conductor to be collected by the probe. If ions are collected from a distance of more than
6 mm where the plasma density is higher (the outer conductor does not really work as a
limiter in this case), then the ion flux to the probe should be further increased.
One should also consider that only ions can be collected which have a velocity vi,k
parallel to the magnetic field low enough to have the parallel displacement of ions during
half of RF period smaller than diameter of the inner conductor dinner :
vi,k
li,k = < dinner (5.7)
f0
Assuming that the parallel velocity equals to the ion acoustic velocity vsi , for f0 =50 MHz
and Te =15 eV one can estimate li,k ≈ 0.54 mm for RF period. Therefore the plasma ions
having vsi can be considered motionless during RF cycle in parallel direction if parallel
electric fields are defined only by electron temperature. Also ions with high parallel
energies can be collected: parallel energies below 100 keV result in the displacement li,k
below 4.4 cm. However at high energies one needs to calculate the exact ion trajectories
to account for collected ions.
Experimental observations show that at high RF voltage the background rectified
current (without intermittent events) grows weakly when the RF voltage is increased. This
is explained by rather low densities of the plasma present in the scrape-off-layer between
the ELMs and the events. In the case of voltages above 10 kV most ions within a certain
radial distance (approximately corresponding to the displacements from Fig. 2.6) from the
plasma are gathered by the limited collecting surface - front surface of the ICRF probe.
Furthermore no steep transition happens when the displacements of ions exceed (according
to Fig. 2.6) the radial distance between the inner and the outer electrodes. This indicates
that the ions from the region in front of the probe have more complicated trajectories
which cannot be quantitatively described in terms of the displacements calculated without
magnetic field.
Thus the existence of the positive ion current can be explained by the fact that the
front surface of the inner electrode collects predominantly ions. Furthermore the ions can
be collected from the regions of higher densities. Therefore the time-averaged current
is higher than the negative contribution of the rectified current expected from the DC
measurements.
5.3. EXPERIMENTS IN ASDEX UPGRADE 109
4 ms
# 15515
a.u.
Dα , outer divertor
Dα , inner divertor
0.4
0.0
3.552 3.553 3.554 3.555
120
80 reflected power
[kW]
40
0 3.552 3.553 3.554 3.555 10 µs
60 60
[kV]
40 40
voltage at the head
20 20
0 0
3.552 3.553 3.554 3.555
-120 0
phase
[degree]
-40
-122
-80
-124
3.552 3.553 3.554 3.555
time, [sec] time to trigger, [sec]
a) b)
no ELMs in time
breakdown voltage, [kV peak]
conditioning
in vacuum
Figure 5.15: Operation in vacuum and plasma. Stars - maximum voltage limited by
plasma, circles - voltage limit in vacuum. a) Voltage limit versus voltage exposure param-
eter U ; b) voltage limit in plasma versus shot number (after torus opening).
In Fig. 5.14 the influence of type I ELM in the discharges with EOC plasma shape
(see chapter 1 and section 5.3.1 for the description) is presented. The positive rectified
current can rise up to few amperes during the type I ELMs. The phase transient of an
arc after the ELM is very similar to that during the cathode spot formation in the semi-
self-sustained RF glow discharge (Fig. 5.5) and develops in a characteristic time of 2 µs.
We consider the first RF breakdown and corresponding maximal voltage before ELM to
be the voltage limit of the system.
In Fig. 5.15a the voltage stand-off in AUG discharges (all breakdowns correlate with
type IR ELMs) and in vacuum is shown. We define the voltage exposure parameter as:
t
U = 0 b V0 · dt, where tb - time of the breakdown counting from the ramp start. The
parameter U is shown on the x-axis of Fig. 5.15a.
The fast power ramp duration of 50 ms used for the voltage limit measurements is of
the order of only a few typical ELM periods. Thus the voltage limit associated with the
ELM appearance can be overestimated, because it is probable that there is no ELM at the
time when the voltage on the probe corresponds to the limit. The second, slower power
ramp is better suited for the measurements of ELM-correlated RF breakdown, though
additional measures need to be applied to control gas desorption in the vacuum lines.
For the majority of points in Fig. 5.15a the early measurements correspond to lower
breakdown voltages. The conditioning effect in vacuum can be observed in the beginning
of each day of operation. Between the pulses of the ICRF probe in AUG plasma, the
voltage stand-off in vacuum was sustained at the level higher than 75 kV by intense high-
voltage conditioning between AUG discharges. One observes the conditioning effect for
the maximal voltage in the presence of the AUG SOL plasma from shot to shot. The
5.3. EXPERIMENTS IN ASDEX UPGRADE 111
45 50 55 60 65 70 75 80
R - R sep , [mm]
Figure 5.16: Dependence of the voltage limit of the ICRF probe conditioned well in plasma
versus radial position R − Rsep .
conditioning effect with the plasma is more emphasised after the AUG torus opening
(Fig. 5.15b). Therefore the surface state is important for the breakdown caused by ELM
activity, and the intense conditioning of the probe in vacuum is a necessary but by far not
sufficient measure to improve the voltage limit in the presence of plasma. Conditioning
with plasma is required.
Variations of the distance between the well-conditioned ICRF probe head and the
separatrix do not affect the voltage stand-off considerably if the breakdowns on the probe
are triggered by ELMs. In Fig. 5.16 the dependence of the voltage limit associated with
ELMs versus the radial position of the probe with respect to the separatrix is shown.
For the specific experiments (see below) only AUG discharges with EOC plasma shape
are used.
Figure 5.17: Affecting the voltage stand-off by controlling the rectified current (the direc-
tion of the current is defined as in Fig. 5.10). Stars - voltage of the breakdown correlated
with an ELM, gray regions - uncertainty of the voltage limit due to the finite ELM fre-
quency, triangles - rectified current averaged between the two last ELMs.
noticeably higher voltages for the same value of electric field, hence an increase of the
power that can be launched by the ICRF antenna (voltage increase from 55 kV to 63 kV
corresponds almost to 30 % increase in a possible power launched by the antenna).
Thus the influence of the DC biasing on the voltage stand-off stresses the role of the
electrode surfaces which face the plasma and are oriented nearly parallel to the magnetic
field. A decrease in the area of these surfaces can lead to a reduction of the number of the
ions displaced by strong radial RF electric fields towards the antenna. If these ions are
not collected by the antenna electrodes they oscillate in the radial electrical fields and can
easily move along magnetic field lines during the time of several RF periods. Therefore
by a reduction of the number of these ions, the density in the gaps between the electrodes
oriented not parallel to the magnetic field can be reduced. The reduction of the plasma
density favors the increase of the voltage stand-off.
The RF breakdown occurs during the injection of ELMs - very unstable, intermittent
events (see below section 5.3.7). Each ELM is unique and introduces a high degree
of uncertainty. Therefore it is still unclear whether the effect of the change of voltage
limit is assigned to some changes of the fluxes during single ELMs that may provoke the
breakdown or the voltage limit changes as a result of a change in total number of charged
particles collected by the surface of the probe during the power ramp.
a) b)
1.0 200
140 # 15430
120
# 15515 150
# 15515 0.8 # 15739
100
80 0.7 100
60 # 15431
0.6
# 15430
40 50
0.5
20 # 15739
0.4 0
0 10 20 30 40 50 60 80 0 50 100 150 200 250 300 350
RF voltage, [kV peak] Net power, [kW]
Figure 5.18: a) Minimal value of time-averaged quality factor during ELMs versus RF
voltage on the probe between ELMs; b) the coupled power versus the net power, the solid
lines – the coupled power normalized to the net power, the dashed lines – the absolute
coupled power.
The increase of the Q-factor (and the decrease of the normalized coupled power) with
voltage can be explained as a consequence of the following effects:
1. The power coupled to the plasma may include a certain fraction corresponding to the
power launched by the waves. The power can have a non-linear relation to the RF
voltage since a variety of wave modes can be launched by the probe (including waves
that can be under some conditions transformed into the non-linear ion Bernstein
mode in front of the probe [109, 110]).
2. The RF probe reduces the particle flux during a part of the ELM phase at high RF
voltages as described in section 3.3.3 for the case of the thick sheath.
Clarification of the role of the excited waves that propagate from the probe inside the
torus chamber and an associated change in the RF resistance of the open end requires
taking into account the variety of the wave types that can be excited by the probes.
Theoretical treatment is not considered in the frameworks of this thesis. The experimental
observations on the propagation of the waves are briefly discussed below in section 5.4.
The effect of the density reduction described in sections 3.3.3 and 3.3.4 gives a con-
sistent explanation of the picture of the coupled power presented in Fig. 5.18. Both the
regime of the thin sheath and the regime of the sheath thickness comparable with the
electrode gap take place during ELMs. The regimes are present during the ELM due to
the intermittent character of the radial transport of plasma during this phase. The thick
sheath regime will lead to a reduction of the charged particle fluxes towards electrodes
and a decrease of the power fraction which is dissipated locally.
Thus the change of the Q-factor can be related to the mechanisms of the local power
dissipation, i.e. the interaction of charged particles with the RF fields and electrodes of
the probe.
5.3. EXPERIMENTS IN ASDEX UPGRADE 115
55 ms
30
# 15430
[kV peak]
voltage,
20
10
1 ms 3.5 ms
1.0 0 1.0
0.2
0.4 0.4
0
3.39 3.40 3.41 3.42 3.43
0.2 0.2
Time, [sec]
0.0 0
-80 3.428 3.429 3.430
Phase, [degree]
Time, [sec]
-84
-88
a) b)
3 2
Ptransmitted 85 -48
80 -53
6 15
75 -58
7 14
400 Dα , outer divertor
[a.u.]
200
ICRF probe
(sender) 8
Pdissipated 13
0
Dα , inner divertor
400
[a.u.]
α
Pnet 9 12
200
0
10 11 3.46 3.48 3.50 3.52 3.54
time [s]
Figure 5.20: a) Scheme of the measurements of the coupled and the transmitted power; b)
measurements of the power transmitted from the ICRF probe to AUG ICRF antenna.
250 500
Q-factor, [a.u.]
150 300
100 200
50 100
0 0
60
8
[kV]
40
[A]
RF voltage RF current 4
20
signal 0 0
saturation 0.8 rectified
current
[A]
0.4
0
3.514 3.515 3.516 3.517 3.518
time, [sec]
Figure 5.21: The coupled power measurements, AUG shot #15515 (see also Fig. 5.20).
mission signal above a noise level of -60 dBm (about the detectable power threshold) can
be detected by a frequency analyzer with a time resolution of 500 µs.
by a mismatch during ELMs. Indeed the net power signal, which is on the same relative
scale in Fig. 5.20b as the transmitted power, only decreases by ≈3 dB (by a factor of
≈2). However, the relevant power, which is the absolute power coupled to the plasma,
as we will see below, actually increases. A mismatch of the directional coupler of the
ICRF antenna used as the receiving probe can also be ruled out. Between the ELMs one
observes a scattering of the value of the transmitted power by more than 10 dBm. This is
larger than the scattering of the power in vacuum (2 dBm). Thus the transmitted power
is very sensitive not only to the ELMs but to any changes at the edge plasma, which can
be qualitatively observed on the Dα signal from the inner divertor tile.
In Fig. 5.21 the Q-factor, the coupled power calculated according to (4.4), the RF
voltage, the RF current and the rectified current are shown during the ELM phase. The
figure shows a strong increase of the coupled power of up to 200 kW during an ELM. The
power increase is confirmed by changes of the RF voltage and of the RF current as well
as by an increase of the rectified current.
The measured reduction of the transmitted power, despite the increase of the coupled
power, indicate that a larger fraction of the power is absorbed during ELMs between the
sender (ICRF probe) and the receiver (ICRF antenna). This absorption can occur either
at the plasma edge or in the plasma center.
Two hypothesis can explain this reduction of the transmitted power:
1. Most of the power is transmitted through the edge and the absorption in the edge
changes (increases during ELMs).
2. Only a fraction of the transmitted power comes through the edge and another
through the center. If now the change of the edge parameters (during ELMs)
increases the fraction of the power traveling through the center (where it can be
efficiently absorbed), then a reduction of the transmitted power can occur even
without an increase in absorption of the power in the edge.
A drastic increase of the rectified current during ELMs supports rather the first hy-
pothesis: the power dissipated in the edge is increased. Furthermore this increase of the
rectified current indicates that the power is absorbed locally, near the head of the ICRF
probe. A rough estimation of the local power absorption at the plasma boundary can be
made. We assume that the RF current carried by charged particles is of the order of the
rectified current, i.e. we neglect the secondary electron emission effects (section 2.6.5) and
the rectified current component due the asymmetry of the surfaces not nearly parallel to
the magnetic field (sections 3.3.2 and 3.3.3). We assume also that ions are accelerated to
energies of the order of eV0 . An estimation for a maximum value of the power dissipated
locally in a time window from 3.514 to 3.515 seconds from Fig. 5.21 results in a value of
' 30 kW. The measured maximal coupled power in this time window is ' 80 kW. Thus
the estimation confirms that a significant fraction of the absorbed power is dissipated
locally, though we cannot make a conclusion that the power does not dissipate in the
central plasma.
The hypothesis of the local power dissipation fits also well to the measurements pre-
sented in section 5.3.7 where the coupled power normalized to the net power was found
to be affected by the high RF voltage on the probe.
5.5. CONCLUSIONS FROM THE ICRF PROBE EXPERIMENT IN AUG 119
As the antenna strap length is comparable to a quarter of the vacuum wavelength, the
RF voltage on the ICRF antenna varies along the strap. The voltage is zero at the
position of the short circuit in the AUG antenna - the most distant position from the
RF generator. The voltage increases along the antenna strap to a maximal value at the
position corresponding to about λ/4 from the short circuit. The strap is positioned in
poloidal direction, therefore the whole range of voltages from zero to the maximal value
is present at approximately the same radial position close to the SOL plasma.
In the ICRF probe only one RF voltage value instantaneously is exposed to the SOL
plasma. The voltage changes according to the standing wave pattern radially away from
the plasma.
The ICRF antenna is designed and operated such as to deliver as much power as possible
to the plasma. A coupling resistance R0 = 2 Ω determines a Q factor of Qant = 12.5.
For this Q-factor a voltage of Vant = 25 kV on the antenna strap corresponds to a power
of Pant = 2 MW (see Fig. 5.4). On the other hand in vacuum the Q-factor is high
(Qvac
ant ≈ 200). Therefore, according to equation (4.4) most part of the input power is
transferred to the plasma.
The probe aims to expose high voltages for a relatively low input power. Therefore it
the has high Q for operation with and without plasma. The value without load Qvac probe =
230 (see section 4.1) is close to Q of the antenna operated in vacuum. Experiments
described in section 5.3.7 showed that the Q-factor of the ICRF probe during operation
with the plasma stays close to 200 and is down to the values of 20-100 during ELMs. This
is considerably higher than Qant .
The relative change of the RF voltage on the probe or antenna during fast changes of
the load (e.g. caused by ELMs) may be critical in terms of comparison of the breakdown
evolution in the RF systems. The relative change of the voltage can be represented as:
r s r
Z0 Pnet Q
δV = δQ + δPnet (5.8)
2 Q Pnet
5.6. HIGH VOLTAGE OPERATION OF THE AUG ICRF ANTENNA 121
The change in voltage is associated with the change in the Q-factor δQ (first term) and
the mismatch in terms of the net power δPnet (second term). Generally δPnet is a function
of the resonant line parameters (the absolute Q value), a change of the load impedance,
matching position and matching circuit characteristics. The value of δQ is only dependent
on the change of the load. One can use experimental data to test the validity of the last
expression.
Let us estimate a difference in δV for a voltage of V =25 kV and changes δPnet and
δQ during the type I ELMs (which do not trigger arcing) in the same discharge, for the
radial position of the ICRF probe close to that of the antenna.
For the ICRF probe: Q ≈ 200, Pnet = 200 kW, δPnet ≈ 100 kW, δQ ≈ 100.
For the ICRF antenna: Q ≈ 12, Pnet = 1500 kW, δPnet ≈ 400 kW, δQ ≈ 5.
As a result one gets a relative change of the voltage on the probe 2.1 times higher than
on the antenna. The difference in δV between the probe and the antenna is confirmed by
the voltage measurements. The ratio δVprobe /δVantenna normally lies in the range 1.5-4 for
high voltages depending on the discharge parameters.
For the EOC plasma shape the voltage on the probe during ELMs can drop from 50
kV to about 25 kV, while on the antenna the voltage decreases from about 23 kV (for the
case of Pnet ≈1.8 MW per antenna) to 15 kV.
Therefore, during ELMs, the ICRF probe may operate at higher voltages than the
AUG ICRF antennas. The Q-factor of the probe is high during ELMs (as shown in
Fig. 5.18, section 5.3.7, Q grows with the RF voltage) compared to Q of the antenna
(Qant tends to 1 due to the transiently increased R0 ). This explains the relatively high
RF voltage on the probe during ELMs though the mismatch δPnet /Pnet of the probe can
be large compared to the mismatch of the antenna.
We conclude that at high RF voltages on the probe between ELMs (>50 kV), the
RF voltage on the probe is higher than 20 kV during type I ELMs. If an RF breakdown
correlates with type I ELMs, the breakdown processes initiated at this RF voltage are
similar to those on the ICRF antenna which operates near its voltage limits (25 kV
between ELMs and 15 kV during ELMs) under the same plasma conditions.
a) breakdown
b)
breakdown
antenna 4 antenna 3 correlates with ELM correlates with ELM
antenna
#15049 quasi-stationary
arc discharge
quasi-stationary
arc discharge
20
voltage, [kV]
15
10
0
stub
tuners
power, [kW]
800
forward
400
splitter
0
dummy 800
reflected
power, [kW]
loads
power for
reflected
400
combiner arc detection
~ ~ 0
2 3 4 5 6
RF generators time, [sec]
Figure 5.22: a) Scheme of the ICRF system of two antennas at AUG, b) time traces od
RF voltage, forward and reflected power for a shot with breakdown associated with ELMs.
in loads between two antennas is rather small, arc usually appears on a single antenna
and should lead to an asymmetry of the loads and substantial increase of the measured
reflected power. However two antennas are not quite equal and one needs to regulate the
critical level of the measured signal to be sure it is due to an arc, but not to the ELMs
that do not trigger an arc.
1. vacuum arc breakdown which happens at the location of the voltage maximum on
the strap and does not correlate with a plasma activity;
3. breakdown at the location of the strap at a certain distance from location of the
voltage maximum, triggered by ELM, poorly detected by the arc detection system
and leading to a quasi-stationary (from 100 µs to 1 s duration) arc discharge.
The type (1) and type (2) breakdown have about the same scales and characteristics
as vacuum arc and ELM-correlated breakdown on the ICRF probe respectively. The type
(1) is the same breakdown as in vacuum without the plasma. The type (2) breakdown
is a RF breakdown caused by a drastic increase of the surface electric field during the
increase of the plasma density in the gap between the antenna strap and the antenna box.
The type (3) breakdown is followed by a quasi-stationary arc discharge. Appearance
of the latter on the antenna is shown in Fig. 5.22a by the gray vertical lines. Our un-
derstanding of this phenomena is the following. The breakdown is associated with the
ignition of an arc at a location far from the voltage maximum since the reflected power is
increased to the values lower than the typical values for the type (1) and type (2) break-
downs. Ionization of the desorbed gas by the external plasma injection during ELMs leads
to formation of the semi-self-sustained discharge (see sections 2.7.4 and 5.1.6). The latter
increases the power absorption (it is a quite non-local effect) which is poorly detected
by the arc detection and a consequent glow-to-arc transition happens at the location not
necessarily coinciding with location of the voltage maximum. The sequence of the events
may be ”transparent” for the arc detection system used and described in the previous
section.
The arc discharge phases which follow type (3) breakdown are shown by the gray
regions in Fig. 5.22a. If the breakdown is not detected at the right time, an arc discharge
occurs which leads to a reduced voltage stand-off of the system. Repeated breakdowns of
type (1) which are detected (the power is switched off often, Fig. 5.22a) follow the quasi-
stationary arc discharge. These breakdowns seem to recondition the system. A good
detection of the ELM-correlated breakdowns is critical since the overall performance of
the RF antennas is affected by the significantly reduced voltage stand-off after ignition
of these quasi-stationary arc discharges. The conditioning of the antenna conductors
124 CHAPTER 5. RESULTS AND DISCUSSION
and transmission lines required after the ignition of the arcs can be partly performed in
vacuum operation.
The ignition of the arc discharges after ELMs is less frequent when the antennas are
well-conditioned. However even few arc discharges may lead to damage of the electrodes
which may set permanent voltage limits of the system lower both for vacuum and plasma
operation.
It is important to understand the movement of the cathode spots of the quasi-stationary
sustained arc discharge. Motion of the cathode spots sustained by RF power should be
treated differently from the motion of spots constituting a DC arc.
2. ignition of the two equivalent other spots to the different sides of the initial spot;
3. the right spot is screened by the left ones and disappears because the power is
dissipated on the left side or is reflected back to the generator before it reaches the
right spot;
1) Z
spot
RF generator antenna
transmission line
~ 2)
Z
arc
3)
4)
Figure 5.23: Explanation of the arc self-screening, different stages of the arc evolution.
The total power which the external circuit can deliver sets the number of spots. The
spots appear and disappear constantly and predominantly the spots to the generator side
survive, so we see the movement of the arc towards the generator.
interelectrode distance;
- sustaining of arc in magnetic field : in a strong external magnetic field the arc tends
to have the current column aligned parallel to the magnetic field;
- movement in magnetic field : in the external magnetic field an arc can move in a
direction perpendicular to the field and to the arc current.
The last point requires a separate discussion. Even for a DC sustained arc the picture
of the arc motion in a transverse magnetic field is not well understood. The motion can
change the direction from the retrograde3 [112, 113] to the amperian [113] when, e.g.
the magnetic field is increased for a given interelectrode distance and a given neutral
pressure (we consider a low-pressure arc). For example, in the experiments [113] the
critical magnetic field for the reversion of the arc motion was found to be about 1 Tesla.
If the magnetic field has an acute angle to the surface, the arc velocity can sometimes
also have a component in the direction of the angle. Sometimes the movement associated
with this velocity component is called Robson drift [114].
For an arc sustained by the RF power one can initially consider both the RF and the
DC components of the arc current when the arc motion in the magnetic field is discussed.
However, if we use the results of the experiments (e.g. [111, 112, 113]), the maximal
velocities (perpendicular to the magnetic field) a DC arc can have, are of the order of
100 m/sec. If we take the values of the velocity from these DC measurements for an arc
with a pure RF current, then the arc does not move to a distance larger than the width
of the current channel in half of a RF period. Therefore if a mechanism of the RF arc
motion is assumed to be similar to that of the DC arc motion, one should consider the
DC component of the arc current as the first candidate which can explain the motion.
A counteraction of the properties of the arc motion described in the list above and
the self-screening effect described in the previous section, can lead to an arc tied to
the particular regions of the electrodes due to, for example, a specific geometry of the
conductors and/or due to a specific distribution of the magnetic field.
In addition to the poorly understood picture of the currents in a RF sustained arc,
other factors, such as complex geometry, uncertain finite ferromagnetic properties of the
stainless steel electrodes (even more so in the regions of welding of the coaxial tube to the
plate of the antenna return conductor) make the interpretation of the arc tracks observed
in the antennas an outstanding issue.
3
The term ”retrograde” is applied to the arc motion in a direction opposite to the amperian (parallel
to (j × B)) direction.
5.7. CONCLUSIONS FROM STUDIES ON AUG ICRF ANTENNAS 127
- arc initiated due to a high surface electric field when plasma presence in an electrode
gap amplifies external electric field on the surface of the RF electrodes (both for the
AUG ICRF antenna and the ICRF probe in AUG during plasma injection which
accompanies ELMs);
- semi-self sustained glow discharge with consequent glow to arc transition (ICRF
probe in test facility and sometimes in the AUG ICRF antenna during type I ELMs);
- more frequent vacuum arc ignition after electrode damage by preceding quasi-
stationary phases of arc discharges (AUG ICRF antenna).
One of the most frequently observed phenomena is the breakdown following the tran-
sient change of the plasma parameters in the scrape-off layer of the tokamak, e.g. ELMs.
ELMs and more local intermittent events which can be seen by the RF probe are ac-
companied by the drastic increase of ion flux to the probe. Therefore the plasma density
is increased during the events in the electrode gaps of the ICRF probe and ICRF an-
tennas. This leads to amplification of the surface electric field and consequent arcing.
For ICRF antenna operation, ELMs can sometimes lead to the formation of a semi-self
sustained discharge involving gas ionization. Therefore ELMs serve as a trigger for the
self-sustained discharge arc discharge both for the ICRF probe and the ICRF antennas.
The arc discharge in AUG ICRF antennas can sometimes be poorly detected and leads
to a need of reconditioning of the electrodes.
A high increase of the power coupled to the plasma and a good absorption of this power
(the power transmitted from the ICRF probe to the ICRF antenna functioning as a pick-
up probe decreases essentially) is observed during the ELM phase. Simple estimation
128
6.2. MEASURES TO IMPROVE THE VOLTAGE STAND-OFF 129
shows that the power can be efficiently absorbed at the plasma edge. Therefore the role
of ELMs and intermittent transport in the SOL should be considered when parasitic
absorption of RF power coupled to the plasma by the ICRF antennas is studied. The
measurements of the averaged Q-factor of the probe during ELMs show that the fraction
of coupled energy normalized to the net power is significantly reduced with an increase
of the RF voltage. This is a further confirmation that the RF voltage affects the charge
particle flux which is essential for a RF breakdown.
The maximal voltage on the ICRF probe in the presence of plasma increases by con-
ditioning the probe with plasma. Despite the intensive conditioning in vacuum, plasma
operation is a necessary measure to bring up the voltage limit in the presence of the
plasma. The additional mechanisms connected to the presence of plasma and accompany-
ing charge particle fluxes, such as gamma-processes on the electrodes, particle stimulated
desorption, and an eventual increase of the electrical fields on the electrode surface, make
the operation with plasma very critical if the electrodes are not trained in these condi-
tions. The operation of AUG ICRF antennas becomes also much more reliable when the
antennas are well conditioned in the presence of the AUG plasma. The existence of the
conditioning effect in plasma stresses the role of the surface state of the electrodes. The
surface state has consequences both for vacuum and plasma operation of the probe. The
same holds true for the antennas.
It was shown experimentally that the voltage stand-off of the well-conditioned ICRF
probe with the plasma presence can be affected by regulating the time-averaged ion flux
into the probe.
In parallel to studies related to the voltage limits, measurements of the rectified current
and of the power transmission and dissipation provided a better understanding of the
interaction of high voltage RF fields with the plasma boundary. The measurements are
also a monitor for the transient changes in plasma edge parameters.
- screening of the antenna from the external plasma: use of an optically closed Faraday
shield;
- regulate particle fluxes to the antenna electrodes facing the plasma: improvement
of the geometrical shape of the antenna (more symmetrical electrodes with surfaces
intersecting magnetic field at high angles, electrodes of less area intersecting the
surface at low angles) and design a DC biased ICRF antenna;
130 CHAPTER 6. SUMMARY AND CONCLUSIONS
- reduce number of arc detection failures: improve the reliability of the RF breakdown
detection system.
time-averaged potential of plasma existing between the surfaces of the antenna strap and
the grounded vacuum vessel can be rather low (see table 3.3.2). The strap is charged to a
potential which is the result of a trade-off of the charge brought by the sheath rectification
(see 3.2) and the ion charge brought perpendicular to the magnetic field (see 5.3.3). Since
the strap is equipotential the time-averaged potential exists on any location of the strap.
In the case with the closed DC circuit the plasma between the strap and the vacuum
vessel is charged to a relatively high potential. The surfaces which are intersected by
magnetic field lines at small angles collect ion current. The plasma is strongly charged
locally only due to the local sheath rectification. Therefore the distribution of DC electric
field is completely different.
To our knowledge there exists no data about these two configurations working in the
similar conditions. The ICRF probe did not allow for experiments with the opened circuit
at AUG because of the technical restrictions connected with very high self-bias voltages.
However an effort should be made to study effects of DC boundary condition and electrode
geometry on plasma convection inside and in front of the ICRF antenna.
Appendix A
Abbreviations
RF Radio Frequency
DC Direct Current
ICRF Ion Cyclotron Range of Frequencies
SOL Scrape-Off-Layer
ICRH Ion Cyclotron Resonance Heating
AUG ASDEX (Axially Symmetric Divertor EXperiment) Upgrade
EOC Edge Optimized Configuration (describes a type of the plasma shape in AUG)
ELM Edge Localized Mode
ESD Electron Stimulated Desorption
ISD Ion Stimulated Desorption
FFT Fast Fourier Transform
132
Appendix B
A heat conduction equation is used to describe the temperature rise on the surface of the
antenna strap due to Joule heating by RF currents. Let us assume that the strap surface
corresponds to (x = l) and a value of (0 < x < l) corresponds to a position inside the
antenna strap (Fig. B.1). The RF current with peak value I0 (peak) flows perpendicular
to x on the surface of the strap and decays exponentially inside the conductor (when l − x
is increased). The heat conduction equation with the volume heat sources is written for
temperature T as:
∂T ∂2T
=κ , (B.1)
∂t ∂x2
where κ = λT /(ρm · c), λT is specific heat conductivity, c is heat capacitance and ρm is the
mass density of the strap material. RF current flows in the skin layer δ. If we assume that
the skin layer is small (δ l), the heating of the strap can be assumed to be provided by
a heat flux to the strap surface:
Q0 I02
q0 = , Q0 = (B.2)
κ 2h2 δ σ
where δ is the skin layer, σ is the electrical conductivity, h is the strap width. We do not
account for losses by a radiation from the surface of the strap nor for the temperature
dependence of the electrical conductivity (these two effects partly compensate each other).
A boundary condition at the surface (x = l) is therefore:
∂T
= Q0 (B.3)
∂x x=l
If the strap has a thickness of 2l and is symmetrically heated from both sides, the boundary
condition of heat insulation can be applied at the middle of the strap thickness:
∂T
=0 (B.4)
∂x x=0
133
134 APPENDIX B. HEAT CONDUCTION IN THE SKIN LAYER
x
q0
T
0
2 ∂T
=0
∂x
q0
h
Figure B.1: Model of a symmetricaly heated antenna strap. Dashed curve - qualitative
temperature distribution inside the strap.
Table B.1: Properties of materials used for the calculation of the temperature rise in the
antenna strap due to the Joule heating.
Let us set the temperature T0 constant across the strap as the initial condition:
T (x) = T0 (B.5)
t=0
The equation B.1 with the boundary conditions B.3, B.4 and initial condition B.5 is
solved by the method of Laplace transforms [116]. The solution is written:
∞ nπ x
Q0 t Q0 l 3x2 − l2 2 X (−1)n −κn2 π2 t/l2
T (x, t) = T0 + + − 2 e cos (B.6)
ρm c l λ 6 l2 π n=1 n2 l
Since the real heat source is distributed in the skin later on the surface, the solution B.6
does not reflect the real temperature profile in the surface layer. This is even more
important as we are mainly interested in the surface temperature. However the error in
the solution is not big and can be estimated.
If one accounts for the volume heat sources instead of using the external flux Q0 on
the surface, the surface temperature should decrease. The estimation of a maximal error
can be done assuming a flat temperature profile in the ”real” skin layer. In this case the
maximal error is:
1 ∂T Q0 δ I2
|∆T | ≈ · δ = = 2 0 (B.7)
λ ∂x x=l λT h σ λT
The maximal error is for the material with the lowest thermal and heat conductivity.
For the working parameters I0 =1 kA, h=0.2 m, stainless steel ∆T ≈ 0.83 K which can
be considered negligible.
The material properties from the table B are used for calculations. The results of the
calculation of the surface temperature are presented in chapter 2 for the ICRF antenna
and in chapter 4 for the ICRF probe.
Appendix C
1 kTn
len = = (C.1)
nn σen p σen
4.18 · 10−21
le , [m] = , p in Pa, σ in m−2 . (C.2)
p σen
For very low electron energies (εe <30 eV) the experimental data for total cross sections
of collisions of electrons with neutrals is available [64, 65]. The table C presents the
electron mean free pass in gases.
Table C.1: Ionization potentials and mean free pass of electrons colliding with neutrals at
1 Pa at low energies (<30 eV). a) - taken from elastic collisions, other - taken from total
collision probability.
135
136 APPENDIX C. DATA FOR COLLISIONAL AND IONIZATION PROCESSES
a) b) c)
-19
10
Ar
cross-section, [m ]
-20 Ar
2
10
H He
p+H → 2p+e
-21
10 H2
plus ionization
He through 2p excitation
-22
10 + p+He →
p+H2→ p+H2+ e +
H p+He + e
-23
10
1 10 100 1000 10000 1 10 100 1000 10000 1 10 100 1000 10000
Figure C.1: Cross sections of charge particle - neutral collisions: a) elastic cross-sections
estimated from eq. C.3; b) total electron impact ionization cross-sections, crosses and
triangles - averaged experimental data for H2 , He, Ar from [66], curves for H, H2 , He -
numerical fit from [118], Ar - from Lotz fitting equation C.4 with N = 2; c) proton impact
ionization for H, H2 and He from [118].
At the low electron energies the total collision cross-section is determined by elastic
collisions. For an electron energy εe >30 eV one can estimate the cross-section of elastic
scattering of electrons on an atom with nuclei charge Zn by the approximate formula:
4/3
−2 −19 Zn
σen , [m ] ≈ 6 · 10 (C.3)
εe
The last equation can be used to estimate elastic collision cross-section for the high
electron energy range.
For the estimation of the ionization cross-section of Ar an empirical fitting formula is
used [117]:
N
X
i ln(εe /Viion ) ion
σen = a i qi 1 − bi exp − ci (εe /Vi − 1) (C.4)
i=1
εe Viion
where i – number of the electron shell counting from the outermost one, ai , bi , ci – tabu-
lated empirical constants from [117], Viion – ionization potential, qi – number of electrons
in the (sub)-shell. The last formula gives a reasonable agreement with experimental data
and the gases listed in the first column of the table C and energies less than 10 keV.
However we use it to extrapolate the cross-section to 20 keV.
Cross sections for elastic collisions estimated by the equation C.3 are shown in Fig C.1a.
Electron impact ionization cross sections are presented in Fig C.1b using the data
from [66, 118] and the fitting formula C.4.
Proton impact ionization cross sections from [118] are found in Fig C.1c. In addition
to the information on the fugure we note that the cross section of ionization of atomic
hydrogen by proton impact has a maximum of 1.5·10−16 m2 at an energy of 50 keV [119].
Appendix D
Fast Fourier transform (FFT) analysis was performed for two signals: signal of the phase
between RF voltage and current, and the rectified current. The analysis allows to look
for the frequency spectrum of the signals.
The phase and the rectified current were averaged for the shots #15515, #15538,
#15539 for the last 4 ms before the breakdown. Since the breakdown correlates with an
ELM, the signals reflect the frequency spectrum of the ICRF probe reaction to the ELM
pre-cursors and intermittent events of higher frequency.
The results are presented in Fig. D.1. An amplitude of the FFT analysis for the phase
signal is shown in Fig. D.1a. To get a noise level for the phase signal three RF pulses
were made in vacuum and the FFT spectrum was averaged. The noise level of the phase
deviations is shown by a gray curve in Fig. D.1a. A power spectrum of the rectified current
is presented in Fig. D.1b. Both spectra of Fig. D.1 show a flat signal distribution over the
range of frequencies presented in the figure. The measurements of the rectified current
are limited in speed by the fiber-optics based signal transfer system (see the section 4.3.2)
with the 128 kHz bandwidth. Therefore artifacts are observed for the higher frequencies.
A spectrum of the rectified current averaged for shots #15430, #15431 is shown for a
wide range of frequencies in Fig. D.2. The frequency range includes also ELM frequencies
(150-300) Hz for the discharges.
The observations by the ICRF probe confirm the intermittent character of the radial
transport in the scrape-off-layer. Similar frequency characteristics are observed by many
diagnostics at the plasma edge of different machines (see discussions in [94, 95, 96, 97, 98,
100, 103]). Therefore the ICRF probe can be used as a tool to monitor the plasma edge
activity.
137
138 APPENDIX D. FREQUENCY SPECTRA OF THE ICRF PROBE SIGNALS
a) b)
bandwidth limit
0.1000
Amplitude, [a.u.]
Amplitude, [a.u.]
0.0100
0.0010
0.0010
noise
0.0001
0.0001
3 4 5 6 3 4 5 6
10 10 10 10 10 10 10 10
Frequency, [Hz] Frequency, [Hz]
Figure D.1: Frequency characteristics of the phase signal (a) and the rectified current (b)
collected for 4 ms before the breakdown averaged for the shots #15515, #15538, #15539.
bandwidth limit
0.0100
Amplitude, [a.u.]
0.0010
0.0001
2 3 4 5 6
10 10 10 10 10
Frequency, [Hz]
Figure D.2: Frequency characteristics of the rectified current collected for 110 ms during
the power ramp averaged for shots #15430, #15431.
Bibliography
[1] Lawson J. D., Some criteria for a power producing thermonuclear reactor, Proc. Phys.
Soc. B70, 6-10, 1957.
[2] Dolan T. J., Fusion research, Pergamon Press, New York, 1982.
[3] Wesson J., Tokamaks (second edition), Clarendon Press, Oxford, 1997.
[4] http://fusedweb.pppl.gov
[5] Rosen M. D. The physics issues that determine inertial confinement fusion target
gain and driver requirements: A tutorial, Phys. Plasmas 6 (5), 1690-1699.
[6] Lindl J., Development of the indirect-drive approach to inertial confinement fusion
and the target physics basis for ignition and gain, Phys. Plasmas 2 (11), 3933, 1995.
[7] Artsimovich L. A., Tokamak devices, Nucl. Fusion 12, 215-252, 1972.
[8] Spitzer L., The stellarator concept, Phys. Fluids 1, 253-264, 1958.
[9] Volkov E. D., Suprunenko V. A., Shishkin A.A., Stellarator, Naukova Dumka, Kiev,
1983, in Russian.
[10] Greenwald M. et. al., A new look at density limits in tokamaks, Nucl. Fusion 28
(12), 2188-2207, 1988.
[11] Greenwald M., Density limits in toroidal plasmas, Plasma Phys. Control. Fusion 44,
R27-R80, 2002.
[12] Jacquinot J., and the JET Team, Deuterium-tritium operation in magnetic confine-
ment experiments: results and underlying physics,Plasma Phys. Control. Fusion 41
(3A), A13-A46, 1999.
[13] Aymar R., Chuyanov V. A., Huguet M., Shimomura Y., ITER Joint central Team,
ITER Home Teams, Overview of ITER-FEAT - The future international burning
plasma experiment, Nucl. Fusion 41 (10), 1301-1310, 2001.
[14] Schumacher U., Status and problems of fusion reactor development, Naturwis-
senschaften 88 (3), 102-112, 2001.
139
140 BIBLIOGRAPHY
[15] Hartmann D., Plasma heating, in Summer University for Plasma Physics, September
17-21, 2001, Garching, ed. Könies A., Krieger K. 185-200, 2001.
[16] Nedospasov A. V., Petrov V. G., Fidel’man G. N., Plasma convection in the poloidal
limiter shadow of a tokamak, Nucl. Fusion 25, 21-27, 1985.
[17] Endler M., Turbulent SOL transport in stellarators and tokamaks, J. Nucl. Mater.
266-269, 84-90, 1999.
[18] Nedospasov A. V., Petrov V. G., Unipolar arcs as impurity source in tokamaks, J.
Nucl. Mater., 93 & 94, 775-779, 1980.
[19] Behrisch R., Surface erosion by electrical arcs, Physics of Plasma-Wall Interaction
in controlled Fusion, NATO asi Series, Series B: Physics 131, 1984.
[20] Wolff H., Arcing in magnetic fusion devices, Atomic and plasma-material interaction
data for fusion, Suppl. to Nuclear Fusion 1, 93-107, 1991.
[22] Höhn F., Multipactor discharges and their interaction with surfaces, PhD Thesis,
1997, in German.
[23] Latham R.V., High voltage insulation: The physical basis, Academic Press, 1981.
[24] Noterdaeme J.-M., Wesner F., Brambilla M., Fritsch R., Kutsch H.-J., Soell M., The
ASDEX Upgrade ICRH antenna, Fusion Eng. Des. 24, 65-74, 1994.
[25] Hofmeister F., Braun F., Wesner F., The RF system and matching procedure for
ASDEX and ASDEX Upgrade, Fusion Eng. Des. 24, 83-89, 1994.
[26] Wedler H., Wesner F., Becker W., Fritsch R., Vacuum insulated antenna feeding lines
for ICRH at ASDEX Upgrade, Fusion Eng. Des. 24, 75-81, 1994.
[27] Braun F., Hofmeister F., Wesner F., Becker W., Faugel H., Hartmann D., Noter-
daeme J.-M., ICRF system enhancements at ASDEX Upgrade, Fusion Eng. Des. 56
(7), 551-555, 2001.
[28] Noterdaeme J.-M., Wukitch S., Hartmann D. A., Brambilla M., Braun F. et al.,
ICRF heating results in ASDEX Upgrade and W7-AS, Fusion Energy (Montreal,
1996) 3, IAEA 335-342.
[29] Wesner F., Becker W., Braun F., Faugel H., Fritsch R. et al., ICRF operation during
H-mode with ELMs: Development and Status at ASDEX Upgrade, Fusion Techonl-
ogy (Lisbon, 1996), (C. Varandas and F. Serra ed.), North-Holland Publ., Amsterdam
537-540.
BIBLIOGRAPHY 141
[30] Noterdaeme J.-M., Hartmann D. A., Stäbler A., Brambilla M. et al., Comparing high
power ion cyclotron resonance frequency heating with neutral injection in ASDEX
Upgrade: differences, similarities and synergies, Fusion Energy (Yokohama, 1998),
IAEA.
[31] Arnold W., De Grassie J.S., Remsen D.B., Pivit E., Russel R., Braun F., Ferrite
hybrid system for generator matching in the ICRF, Fusion Techonlogy (Marseille,
1998), (B. Beaumont, P. Libeyre, B. de Gentile, G. Tonon ed.), EURATOM-CEA,
Cadarache, 427.
[32] Wesner F., Becker W., Brambilla M., Braun F., Faugel H. et al., Status and planning
for ICRF at ASDEX Upgrade, Radio Frequency Heating and Current Drive of Fusion
Devices (Brussels, 1998; J. Jacquinot, G. Van Oost and R.R. Weynants ed.) 22A,
Europ. Phys Society 13-16.
[34] Lieberman M.A., A.J. Lichtenberg, Principles of plasma discharges and materials
processing, John Wiley & Sons, Inc., 1994.
[35] Godyak V.A., Khanneh A.S., IEEE Trans. Plasma Sci. 14, 112, 1986.
[36] Godyak V.A., Soviet Radio Frequency Discharge Research, Delphic Associates Inc.
Falls Church VA, 1986.
[37] Godyak V.A., Piejak R. B., Alexandrovich B. M., IEEE Trans. Plasma Sci. PS-19,
660, 1991.
[38] Lisovskii V.A., Features of the α-γ transition in a low-pressure rf argon discharge,
Techn. Phys. 43 (5), 526 - 534, 1998.
[39] Polozhiy K., α–γ transition in a high Q circuit, private communication, September
2002.
[41] Loew G. A. and Wang J.W., RF breakdown studies in room temperature electron
LINAC structures, SLAC-PUB-4647, 1988.
[42] Loew G. A. and Wang J.W., RF breakdown and field emission, SLAC-PUB-4845,
1989.
[43] Tan J., PhD Thesis, Université Pierre et Marie Curie (Paris), 1989.
[44] Luong M., Safa H., Bonin B., et al., A comparison of enchanced field emission from
broad surfaces in DC and RF regimes, J. Phys. D: Appl. Phys. 30 (8), 1248-1251,
April 1997.
142 BIBLIOGRAPHY
[45] Paulini J., Klein T., Simon G., Thermo-field emission and the Nottingham effect, J.
Phys. D: Appl. Phys. 26, 1310-1315, 1993.
[46] Fomenko V. S., Podchernyaeva I. A., Emission and adsorption properties of materials.
Handbook, Atomizdat, Moscow 1975.
[47] Padamsee H., Knobloch J., The nature of field emission from microparticles and the
ensuing voltage breakdown, High Energy Density Microwaves, 1998 (AIP Conference
Proceedings, edited by R. M. Phillips) 474, 212, AIP, NY, 1999.
[48] Werner G. R., Betzwieser J. C., Knobloch J., Padamsee H., Qureshi M., Shipman J.
E., Investigation of voltage breakdown caused by microparticles, Particle Accelerator
Conference, 2001 Chicago, 1071-1073, IEEE, 2001.
[49] Latham R. V., Prebreakdown electron emission, IEEE Trans, Elec. Insul. EI-18 (3),
1983.
[50] Halbritter J., On contamination on electrode surfaces and electric field limitations,
EI-20 (4), 1985.
[51] Cuneo M.E., The effect of electrode contamination, cleaning and conditioning on
high-energy pulsed-power device performances, IEEE Trans, Elec. Insul. 6 (4), 469-
485, 1999.
[52] Noer R. J., Electron field emission from broad-area electrodes, Appl. Phys. A 28,
1-24, 1982.
[53] Mesyats G.A., Proskurovsky D.I., Pulsed electrical discharge in vacuum, Springer
Series on Atoms and Plasmas, 1989.
[55] Charbonnier F. M., Strayer R. W., Swanson L. W., Martin E. E., Nottingham effect
in field and T-F emission: heating and cooling domains, and inversion temperatures,
Phys. Rev. Lett. 13 (13), 397-401, 1964.
[56] Niayesh K., Influence of electrode surface microstructure on the state of short vacuum
gaps after interuption of high frequency currents, J. Phys. D: Appl. Phys. 33, 2189-
2191, 2000.
[58] McCracken G.M., Stott P.E., Plasma-surface interactions in tokamaks, Nucl. Fusion
19 (7), 889-981, 1979.
[60] Bruining H., Physics and applications of secondary electron emission, Pergamon Press
Ltd, London, 1954.
[61] Thomas E. W., Particle induced electron emission, Atomic and plasma-material in-
teraction data for fusion, Suppl. to Nuclear Fusion 1, 79-91, 1991.
[62] Rusic D., Moore D., Manos D., Cohen S., Secondary electron yeilds of carbon-coated
and polished stainless steel, J. Vac. Sci. technol. (20) (4), 1313-1316, 1982.
[63] Höhn F., Jacob W., Beckmann R., Wilhelm R., The transition of a multipactor to a
low-pressure gas discharge, Phys. Plasmas 4 (4), 940-944, 1997.
[64] McDaniel E. W., Collision phenomena in ionized gases, edited by S. C. Brown, John
Wiley & Sons, Inc., New York, London, Sydney, 1964.
[65] Brown S. C., Basic data of plasma physics, sec. edition, the M.I.T. press, 1966.
[66] Märk T. D., Dunn G. H. ed., Electron impact ionization, Springer-Verlag, Wien New
York, 1985.
[67] Henrikse B. B., Keefer D. R., Clarkson M. H., Electromagnetic field in electrodeless
discharge, J. Appl. Phys. 42 (13), 5460-5464, 1971.
[68] Raizer Yu.P., Shneider M. N., Yatsenko N. A., Capacitively coupled RF discharges,
CRC Press, Boca Raton, London, Tokyo, 1995.
[69] Baity F. W., Barber G. C., Bobkov V., Goulding R. H., Noterdaeme J.-M., Swain
D. W., 14th Topical Conference on Radiofrequency Power in Plasmas, Oxnard 2001,
AIP Conference proceedings 595, AIP Press, Melville, NY 2001, 410-513.
[70] Shoji M., Sato M., Breakdown measurement of a capacitive radio frequency discharge
with insulated electrode, Appl. Phys. Lett. 71 (23), 3362-3363, 1997.
[71] Lisovskiy V. A., Yegorenkov V. D., RF breakdown of low-pressure gas and a novel
method for determination of electron-drift velocities in gases, J. Phys. D: Appl. Phys.
31, 3349-3357, 1998.
[72] Raizer Yu.P., Shneider M.N., Near-electrode sheath dynamics, current anharmonicity
and battery effect in asymmetrical, low-pressure radio-frequency discharges, Plasma
Sources Sci. Technol. 1, 102-108, 1992.
[73] Verplancke P., Langmuir probes at high frequencies in a magnetized plasma: theoret-
ical and experimental plasma: a theoretical and experimental study, Ph. D. Thesis,
Max-Planck-Institut für Plasmaphysik, 1997.
[74] Riemann K.-U., The Bohm criterion and sheath formation, J. Phys. D: Appl. Phys.,
24, 493-518, 1991.
144 BIBLIOGRAPHY
[75] Riemann K.-U., The validity of Bohm’s sheath criterion in RF discharges, Phys.
Fluids B, 4 (9), 2693-2695, 1992.
[76] Lieberman M. A., Analytical Solution for Capacitive RF sheath, IEEE Trans. Plasma
Sci., 16 (6), 638-644, 1988.
[77] Chodura R., Neuhauser J., Energy gain of plasma ions in a strong high frequency
electric field between two target plates, 16th EPS Conference on Controlled Fusion
and Plasma Physics (Venice, 1989) 13B III, 1089.
[78] Chen F. F., Introduction to Plasma Physics, Plenum, New York, London, 1984.
[79] Van Nieuwenhove R., Koch R., Van Oost G., Observation of plasma expulsion from
a powered screenless ICRF antenna, 21st EPS Conference on Controlled Fusion and
Plasma Physics (Montpellier, 1994) 18B II, 976.
[80] Bohm D., Characteristics of Electrical Discharges in Magnetic Fields (ed. Cuthrie
A., Wakerling R.K.), Mc Graw-Hill, N.Y., 1949.
[81] D’Ippolito D. A., Myra J. R., Scrape-off layer profile modifications by convective
cells, Phys. Plasmas 3 (2), 699-701, 1996.
[82] D’Ippolito D. A., Myra J. R., Rogers J. H., et al., Analysis of rf sheath interactions
in TFTR, Nuclear Fusion 38 (10), 1543-1563, 1998.
[83] Bécoulet M., Colas L., Pécoul S. et al., Edge plasma density convection during ion
cyclotron resonance heating on Tore Supra, Phys. Plasmas, 9 (6), 2619-2632, 2002.
[84] Riyopoulos S., Chernin D., Dialetis D.,Theory of electron multipactor in crossed
fields, Phys. Plasmas 2 (8), 1995.
[85] Wesner F., Becker W., Multipacting in magnetic field, private communication, May,
2002.
[86] Bobkov V., Azarenkov M., Bizyukov A., Bobkov Vl., Noterdaeme J.-M., Wilhelm R.,
Development of coatings to improve RF voltage stand-off, 8th International Confer-
ence on Plasma Surface Engineering, Garmisch-Partenkirchen, 2002.
[89] Zhurin V.V., Kaufmann H.R., Robinson R.S., Physics of closed drift thrusters,
Plasma Sources Sci. Technol. 8 (1), R1-R20, 1999.
[92] Terman F. E., Electronic and radio engineering, McGraw-Hill Book Company, Inc.,
1955.
[93] Guo W., DeJoseph C. A. Jr, Time-resolved current and voltage measurements on a
pulsed rf inductively coupled plasma, Plasma Sources Sci. Technol. 10, 43-51, 2001.
[94] Endler M., Niedermeyer H., Giannone L., Holzhauer E., Rudyj A., Theimer G., Tsois
N., ASDEX TEAM, Measurements and modelling of electrostatic fluctuations in the
scrape-off layer of ASDEX, Nucl. Fusion 35 (11), 1307-1339, 1995.
[96] Antar G. Y., Devynck P., Garbet X., Luckhardt S. C., Turbulence intermittency and
burst properties in tokamak scrape-off layer, Phys. Plasmas 8 (5), 1612-1624, 2001.
[97] Boedo J.A., Rudakov D., Moyer R., Krasheninnikov S. et al., Transport by inter-
mittent convection in the boundary of the DII-D tokamak, Phys. Plasmas 8 (11),
4826-4833, 2001.
[98] Neuhauser J.et al., Transport into and across the scrape-off layer in the ASDEX
Upgrade divertor tokamak, Plasma Phys. Control. Fusion 44, 855-869, 2002.
[99] Mueller H.W., Bobkov V., Haas G., Jakobi M., Laux M., Maraschek M., Neuhauser
J., Reich M., Rohde V., Schweinzer J., Wolfrum E., and ASDEX Upgrade Team,
Profile and transport studies in the outer scrape-off-layer at ASDEX Upgrade, 29th
conf. on Controlled Fusion and Plasma Physics
[100] Sarazin Y., Ghendrih Ph., Intermittent particle transport in two-dimensional edge
turbulence, Phys. Plasmas 5 (12), 4214-4228, 1998.
[101] Maslov V.I., Coherent structures at the plasma edge, private communication, May
2002.
[102] Krasheninnikov S., On scrape off layer plasma transport, Phys. Lett. A 283 (5-6),
368-370, 2001.
[103] D’Ippolito D. A., Myra J. R., Krasheninnikov S. I., Cross-field blob transport in
tokamak scrape-off layer plasmas, Phys. Plasmas 9 (1), 222-233, 2002.
[104] Zohm H., Edge localized modes (ELMs), Plasma Phys. Control. Fusion 38, 105-128,
1996.
[105] Suttrop W., The physics of large and small edge localized modes, Plasma Phys.
Control. Fusion 42, A1-A13, 1996.
[106] Porter G. D., Casper T. A., Moller J. M., Edge localized mode particle losses from
the DII-D tokamak, Phys. Plasmas 8 (12), 5140-5150, 2001.
146 BIBLIOGRAPHY
[107] Murmann H., Goetsch S., Salzmann H., and Steuer K.-H., The Thomson scattering
systems of ASDEX Upgrade tokamak, Rev. Sci. Instrum. 63, 4941-4943, 1992.
[108] Rohde V., Laux M., Bachmann P., Herrmann A., Weinlich M., and ASDEX Upgrade
team, Direct measurements of the plasma potential in the edge of ASDEX Upgrade
tokamak, J. Nucl. Mater 241-243, 712-715, (1997).
[109] Golant V. E., Fedorov V. I., Radiofrequency methods of plasma heating in toroidal
thermonuclear devices, Energoatomizdat, Moscow, 1986, in Russian.
[110] Stix T. H., Waves in plasmas, American Institute of Physics, New York, 1992.
[111] Juettner B., Cathode spots of electric arcs, J. Phys. D: Appl. Phys. 34, R103-R123,
2001.
[112] Juettner B., and Kleberg I., The retrograde motion of arc cathode spots in vacuum,
J. Phys. D: Appl. Phys. 33, 2025-2036, 2001.
[113] Robson A. E., The motion of a low-pressure arc in a strong magnetic field, J. Phys.
D: Appl. Phys. 11, 1917-1923, 1978.
[114] Laux M., Drift of the cathode spots in a magnetic field, private communication,
February 2003.
[116] Carslaw H.S., Jaeger J. C., Conduction of heat in solids (second edition), Oxford
University press, 1959.
[117] Lotz W., An empirical formula for the electron impact ionization cross section, Z.
Physik 206, 205-211, 1967.
[118] Janev R. K., Langer W. D., Evans K. Jr., Post D. E. Jr., Elementary processes in
Hydrogen-Helium Plasmas, cross sections and reaction rate coefficients, Springer -
Verlag Berlin Heidelberg New York London Paris Tokyo, 1987.
[119] Janev R. K., Smith J. J., Cross sections for collision processes of hydrogen atoms
with electrons, protons and multiply charged ions, Atomic and plasma-material in-
teraction data for fusion, Suppl. to Nuclear Fusion 4, 1993.
Acknowledgements
Here I would like to thank everybody who made it possible to conduct and to complete
this work. Especially I am grateful to Prof. Dr. R. Wilhelm and Dr. J.-M. Noterdaeme
who gave me a chance to carry out this PhD project and supervised me during the work.
I thank them for the scientific and the organizational discussions which helped me to find
solutions in the situations where I could not see any way out at first.
I thank kindly Mrs. W. Sinz, Mrs. H. Martin and Mrs. C. Alberter for helping to
solve a huge number of problems including organizational and personal ones.
I thank Dr. F. Wesner, Dr. D. Hartmann, F. Braun, J. Kneidl, H. Faugel, W. Becker
and F. Hofmeister for the enormous help and support during design of the electrical part
of the experimental device and during tests and operation of the device in the test facility
and on ASDEX Upgrade.
I would like to thank separately G. Heilmaier and P. Leitenstern for working on the
design of the experimental device. I thank G. Heilmaier, P. Leitenstern, G. Siegl and P.
Angene for assembling the device and mounting it in the test facility and on the midplane
manipulator of ASDEX Upgrade. I would also like to thank them for their patience when
we found out that assembling and mounting the device more than once is required.
I thank M. Münich for the help with the capacitors which eventually solved the problem
of the RF screening of the diagnostic cables in the experimental device. I would also like
to thank other members of ECRH and NBI groups of the Technology division of IPP in
Garching for their help in many situations.
I am grateful to Dr. J. Neuhauser for numerous discussions on physics and the ways of
doing the things. I thank Dr. Rohde, Dr. H. Maier, Dr. H.-W. Müller and Dr. M. Laux
for the discussions which helped to conduct the experiments on ASDEX Upgrade. I feel
grateful to Dr. O. Gruber and the experimental leaders of ASDEX Upgrade for giving
me the experimental time and coordinating the work with the diagnostics involved in the
experiments. I thank the ASDEX Upgrade Team for technical support and for providing
the data from different diagnostics as well as giving me answers on different questions
about these diagnostics.
I thank Dr. D. Swain, Dr. R. Goulding and Dr. F. Baity from the Oak Ridge National
Laboratory for cooperation and discussions about physics and the RF technology.
I would like to thank Dr. A. Lyssoivan for the discussion about the aspects of inter-
action of the RF power with the plasma.
I thank Dr. A. Bizyukov and Dr. V. Maslov from Kharkov National University for
many helpful discussions about non-linear plasma physics. A lot of ideas which raised
from these discussions were used in this thesis. I want to thank Prof. Dr. M. Azarenkov
and Dr. I. Girka for promoting me as a promising post-graduate student and for the
administrative support.
Especially I want to thank my father Dr. Valentyn Bobkov for his parental support
and supervision, and also for the successful scientific cooperation. I want also to thank
all my family, in particular my wife, for their support though the geographical distance
between us was large most of the time.