Doctoral Thesis
Refine
Year of publication
Document Type
- Doctoral Thesis (13) (remove)
Language
- English (13)
Has Fulltext
- yes (13)
Is part of the Bibliography
- no (13)
Keywords
- Stellarator (13) (remove)
Institute
- Institut für Physik (13)
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.
This work investigates turbulence in the core plasma of the optimised stellarator
Wendelstein 7-X. It focuses on experimental characterisation and
evaluation of the electrostatic micro-instabilities, which drive turbulent fluctuations,
and the saturation of turbulence by zonal flows. Expectations for
Wendelstein 7-X are formulated by reviewing theoretical work and with
the help of gyrokinetic simulations. The experimental analysis centres on
line-integrated density fluctuation measurements with the phase contrast
imagining diagnostic in electron cyclotron heated hydrogen discharges. An
absolute amplitude calibration was implemented, and a method for reliable
determination of dominant phase velocities in wavenumber-frequency
spectra of density fluctuations has been developed. Line-averaged density
fluctuation levels are observed to vary between magnetic configurations.
The wavenumber spectra exhibit a dual cascade structure, indicating fully
developed turbulence. The dominant instability driving turbulent density
fluctuations on transport relevant scales is identified as ion-temperaturegradient-
driven modes, which are mainly localised in the edge region of the
confined plasma. Despite the line-integrated nature of the measurement, the
localisation of density fluctuations is shown by comparing their dominant
phase velocity with the radial profile of the E × B rotation velocity due to
the ambipolar neoclassical electric field. Nonlinear gyrokinetic simulations
and a simplified plasma rotation model within a synthetic diagnostic confirm
the localisation. Oscillations of the dominant phase velocity indicate
the existence of zonal flows as a saturation mechanism of ion-temperaturegradient-
driven turbulence. A direct effect on turbulent density fluctuation
amplitudes and radial transport is observed.
In the present thesis, a systematic study of beam driven Alfvén eigenmodes in high-density and low-temperature plasmas of the W7-AS stellarator is performed. The device went out of operation in 2002 and the study is based on stored experimental data. Alfvén instabilities can roughly be divided into ideal MHD Alfvén eigenmodes and those existing due to kinetic effects. The spectrum of ideal MHD Alfvén waves in toroidal fusion devices consists of a continuum of stable waves that are strongly localized. Weakly damped, discrete eigenmodes can exist in gaps of the continuous spectrum which are formed by plasma inhomogeneities and the coupling of Alfvén continua. This allows an identification of ideal MHD Alfvén eigenmodes in terms of their frequency and mode numbers. Kinetic effects can modify this spectrum and cause additional types of eigenmodes, the kinetic Alfvén eigenmodes (KAE) and energetic particle modes (EPM). The goal of this thesis is twofold: (I) identification and description of fast particle driven Alfvén instabilities in W7-AS, and (II) study of energetic particle losses induced by Alfvén instabilities. The reconstruction of the ideal MHD plasma equilibrium for each discharge with sufficient accuracy is the very foundation of all subsequent steps. This is achieved, based on measured plasma parameter profiles that are further refined by validating them to the measurements of other, independent plasma diagnostics. The applied scheme is inspired by an approach of Integrated Data Analysis (IDA) to combine different diagnostic data and provide combined uncertainties. After mode number analysis and eigenmode identification, the theoretically expected, linear growth rate of the instability is calculated where possible, and the various contributions of the fast particle drive to the instability of the mode are identified. Alfvénic activity recorded by the Mirnov diagnostic is analyzed, which consists of a set of spatially distributed coils that measure magnetic fluctuations. On W7-AS, the probes are arranged in three poloidal arrays at different toroidal positions. The spacing between the probes is non-equidistant. In addition, the signals of one probe array are digitized with a different sample rate. These characteristics prohibit the straight-forward use of standard tools available for harmonic analysis. Instead, a new tool has been developed and thoroughly tested. It is a multi-dimensional extension of the Lomb periodogram, able to provide reliable time-resolved frequency and mode number spectra in the case of uneven datapoint spacing. Numerical studies of this periodogram show a good performance with respect to mode number resolution given the low number of available probes, and robustness against perturbations of the signal. Only two of the probe arrays can be used for the analysis of eigenmodes with frequencies >70 kHz, such that for high-frequency phenomena insufficient information about the mode numbers is available. A total of 133 different Alfvén eigenmodes is studied in discharges from different experimental campaigns. A restriction to discharges from various high-beta campaigns with neutral beam heating is required to allow for a realistic reconstruction of plasma equilibrium and velocity distribution functions of energetic particles. The discharges are characterized by high density, ne = 5 x 1019 m-3 to 2.5 x 1020 m-3 at relatively low temperatures of Te = Ti = 150 ... 600 eV. Alfvén eigenmodes often appear transiently in the startup phase of these discharges, where density and heating power are being ramped up. Occasionally, Alfvén eigenmodes are seen in the stationary, high-beta phase in the presence of considerable neutral beam heating. Most of the Alfvén eigenmodes are successfully classified as ideal MHD eigenmodes. 19 global, 47 toroidicity-induced and 8 ellipticity-induced Alfvén eigenmodes (GAEs, TAEs, and EAEs, respectively) are unambiguously identified by their mode numbers and frequencies. Excellent agreement between experimentally observed mode number spectra and theoretically calculated eigenmode structure is shown for a TAE example. Additional 13 events are found to have frequencies inside the EAE gap and could possibly be EAEs. Evidence for high-frequency Alfvén eigenmodes (mirror- and helicity-induced Alfvén eigenmodes) is seen, but can not be proven rigorously due to uncertain mode numbers and the complexity of the Alfvén continuum. The remaining 41 Alfvén eigenmodes can not be classified to be one of the above cases. Reasons are either high frequencies, mode numbers obscured by far-field effects, or mode numbers that could not be related to ideal MHD Alfvén eigenmodes. A selection of these shows indications of strong non-linear wave-particle interactions and are assumed to be EPMs. Kinetic Alfvén eigenmodes are not expected to exist in the experimental conditions that were studied. The radially resolved velocity distribution function is used to describe the parameter regimes in which the modes are observed in terms of the dimensionless parameters vb/vA (beam velocity normalized to the Alfvén velocity) and ßfast/ßth, where beta is the ratio of plasma pressure to magnetic pressure. The first parameter describes through which of the possible resonance velocities particles can interact with the eigenmode. A peculiarity of the fast particle dynamics in fusion devices is that they can resonantly interact with Alfvén eigenmodes through sideband resonances even if v < vA. The second parameter describes the energy content of the destabilizing fast particle population compared to the potentially stabilizing thermal plasma component. These parameters contain relevant information about the instability of an eigenmode and such diagrams are given for all observed modes. In addition to that, the expected linear growth rate of gap modes is calculated based on a theoretical model that extends the ideal MHD by a perturbative, drift-kinetic description of the energy exchange between waves and circulating particles, neglecting the effects of trapped particles. For the discharges under consideration the thermal electron speed is comparable to vA and the electrons provide a significant Landau damping contribution. Due to strong density gradients near the plasma boundary in most of the discharges, the thermal ions can provide a small drive via the spatial inhomogeneity which does not overcome the electron damping, however. The drive by spatial inhomogeneity of thermal ions requires a certain propagation direction of the mode and is equally stabilizing for opposite mode numbers. The fast particles also contribute to the growth rate via spatial inhomogeneity, velocity gradients and velocity anisotropy terms are negligible in W7-AS. Most of the observed GAE or EAE modes have negative mode numbers, which correspond to a propagation direction for which the spatial inhomogeneity of thermal and beam ions is predicted to be stabilizing. A fast particle drive of these modes is not confirmed, whereas the TAEs are found to be strongly destabilized by neutral beam injection. The distribution of plasma parameters for discharges showing TAEs in terms of the dimensionless stability parameters suggests an instability threshold that is qualitatively confirmed by an exploration of the parameter space with the theoretical model. Wave-induced, resonant losses of energetic ions scale linearly with the wave amplitude. To identify them, correlations between ion loss probe signals and wave amplitudes are searched, where correlation times in the order of the slowing-down time of energetic particles are expected. Significant correlations can be established only exceptionally for 3 of the identified ideal MHD Alfvén eigenmodes. Those Alfvén eigenmodes, however, which are assumed to be EPMs frequently show severe losses of energetic ions that are visible in the time traces of the plasma energy as well.
The WEGA stellarator is used to confine low temperature, overdense (densities exceeding the cut-off density of the heating wave) plasmas by magnetic fields in the range of B=50-500 mT. Microwave heating systems are used to ignite gas discharges using hydrogen, helium, neon or argon as working gases. The produced plasmas have been analyzed using Langmuir and emissive probes, a single-channel interferometer and ultra-high resolution Doppler spectroscopy. For a typical argon discharge in the low field operation, B=56 mT, the maximum electron density is n_e~10^18m^{-3} with temperatures in the range of T=4-12 eV. The plasma parameters are determined by using Langmuir probes and are cross-checked with interferometry. It is demonstrated within this work that the joint use of emissive probes and ultra-high resolution Doppler spectroscopy allows a precise measurement of the radial electric field. Here the floating potential measurements using emissive probes have been compared to measurements of the poloidal rotation of the plasma which is also linked to the radial electric field. In order to alter the plasma parameters a biasing probe setup has been used during this work. The focus of this work is on demonstrating the ability to modify the existing radial electric field in a plasma by using the biasing probe. This technique is in principle not new, as it has been around for decades. Looking at details, it turns out that describing low field operation WEGA argon plasmas in connection with biasing is not covered by the present set of theoretical approaches and experimental cognition. This work will commence with a basic approach and first establishes the diagnostic tools in a well-known discharge. Then the perturbation caused by the biasing probe is assessed. Following the characterization of the unperturbed plasmas, plasma states altered by the operation of the energized biasing probe will be characterized. It is demonstrated that modifying the existing radial electric field can be achieved and reliably diagnosed using spectroscopy and probe measurements. In order to verify the different approaches for determining the radial electric field the diagnostics are cross-checked against another whenever possible. During biasing the plasma two different stable plasma states have been found. Stable here refers to the state existing much longer than the confinement time for WEGA. The presence of a calorimetric limiter placed in the scrape-off layer has an impact on the type of the plasma state. The two observed plasma states differ in plasma parameter profiles, such as density, temperature, electric field and confined energy. The results are compared to two simple models. One model relies on the relevant atomic processes and a second one is based on neoclassical theory. Both models can be used to derive the particle and power flux from the plasma. The losses predicted by the atomic models can be tested using bolometry. It can be shown that both models agree well in the description of the particle balance of the electrons for large regions of the plasma. By comparing the models the neoclassical heat flux turns out to be small compared to the energy fluxes caused by atomic processes. For the reference discharge taking the energy flux due to the atomic processes and balancing it by the input microwave power is satisfying the energy balance, without the need for transport. For the biased discharges it turns out that neoclassical transport can be neglected as well, but the additional biasing power has to be taken into account. A simple model for the biasing power is motivated and tested. An agreement in the energy balance can be reached in this way as far as the models are applicable. The models also allow drawing conclusions on the amount of absorbed microwave power.
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.
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 dissertation focusses on the numerical modelling of resonant destabilization of Alfvén eigenmodes by fast ions in fusion plasmas. It especially addresses non-linear simulations of stellarator plasmas in which particle collisions are retained. It is shown that collisions are required for a realistic description of Alfvén waves in plasmas relevant to nuclear fusion.
We start by carefully verifying the implementation of the collision operators into the electromagnetic version of the gyro-kinetic delta-f particle-in-cell code EUTERPE. After these initial benchmarks are completed successfully, the code is in a position to be applied to realistic tokamak and stellarator scenarios.
Since every collision operator needs to fulfil conservation laws, a momentum-conserving version of the pitch-angle scattering operator is implemented. This is in particular important for neoclassical transport simulations aimed at computing flux-surface variations of the electrostatic potential in stellarators.
Using the simplified CKA-EUTERPE model (employing a fixed-mode-structure approximation), we perform non-linear simulations in tokamaks and stellarators. We show that the non-linear dynamics of fast-ion-driven Alfvén eigenmodes is significantly influenced by collisions. They have the potential to enhance the saturation level and to affect the frequency chirping of the modes.
It is thus concluded that collisions play an essential role in determining Alfvén-eigenmode-induced fast-ion transport - an important issue for future fusion devices. In order to address this issue the CKA-EUTERPE model is extended to evolve multiple modes at the same time. First results of this multi-mode version (which enhances the level of realism of the simulations) are shown in the Appendix of the thesis.
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 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.
There is a wide variety of Alfvén waves in tokamak and stellarator plasmas. While most of them are damped, some of the global eigenmodes can be driven unstable when they interact with energetic particles. By coupling the MHD code CKA with the gyrokinetic code EUTERPE, a hybrid kinetic-MHD model is created to describe this wave–particle interaction in stellarator geometry. In this thesis, the CKA-EUTERPE code package is presented. This numerical tool can be used for linear perturbative stability analysis of Alfvén waves in the presence of energetic particles. The equations for the hybrid model are based on the gyrokinetic equations. The fast particles are described with linearized gyrokinetic equations. The reduced MHD equations are derived by taking velocity moments of the gyrokinetic equations. An equation for describing the Alfvén waves is derived by combining the reduced MHD equations. The Alfvén wave equation can retain kinetic corrections. Considering the energy transfer between the particles and the waves, the stability of the waves can be calculated. Numerically, the Alfvén waves are calculated using the CKA code. The equations are solved as an eigenvalue problem to determine the frequency spectrum and the mode structure of the waves. The results of the MHD model are in good agreement with other sophisticated MHD codes. CKA results are shown for a JET and a W7-AS example. The linear version of the EUTERPE code is used to study the motion of energetic particles in the wavefield with fixed spatial structure, and harmonic oscillations in time. In EUTERPE, the gyrokinetic equations are discretized with a PIC scheme using the delta-f method, and both full orbit width and finite Larmor radius effects are included. The code is modified to be able to use the wavefield calculated externally by CKA. Different slowing-down distribution functions are also implemented. The work done by the electric field on the particles is measured to calculate the energy transfer between the particles and the wave and from that the growth rate is determined. The advantage of this approach is that the full magnetic geometry is retained without any limiting assumptions on guiding center orbits. Extensive benchmarks have been performed to test the new CKA-EUTERPE code. Three tokamak benchmarks are presented, where the stability of TAE modes are studied as a function of fast particle energy, or in one case as a function of the fast particle charge. The benchmarks show good agreement with other codes. Stellarator calculations were performed for Wendelstein 7-AS and the results demonstrate that the finite orbit width effects tend to be strongly stabilizing.