The presence of a large pinch velocity in the edge pedestal of high confinement (H-mode) tokamak plasmas implies that particle transport in the plasma edge must be treated by a pinch-diffusion theory, rather than a pure diffusion theory. Momentum balance also requires the inclusion of a pinch term in descriptions of edge particle transport. A numerical investigation of solving generalized pinch-diffusion theory using methods extended from the numerical solution methodology of pure diffusion theory has been carried out. The generalized diffusion equation has been numerically integrated using the central finite-difference approximation for the diffusion term and three finite difference approximations of the pinch term, and then solved using Gauss reduction. The pinch-diffusion relation for the radial particle flux was solved directly and used as a benchmark for the finite-difference algorithm solutions to the generalized diffusion equation. Both equations are solved using several mesh spacings, and it is found that a finer mesh spacing will be required in the edge pedestal, where the inward pinch velocity is large in H-mode plasmas, than is necessary for similar accuracy further inward where the pinch velocity diminishes. An expression for the numerical error of various finite-differencing algorithms is presented.

The transport of heat out of tokamak plasmas by turbulence is the dominant mechanism limiting the performance of fusion reactors. Turbulence can be driven by the ion temperature gradient (ITG) and suppressed by toroidal equilibrium scale sheared flows. Numerical simulations attempting to understand, and ultimately reduce, turbulence are crucial for guiding the design and optimisation of future reactors. In this thesis, we investigate ion-scale turbulence by means of local gyrokinetic simulations in the outer core of the Mega Ampere Spherical Tokamak (MAST).We perform a parameter scan in the values of the ITG and the flow shear. We show that nonlinear simulations reproduce the experimental ion heat flux and that the experimentally measured values of the ITG and the flow shear lie close to the turbulence threshold. We demonstrate that the system is subcritical in the presence of flow shear, i.e., the system is formally stable to small perturbations, but transitions to a turbulent state given a large enough initial perturbation. We propose a scenario for the transition to subcritical turbulence previously unreported in tokamak plasmas: close to the threshold, the plasma is dominated by a low number of coherent long-lived structures; as the system is taken away from the threshold into the more unstable regime, the number of these structures increases until they fill the domain and a more conventional turbulence emerges. We make quantitative comparisons of correlation properties between our simulations and experimental measurements of ion-scale density fluctuations from the MAST BES diagnostic. We apply a synthetic diagnostic to our simulation data and find reasonable agreement of the correlation properties of the simulated and experimental turbulence, most notably of the correlation time, for which significant discrepancies were found in previous numerical studies of MAST turbulence. We show that the properties of turbulence are essentially functions of the distance to threshold, as quantified by the ion heat flux. We find that turbulence close to the threshold is strongly affected by flow shear, whereas far from threshold, the turbulence resembles a conventional ITG-driven, zonal-flow damped regime.

The future fusion power plants that are based on magnetic confinement will deal with plasmas that inevitably contain energetic (non-thermal) particles. These particles come, for instance, from fusion reactions or from external heating of the plasma. Ensembles of energetic ions can excite eigenmodes in the Alfvén frequency range to such an extent that the resulting wave fields redistribute the energetic ions, and potentially eject them from the plasma. The redistribution of ions may cause a substantial reduction of heating efficiency. Understanding the dynamics of such instabilities is necessary to optimise the operation of fusion experiments and of future fusion power plants. Two models have been developed to simulate the interaction between energetic ions and Alfvén eigenmodes. One is a bump-on-tail model, of which two versions have been developed: one fully nonlinear and one quasilinear. The quasilinear version has a lower dimensionality of particle phase space than the nonlinear one. Unlike previous similar studies, the bump-on-tail model contains a decorrelation of the wave-particle phase in order to model stochasticity of the system. When the characteristic time scale for macroscopic phase decorrelation is similar to or shorter than the time scale of nonlinear wave-particle dynamics, the nonlinear and the quasilinear descriptions quantitatively agree. A finite phase decorrelation changes the growth rate and the saturation amplitude of the wave mode in systems with an inverted energy distribution around the wave-particle resonance. Analytical expressions for the correction of the growth rate and the saturation amplitude have been derived, which agree well with numerical simulations. A relatively weak phase decorrelation also diminishes frequency chirping events of the eigenmode. The second model is called FOXTAIL, and it has a wider regime of validity than the bump-on-tail model. FOXTAIL is able to simulate systems with multiple eigenmodes, and it includes effects of different individual particle orbits relative to the wave fields. Simulations with FOXTAIL and the nonlinear bump-on-tail model have been compared in order to determine the regimes of validity of the bump-on-tail model quantitatively. Studies of two-mode scenarios confirmed the expected consequences of a fulfillment of the Chirikov criterion for resonance overlap. The influence of ICRH on the eigenmode-energetic ion system has also been studied, showing qualitatively similar effects as seen by the presence of phase decorrelation. Another model, describing the efficiency of fast wave current drive, has been developed in order to study the influence of passive components close to the antenna, in which currents can be induced by the antenna generated wave field. It was found that the directivity of the launched wave, averaged over model parameters, was lowered by the presence of passive components in general, except for low values of the single pass damping of the wave, where the directivity was slightly increased, but reversed in the toroidal direction. / De framtida fusionskraftverken baserade på magnetisk inneslutning kommer att hantera plasmor som oundvikligen innehåller energetiska (icke-termiska) partiklar. Dessa partiklar kommer exempelvis från fusionsreaktioner eller från externa uppvärmningsmekanismer av plasmat. Ensembler av energetiska joner kan excitera egenmoder i Alfvén-frekvensområdet i en sådan utsträckning att de resulterande vågfälten omfördelar de energetiska jonerna i rummet, och potentiellt slungar ut jonerna ur plasmat. Omfördelningen av joner kan orsaka en väsentligen minskad uppvärmningseffekt. Det är nödvändigt att förstå dynamiken hos denna typ av instabilitet för att kunna optimera verkningsgraden hos experiment och hos framtida fusionskraftverk. Två modeller har utvecklats för att simulera interaktionen mellan energetiska joner och Alfvén-egenmoder. Den första är en bump-on-tail-modell, av vilken två versioner har utvecklats: en fullt icke-linjär och en kvasi-linjär. I den kvasi-linjära versionen har partiklarnas fasrum en lägre dimensionalitet än i den icke-linjära versionen. Till skillnad från tidigare liknande studier innehåller denna bump-on-tail-modell en dekorrelation av våg-partikelfasen för att modellera stokasticitet hos systemet. När den karakteristiska tidsskalan för makroskopisk fasdekorrelation är ungefär samma som eller kortare än tidsskalan för icke-linjär våg-partikeldynamik så stämmer den icke-linjära och den kvasi-linjära beskrivningen överens kvantitativt. En ändlig fasdekorrelation förändrar vågmodens tillväxthastighet och satureringsamplitud i system med en inverterad energifördelning omkring våg-partikelresonansen. Analytiska uttryck för korrektionen av tillväxthastigheten och satureringsamplituden har härletts, vilka stämmer väl överens med numeriska simuleringar. En relativt svag fasdekorrelation försvagar även "frequency chirping events" (snabba frekvensskiftningar i korttids-Fourier-transformen av egenmodens amplitudutveckling) hos egenmoden. Den andra modellen, kallad FOXTAIL, har ett mycket bredare giltighetsområde än bump-on-tail-modellen. FOXTAIL kan simulera system med flera egenmoder, och den inkluderar effekter av olika enskilda partikelbanor relativt vågfälten. Simuleringar med FOXTAIL och med bump-on-tail-modellen har jämförts för att kvantitativt bestämma bump-on-tail-modellens giltighetsområde. Studier av scenarier med två egenmoder bekräftar de förväntade effekterna av när Chirikov-kriteriet för resonansöverlapp uppfylls. Även inflytandet av ICRH på dynamiken mellan egenmoder och energetiska joner har studerats, vilket har visat kvalitativt liknande effekter som har observerats i närvaron av fasdekorrelation. En annan modell, vilken beskriver effektiviteten hos "fast wave current drive" (strömdrivning med snabba magnetosoniska vågor), har utvecklats för att studera inflytandet av passiva komponenter nära antennen, i vilka strömmar kan induceras av vågfälten som genereras av antennen. Det visades att den utskickade vågens direktivitet, medelvärdesbildat över modellparametrar, generellt sett minskade vid närvaron av passiva komponenter, förutom vid låg "sinlge pass damping" (dämpning av vågen vid propagering genom hela plasmat), då direktiviteten istället ökade något, men bytte tecken i toroidal riktning. / <p>QC 20160927</p>

Experiments for the development of fusion power that are based on magnetic confinement deal with plasmas that inevitably contain energetic (non-thermal) particles. These particles come e.g. from fusion reactions or from external heating of the plasma. Ensembles of energetic ions can excite plasma waves in the Alfvén frequency range to such an extent that the resulting wave fields redistribute the energetic ions, and potentially eject them from the plasma. The redistribution of ions may cause a substantial reduction heating efficiency, and it may damage the inner walls and other components of the vessel. Understanding the dynamics of such instabilities is necessary to optimise the operation of fusion experiments and of future fusion power plants. A Monte Carlo model that describes the nonlinear wave-particle dynamics in a toroidal plasma has been developed to study the excitation of the abovementioned instabilities. A decorrelation of the wave-particle phase is added in order to model stochasticity of the system (e.g. due to collisions between particles). Based on the nonlinear description with added phase decorrelation, a quasilinear version of the model has been developed, where the phase decorrelation has been replaced by a quasilinear diffusion coefficient in particle energy. When the characteristic time scale for macroscopic phase decorrelation becomes similar to or shorter than the time scales of nonlinear wave-particle dynamics, the two descriptions quantitatively agree on a macroscopic level. The quasilinear model is typically less computationally demanding than the nonlinear model, since it has a lower dimensionality of phase space. In the presented studies, several effects on the macroscopic wave-particle dynamics by the presence of phase decorrelation have been theoretically and numerically analysed, e.g. effects on the growth and saturation of the wave amplitude, and on the so called frequency chirping events with associated hole-clump pair formation in particle phase space. Several effects coming from structures of the energy distribution of particles around the wave-particle resonance has also been studied. / <p>QC 20150330</p>