Refine
Year of publication
Document Type
- Doctoral Thesis (23)
Language
- English (23) (remove)
Has Fulltext
- yes (23)
Is part of the Bibliography
- no (23)
Keywords
- Plasma (23) (remove)
Institute
- Institut für Physik (23) (remove)
Recent experimental campaigns in the Wendelstein 7-X stellarator, a
plasma-confining device designed to investigate the Magnetic Confinement Fusion
(MCF) approach to generating electrical power, have shown that the injection of
fuelling pellets had an unexpected and considerable impact on the performance of
the plasma. Rather than simply refuelling the device and `diluting' the plasma
energy, pellet injection is followed by a significant increase in the ratio of
the ion temperature to the electron temperature. It has been suggested that this
is not merely due to the improved confinement following the reduction of
turbulent transport after the pellet material has homogenised with the bulk
plasma, but also due to a direct transfer of energy from electrons to ions. The
proposed mechanism for this energy transfer is the ambipolar expansion of the
pellet plasmoid, the localised plasma structure produced by the
ionisation of ablated pellet material, along magnetic field lines.
Early work on pellet plasmoid expansion predicted that half the heating power
deposited in plasmoid electrons by collisions with hot ambient electrons is
transferred to plasmoid ions in the form of flow velocity as the plasmoid
expands. The complicated nature of the system of the pellet plasmoid embedded in
the ambient plasma, particularly the behaviour of electrons, which experience
many collisional and collisionless phenomena on multiple disparate timescales,
means that early models of the expansion were not wholly self-consistent, but
rather made use of strong approximations that apply in some regions of the
plasmoid but not in others. For example, only electrons and ions associated with
the plasmoid were rigorously treated, meaning that the framework was one of
`expansion into vacuum'. Combined with the assumption of Maxwellian electrons,
this led to an electric potential that was unbounded at infinity. Naturally, the
validity of the conclusions of such a model are called into question because the
approximations lose their validity far from the plasmoid and as time advances,
yet predictions about the final state of the plasma are desired. A deeper
investigation is required: careful consideration of the phenomena in question
and the timescales (and lengthscales) on which they act must be made in order to
rigorously construct a model that is valid throughout the entire expansion.
The first two papers presented in this thesis iterate on the model established
in the paper that first predicted the electron-to-ion energy transfer; their aim
was to find out how the character of the expansion changes with a more
sophisticated and accurate description of various phenomena, while remaining
within the existing framework of expansion into vacuum. Ultimately, we find that
the qualitative character is unchanged, and that approximately half the heating
power deposited in plasmoid electrons is transferred to ions.
Two other papers in this thesis address the limitations of the original model.
This is achieved by properly considering the electron kinetic problem in a
plasmoid. One paper considers the electron kinetic problem when electrons are
highly isotropised. In this case the kinetic equation can be integrated to
remove all but two independent variables, which is the maximum possible
reduction considering it is a time-dependent problem. The full nonlinear
integro-differential Landau self-collision operator is integrated exactly and
few approximations are made, leading to a rather general kinetic equation.
However, for fuelling pellets some anisotropy in the electron distribution is
expected. Another paper considers the electron kinetic problem (and the entire
plasmoid expansion) allowing for electron anisotropy. Careful consideration of
the ordering of timescales of electron phenomena in a pellet plasmoid leads to a
steady-state kinetic problem that we call collisional quasi-equilibrium (QE). QE
appears in many ways similar to the collisional steady-state characterising a
true thermal equilibrium. It was found that the time-dependent kinetic problem
of the earlier paper, with isotropic electrons, produces the QE distribution
function, corroborating the existence of the QE state. We then take moments of
the electron kinetic equation that is valid on the expansion timescale, assuming
that the electron distribution is that given as the solution to the QE kinetic
problem. This is completely analogous to what is done to obtain the Braginskii
equations or any Chapman-Enskog theory. The result is a set of equations for the
long-term evolution of the macroscopic quantities that describe the distribution
function existing in a quasi-steady-state at each point in time. It is from this
point that one may feasibly describe the plasmoid expansion with an accurate
picture of the electron kinetics and finally obtain the electron-to-ion energy
transfer so desired in a rigorous model of the expansion.
From a broader point of view, the two frameworks provided by these rigorous
investigations of the electron kinetic problem serve as a basis for the future
study of plasmoids. Such a `first-principles' approach to plasmoid dynamics is
novel and interesting in its own right, but it will be demonstrated that such an
approach is essential for pellet plasmoids owing to the fact that they are
poorly described by the `standard tools' of plasma physics.
Using the QE framework it was found that, once more, about half the heating
power experienced by plasmoid electrons is transferred to plasmoid ions. The
incredible robustness of the prediction of such an energy transfer is, in the
author's opinion, the result of the self-similar nature of the expansion found
as a solution to the original model. As a rule, the profiles of self-similar
solutions tend to be attractors for the `real', more complicated, system, and
the qualitative predictions involving no parameters, of which the
electron-to-ion energy transfer is one, tend to be very sturdy.
Aside from fuelling pellets, composed of hydrogen or deuterium, one paper in
this thesis investigates the physics of high-Z pellets that are designed to
terminate the plasma safely in the event of a `disruption', where much of the
magnetic field energy is channelled into a runaway electron beam with
potentially disastrous consequences if the beam encounters a plasma-facing
component. The paper draws on the work carried out in the paper concerning the
kinetic problem of isotropised electrons in a plasmoid.
This thesis is `cumulative'; the vast majority of the work carried out is
described within a set of Papers, labelled A-E, placed at the back of the text.
There is a preceding `wrapper text' (given in numbered Sections) tasked with
introducing the reader to the topic, guiding the reader through the papers, and
expounding some of their main results. Some amount of material not present in
the papers is also provided in the wrapper text. Naturally, the wrapper text
mainly focusses on the results of the papers which are under my first
authorship. In the course of publishing papers over an extended period of time
the nomenclature is bound to vary. Although it is mostly consistent between the
papers, a few difference do arise, and the section `Common symbols and
subscripts' is provided in the frontmatter to alleviate confusion. Particular
care should be taken with the symbols x and z; both can refer to the
coordinate parallel to the magnetic field line, but in papers where z is used
for this purpose x tends to have another definition. In the wrapper text the
choice of symbols is generally chosen to reflect those in the corresponding
paper.
Kinetic modeling and infrared spectroscopy of charge carriers across the plasma-wall interface
(2022)
In this thesis, charge transport at the plasma-wall interface is investigated theoretically, on a semiclassical, microscopic level. Based on the Boltzmann and Poisson equations a set of equations is derived and numerically solved to model charge carriers both within a semiconducting wall and a gaseous plasma in front of it. While the plasma is considered collision-free, within the solid, phonon collisions, as well as recombination processes between conduction band electrons and valence band holes are considered. This results, for the first time, in a self-consistent modeling of both the gaseous electron-ion plasma and the electron-hole plasma in the solid on the same footing. Utilizing specific approximations for different physical scenarios, numerical solutions are presented both for the floating and the electronically contacted (biased) interface. In the latter case, the current voltage characteristic is calculated and shown to heavily depend on the charge kinetics within the wall.
Furthermore, we present optical methods to measure the wall charge noninvasively. These utilize the influence of the deposited surplus charges on the optical reflection coefficient of the surface. By calculating the optical response of these charges, we show that the magnitude of the surface charge can be inferred from the change in the reflectivity of the surface caused by the presence of the plasma. While nonlocal effects are considered, it is shown analytically and numerically that these can be neglected at the scales of the considered physical systems.
Modern space missions depend more and more on electric propulsion devices for in-space
flights. The superior efficiency by ionizing the feedgas and propelling them using electric
fields with regard to conventional chemical thrusters makes them a great alternative. To
find optimized thruster designs is of high importance for industrial applications. Building
new prototypes is very expensive and takes a lot of time. A cheaper alternative is to rely
on computer simulations to get a deeper understanding of the underlying physics. In order
to gain a realistic simulation the whole system has to be taken into account including the
channel and the plume region. Because numerical models have to resolve the smallest time
and spatial scales, simulations take up an unfeasible amount of time. Usually a self-similarity
scaling scheme is used to greatly speed up these simulations. Until now the limits of this
method have not been thoroughly discussed. Therefore, this thesis investigates the limits
and the influence of the self-similarity scheme on simulations of ion thrusters. The aim
is to validate the self-similarity scaling and to look for application oriented tools to use
for thruster design optimization. As a test system the High-Efficiency-Multistage-Plasma
thruster (HEMP-T) is considered.
To simulate the HEMP-T a fully kinetic method is necessary. For low-temperature plasmas,
as found in the HEMP-T, the Particle-in-Cell (PIC) method has proven to be the best
choice. Unfortunately, PIC requires high spatial and temporal resolution and is hence
computationally costly. This limits the size of the devices PIC is able to simulate as well
as limiting the exploration of a wider design space of different thrusters. The whole system
is physically described using the Boltzmann and Maxwell equations. Using these system
of equations invariants can be derived. In the past, these invariants were used to derive a
self-similarity scaling law, maintaining the exact solution for the plasma volume, which is
applicable to ion thrusters and other plasmas. With the aid of the self-similarity scaling
scheme the computation cost can be reduced drastically. The drawback of the geometrical
scaling of the system is, that the plasma density and therefore the Debye length does not
scale. This expands the length at which charge separation occurs in respect to the system
size. In this thesis the limits of this scaling are investigated and the influence of the scaling
at higher scaling factors is studied. The specific HEMP-T design chosen for these studies is
the DP1.
Because the application of scaling laws is limited by the increasing influence of charge separation with increased scaling, PIC simulations still are computationally costly. Another approach to explore a wider design space is given using Multi-Objective-Design-Optimization
(MDO). MDO uses different tools to generate optimized thruster designs in a comparatively
short amount of time. This new approach is validated using the PIC method. During this
validation the drawback of the MDO surfaces. The MDO calculations are not self-consistent
and are based on empirical values of old thruster designs as input parameters, which not
necessarily match the new optimized thruster design. By simulating the optimized thruster
design with PIC and recalculate the former input parameters, a more realistic thruster design is achieved. This process can be repeated iteratively. The combination of self-consistent
PIC simulations with the performance of MDO is a great way to generate optimized thruster
designs in a comparatively short amount of time. The proof of concept of such a combination
is the pinnacle of this thesis.
In this work, studies with respect to the exhaust problem were performed
in the stellarator experiment Wendelstein 7-X with different target concepts and different magnetic field geometries. Different infrared cameras were used to study the heat flux from the plasma onto the PFC. In the first publication, the limiter set-up was used with a simpler magnetic topology in the plasma edge. The radial fall-off of the parallel heat flux for inboard limiters in W7-X shows, similar to inboard limiters in tokamaks, two different radial fall-off lengths, a short (narrow) one, characterizing the near-SOL, and a long (broad) characterizing the far-SOL. For the far-SOL, the heating power and connection length have been identified as the main scaling parameters, while for the near-SOL, the electron temperature close to the LCFS has been identified as the main scaling parameter. The two fall-off lengths differ by a factor 10, and the found scalings for both regimes differ from known models and experimental scalings in tokamaks. A turbulent-driven feature was discussed in the publication as a possible explanation for the behavior of the fall-off length in W7-X.
The gained information and data have been further used to support many
other publications, covering the symmetry of the heat loads, the
energy balance of the machine, and seeding experiments.
The heat exhaust in W7-X with an island divertor was studied in the second
and third publication. Definitions of parameters such as peaking factor and
wetted area were applied for the heterogeneous heat flux pattern on the
W7-X divertor. It was shown that the island divertor concept is capable
of spreading out the heat efficiently, resulting in large wetted areas of up to 1.5 m2. The reached values for the wetted area are comparable to the ones of the larger tokamak JET but with a much smaller ratio of wetted
area to the area of the last closed flux surface. Furthermore, a positive
scaling of the wetted area with the power in the SOL was observed. This
scaling is beneficial for future reactors but needs further investigation of the involved transport processes. The peaking factor (discussed in the second publication) describes how concentrated the heat load is within the region of the strike line. It was shown that this factor is decreasing for increasing densities without affecting the wetted area. The present work paves the way for further analysis of the transport processes of the heat flux towards the island divertor of Wendelstein 7-X.
The active screen plasma nitrocarburizing (ASPNC) technology is a state-of-the-art plasma-assisted heat treatment for improving surface hardness and wear resistance of metallic workpieces based on thermochemical diffusion. In comparison to conventional plasma nitrocarburizing, the use of an active screen (AS) improves thermal homogeinity at the workload and reduces soot formation. Further it can serve as a chemical source for the plasma processes, e.g. by use of an AS made of carbon-fibre reinforced carbon. This compilation of studies investigates the plasma-chemical composition of industrial- and laboratory-scale ASPNC plasmas, predominantly using in-situ laser absorption spectroscopy with lead-salt tuneable diode lasers, external-cavity quantum cascade lasers, and a frequency comb. In this way the temperatures and concentrations of the dominant stable molecular species HCN, NH3, CH4, C2H2, and CO, as well as of less prevelant species, were recorded as functions of e.g. the pressure, the applied plasma power, the total feed gas flow and its composition. Additionally, the diagnostics were applied to a chemically similar plasma-assisted process for diamond deposition.
Resulting from this thesis are new insights into the practical application of an AS made of CFC, the plasma-chemistry involving hydrogen, nitrogen, and carbon, and the particular role of CO as an indicator for reactor contamination. The effect of the feed gas composition on the resulting nitrogen- and carbon-expanded austenite layers was proven by combination of in-situ laser absorption spectroscopy with post-treatment surface diagnostics. Furthermore this work marks the first use of frequency comb spectroscopy with sub-nominally resolved Michelson interferometry for investigation of a low-pressure molecular discharge. This way the rotational bands of multiple species were simultaneously measured, resulting in temperature information at a precision hitherto not reached in the field of nitrocarburizing plasmas.
The importance of ion propulsion devices as an option for in-space propulsion of space
crafts and satellites continues to grow. They are more efficient than conventional chemi-
cal thrusters, which rely on burning their propellant, by ionizing the propellant gas in a
discharge channel and emitting the heavy ions at very high velocities. The ion emission
region of a thruster is called the plume and extends several meters axially and radially
downstream from the exit of a thruster. This region is particularly important for the effi-
ciency of a thruster, because it determines energy and angular distribution of the emitted
ions. It also determines the interaction with the carrier space craft by defining the electric
potential shape and the fluxes and energies of the emitted high energy ions, which are the
key parameters for sputter erosion of satellite components such as solar panels. Developing
new ion thrusters is expensive because of the high number of prototypes and testing cycles
required. Numerical modeling can help to reduce the costs in thruster development, but
the vastly differing length and time scales of the system, particularly the large differences of
scales between the discharge chamber and the plume, make a simulation challenging. Often
both regions are considered to be decoupled and are treated with different models to make
their simulation technically feasible. The coupling between channel and plume plasmas and
its influence on each other is disregarded, because there is no interaction between the two
regions. Therefore, this thesis investigates the physical effects which arise from this cou-
pling as well as models suitable for an integrated simulation of the whole coupled problem
of channel and plume plasmas. For this purpose the High Efficiency Multistage Plasma
Thruster (HEMP-T) ion thruster is considered.
For the discharge channel plasma, a fully kinetic model is required and the Particle-in-Cell
(PIC) method is applied. The PIC method requires very high spatial and temporal resolu-
tions which makes it computationally costly. As a result, only the discharge channel and the
near-field plume close to the channel exit can be simulated. In the channel, the results show
that electrons are magnetized and follow the magnetic field lines. The orientation of the
magnetic field there is mostly parallel to the symmetry axis and the channel walls which re-
sults in a high parallel electron transport and leads to a flat electric potential and a reduced
plasma-wall sheath. Only at the magnetic cusps, which are characteristic of HEMP-Ts the
electrons are guided towards the wall, with ions following due to quasineutrality, where a
classical plasma-wall sheath develops. The ion-wall contact is thus limited to the cusp re-
gion. The small radial drop of the potential towards the wall gives rather low energies of
ions impinging at the wall and minimizes erosion in the HEMP-T.
In the near-field plume, which extends from the thruster exit plane to some centimeters
downstream, the ion emission characteristics is defined. The ratio of radial and axial elec-
tric field components in this region determines the ion emission angle which should be
minimized for maximum thruster efficiency. The plasma discharge in the channel produces
high plasma densities and the subsequent drop from plasma to vacuum potential occurs
further downstream for higher densities. This increases the ratio of radial and axial electric
field components because the plasma expands radially outside of the confinement from the
dielectric discharge channel walls. The potential structure in the near-field plume impacts
also the supply of electrons for the channel discharge because the electrons enter the channel
from the plume. An effect which arises from this coupling is the breathing mode oscilla-
tion. It is an oscillation which is observed in all plasma quantities and is located near the
thruster exit. The oscillation frequency measured in the simulation is in good agreement
with a predator-prey estimate which validates this ansatz. However, the electron tempera-
ture, assumed constant in the predator-prey model, correlates inversely with the oscillation,
i.e. it is minimal at the current maximum and vice versa, which contributes to the observed
oscillations. Because of the oscillation of the plasma number density, the potential drop also
oscillates in the exit region and thus the ratio of radial to axial electric field components,
which results in the oscillation of the mean ion emission angle.
Regarding suitable models for a combined simulation of channel and plume plasmas, the
PIC model for channel and near-field plume is explicitly coupled to a hybrid fluid-PIC
model for the plume. The latter treats the electrons as a fluid, hence increasing the effective
spatial and temporal resolutions which can be applied in the plume simulations at the cost
of reduced accuracy of the electron model. Plasma densities decrease by two orders of
magnitude two meters downstream from the channel exit. The explicitly coupled kinetic
and hybrid PIC models are well suited for the computation of a HEMP-T and its plume
expansion, but they disregard the coupling of channel and plume plasmas for which other
methods are necessary. For this purpose a new approach is presented with a proof-of-
principle validation. The limited spatial resolution in the plume can be overcome with the
mesh-coarsening method, which increases the resolution in regions of low plasma density
without numerical artifacts. Sub-cycling for the electrons in the plume can then be used
to increase the temporal resolution in the plume. The combination of both methods, called
the sub-cycling mesh-coarsening (SMC) algorithm in the scope of this work, promises high
savings in computational cost which can make a combined simulation of plume and channel
plasmas feasible.
Manipulating and utilizing plasmas becomes a more and more important task in various research fields of physics and in industrial developments. Especially in nowadays spacerelevant applications there are different ideas to modify plasmas concerning particular tasks.
One major point of interest is the ability to influence plasmas using magnetic fields. To study the underlying physical effects that were achieved by these magnetic fields for both scenarios Particle-in-Cell simulations were done. Two examples are discussed in this thesis.
The first example originates from an experiment performed by the European Space Agency ESA in collaboration with the German Space Agency DLR. To verify the possibility of heat-flux reduction by magnetic fields onto the thermal protection system of a space vehicle a simplified experiment on earth was developed. Most of the heat that is created during re-entry comes from compression of the air ahead of the hypersonic vehicle, as a result of the basic thermodynamic relation between temperature and pressure. The shock front, which builds up in front of the vehicle deflects most of the heat and prohibits the surface of the space vehicle from direct contact with the maximum flux. State of the art spacecrafts use highly developed materials like ceramics to handle the enormous heat. An attractive approach to reduce costs is to use magnetic fields for heat-flux reduction. This would allow the use of cheaper materials and thus reduce costs for the whole space mission. A partially-ionized Argon beam was used to create a certain heat-flux onto a target. The main finding of the experimental campaign was a large mitigation of heat-flux by applying a dipole-like magnetic field. The Particle-in-Cell method was able to reproduce experimental observations like the heat-flux reduction. An additionally implemented optical diagnostics module allowed to confirm the results of the spectroscopy done during the experiment. The underlying effect that is responsible for the heat-flux reduction was identified as a coupling between the modified plasma and the dominating neutral flux component. The plasma, that is guided towards the target, act as a shield in front of the target surface for arriving neutrals. These neutrals are slowed down by charge-exchange collisions. Furthermore the magnetic field induces an increased turbulent transport that is also needed to reach a reduction in heat-ux. The turbulent transport was also obtained by three-dimensional Direct Simulation Monte Carlo simulations. Unfortunately, such source driven turbulence can not be expected in space, so that a heat flux reduction in real space applications is questionable. Nevertheless, other effects like the induced turbulence by the rotating vehicle can compensate the missing source driven effect.
The second scenario in which a magnetic field is used to modify the heat flux of a plasma is the operation of the pulsed cathodic arc thruster. The same Particle-in-Cell code was used to simulate a typical pulse of this newly developed thruster of Neumann Space Pty Ltd. The typical behavior of the thruster could be reproduced numerically. The thrust is mainly produced by fast electrons. These electrons are accelerated by electric fields as a result of a plasma-beam instability. This plasma-beam instability was verified by a phase space diagnostics for the electrons. To demonstrate the influence of the magnetic field a simulation of the cathodic arc thruster without magnetic field and one with magnetic field were compared. It was shown that the use of a magnetic field leads to a ten times larger thrust by directing the heat ux. The resulting narrow plume is an additional Advantage of the particle guiding magnetic field. This narrowness of the plume reduces the danger of interaction with other components of the space vehicle.
Both scenarios demonstrate the different capabilities for electromagnetic fields to manipulate plasmas and especially the corresponding heat-flux with respect to certain tasks. The possibilities range from reducing the heat-flux onto a target to maximizing the thrust by directing the heat-ux. This thesis demonstrates that simulations are a great tool to support experiments and to deliver an improved physics understanding. They help to identify the basic physics principles in the different systems, because they can deliver information not accessible to experiments.
In particular, a better understanding of the influence of electromagnetic fields on the heat-flux distribution in space-relevant applications was obtained. This can be the basis for further simulation-guided optimization, e.g. for the design of more effective cathodic arc thrusters. Here, the goal is to minimize costs for prototypes by replacing the hardware by virtual prototypes in the simulations. This allows to test basic design ideas in advance and get more highly-optimized designs at a fraction of time and costs.
The content of this thesis can be summarized as follows: (i) The deposition processes of SiOx and SiOxCyHz coatings were investigated in a low-pressure, low temperature HMDSO-O2-N2 plasmas. Infrared laser absorption spectroscopy (IRLAS) and optical emission spectroscopy (OES) were combined to measure the gas temperatures in the hot and colder zones of the plasma as well as to monitor the concentration of the methyl radical, CH3, and of seven stable molecules, HMDSO, CH4, C2H2, C2H4, C2H6, CO and CO2. Tunable lead salt diode lasers (TDLs) and an external-cavity quantum cascade laser (EC-QCL) were simultaneously employed as radiation sources to perform the IRLAS measurements. They were found to be in the range between 10^{11} to 10^{15} cm^{−3}. The influence of the discharge parameters of power, pressure and gas mixture on the molecular concentrations was studied. The plasma generation is characterized by a certain degree of inhomogeneity with different temperature zones, i.e., hottest, hot and colder zones depending on the construction of the reactor. This complexity is characterized by the multiple molecular species including the HMDSO precursor and products in ground and excited states existing in the plasma. (ii) Employing similarly IRLAS and OES techniques, the deposition of nanocrystalline diamond at relatively low temperature in low-pressure MW H2 plasmas with small ad-mixtures of methane and carbon dioxide was investigated. Five methods were applied for an extensive temperature analysis, providing new insights into energetic aspects of the multi-component non-equilibrium plasma. The OES method provided information about the gas temperature of H2 inside the MW plasma. Using lead salt diode lasers, the rotational temperature of the methyl radical, CH3 , and gas temperature of methane molecule, CH4 , was measured. A variety of CO lines in the ground and in three excited states have been analysed using an EC-QCL with a relatively wide spectral range. These methods have shown that based on the construction of the DAA reactor using 16 single plasma sources the plasma generation is characterized by a variety of hottest, hot and colder zones. Extensive measurement of these various species temperatures in the complex plasma enabled the concentration determination of the various stable and unstable plasma species, which were found to be in the range between 10 11 to 10 15 cm −3 . The influence of the discharge parameters, power and pressure, on the molecular concentrations has been studied. To achieve insight into general plasma chemical aspects, the dissociation of the carbon precursor gases including their fragmentation and conversion to the reaction products was analysed in detail. The evolution of the concentration of the methyl radical, CH 3 , of five stable molecules, CH4, CO2, CO, C2H2 and C2H4, and of vibrationally excited CO in the first and second hot band was monitored in the plasma processes by in situ infrared laser absorption spectroscopy using lead salt diode lasers (TDL) and an external-cavity quantum cascade laser (EC-QCL) as radiation sources. OES was applied simultaneously to obtain complementary information about the degree of dissociation of the H2 precursor gas. The analysis of the carbon and oxygen mass balances shows clearly, that the deposition on the reactor walls and the production of other hydrocarbons species may act as sinks for carbon and oxygen. (iii) The absolute line strengths of many P-branch transitions of the ν3 fundamental of {28}^SiH4 were determined using the wide tuning range and the narrow line width of a cw EC-QCL between 2096 and 2178 cm^{−1}. The line positions and line strengths of transitions of the stretching dyad within the P-branch of {28}^SiH4 were determined with an estimated experimental measurement accuracy of 10%. The high spectral resolution available has enabled us to resolve and measure representative examples of the tetrahedral splittings associated with each component of the P-branch. The positions of these components are in excellent agreement with spherical top data system (STDS) predictions and theoretical transitions from the TDS spectroscopic database for spherical top molecules. To our knowledge, this is the first reported measurement of these line strengths in this band and is an example of the applicability of high-powered, widely tunable EC-QCLs to high resolution spectroscopy in the MIR. (iv) Similarly, the determination of the silyl radicals, ν3 band, line strengths is ongoing using the same cw EC-QCL. This effort was impaired by silane and other unknown species lines overlap; however, the silyl radicals was successfully detected in a SiH4/H2 plasma. A method to determine the silyl line strengths has been presented through its iterative decay measurements which relied on the value of the silyl radical self reaction constant. There was a consensus of its value in the literature.
The aim of this thesis is to concentrate on the investigation of these ROS&RNS composition distribution and their production pathways in the gas phase produced by a plasma jet. By understanding the physical mechanisms behind the generation of the ROS&RNS a precise tuning and design of the composition distribution in the gas phase can be achieved. One crucial physical parameter is the dissipated power inside the plasma. Only if this parameter is known a meaningful comparison of different feed gas settings is possible. Therefore, a concept for measuring the dissipated power inside the plasma for the modified micro-scaled atmospheric pressure plasma jet( µAPPJ) is designed. Additionally, due to achievements within this thesis it is now possible to ignite a homogeneous discharge in argon and helium within the geometry of the µAPPJ. The used feed gas is a determining factor concerning the electron energy distribution function and consequently influencing the production mechanism of the ROS&RNS. First of all, the electrical characterisation of the modified µAPPJ was performed including the alpha-to-gamma transition. It is shown that the alpha-to-gamma transition power is increasing with increasing frequency. For the first time it is now feasible to investigate the influence of the dissipated power on the neutral gas temperature, the metastable atom densities and the ROS&RNS production for the modified µAPPJ with argon and helium as feed gas. Due to the possibility of changing the feed gas and controlling the dissipated power a fundamental insight into the production mechanism of the ROS&RNS generated by the plasma jet is achieved. With rising dissipated power the temperature and the metastable densities as well as the ozone and nitrogen dioxide concentrations are increasing. By adding molecular oxygen and nitrogen to the feed gas of a plasma jet the ROS&RNS composition can be tuned. However, also the dissipated power is changed by the small amount of admixtures. Due to the developed dissipated power measurements within this thesis it was possible to disentangle the influence of the admixture on the power and on the ROS&RNS production. If the dissipated power is fixed for the µAPPJ with argon and helium feed gas, respectively, the highest amount of ozone was measured with oxygen admixture in an argon discharge, the highest amount of dinitrogen pentoxide with nitrogen admixture in an argon discharge and the highest amount of nitrogen dioxide with nitrogen admixture in a helium discharge. Beyond the influence of the dissipated power and the molecular admixture on the ROS&RNS production the feed gas temperature is a crucial parameter for the corresponding chemical reactions. By changing this parameter the distribution of ozone and nitrogen dioxide can be tuned precisely in such a way that with increasing temperature the ozone density goes down and the nitrogen dioxide density rises. Another determinant for the ROS&RNS composition produced by an atmospheric pressure plasma jet is the influence of ambient air. If the ambient air is changing from pure nitrogen to pure oxygen atmosphere the ozone density produced by the plasma jet is increasing. For the same conditions the nitrogen dioxide has a maximum at an oxygen-to-nitrogen ratio of 1:1. To avoid the influence of the ambient air on the reactive species production the afterglow of the µAPPJ was prolonged with a glass tube. By increasing the amount of molecular admixtures to the feed gas with each in equal quantities a totally different ROS&RNS composition can be obtained compared without the glass tube. It figures out that for small molecular admixtures the reactive species composition is nitrogen dominated and for higher admixtures it is oxygen dominated. Consequently, by shielding the ambient air from the active effluent and by admixing molecular oxygen and nitrogen the ROS&RNS composition can be designed.
The present thesis deals with dynamic structures that form during the expansion of plasma into an environment of much lower plasma density. The electron expansion, driven by their pressure, occurs on a much faster time scale than the ion expansion, owed to their mobility. The high inertia of the ions causes the generation of an ambipolar electric field that decelerates the escaping electrons while accelerating the ions. The ambipolar boundary propagates outwards and forms a plasma density front. For a small density differences, the propagation of the front can be described with the linear ansatz for ion acoustic waves. For a large density differences, experiments have shown that the propagation velocity of such a front is still related to the ion sound velocity. However, the reported proportionality factors are scattered over a wide range of values, depending on the considered initial and boundary conditions. In this thesis, the dynamics during plasma expansion are studied with the use of experiments and a versatile particle-in-cell simulation. The experimental investigations are performed in the linear helicon device Piglet. The experiment features a fast valve, which is used to shape the neutral gas density profile. During the pulsed rf-discharges, plasma is generated in the source region and expands collisionless into the expansion chamber. The computer simulation is tailored very close to the experiment and provides a deeper insight in the particle kinetics. The experimental results show the existence of a propagating ion front. Its velocity is typically supersonic and depends on the density ratio of the two plasmas. The ion front features a strong electric field. The front can have similar properties to a double layer is not necessarily a double layer by definition. The computer simulation reveals that the propagating electric field repels the downstream ambient ions. These ions form a stream with velocities up to twice as high as the front velocity. The observed ion density peak is due to the accumulation of the repelled ions and is located at their turning point. The ion front formation depends strongly on the initial ion density profile and is part of a wave-breaking phenomenon. The observed front is followed by a plateau of little plasma density variation. This could be confirmed for the expansion experiment by a comparison with virtual diagnostics in the computer simulation. The plateau has a plasma density determined by the ratio between the high and low plasma density. It consists of streaming ions that have been accelerated in the edge of the main plasma. The presented results confirm and extend findings obtained by independent numerical models and simulations.
Magnetic reconnection is a fundamental plasma process where a change in field line connectivity occurs in a current sheet at the boundary between regions of opposing magnetic fields. In this process, energy stored in the magnetic field is converted into kinetic and thermal energy, which provides a source of plasma heating and energetic particles. Magnetic reconnection plays a key role in many space and laboratory plasma phenomena, e.g. solar flares, Earth’s magnetopause dynamics and instabilities in tokamaks. A new linear device (VINETAII) has been designed for the study of the fundamental physical processes involved in magnetic reconnection. The plasma parameters are such that magnetic reconnection occurs in a collision-dominated regime. A plasma gun creates a localized current sheet, and magnetic reconnection is driven by modulating the plasma current and the magnetic field structure. The plasma current is shown to flow in response to a combination of an externally induced electric field and electrostatic fields in the plasma, and is highly affected by axial sheath boundary conditions. Further, the current is changed by an additional axial magnetic field (guide field), and the current sheet geometry was demonstrated to be set by a combination of magnetic mapping and cross-field plasma diffusion. With increasing distance from the plasma gun, magnetic mapping results in an increase of the current sheet length and a decrease of the width. The control parameter is the ratio of the guide field to the reconnection magnetic field strength. Cross-field plasma diffusion leads to a radial expansion of the current sheet at low guide fields. Plasma currents are also observed in the azimuthal plane and were found to originate from a combination of the field-aligned current component and the diamagnetic current generated by steep in-plane pressure gradients in combination with the guide field. The reconnection rate, defined via the inductive electric field, is shown to be directly linked to the time-derivative of the plasma current. The reconnection rate decreases with increasing ratio of the guide field to the reconnection magnetic field strength, which is attributed to the plasma current dependency on axial boundary conditions and the plasma gun discharge. The above outlined results offer insights into the complex interaction between magnetic fields, electric fields, and the localized current flows during reconnection.
The laser-matter interaction is a topic of current research. In this context, the interaction of intensive laser radiation with atomic clusters is of special interest. Du to the small cluster size, the laser field can penetrate the whole cluster volume, which leads to a high absorption of energy in the cluster. As a result, plasmas with high density and high temperature are produced. In the early phase of the laser-cluster interaction, free electrons are initially created in the cluster due to tunnel ionization or photoionization. Via collisions of these electrons with the cluster atoms, the ionization is increased and thus a dense nanoplasma is produced, which is heated by the laser. If free electrons leave the cluster during the laser-cluster interaction (outer ionization), a positive charge buildup is created. The associated charge repulsion finally can lead to the fragmentation of the cluster due to Coulomb explosion. Experimentally, interesting phenomena emerging from laser-excited clusters are observed, e.g., the creation of fast electrons, the production of highly charged ions, and X-ray emission. In this dissertation, the interaction of Gaussian laser pulses in the infrared regime with argon and xenon clusters is simulated by means of a nanoplasma model. Considering laser intensities in the non-relativistic regime, the relevant processes such as ionization, heating and expansion are theoretically described in this model with a set of coupled rate equations and hydrodynamic equations. One focus of the thesis is on the heating of the nanoplasma via inverse bremsstrahlung (IB), which is due to the absorption of laser photons in electron-ion collisions. In particular, the important question is investigated whether the consideration of the ionic structure – that means, the nuclear charge and the bound electrons – modifies the electron-ion collisions and thus the IB heating rate. Starting from a quantum statistical description, effective electron-ion potentials are used which account for both the screening due to the dense plasma and the inner ionic structure. Within the quantum mechanical first Born approximation, the consideration of the ionic structure leads to a drastic increase of the IB heating rate, in particular for high nuclear charges and low ionic charge states. However, for the parameters relevant in experiments, the applicability of the first Born approximation is questionable. Therefore, quantum mechanical calculations going beyond the first-order perturbation theory are performed. In addition, the IB heating rate is investigated with different classical methods. These are based either on transport cross sections for elastic electron-ion scattering or on classical simulations of inelastic scattering processes. Also within the classical approaches, the consideration of the ionic structure leads to an increase of the heating rate. However, this increase is shown to be only moderate. In a further part, the thesis focuses on the question how the dynamics of the laser-cluster interaction is influenced by the consideration of excited states. This is explored exemplarily for argon clusters excited by single or double laser pulses. The consideration of excitation processes in the nanoplasma leads to a decrease of the electron temperature and to an increase of the density of free electrons. Moreover, it is shown that the consideration of excitation processes results in an essential acceleration of the ionization dynamics. As a consequence, the mean ionic charge state in the plasma as well as the number of highly charged ions is significantly increased. For the population of ground states and excited states within an ionic charge state Z, collisional deexcitation processes play an important role. By means of an analytical relation between excitation and deexcitation cross sections, the rates for the respective processes in the presence of the laser field are calculated. The role of deexcitation processes is studied in detail, showing that the inclusion of these processes is essential for the correct theoretical description of the photon emission from laser-excited clusters. Based on these results, the photon yield is calculated for selected radiative transitions resulting from highly charged argon ions in the UV and X-ray regime.
During the past decade, various physical properties of the Yukawa ball, like structure and energy states, were unraveled using experiments. However, the dynamical features served further attention. Therefore, the main aim of my thesis was to investigate and understand how a finite system-represented by Yukawa clusters-evolves from a solid, crystalline structure to a liquid-like system, how it behaves in this phase and in what manner the reordering back into the solid state can be described. As a method of choice to reach this goal, laser heating has been proven successful. Moreover, the special importance of wakefields for dust clusters confined at low neutral gas pressure was addressed. Melting of finite dust clouds can be induced in two ways, either by altering the properties of the ambient plasma or by laser heating. The latter was shown to be a generic melting scenario, allowing to estimate a critical coupling parameter at the melting point. Moreover, the melting transition of finite 3D dust systems was found to be a two-step process where angular order is lost before the radial order starts to diminish at higher energies. Next, the mode dynamics of finite 3D dust ensembles in the solid and the liquid phase was studied. Crystal and fluid modes revealed the main spectral properties of the system. The normal modes are mainly suited to describe crystalline states. Fluid modes were excited naturally and via laser heating, with excitation frequencies almost independent of the coupling parameter in the solid and the liquid-like regime. Tuning the plasma parameters can be used to vary the particle-particle interaction via the ion focus. Both methods, even though assuming equilibrium situations, allowed to hint at these wakefields. The corresponding peaks in the fluid and normal mode spectra were no eigenmodes, confirming the nonequilibrium character of the ion focusing effect. First steps to extend the normal mode theory to achieve the dynamics of wake-affected nonequilibrium dust clusters were presented. Statistical quantities were obtained evaluating long-run experiments and transport coeffcients for finite dust systems were calculated via the instantaneous normal mode technique. Diffusion was found considerably higher for 3D than for 2D dust clusters. Using the configurational entropy, we have shown that in 2D and 3D disorder increases with increasing size of the system, in agreement with simulations. The temperature dependence of the configurational entropy differs for 2D and 3D dust clouds, with a threshold behavior found for finite 2D ensembles only. Finally, using instantaneous normal modes to reveal the total fraction of unstable modes, the predictive connection of Keyes (Phys Rev E 62, p7905, 2000), between transport and disorder was tested and verified for 2D, but not for 3D clusters. The reason for this has to be left open. Finally, laser-mediated recrystallization processes of finite 3D dust clouds were investigated. First, the temporal evolution of the Coulomb coupling parameter was traced during heating and recrystallization. A cooling rate has been determined from the initial phase of recrystallization. This cooling rate is lower than damping by the neutral gas, in agreement with simulations. We have observed a large fraction of metastable states for the final cluster configurations. Further, we have revealed that the time scale for the correlation buildup in the finite 3D ensemble was on even slower scales than cooling. Thus, different time scales can be attributed to the fast emergence of the shells and to the slower individual ordering within the shells.
The confinement of energy has always been a challenge in magnetic confinement fusion devices. Due to their toroidal shape there exist regions of high and low magnetic field, so that the particles are divided into two classes - trapped ones that are periodically reflected in regions of high magnetic field with a characteristic frequency, and passing particles, whose parallel velocity is high enough that they largely follow a magnetic field line around the torus without being reflected. The radial drift that a particle experiences due to the field inhomogeneity depends strongly on its position, and the net drift therefore depends on the path taken by the particle. While the radial drift is close to zero for passing particles, trapped particles experience a finite radial net drift and are therefore lost in classical stellarators. These losses are described by the so-called neoclassical transport theory. Recent optimised stellarator geometries, however, in which the trapped particles precess around the torus poloidally and do not experience any net drift, promise to reduce the neoclassical transport down to the level of tokamaks. In these optimised stellarators, the neoclassical transport becomes small enough so that turbulent transport may limit the confinement instead. The turbulence is driven by small-scale-instabilities, which tap the free energy of density or temperature gradients in the plasma. Some of these instabilities are driven by the trapped particles and therefore depend strongly on the magnetic geometry, so the question arises how the optimisation affects the stability. In this thesis, collisionless electrostatic microinstabilities are studied both analytically and numerically. Magnetic configurations where the action integral of trapped-particle bounce motion, J, only depends on the radial position in the plasma and where its maximum is in the plasma centre, so-called maximum-J configurations, are of special interest. This condition can be achieved approximately in quasi-isodynamic stellarators, for example Wendelstein 7-X. In such configurations the precessional drift of the trapped particles is in the opposite direction from the direction of propagation of drift waves. Instabilities that are driven by the trapped particles usually rely on a resonance between these two frequencies. Here it is shown analytically by analysing the electrostatic energy transfer between the particles and the instability that, thanks to the absence of the resonance, a particle species draws energy from the mode if the frequency of the mode is well below the charateristic bounce frequency. Due to the low electron mass and the fast bounce motion, electrons are almost always found to be stabilising. Most of the trapped-particle instabilities are therefore predicted to be absent in maximum- J configurations in large parts of parameter space. Analytical theory thus predicts enhanced linear stability of trapped-particle modes in quasi-isodynamic stellarators compared with tokamaks. Moreover, since the electrons are expected to be stabilising, or at least less destabilising, for all instabilities whose frequency lies below the trapped-electron bounce frequency, other modes might benefit from the enhanced stability as well. In reality, however, stellarators are never perfectly quasi-isodynamic, and the question thus arises whether they still benefit from enhanced stability. Here the stability properties of Wendelstein 7-X and a more quasi-isodynamic configuration, QIPC, are investigated numerically and compared with another, non-quasiisodynamic stellarator, the National Compact Stellarator Experiment (NCSX) and a typical tokamak. In gyrokinetic simulations, performed with the gyrokinetic code GENE in the electrostatic and collisionless approximation, several microinstabilities, driven by the density as well as both ion and electron temperature gradients, are studied. Wendelstein 7-X and QIPC exhibit significantly reduced growth rates for all simulations that include kinetic electrons, and the latter are indeed found to be stabilising when the electrostatic energy transfer is analysed. In contrast, if only the ions are treated kinetically but the electrons are taken to be in thermodynamic equilibrium, no such stabilising effect is observed. These results suggest that imperfectly optimised stellarators can retain most of the stabilising properties predicted for perfect maximum-J configurations. Quasi-isodynamic stellarators, in addition to having reduced neoclassical transport, might therefore also show reduced turbulent transport, at least in certain regions of parameter space.
The concept of the electron surface layer introduced in this thesis provides a framework for the description of the microphysics of the surplus electrons immediately at the wall and thereby complements the modelling of the plasma sheath. In this work we have considered from a surface physics perspective the distribution and build-up of an electron adsorbate on the wall as well as the effect of the negative charge on the scattering of light by a spherical particle immersed in a plasma. In our electron surface layer model we treat the wall-bound electrons as a wall-thermalised electron distribution minimising the grand canonical potential and satisfying Poissons equation. The boundary between the electron surface layer and the plasma sheath is determined by a force balance between the attractive image potential and the repulsive sheath potential and lies in front of the crystallographic interface. Depending on the electron affinity x, that is the offset of the conduction band minimum to the potential in front of the surface, two scenarios for the wall-bound electrons are realised. For x<0 electrons do not penetrate into the solid but are trapped in the image states in front of the surface where they form a quasi two-dimensional electron gas. For x>0 electrons penetrate into the conduction band where they form an extended space charge. These different scenarios are also reflected in the electron kinetics at the wall which control the sticking coefficient and the desorption time. If x<0 electrons from the plasma cannot penetrate into the solid. They are trapped in the image states in front of the surface. The transitions between unbound and bound states are due to surface vibrations. Trapping of electrons is mediated by one-phonon transitions and takes place in the upper bound states. Owing to the large binding energy of the lowest bound state transitions from the upper bound states to the lowest bound state are due to multi-phonon processes. For low surface temperatures relaxation to the lowest bound state takes place while for higher temperature a relaxation bottleneck emerges. Desorption occurs in cascades for systems without relaxation bottleneck and as a one-way process in systems with a relaxation bottleneck. From the perspective of plasma physics the most important result is that the sticking coefficient for electrons is relatively small, typically on the order of 0.001. For x>0 electron physisorption takes place in the conduction band. For this case sticking coefficients and desorption times have not been calculated yet but in view of the more efficient scattering with bulk phonons, responsible for electron energy relaxation in this case, we expect them to be larger than for the case of x<0. Finally, we have studied the effects of surplus electrons on the scattering of light by a spherical particle. For x<0 the electrons form a spherical electron gas around the particle and their electrical conductivity modifies the boundary condition for the magnetic field. For x>0 the electrons in the bulk of the particle modify the refractive index through their bulk electrical conductivity. In both cases the conductivity is limited by scattering with surface or bulk phonons. Surplus electrons lead to an increase of absorption at low frequencies and, most notably, to a blue-shift of an extinction resonance in the infrared. This shift is proportional to the charge and is strongest for submicron-sized particles. The particle charge is also revealed in a blue-shift of the rapid variation of one of the two polarisation angles of the reflected light. From our work we conclude that the electron affinity is an important parameter of the surface which should affect the charge distribution as well as the charge-up. Therefore, we encourage experimentalists to study the charging of surfaces or dust particles as a function of x. Interesting in this respect is also if or under what conditions the electron affinity of a surface exposed to a plasma remains stable. Moreover, we suggest to use the charge signatures in Mie scattering to measure the particle charge optically. This would allow a charge measurement independent of the plasma parameters and could be applied to nano-dust where conventional methods cannot be applied.
In this work, various aspects of fundamental physics and chemistry of molecular gas discharges are presented with emphasis on the interaction between species, activated by low-pressure plasmas, and surfaces. As already known, synergistic effects of multiple plasma-generated species are responsible for surface modification. However, due to the large number of internal parameters of a discharge and the complex plasma processes the identification of correlations between plasma characteristics and their effects on surfaces are complicated. Therefore, the aim of this thesis is to improve the understanding of several phenomena associated with plasma–surface interactions by measuring or calculating fundamental kinetic, transport or spectroscopic data needed to interpret measurements and hereby, to support some future applications of plasmas.
Beams of ions and electrons are a source of free energy which can be transferred to waves via an instability. Beams exist in almost all plasma environments, but their instabilities are particularly important for the dynamics of space plasmas. In the absence of collisions, the instability drives waves to large amplitudes and forms nonlinear structures such as solitary waves. The electric fields in these waves can scatter particles in the background plasma, or disrupt currents. Both of these effects are important for the overall dynamics of the plasma. In this thesis, both electron and ion beam plasma instabilities have been investigated in the linear plasma device VINETA and using a Particle-in-Cell simulation. The electron beam instability has been demonstrated by previous authors to be a useful diagnostic for the plasma density. The spatial resolution of previous results was confirmed at a few millimetres, and a temporal resolution of 1ms was shown for the first time. An ion beam was generated with a double plasma discharge. Compared to space, this environment and indeed most laboratory plasmas have considerably higher collisionality and a limited spatial extent which introduces gradients in the plasma. Gradients perpendicular to the beam propagation direction are linked to a decrease of both the wavelength and amplitude of the instability. It was observed in both experiment and simulation that gradients in sheaths at the boundaries of the plasma not only affect the time averaged plasma parameters, but also excite instabilities. Fluctuations within the sheath spread the beam in velocity space, effectively increasing its temperature. Warmer beams require a higher drift velocity to excite an instability. This was also confirmed by experimental and numerical results. Collisions are shown to be the dominant damping force for the electron beam instability. For ions, collisions play an important role in the simulation, but appear to be overshadowed by Landau damping from impurities in the experiment. When boundary conditions are removed from the simulation, wave amplitudes increase and nonlinear effects become important. Saturation by particle trapping and coalescence of phase space holes is observed, which could eventually lead to the solitary waves as they are observed in space plasmas.
Electromagnetic Drift Waves
(2010)
In the rf-plasma of the linear magnetized VINETA experiment, different types of low-frequency waves are observed. The emphasis in this work is on the interaction mechanism between drift waves on the one and kinetic Alfven waves on the other hand. In the peaked density profile of the plasma column drift waves occur as modulation of the plasma density. As gradient driven instability, they draw their energy from the radial density gradients. Alfven waves as magnetic field fluctuations are stable in the present configuration. They are launched by a magnetic excitation antenna. Parallel conduction currents in the plasma are common to both wave phenoma. A B-dot probe as standard diagnostic tool is used to detect the fluctuating magnetic fields of both wave types. The challenge are the small induced voltages due to the low wave frequency. The probe design with an integrated amplifier close to the probe head takes this into acount. The developed B-dot probe is mounted to different positioning systems to characterize both wave phenomena. For Alfven waves, the dispersion relation is recorded experimentally. It is found to be in good agreement with the prediction of the Hall-MHD theory with included resistive term, accounting for the cold collisional plasma. The fluctuating magnetic field pattern is recorded with azimuthal scans. The current density is obained by Amperes law. It is concentrated in helically twisted current filaments. For the unstable drift waves, similar investigations are done with simultaneously recorded density fluctuations. In the azimuthal plane, the locations of the parallel current filaments and the fluctuating density are found to be in phase, supporting the predicted drive of parallel currents by pressure gradients. A mutual influence of the two wave types is observed in an interaction experiment. Assuming parallel currents as coupling quantity, an interpretation of the experimental findings is given based on the linear theory of drift waves.
This thesis constitutes a computational study of charge and ion drag force on micron-sized dust particles immersed in rf discharges. Knowledge of dust parameters like dust charge, floating potential, shielding and ion drag force is very crucial for explaining complex laboratory dusty plasma phenomena, such as void formation in microgravity experiments and wakefield formation in the sheaths. Existing theoretical models assume standard distribution functions for plasma species and are applicable over a limited range of flow velocities and collisionality. Kinetic simulations are suitable tools for studying dust charging and drag force computation. The main aim of this thesis is to perform three dimensional simulations using a Particle-Particle-Particle-Mesh ($P^3M$) model to understand how the dust parameters vary for different positions of dust in rf discharges and how these parameters on a dust evolve in the presence of neighboring dust particles. At first, rf discharges in argon have been modelled using a three-dimensional PIC-MCC code for the discharge conditions relevant to the dusty plasma experiments. All necessary elastic and inelastic collisions have been considered. The plasma background is found collisional, charge-exchange collisions between ions and neutrals being dominant. Electron and ion distributions are non-Maxwellian. The dominant heating mechanism is Ohmic. Then, simulations have been done to compute the dust parameters for various sizes of dust located at different positions in the rf discharges. Dust charge and floating potential in the presheath are slightly larger than the values in the bulk due to the higher electron flux to the dust particle in the presheath. From presheath to the sheath the charge and floating potential values decrease due to the decrease of the electron current to the dust. A linear dependence of dust potential on dust size has been found, which results in a nonlinear dependence of the dust charge with the dust size when the particle is assumed to be a spherical capacitor. This has been verified by independently counting the charges collected by the dust. %where indeed it has been noted that the dust charge %scales nonlinearly with the dust size. The computed dust parameters are also compared with theoretical models. Simulated dust floating potentials are comparable to values obtained from Allen-Boyd-Reynolds (ABR) and Khrapak models, but much smaller than the values obtained from Orbit Motion Limited (OML) model. The dust potential distribution behaves Debye-H\"{u}ckel-like. The shielding lengths are in between ion and electron Debye lengths. % indicating shielding by both ions and electrons. Further, the orbital drag force is typically larger than the collection drag force. The total drag force for the collisional case is larger than for the collisionless case and it scales nonlinearly with the dust size. The collection drag values and size-scaling agrees with Zobnin's model. The charging and drag force computation is then extended to two and multiple static dust particles in the rf discharge to study the influence of neighboring dust particles on the dust parameters. Initially, the dust parameters on two dust particles are computed for various interparticle separation distances and for dust particles placed at different locations in the rf discharge. It is observed that for dust separations larger than the shielding length the dust parameters for the two dust particles match with the single dust particle values. As the dust separation is equal to or less than the shielding length the ion drag force increases due to the buildup of a parallel drag force component. However, the main dust properties like charge, potential, vertical component of ion drag are not affected considerably. This is attributed to the smaller collection impact parameter values compared to the dust separation. %This is because the %collection impact parameter values in the sheath and the presheath are smaller %than the smallest dust separation and in case of the dust in the bulk, the %collection impact parameter is comparable with the dust separation. Then the dust charges on multiple dust particles located at different positions in the discharge and arranged along the discharge axis are also computed. It is found that the charges of the multiple dust particles in the bulk or presheath do not differ much from the single particle values at that location. But the dust charges of multiple dust particles located in the sheath drastically differ from the single dust parameter values. Due to ion focusing from dust particles in the upper layers, the ion current increases to dust particles in the lower layers resulting in smaller charge values. This is as well the case where dust particles are vertically aligned as in the standard experiments of dusty plasmas. In conclusion, this work used a fully kinetic (PIC and MD or $P^3M$) model to study the physics of dust charging in rf plasmas. Our simulations revealed that the dust parameters vary considerably from the bulk to the sheath. The CX collisions increase flux to the dust thereby affecting the dust parameters and their scaling with dust size. Also, a dust particle affects the charging dynamics of its neighbor only when their separation is within the shielding length. In the plasma sheath, ion focussing can cause great reduction in dust charges.
The present work is the first work dealing with turbulence in the WEGA stellarator. The main object of this work is to provide a detailed characterisation of electrostatic turbulence in WEGA and to identify the underlying instability mechanism driving turbulence. The spatio-temporal structure of turbulence is studied using multiple Langmuir probes providing a sufficiently high spatial and temporal resolution. Turbulence in WEGA is dominated by drift wave dynamics. Evidence for this finding is given by several individual indicators which are typical features of drift waves. The phase shift between density and potential fluctuations is close to zero, fluctuations are mainly driven by the density gradient, and the phase velocity of turbulent structures points in the direction of the electron diamagnetic drift. The structure of turbulence is studied mainly in the plasma edge region inside the last closed flux surface. WEGA can be operated in two regimes differing in the magnetic field strength by almost one order of magnitude (57mT and 500mT, respectively). The two regimes turned out to show a strong difference in the turbulence dynamics. At 57mT large structures with a poloidal extent comparable to the machine dimensions are observed, whereas at 500mT turbulent structures are much smaller. The poloidal structure size scales nearly linearly with the inverse magnetic field strength. This scaling may be argued to be related to the drift wave dispersion scale. However, the structure size remains unchanged when the ion mass is changed by using different discharge gases. Inside the last closed flux surface the poloidal ExB drift in WEGA is negligible. The observed phase velocity is in good agreement with the electron diamagnetic drift velocity. The energy in the wavenumber-frequency spectrum is distributed in the vicinity of the drift wave dispersion relation. The three-dimensional structure is studied in detail using probes which are toroidally separated but aligned along connecting magnetic field lines. As expected for drift waves a small but finite parallel wavenumber is found. The ratio between the average parallel and perpendicular wavenumber is in the order of 10^-2. The parallel phase velocity of turbulent structures is in-between the ion sound velocity and the Alfvènvelocity. In the parallel dynamics a fundamental difference between the two operational regimes at different magnetic field strength is found. At 500mT turbulent structures can be described as an interaction of wave contributions with parallel wavefronts. At 57mT the energy in the parallel wavenumber spectrum is distributed among wavenumber components pointing both parallel and antiparallel to the magnetic field vector. In both cases turbulent structures arise preferable on the low field side of the torus. Some results on a novel field in plasma turbulence are given, i.e. the study of turbulence as a function of resonant magnetic field perturbations leading to the formation of magnetic islands. Magnetic islands in WEGA can be manipulated by external perturbation coils. A significant influence of field perturbations on the turbulence dynamics is found. A distinct local increase of the fluctuation amplitude and the associated turbulent particle flux is found in the region of magnetic islands.
The present experimental work investigates plasma turbulence in the edge region of magnetized high-temperature plasmas. A main topic is the turbulent dynamics parallel to the magnetic field, where hitherto only a small data basis existed, especially for very long scale lengths in the order of ten of meters. A second point of special interest is the coupling of the dynamics parallel and perpendicular to the magnetic field. This anisotropic turbulent dynamics is investigated by two different approaches. Firstly, spatially and temporally high-resolution measurements of fluctuating plasma parameters are investigated by means of two-point correlation analysis. Secondly, the propagation of signals externally imposed into the turbulent plasma background is studied. For both approaches, Langmuir probe arrays were utilized for diagnostic purposes. The main findings can be summarized as follows: Greatly elongated fluctuation structures exist in plasma edge turbulence. The structures are aligned along the confining magnetic field (k|| = 0). The correlation degree of fluctuations for a short connection length of 0.75m is greater than 80%. For much longer connection lengths of 23m and 66m, the correlation degree is reduced to approximately 40%. A conceptual interpretation of these observations is the coexistence of two different fluctuation components. One component has a correlation length parallel to the magnetic field below 20m and the other component a correlation length greater than 70m. Sine signals in the frequency range 1-100 kHz were injected into the turbulent plasma background. The propagation parallel and perpendicular to the magnetic field of the signals was studied. In poloidal direction, an asymmetry is observed, that can be explained by a copropagation of the signal with the background E × B-rotation of the plasma. The signal propagation parallel to the magnetic field shows no such asymmetry. As an advanced approach, spatio-temporal wave patters were injected into the edge plasma. The waves launched that way can be seen as test waves' in a turbulent background. The coupling strength of the imposed wave patterns to the background turbulence relies on the match of the imposed waves to the dynamics of turbulent structures. If the propagation direction of the imposed waves is parallel to the propagation direction of the background plasma, improved coupling is observed. This finding underlines the importance of the background plasma rotation for future attempts of controlling the plasma edge turbulence. Further optimization of frequency and wave vector of the imposed waves is probably a promising approach for achieving a significant and systematic influence of turbulence. Taking into account the present experimental state-of-the-art, for a deeper insight into the mechanism of the plasma edge turbulence of magnetized high-temperature plasmas a joint effort of numerical modeling and experimental results is a valuable approach. Such a cooperation should cover the explanation of the correlation observations as well as the experiments on signal injection into background turbulence. A quantitative comparison between the results presented in this work and a dedicated numerical drift wave simulation would be a significant step forward to a better understanding of plasma edge turbulence.
Two main aspects concerning drift wave dynamics in linear, magnetized plasma devices are addressed in the work: In part I of the thesis, drift waves are studied in a helicon plasma. The plasma parameter regime is characterized by comparably high collision frequencies and comparably high plasma-p exceeding the electron-ion mass ratio. Single Langmuir probes and a poloidal probe array are used for spatiotemporal studies of drift waves as well as for characterization of background plasma parameters. The main goals are the identification of a low-frequency instability and its major destabilization mechanisms. All experimentally observed features of the instability were found to be consistent with drift waves. A new code, based on a non-local cylindrical linear model for the drift wave dispersion, was used to gain more insight into the dominating destabilzation mechanisms, and also into dependencies of mode frequencies and growth rates on different parameters. In the experiment and in the numerical model, poloidal mode structures were found to be sheared. Part II of the thesis reports about mode-selective spatiotemporal synchronization of drift wave dynamics in a low-P plasma. Active control of the fluctuations is achieved by driving a preselected drift mode to the expense of other modes and broadband turbulence. It is demonstrated that only if a resonance between the driver signal and the drift waves in both space and time is reached, the driver has a strong influence on the drift wave dynamics. The synchronization effect is qualitatively well reproduced in a numerical simulation based on a Hasegawa-Wakatani model.