Doctoral Thesis
Refine
Year of publication
Document Type
- Doctoral Thesis (26) (remove)
Has Fulltext
- yes (26)
Is part of the Bibliography
- no (26)
Keywords
- Plasmaphysik (26) (remove)
Institute
In course of the recent results from Wendelstein 7-X, stellarators are on the brink for assessing their maturity as a fusion reactor. To this end, stellarator specific transport regimes need detailed exploration both with appropriate systematic experimental investigations and models. A way to enhance the efficiency of this process is seen in an systematic evaluation of existing experimental data. We propose appropriate tools developed in information theory for examining large datasets. Information entropy calculations, that have proven to assist the systematic assessment of datasets in many other scientific fields, are used for novelty detection.
Potentially, as a first use-case of this holistic process, this thesis attempts to link and to develop approaches to examine the stellarator specific core-electron-root-confinement (CERC) regime. The specific interest for CERC emerges from the behavior of the radial electric field. While ion-root conditions exhibit negative radial electric fields, CERC’s positive field in the very core of fusion grade plasmas adds an outward thermodynamic force to high-Z impurities and could add to potential actuators to control impurity influx as to be examined for full-metal wall operation in large stellarators. Recently, this feature received revived intent for reactor scale stellarators.
Also, in this work, parameter regions close to the transition from ion-root to CERC are
examined. At lower rotational transform (a characteristic feature of the magnetic field confining fusion grade plasmas), transitions were detected when the plasma current evolved. As in smaller stellarators, it is concluded that low-order rationals and magnetic islands are related to the transitions. This is widely supported by extensive MHD simulations which finally provide indications for the role of zonal flow oscillations. As one of the outcomes, gyrokinetic instabilities are seen interacting for the first time with the neoclassical mechanisms in experiments.
In order to cope with the vast number of highly sampled spatio-temporal plasma data, new
techniques for novelty detection are required. Fundamental prerequisites for the detailed
physics investigations were the feasibility study of entropy-based data analysis techniques, and their adaptation to detect previously unrevealed transition mechanisms. These tools were applied to multivariate bulk plasma emissivity data, which allowed the exploration of large parameter spaces and provided insights in the spatio-temporal dynamics of CERC transitions.
In this manner, this research highlights the feasibility of information flow measure analysis in fusion studies. Applications of different entropy-based complexity measures are explored and this work sheds light on the capabilities, added value and limitations of these techniques. This investigation presents the integration of information flow measures to gain deeper understanding of plasma transport phenomena, by providing an approach to fast systematic data mining suited for real-time analysis. This work paves the way for further development and implementation of information-theoretic methods for plasma data analysis.
In summary, this research highlights the gained insight on CERC transitions, while showcasing the feasibility, added values and limitations of information flow measure analysis for fusion studies, to induce theory based analysis revealing new insights in fundamental, stellarator-specific transport mechanisms.
This thesis describes how the data of the Langmuir probes in the Wendelstein 7-X (W7X) Test Divertor Unit (TDU) were evaluated, checked for consistency with other diagnostics and used to analyse plasma detachment.
Langmuir probes are an electronic diagnostic, and were among the first to be used in plasma physics to determine particle fluxes, potentials, temperatures and densities.
W7X is a large, advanced stellarator, magnetic confinement fusion experiment, operated at the Max-Planck-Institut for Plasma Physics(IPP) in Greifswald, Germany.
Its TDU is an uncooled graphite component, shaped and positioned to intercept the convective heat load of the plasma.
Detachment describes a desirable operation state of strongly reduced loads on this component.
The evaluation of Langmuir probe data relies heavily on models of the sheath, formed at the interface between plasma and a solid surface, to infer plasma parameters from the directly measured quantities.
Multiple such models are analysed, generalised, and adapted to our use case.
A detailed comparison is made to determine the most suitable model, as this choice strongly affects the predicted parameters.
Special attention is paid to uncertainties on the parameters, which are determined using a Bayesian framework.
From the inferred parameters, heat and particle fluxes are calculated.
These are also indirectly measured by two other, camera-based diagnostic systems.
Observations are compared to test the validity of assumptions and calculations in the evaluation of all three diagnostics by checking their results for consistency.
The first comparison, with the infrared emission camera system, shows good agreement with theoretical predictions and reported measurements of the sheath transmission factor, for which we derive and measure a value in W7X.
Parameter dependencies in the quality of this agreement hint at remaining issues.
The second comparison, with the Hydrogen alpha photon flux camera system, shows significant discrepancy with expectations.
These are argued to originate from systematic differences in the measurement locations, which are quantified and related to the magnetic topology.
Langmuir probe observations of individual discharges are analysed to discuss conditions under which detachment occurs, transition into that state and fluctuations observed prior to and during it.
A spatial parametrisation of the data is developed and used to facilitate this.
These observations contribute to the larger aim of understanding particle balance control and fusion plasma edge processes.
In der vorliegenden Arbeit wurde die Argon Plasma Quelle „MiniJet-R“ von HHF-Elektronik, Aachen, auf ihre Eignung als medizinischer CAP-Generator und ihre Emission von UVC-Strahlung und NO2 untersucht. Dabei wurde die Emission von UVC-Strahlung auf ihre Abstandsabhängigkeit und ihre Winkelverteilung vermessen. Die UVC-Intensität nimmt im für Entfernungen bis ca. einer Plasmaflammenlänge weniger schnell als 1/r² ab. Erst ab Entfernungen die ca. zwei Plasmaflammenlängen entsprechen nimmt die UVC-Intensität mit 1/r² ab. Die Intensität ist über alle Winkel gleich verteilt, bis eine durch den apparativen Aufbau der Quelle bedingte Abschattung bei Winkeln ab 85° einsetzt. Weiter wurde die Abhängigkeit der UVC-Intensität von den Betriebsparametern Argon Gasdurchfluss und Power Level untersucht. Dabei wurde gezeigt, dass die UVC-Intensität mit steigendem Gasdurchfluss abnimmt. Bei der Charakterisierung der Power Level – Abhängigkeit zeigte sich, dass die UVC-Intensität bei Power Level 3 ein Minimum hat. Zur näheren Bestimmung der UVC-Strahlung des Plasmas wurde das Spektrum des Plasmas von 168nm bis 275nm aufgenommen. Durch den Vergleich des gemessenen Spektrums und eines berechneten NO-Spektrums konnte NO als Hauptquelle der UVC-Strahlung nachgewiesen werden. Mittels Chemilumineszens-Messung konnte außerhalb des Plasmas NO2 als verbleibende Komponente identifiziert werden, während NO nur in vernachlässigbaren Konzentrationen außerhalb des Plasmas nachweisbar war. In weiteren Messungen wurde die NO2-Erzeugung des Plasmas in Abhängigkeit der Betriebsparameter Gasdurchfluss und Powerlevel sowie die NO2-Konzentration in der Raumluft in Abhängigkeit vom Abstand und von der Richtung zum Plasma bestimmt. Dabei wurde nachgewiesen, dass bei ausreichendem Abstand zum Plasma die NO2-Konzentration unterhalb des Arbeitsplatzgrenzwertes liegt. Unter ungünstigen Betriebsbedingungen und in unmittelbarer Umgebung konnten allerdings auch erheblich höhere Konzentrationen festgestellt werden. Die gemessenen UVC-Intensitäten und NO2-Konzentrationen werden mit den geltenden Maximalwerten unter dem Aspekt der Arbeitsplatzsicherheit verglichen. Abschließend erfolgt eine Beurteilung der CAP-Quelle „Minijet-R“ und eine Beschreibung einer idealen CAP-Quelle.
An experimental investigation of particle parallel flows has been carried out at Wendelstein 7-X (W7-X), one of the most advanced stellarators in the world. The studies are restricted to the outermost plasma region, the scrape-off layer (SOL), which is shaped to tackle the exhaust problem in vision of future fusion reactors based on plasma magnetic confinement. The aim of the measurements is to set the basis for a physics analysis of the SOL dynamics by obtaining direct information on convective heat transport, together with the assessment of the predominant flow directions of the main plasma ions and of fusion-products or wall-released impurities. In this way, a better comprehension of the interplay between the transport parallel and perpendicular to the SOL field lines can be achieved, contributing to the understanding of the effectiveness of the island divertor configuration.
The chosen instrument for the experimental studies is the Coherence Imaging Spectroscopy (CIS) diagnostic, a camera-based interferometer capable of measuring 2D Doppler particle flows associated with a selected visible line from the plasma. The diagnostic is distinguished by its high time resolution and spatial coverage, allowing the visualisation and measurements of flow velocities for a full module of W7-X simultaneously. A CIS diagnostic has been fully designed for W7-X with an improved level of accuracy achieved thanks to the implementation of a new calibration source, a continuous-wave-emission tunable laser. The laser allowed a full characterization of the diagnostic and a frequent precise calibration, making the CIS system reliable for parallel flow investigations during the operational campaign OP1.2. The validity and importance of the CIS measurements have been further confirmed with dedicated simulation of the SOL plasma parameters by the EMC3-EIRENE code, and by comparisons with other edge diagnostics. The CIS results show the effects related to dynamical changes in the SOL due to impurity gas puffs or the development of a plasma current. Moreover, CIS can be used as a powerful tool to test the limits of the current theoretical models, for example in the case of forward and reversed field experiments.
With this thesis, studies which form the bedrock for the long term goal of first wall heat load control and optimization for the advanced stellarator Wendelstein 7-X are developed, described and put into context. It is laid out how reconstruction of features of the edge magnetic field from plasma facing component heat loads is an important first step and can successfully be achieved by artificial neural networks. A detailed study of plasma facing component heat load distribution, potential overloads and overload mitigation possibilities is made in first order approximation of the impact of the main plasma dynamic effects.
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.
Ion thrusters are Electric Propulsion systems used for satellites and space missions. Within
this work, the High Efficient Multistage Plasma Thruster (HEMP-T), patented by the
THALES group, is investigated. It relies on plasma production by magnetised electrons.
Since the confined plasma in the thruster channel is non-Maxwellian, the near-field plume
plasma is as well. Therefore, the Particle-In-Cell method combined with a Monte-Carlo
Collision model (PIC-MCC) is used to model both regions. In order to increase the sim-
ulated near-field plume region, a non-equidistant grid is utilised, motivated by the lower
plasma density in the plume. To minimise artificial self-forces at grid points bordered by
cells of different size a modified method for the electric field calculation was developed in
this thesis. In order to investigate the outer plume region, where electric field and collisions
are negligible, a ray-tracing Monte-Carlo model is used. With these simulation methods,
two main questions are addressed in this work.
What are the basic mechanisms for plasma confinement, plasma-wall-interaction
and thrust generation?
For the HEMP-T the plasma is confined by magnetic fields in the thruster channel, generated
by cylindrical permanent magnets with opposite polarity. Due to different Hall parameters,
electrons are magnetised, while ions are not. Therefore, the dominating electron transport
is parallel to the magnetic field lines. In the narrow cusp regions, the magnetic mirror effect
reduces the electron flux towards the wall and confines the electrons like in a magnetic
bottle. At the anode, propellant gas streams into the thruster channel, which gets ionised
by the electrons creating the plasma. As a result of the electron oscillation between the two
cusp regions, ionisation of the propellant gas is efficient.
The magnetic field configuration of the HEMP-T also influences the plasma potential inside
the thruster channel. Close to the symmetry axis, the mainly axial magnetic field results in
a flat potential. At the inner wall, the field configuration reduces the plasma wall interaction
to only the narrow cusp regions. Here, the floating potential of the dielectric channel wall
and its plasma sheath result in a rather low radial potential drop compared to the applied
anode potential. As a result, the electric potential is rather flat and impinging ions at the
thruster channel wall have energies below the sputter threshold energy of the wall material.
Therefore, no sputtering appears at the dielectric wall. At the thruster exit the confinement
by the magnetic field is weakened and the potential drops with nearly the full anode voltage.
The resulting electric field accelerates the generated ions into the plume and generate the
thrust, but they are also able to sputter surfaces. During terrestrial testing, sputteringat vacuum vessel walls leads to the production of impurities. The amount of back-flux
towards the channel exit is determined by the sputter yield of the vacuum chamber wall. A
large distance between thruster exit and vessel wall reduces the back-flux and smooths the
pattern of deposition inside the thruster channel. Dependent on their material, the evolving
deposited layers can get conductive, modify by this the potential distribution and reduce
the thrust.
For the HEMP-T, ions are mainly generated at high potential close to the applied anode
potential. Therefore, the accelerated ions producing the thrust gain the maximum energy
as observed in experiment. Ions emitted from the thruster into different angles in the
plume contribute mainly to the ion current at angles between 30 ◦ and 90 ◦ . They mainly
originate from ionisation at the thruster exit. The resulting angular distribution of the
ejected ion current is close to the one of the experiment, slightly shifted by about ten
degrees to higher emission angles. In front of the thruster exit, electrons are trapped by
electrostatics forces. This enhanced density allows ionisation and an additional electron
density structure establishes.
What are possible physics based ideas for optimisation of an ion thruster?
An optimised thruster should have a high ionisation rate inside the thruster channel, low
erosion and an ion angular distribution with small contributions at high angles for min-
imised thruster satellite interactions. In experiments, the HEMP-T satisfies already quite
nicely these requests. In the simulations, low erosion inside the thruster channel and angular
ion distributions close to the experimental data are demonstrated. However, the ionisation
efficiency is lower and radial ion losses are larger than in experiment. A possible explanation
of these differences is an underestimated transport perpendicular to the magnetic field lines,
well known for magnetised plasmas.
A successful example for an optimisation using numerical simulations is the reduction of
back-flux of sputtered impurities during terrestrial experiments by an improved set-up of
the vacuum vessel. The implementation of baffles reduces the back-flux towards the thruster
exit and therefore deposition inside the channel. These improvements were successfully im-
plemented in the experiment and showed a reduction of artefacts during long time measure-
ments. This leads to a stable performance, as it is expected in space.
In this thesis, size-sensitive phenomena of three-dimensional dust crystals emerged in a low temperature plasma are presented. Depending on the number of particles in the system phase transitions, collective vortex motions and large-scaled expansions can be observed. To investigate these fascinating effects an advanced experimental setup as well as new evaluation methods have been developed. This thesis will present these new techniques and the gained insights.
Achieving commercial production of electricity by magnetic confinement fusion requires improvements in energy and particle confinement. In order to better understand and optimise confinement, numerical simulations of plasma phenomena are useful. One particularly challenging regime is that in which long wavelength MHD phenomena interact with kinetic phenomena. In such a regime, global electromagnetic gyrokinetic simulations are necessary. In this regime, computational requirements have been excessive for Eulerian methods, while Particle-in-Cell (PIC) methods have been particularly badly affected by the "cancellation problem", a numerical problem resulting from the structure of the electromagnetic gyrokinetic equations. A number of researchers have been working on mitigating this problem with some significant successes. Another alternative to mitigating the problem is to move to a hybrid system of fluid and gyrokinetic equations. At the expense of reducing the physical content of the numerical model, particularly electron kinetic physics, it is possible in this way to perform global electromagnetic PIC simulations retaining ion gyrokinetic effects but eliminating the cancellation problem. The focus of this work has been the implementation of two such hybrid models into the gyrokinetic code EUTERPE. The two models treat electrons and the entire bulk plasma respectively as a fluid. Both models are additionally capable of considering the self-consistent interaction of an energetic ion species, described gyrokinetically, with the perturbed fields. These two models have been successfully benchmarked in linear growth rate and frequency against other codes for a Toroidal Alfvén Eigenmode (TAE) case. The m=1 internal kink mode, which is particularly challenging in terms of the fully gyrokinetic cancellation problem, has also been successfully benchmarked using the hybrid models with the MHD eigenvalue code CKA. Non-linear simulations in this TAE case have been performed confirming the analytical prediction of a quadratic relationship between the linear growth rate of the TAE and the saturated amplitude of the TAE for a range of moderate values of the linear growth rate. At higher linear growth rate, a slower scaling of saturated amplitude with linear growth rate is observed. This analysis has been extended to include the non-linear wave-wave coupling between multiple TAE modes. It has been shown that wave-wave coupling results in a significant reduction in the saturated amplitude. It has been demonstrated that both plasma elongation and ion kinetic effects can exert a stabilising influence on the internal kink mode. A population of energetic particles can also exert a stabilising influence at low normalised pressure. At high normalised fast particle pressure the stabilised kink mode has been shown to give way to the m=1 EPM, which has been simulated both linearly and non-linearly (the "fishbone" mode). The first self-consistent simulations of global modes in the magnetic geometry of the optimised stellarator Wendelstein 7-X have been performed both linearly and non-linearly. Limitations have been encountered in performing simulations in 3D geometry. A hypothesis for the cause of these problems is outlined and ideas for mitigation are briefly described. In addition to the hybrid model simulations, some of the first utilisations of a new scheme for mitigating the cancellation problem in the fully gyrokinetic regime have been carried out in the framework of this thesis. This scheme, which was developed separately, is concisely described in this work. The new scheme has been benchmarked with existing gyrokinetic and hybrid results. The linear Wendelstein 7-X simulations and linear and single mode non-linear TAE simulations have been repeated with the new model. It is shown that bulk plasma kinetics can suppress the growth rate of global modes in Wendelstein 7-X. The results of fully gyrokinetic TAE simulations, the first to have been performed to our knowledge, are shown to be in close agreement with those results obtained using hybrid models. In the TAE case, the hybrid models are an order of magnitude less computationally demanding than the new gyrokinetic scheme, which is in turn at least an order of magnitude less computationally demanding than the previous gyrokinetic scheme.
This thesis is devoted to experiments on three-dimensional dust clouds which are confined in low temperature plasmas. Such ensembles of highly electrically charged micrometer-sized particles reveal fascinating physics, such as self-excited density waves and vortices. At the same time, these systems are challenging for experimental approaches due to their three-dimensional character. In this thesis, new optical diagnostics for dusty plasmas have been developed and, in combination with existing techniques, have been used to study these 3D dusty plasmas on different size and time scales.