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
The collisionless tearing mode is investigated by means of the delta-f PIC code EUTERPE solving the gyrokinetic equation. In this thesis the first simulations of electromagnetic non-ideal MHD modes in a slab geometry with EUTERPE are presented. Linear simulations are carried out in the cases of vanishing and finite temperature gradients. Both cases are benchmarked using a shooting method showing that EUTERPE simulates the linearly unstable tearing mode to a very high accuracy. In the case of finite diamagnetic effects and values of the linear stability parameter Delta of order unity analytic predictions of the linear dispersion relation are compared with simulation results. The comparison validates the analytic results in this parameter range. Nonlinear single-mode simulations are performed in the small- to medium-Delta range measuring the dependency of the saturated island half width on the equilibrium current width. The results are compared with an analytic prediction obtained with a kinetic electromagnetic model. In this thesis the first simulation results in the regime of fast nonlinear reconnection~(medium- to high-Delta range) are presented using the standard gyrokinetic equation. In this regime a nonlinear critical threshold has been found dividing the saturated mode from the super-exponential phase for medium-Delta values. This critical threshold has been proven to occur in two slab equilibria frequently used for reconnection scenarios. Either changing the width of the equilibrium current or the wave number of the most unstable mode makes the threshold apparent. Extensive parameter studies including the variation of the domain extensions as well as the equilibrium current width are dedicated to a comprehensive overview of the critical threshold in a wide range of parameters. Additionally, a second critical threshold for high-Delta equilibria has been observed. A detailed comparison between a compressible gyrofluid code and EUTERPE is carried out. The two models are compared with each other in the linear regime by measuring growth rates over wave numbers of the most unstable mode for two setups of parameters. Analytical scaling predictions of the dispersion relation relevant to the low-Delta regime are discussed. Employing nonlinear simulations of both codes the saturated island half width and oscillation frequency of the magnetic islands are compared in the small-Delta range. Both models agree very well in the limit of marginal instability and differ slightly with decreasing wave vector. Recently, the full polarisation response in the quasi-neutrality equation was implemented in EUTERPE using the Padé approximation of the full gyrokinetic polarisation term. Linear simulation results including finite ratios of ion to electron temperature are benchmarked with the dispersion relation obtained from a hybrid model. Finite temperature effects influence the saturated island width slightly with increasing ion to electron temperature ratio which has been verified by both models.
Im ersten Teil der Arbeit wird der erfolgreiche Aufbau einer Diagnostik zur quantitativen Bestimmung von Oberflächenladungsdichten beschrieben. Das Messprinzip bedient sich des elektro-optischen Pockelseffekts eines BSO-Kristalls, der in der Entladungszelle als Dielektrikum eingesetzt ist. Diese Methode arbeitet zeitlich und lateral aufgelöst, was die Untersuchung der Dynamik von Oberflächenladungen auf drei verschiedenen Zeitskalen ermöglicht. Die erste Zeitskala liegt in der Größenordnung von einigen 100 ns. Damit kann erstmals die Deposition von elektrischer Ladung auf einer dielektrischen Oberfläche während eines Entladungsdurchbruchs beobachtet werden. Die Deposition beginnt im Zentrum eines zuvor deponierten Ladungsspots. Die Polarität der neudeponierten Ladung ist der des ursprünglichen Ladungsspots entgegengesetzt. Die Folge ist, dass die absolute Ladungsdichte im Zentrum im Verlauf einiger hundert Nanosekunden kleiner wird als in den Randbereichen. Der Umladungsprozess wird so lange fortgesetzt, bis das elektrische Feld der neu deponierten Ladungen dem äußeren Feld so stark entgegenwirkt, dass die Spannung zur Aufrechterhaltung der Entladung unterschritten wird und die Entladung erlischt. Die zweite untersuchte Zeitskala liegt in der Größenordnung der Periodendauer der externen Spannung. Im Nulldurchgang der Spannung liegen zeitlich stationäre Ladungsdichteverteilungen auf dem Dielektrikum vor. Die Geometrie eines mittleren Ladungsspots wird in Abhängigkeit der anliegenden Spannungen und des Gasdrucks untersucht. Einerseits ist der Spotradius abhängig von den Ionisationsprozessen im Volumen, weil die Dichte der Raumladungen die Stärke des Elektronenfokus in das Innere der Entladung steuert. Andererseits wird die Spotbildung durch eine laterale Drift von Ladungsträgern kurz vor der Oberfläche aufgrund des elektrischen Feldes deponierter Ladungsträger beeinflusst. Die dritte untersuchte Zeitskala liegt in einer Größenordnung von Sekunden. Im Fall einer initial homogenen Oberflächenladungsverteilung nimmt die mittlere Ladungsdichte in einer Größenordnung von Sekunden monoton ab. Dieser Prozess stellt einen Ladungsabbau dar, dessen zeitliches Verhalten durch zwei überlagerte Exponentialfunktionen beschreiben ließ. Dadurch werden zwei Ladungsträgerpopulationen im BSO angenommen, die verschieden abgebaut werden. Im Fall einer initial inhomogenen Ladungsdichteverteilung wird ein Transport elektrischer Ladung auf der BSO-Oberfläche in einer Größenordnung von Sekunden beobachtet. Es wird weiterhin erstmals die durch einen Atmosphärendruck-Plasmajet deponierten Ladungen auf BSO zeitaufgelöst gemessen. Die zeitliche Entwicklung der Oberflächenladungen kann mit der Messung des elektrischen Stroms an einer der Ringelektroden des Jets korreliert werden. Dadurch wird geschlossen, dass der Ladungsaustauch nicht direkt durch einen Bullet verursacht wird. Er erzeugt stattdessen einen elektrisch leitfähigen Kanal zwischen der Düse des Jets zur BSO-Oberfläche. Infolgedessen kann Ladung, die sich auf der Innenseite der Jetkapillare befindet, auf den BSO-Kristall transportiert werden. Im zweiten Teil der Arbeit werden Kenngrößen entwickelt, die den Ordnungszustand einer aus Einzelobjekten zusammengesetzten Entladungsstruktur quantitativ beschreiben. Die Kenngrößen werten dabei die laterale Leuchtdichteverteilung der Entladungsemisssion, u.a. auf Basis der Tripel-Korrelationsfunktion. Dabei werden zwei separate Bifurkationsspannungen zwischen einer hexagonalen und einer ungeordneten Anordnung beobachtet: Bei der Verringerung der Spannung wird zunächst der Bifurkationspunkt der azimutalen Ordnung durchlaufen und anschließend der Bifurkationspunkt der radialen Ordnung. Die Systeme gehen jeweils in einen Zustand geringerer Ordnung über. Die Ursache des Ordnungsverlusts ist das zunehmende Fehlen von Entladungsspots, was im Mittel zu einer geringeren Wechselwirkung der Spots untereinander führt und das System an Freiheitsgraden gewinnt. Im dritten Teil dieser Arbeit wird erstmals ein Ansatz verfolgt, der die Steuerung lateral strukturierter Entladungen ermöglicht. Dafür wurde ein Aufbau konstruiert, bei dem ein gekühlter Halbleiter als Dielektrikum in der Entladungszelle dient. Dessen externe Beleuchtung führt bei einer anliegenden Spannung zu einer Änderung des Spannungsteilerverhältnisses der kapazitiven Elemente und schließlich zu einer lokalen Erhöhung der Spannung über dem Entladungsraum. Die Größe und Leuchtintensität der durch die Beleuchtung gezündeten Entladung ist stark abhängig von der beleuchteten Fläche, der Leistungsdichte der Beleuchtung und der anliegenden Spannung.
Diese Dissertation beschäftigt sich mit der Erzeugung von edelmetallfreien Katalysatoren für die Sauerstoffreduktion in Brennstoffzellen. Dabei wird ein neuartiger, dualer Plasmaprozess entwickelt, aufgebaut und die so-erzeugten Schichten mit verschiedenen elektrochemischen (CV, RDE und RRDE) und strukturanalytischen Methoden (SEM, EDX, IR, XPS, Leitfähigkeit, XRD, NEXAFS, EXAFS und TEM) untersucht. Auf diese Weise ist es erstmalig gelungen edelmetallfreie Katalysatoren mit einem Plasmaprozess herzustellen, ohne dass eine zusätzliche Pyrolyse benötigt wird. Die katalytische Aktivität der Schichten ist außerdem deutlich höher als die von rein chemisch hergestellten Metall–Polypyrrol-Schichten.
Magnetic reconnection is a ubiquitous phenomenon observed in a wide range of magnetized plasmas from magnetic confinement fusion devices to space plasmas in the magnetotail. The process enables the release of accumulated magnetic energy by rapid changes in magnetic topology, heating the plasma in the vicinity of the reconnection site, generating fast particles and allowing a wealth of instabilities to grow. This thesis reports on the results from a newly constructed linear, cylindrical and modular guide field reconnection experiment with highly reproducible events, VINETA.II. A detailed analysis of the reconnecting current sheet properties on a macroscopic and microscopic scale in time and space is presented. In the experiment, four parallel axial wires create a figure-eight in-plane magnetic field with an X-line along the central axis, as well as an axial inductive field that drives magnetic reconnection. Particle-in-cell simulations show that the axial current is limited by sheaths at the boundaries and that electrostatic fields along the device axis always set up in response to the induced electric field. Current sheet formation requires an additional electron current source, realized as a plasma gun, which discharges into a homogeneous background plasma created by a rf antenna. The evolution of the plasma current is found to be dominantly set by its electrical circuit. The current response to the applied electric field is mainly inductive, which in turn strongly influences the reconnection rate. The three-dimensional distribution of the current sheet is determined by the magnetic mapping of the plasma gun along the sheared magnetic field lines, as well as by radial cross-field expansion. This expansion is due to a lack of equilibrium in the in-plane force balance. Resistive diffusion of the magnetic field by E=η j is found to be by far insufficient to account for the high reconnection rate E=-dΨ/dt at the X-line, indicating the presence of large electrostatic fields which do not contribute to dissipative reconnection. High-frequency magnetic fluctuations are observed throughout the current sheet which are compared to qualitatively similar observations in the Magnetic Reconnection Experiment (MRX, Princeton). The turbulent fluctuation spectra in both experiments display a spectral kink near the lower hybrid frequency, indicating the presence of lower hybrid type instabilities. In contrast to the expected perpendicular propagation of mainly electrostatic waves, an electromagnetic wave is found in VINETA.II that propagates along the guide field and matches the whistler wave dispersion. Good correlation is observed between the local axial current density and the fluctuation amplitude across the azimuthal plane. Instabilities driven by parallel drifts can be excluded due to the large required drift velocities or low resulting phase velocities that are not observed. It is instead suggested that a perpendicular, electrostatic lower hybrid mode indeed exists that resonantly excites a parallel, electromagnetic whistler wave through linear mode conversion. The resulting fluctuations are found to be intrinsic to the localized current sheet and are independent of the slower reconnection dynamics. Their amplitude is small compared to the in-plane fields, and have a negligible contribution to anomalous resistivity through momentum transport in the present parameter regime.
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.
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.
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.
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.
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.
The main objective of this work is to contribute to the understanding of the grafting of nitrogen and amino surface functional groups on polymers by means of plasmas containing nitrogen and hydrogen. For this purpose, many aspects of plasma surface modification were studied. In the frame of this work, a new, UHV-sealed plasma reactor system was put into operation. The system is special for its clean reaction environment and the possibility to perform quasi in situ XPS measurements. A comparison of the UHV system to a fine vacuum reactor showed that a clean reaction environment is mandatory for reproducible plasma processing and efficient nitrogen and amino functionalisation. A key motivation for the present work was the observation that the non-coating plasma processes reported in literature fail to graft primary amino groups on polymer surfaces with densities that significantly exceed 3 - 4% NH2/C. In order to investigate this phenomenon in detail, this work followed two experimental tracks: On the one hand, a broad systematic study of plasma processing parameters was performed. On the other, the surface diagnostics methods used for the quantification of amino groups were critically reviewed. For this, a numerical algorithm was developed to reconstruct the element depth profile from angle-resolved XPS data. In the scope of the process parameter study, cw and pulsed microwave (MW) plasma excitation was compared to radio-frequency (RF) excitation. The home-built MW source was studied and optimised with respect to ignition behaviour and power efficiency. The performance of the MW and RF plasmas in polymer surface modifications was studied in various gas mixtures containing NH3 and H,, or N2 and H,. Also the differences of glow and afterglow processing of polymers were investigated. Large variations of the nitrogen and primary amino grafting efficiencies were obtained. They triggered a number of new ideas for the underlying reaction mechanisms. Special attendance was devoted to the selectivity of the functionalisation processes for primary amino groups. Nitrogen-containing discharges that were rich in hydrogen achieved selectivities up to 100%. The upper limit of 3 - 4% amino groups on the surface, however, was not passed. Angle-resolved XPS measurements revealed a systematic problem for the definition of a surface density, which is capable of explaining the upper limit for amino groups. It is either due to a limited labelling depth of amino groups by the applied TFBA derivatisation reaction, or to a limited functionalisation depth of the plasma process. One very efficient nitrogen-grafting plasma process that was developed on polystyrene was applied to seven other unfluorinated polymers. The similarity of the resulting functionalisation demonstrated a good transfer-ability of plasma surface functionalisation processes. Plasma treatments of polymer surfaces, especially in hydrogen-containing gases, are known to be generally followed by uncontrollable oxidation phenomena. The properties of plasma-functionalised polymer surfaces were therefore studied in conjunction with ageing effects. Quasi in situ XPS analysis allowed to distinguish the influence of oxygen contamination during the plasma process from post-process oxidation due to contact of plasma-treated samples to atmospheric oxygen. The surface modification experiments were accompanied by several gas phase diagnostic techniques. In the scope of this work, the UHV reactor system was equipped with optical emission spectroscopy (OES), two-photon absorption laser-induced fluorescence (TALIF), and tunable diode laser absorption spectroscopy (TDLAS). A separate plasma source was setup to perform an absolute quantification of the vacuum-ultra-violet (VUV) emission intensity of hydrogen-containing MW-excited plasmas. The techniques were evaluated with respect to their contribution to an understanding of the plasma processing of polymers. The rich experimental data allowed to suggest new reaction mechanisms for the grafting of nitrogen- and amino functional groups. Surface passivation experiments in H, plasmas of nitrogen-functionalised surfaces initiated a re-evaluation and an extension of the mechanism of selective etching [1]. Together with two other new reaction mechanisms, a hypothetical reaction scheme was suggested. It was studied by the help of two numerical models for heterogenous reactions of radicals with the surface. In order to avoid the complexity of the fragmentation process of NH,, the models were restricted to discharges in N, and H9. Despite the sparse information on the composition of the gas phase, the data of two experimental series showed a very particular phenomenology that allowed a first test of the model. The test supports the newly-suggested reaction mechanisms. Especially the role of NH2 attachment to open reaction sites for the grafting of amino groups was emphasised. A more stringent test of the model is left to future experiments with extended gas phase diagnostic means.