Refine
Year of publication
Document Type
- Doctoral Thesis (22)
Language
- English (22) (remove)
Has Fulltext
- yes (22)
Is part of the Bibliography
- no (22)
Keywords
- Plasmaphysik (22) (remove)
Institute
- Institut für Physik (22) (remove)
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.
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.
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.
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.
The development of innovative coatings with multifunctional properties is an ambitious task in modification of material surfaces. A novel approach is a hybrid method combining the non-thermal plasma processing with nanotechnology for the development of multifunctional surface coatings. The conception of the hybrid coating process is based on three steps: the preparation of a suspension consisting of an organic liquid and functional nanoparticles, the deposition of the suspension as a thin liquid film on the material surface, and the plasma modification of the liquid organic film to achieve a thin solid composite film with embedded nanoparticles demonstrating multifunctional properties and good adherence on the substrate material. In this work the liquid polydimethylsiloxane (PDMS) was applied as a model system, and the experimental investigations were focused on the PDMS plasma modification. In particular, the specific role of the different plasma components and the influence of the plasma and processing parameters on the PDMS modification were studied. The applied capacitively coupled radio frequency (CCRF) plasma was analyzed by electric probe measurements and optical emission spectroscopy, whereas the molecular changes in PDMS due to plasma-induced chemical reactions were studied by the Fourier transform infrared reflection absorption spectroscopy. Additionally, the photocatalytic activity of thin composite films consisting of plasma cross-linked PDMS with embedded TiO2 nanoparticles was demonstrated. During the investigation it was found that the CCRF discharge modifies efficiently thin liquid PDMS films to solid coatings. The samples were positioned in the plasma bulk at floating potential. The penetration depth of particles like neutrals, ions, electrons and radicals in the film is strongly limited. The heating of samples in the CCRF discharge is weak to modify PDMS by itself and only the plasma radiation is able to transform the liquid bulk to solid one. It is known that the absorption onset of PDMS lies in the VUV region (below 200 nm). The energetic VUV radiation penetrates into the PDMS film on a thickness from several hundred nanometers to few micrometers and initiates photochemical reactions there. Thus, different gases like Ar, Xe, O2, H2O, air and H2 were tested to provide the strongest VUV emission intensity of the CCRF discharge. Discharge pressure and power were varied for all these gases and it was found that at all conditions the H2 plasma demonstrates drastically stronger emission. Thus, H2 gas was selected for the plasma treatment of liquid PDMS films. The IRRAS analysis revealed the transformation process of PDMS with the degradation of CH3 groups, the formation of new groups like SiOH, CH2 and SiH, the formation of the SiOx material and crosslinking. It was found that the modification effect is not uniform across the film thickness. The top region with an initial thickness up to 100 nm loses all CH3 groups, in the underlying region the CH3 concentration increases gradually from zero to the value for PDMS, if the film was thick enough. The methyl-free SiOx top layer contains also SiOH and SiH groups. Furthermore, the SiH groups are concentrated only in a very thin layer with a thickness below 10 nm. The presence of the unscreened polar SiOSi and SiOH groups on the surface causes the adsorption of H2O from the atmosphere, which was also observed by IRRAS. By means of the spectroscopic ellipsometry it was found out that all above described regions experience a shrinking. The reason is the crosslinking and loss of material. The most shrunken layer is the top SiOx layer with the shrinking ratio (final thickness/initial thickness) of 0.55 - 0.60. Further, this ratio gradually rise up to the value of 0.95 in the deeper region, which has the concentration of CH3 groups of about that for PDMS. After the analysis of all results the depth of effective modification was estimated at 300 400 nm for the most optimal conditions. The optimization of the plasma VUV intensity was realized by variation of discharge pressure and power. The strongest plasma emission at studied conditions provided the irradiance of the sample of ca. 13 mW/cm2. However, such strong radiation causes very strong production rate of the gases. These products leave the modifying film slower as they are produced, what causes their accumulation in there. Their pressure grows up leading to formation of bubbles, which later explode. Finally, the film becomes heavily damaged. To avoid this effect the pressure and the RF power were changed to reduce the irradiance to 6 - 7 mW/cm2. This resulted in the absence of any damages.
In the last decade a new domain has developed in plasma physics: plasma medicine. Despite the successes that have already been achieved in this exciting new field, the interaction of plasmas with “biological materials” is not yet fully understood. Further investigations in particular with respect to the properties of the applied plasmas sources are therefore essential in order to decode this complex interaction process. Currently, a great variety of different discharge types are used in plasma medical investigation which are generally are operated in noble gases like helium and argon or with dry air. In the present work, the main focuses is on the diagnostics of reactive oxygen and nitrogen species (RONS) resulting from the plasma chemistry of an argon radio-frequency (RF) atmospheric pressure plasma jet (APPJ) and its interaction with the ambient atmosphere. To conduct this study, a commercially available plasma device, so-called kinpen is used due to its technical development maturity and its accessibility on the market. As a method of choice, diagnostic techniques are based on optical spectroscopy known to be a reliable tool to investigate plasmas. Consequently, three complementary optical laser diagnostics, namely quantum cascade laser absorption spectroscopy (QCLAS), laser induced fluorescence (LIF) and planar single shot LIF (PLIF), have been successfully applied to the plasma jet itself or its effluent. All of these diagnostics offer a high species selectivity and an excellent spatial and temporal resolution. They are used in this work for i) the characterization of the plasma chemical dynamics with respect to the generation of biological active RONS – in particular for the case of N2 and O2 admixtures. ii) the measurement of the NO density profile in the plasma effluent iii) the investigation of the flow characteristics of the neutral gas component (laminar vs. turbulent) and its influence on the plasma chemistry. Numerical analysis have been carried out in collaboration with PLASMANT (University of Antwerp) via kinetic simulations of the entire plasma chemistry. Expectingly, atomic oxygen (O) and nitric oxide (NO) turn out to be precursors of ozone (O3) and nitric dioxide (NO2). However, it was intriguing to unveil that atomic oxygen and nitrogen metastable (N2(A)) play together a key part --as intermediate species-- in the generation of more stable RONS, e.g. NO. The absolute density of NO space resolved was measured by LIF and absolutely calibrated molecular beam mass spectrometer. LIF was used to determine relative density of OH radical in the plasma plume. 2D-LIF was used to investigate the gas flow pattern with OH as a flow tracer. The results are discussed in details and show different operating mode of the jet, e.g. laminar or turbulent and that the plasma influences these regimes. The first detection and relative measurement by LIF of nitrogen metastable (N2(A)) produced by an argon APPJ is also shortly reported in this work. The outcome of this thesis will bring new insights in the field of argon APPJs chemistry and its interaction with the ambient atmosphere which can be valuable to support plasma modelling and to consider for the applications in plasma medicine.
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.
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.
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.
Particle and heat transport in fusion devices often exceed the neoclassical prediction. This anomalous transport is thought to be produced by turbulence caused by microinstabilities such as ion and electron-temperature-gradient (ITG/ETG) and trapped-electron-mode (TEM) instabilities, the latter ones known for being strongly influenced by collisions. Additionally, in stellarators, the neoclassical transport can be important in the core, and therefore investigation of the effects of collisions is an important field of study. Prior to this thesis, however, no gyrokinetic simulations retaining collisions had been performed in stellarator geometry. In this work, collisional effects were added to EUTERPE, a previously collisionless gyrokinetic code which utilizes the δ f method. To simulate the collisions, a pitch-angle scattering operator was employed, and its implementation was carried out following the methods proposed in [Takizuka & Abe 1977, Vernay Master's thesis 2008]. To test this implementation, the evolution of the distribution function in a homogeneous plasma was first simulated, where Legendre polynomials constitute eigenfunctions of the collision operator. Also, the solution of the Spitzer problem was reproduced for a cylinder and a tokamak. Both these tests showed that collisions were correctly implemented and that the code is suited for more complex simulations. As a next step, the code was used to calculate the neoclassical radial particle flux by neglecting any turbulent fluctuations in the distribution function and the electric field. Particle fluxes in the neoclassical analytical regimes were simulated for tokamak and stellarator (LHD) configurations. In addition to the comparison with analytical fluxes, a successful benchmark with the DKES code was presented for the tokamak case, which further validates the code for neoclassical simulations. In the final part of the work, the effects of collisions were investigated for slab and toroidal ITGs and TEMs in a tokamak configuration. The results show that collisions reduce the growth rate of slab ITGs in cylinder geometry, whereas they do not affect ITGs in a tokamak, which are mainly curvature-driven. However it is important to note that the pitch-angle scattering operator does not conserve momentum, which is most critical in the parallel direction. Therefore, the damping found in a cylinder could be the consequence of this missing feature and not a physical result [Dimits & Cohen 1994]. Nonetheless, the results are useful to determine whether the instability is mainly being driven by a slab or toroidal ITG mode. EUTERPE also has the feature of including kinetic electrons, which made simulations of TEMs with collisions possible. The combination of collisions and kinetic electrons made the numerical calculations extremely time-consuming, since the time step had to be small enough to resolve the fast electron motion. In contrast to the ITG results, it was observed that collisions are extremely important for TEMs in a tokamak, and in some special cases, depending on whether they were mainly driven by density or temperature gradients, collisions could even suppress the mode (in agreement with [Angioni et al. 2005, Connor et al. 2006]). In the case of stellarators it was found that ITGs are highly dependent on the device configuration. For LHD it was shown that collisions slightly reduce the growth rate of the instability, but for Wendelstein 7-X they do not affect it and the growth rate showed a similar trend with collisionality to that of the tokamak case. Collisions also tend to make the ballooning structure of the modes less pronounced.