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.

Incoherent scatter radar (ISR) measurements were used in conjunction with plasma simulations to study two micro-scale plasma processes that commonly occur in the auroral ionosphere. These are 1) ion acoustic turbulence and 2) Langmuir turbulence. Through an ISR experiment we investigated the dependence of ion acoustic turbulence on magnetic aspect angle. The results showed a very strong aspect angle sensitivity which could be utilized to classify the turbulence according to allowable generation mechanisms and sources of free energy. In addition, this work presents results that led to the discovery of a new type of ISR echo, explained as a signature of cavitating Langmuir turbulence. A number of incoherent scatter radar experiments, exploiting a variety of beam and pulse patterns, were designed or revisited to investigate the Langmuir turbulence underlying the radar echoes. The experimental results revealed that Langmuir turbulence is a common feature of the auroral ionosphere. The experimental efforts also led to uncovering a relationship between Langmuir turbulence and one type of natural electromagnetic emission that is sometimes detected on the ground, so-called “medium frequency burst”, providing an explanation for the generation mechanism of these emissions. In an attempt to gain insights into the source mechanism underlying Langmuir turbulence, 1-dimensional Zakharov simulations were employed to study the interactions of ionospheric electron beams with a broad range of parameters with the background plasma at the F region peak. A variety of processes were observed, ranging from a cascade of parametric decays, to formation of stationary wave packets and density cavities in the condensate region, and to direct nucleation and collapse at the initial stage of the turbulence. The simulation results were then compared with the ISR measurements where inconsistencies were found in the spectral details and intensity of the simulated and measured Langmuir turbulence echoes, suggesting the possibility that the direct energy for the turbulence was provided by unstable low-energy (5 − 20 eV) electron populations produced locally in the F region of the ionosphere rather than by electron beams originating from the magnetosphere.

(Hamiltonian chaos and its application to anomalous dynamics in turbulent environment) RMP-induced ELM control has been tested on several tokamaks. It is believed that increase of electron transport across the magnetic field plays an important role. Edge plasma turbulence also affects dynamics in the edge region of tokamak. We study the simultaneous effect of plasma turbulence and RMP-induced stochastic magnetic field within the single-particle framework. We find out that the plasma turbulence is an important element of dynamics that should be taken into account in further (especially single-particle) studies.

La combustion assistée par plasma a reçu une attention croissante dans les deux communautés de plasma et de combustion. Les décharges Nanoseconde Répétitive Pulsée (NRP) sont des techniques prometteuse et efficaces pour initier et contrôler les processus de la combustion en particulier quand les systèmes d’allumage conventionnels sont inefficaces ou trop coûteux en énergie. Néanmoins, les phénomènes rencontrés dans la combustion assistée par plasma sont encore mal connus. Les études numériques présentées dans la littérature sont limitées à des simulations 1-D et 2-D dans des conditions au repos. La complexité du problème augmente dans les configurations pratiques où le phénomène d’allumage est contrôlé par le mouvement du fluide ainsi que le mélange autour de la zone de décharge. La simulation numérique directe (DNS) est un outil de recherche puissant pour la compréhension des interactions plasma/combustion/écoulement. Toutefois, le coût de calcul de la combustion turbulente avec un nombre de Reynolds élevé et la cinétique chimique détaillée couplée avec le plasma hors-équilibre est prohibitif. Cette thèse présente un nouveau modèle de couplageplasma-combustion pour introduire les effets des décharges de plasma hors-équilibre dans le système d’équations qui décrit le phénomène de la combustion. Le modèle est construit en analysant les chemins par lesquels l’énergie électrique est transférée au gaz. Ce modèle de décharges NRP permet des simulations multidimensionalesDNS de la combustion et l’allumage assistés par plasma. Les phénomènes physiques complexes de l’allumage assisté par décharges multiples de plasma dans des mélanges au repos et en régime d’écoulement turbulent sont analysés dans cette thèse. / Plasma-assisted combustion has received increasing attention in both plasma and combustion communities. Nanosecond Repetitively Pulsed (NRP) discharges are a promising and efficient technique to initiate and control combustion processes particularly when conventional ignition systems are rather ineffective or too energy costly. Even though a promising technique, the phenomena occurring in NRP discharges-assisted combustion are still poorly understood. The numerical studies presented in the literature are limited to 1-D and 2-D simulations in quiescent conditions. The problem complexity increases in practical configurations as ignition phenomena are also controlled by the flow and mixing field characteristics in and around the discharge channel. Direct Numerical Simulations (DNS) is a powerful research tool to understand these plasma/combustion/flow interactions. However, the computational cost of fully coupled detailed non-equilibrium plasma and combustion chemistry, and high Reynolds number simulations is prohibitive. This thesis presents a model to describe the effects of non-equilibrium plasma discharges in the set of equations governing the combustion phenomena. Based on the results reported in the literature, the model is constructed by analyzing the channels through which the electric energy is deposited. The two main channels by which the electrons produced during the discharge impact the reactive mixture are considered: 1) the excitation and the subsequent relaxation of the electronic states of nitrogen molecules, which leads to an ultrafast increase of the gas temperature and dissociation of species; and 2) the excitation and relaxation of vibrational states of nitrogen molecules which causes a much slower gas heating. This high level model of NRP discharges allows DNS studies of plasma-assisted combustion / ignition in high turbulent Reynolds number. The complex physics underlying plasma-assisted ignition by multiple discharges in both quiescent and turbulent flow conditions are discussed in the present thesis.

Mémoire (M.Eng.)--Université du Québec à Chicoutimi, 1996. / Document électronique également accessible en format PDF. CaQCU

Title: Anomalous diffusion of plasma in tokamak edge region Author: Jakub Seidl Tutoring institution: Institute of Plasma Physics AS CR, v.v.i. Supervisor: Doc. Ing. Ladislav Krlín, DrSc., IPP AS CR Consultant: RNDr. Radomír Pánek, Ph.D., IPP AS CR Abstract: This work provides brief introduction to the topic of plasma turbu- lence in tokamak edge region and several aspects of anomalous plasma diffusion are discussed. We use numerical code ESEL to model interchange turbulence and investigate properties of turbulent structures in different regimes of parallel trans- port. Means of experimental verification of the results are discussed. Results of the modelling are then used to interpret unexpected results experimentally obtained by electrostatic ball-pen probes on tokamak ASDEX Upgrade, mainly appearance of a 'bump' in power spectra of measured plasma potential. Next, we explain be- haviour of cross-correlation function of density signals measured by two spatially separated Langmuir probes in the vicinity of magnetic separatrix and we point out an ambiguity in interpretation of results of vorticity measurement made by set of floating Langmuir probes. In the last part, transport of plasma impurities by electrostatic turbulent potential is modelled. We identify reversal of radial particle velocity for particles with...

Electromagnetic waves can be used to transmit information over long distances and are therefore often employed for communication purposes. The electromagnetic waves are reflected off material objects on their paths and interact with the medium through which they propagate. For instance, the plasma in the ionosphere can refract and even reflect radio waves propagating through it. By increasing the power of radio waves injected into the ionosphere, the waves start to modify the plasma, resulting in the generation of a wide range of nonlinear processes, including turbulence, in particular near the reflection region. By systematically varying the injected radio waves in terms of frequency, power, polarisation, duty cycle, inclination, etc. the ionosphere can be used as an outdoor laboratory for investigating fundamental properties of the near-Earth space environment as well as of plasma turbulence. In such ionospheric modification experiments, it has been discovered that the irradiation of the ionosphere by powerful radio waves leads to the formation of plasma density structures and to the emission of secondary electromagnetic radiation, a phenomenon known as stimulated electromagnetic emission. These processes are highly repeatable and have enabled systematic investigations of the nonlinear properties of the ionospheric plasma. In this thesis we investigate features of the plasma density structures and the secondary electromagnetic radiation. In a theoretical study we analyse a certain aspect of the formation of the plasma structures. The transient dynamics of the secondary radiation is investigated experimentally in a series of papers, focussing on the initial stage as well as on the decay. In one of the papers we use the transient dynamics of the secondary radiation to reveal the intimate relation between certain features of the radiation and structures of certain scales. Further, we present measurements of unprecedentedly strong secondary radiation, attributed to stimulated Brillouin scattering, and report measurements of the secondary radiation using a novel technique imposed on the transmitted radio waves.

Magnetic nozzles, like Laval nozzles, are observed in several natural systems and have application in areas such as electric propulsion and plasma processing. Plasma flowing through these nozzles is inherently tied to the field lines and must separate for momentum redirection or particle transport to occur. Plasma detachment and associated mechanisms from a magnetic nozzle are investigated. Experimental results are presented from the plume of the VASIMR® VX-200 device flowing along an axisymmetric magnetic nozzle and operated at two ion energies to explore momentum dependent detachment. The argon plume expanded into a 150m3 vacuum chamber where the background pressure was low enough that charge-exchange mean-free-paths were longer than experiment scale lengths. This magnetic nozzle system is demonstrated to hydrodynamically scale up to astrophysical plasmas, particularly the solar chromosphere, implying general relevance to all systems. Plasma parameters were mapped over a large spatial range using measurements from multiple plasma diagnostics. The data show that the plume does not follow the magnetic field lines. A mapped integration of the ion flux shows the plume may be divided into three regions where 1) the plume briefly follows the magnetic flux, 2) diverges quadratically before 3) expanding with linear trajectories. Transitioning from region 1→2, the ion flux departs from the magnetic flux suggesting ion detachment. An instability forms in region 2 driving an oscillating electric field that causes ions to expand before enhancing electron cross-field transport through anomalous resistivity. Transitioning from region 2→3 the electric field dissipates, the trajectories linearize, and the plume effectively detaches. A delineation of sub-to-super Alfvénic flow aligns well with the inflection points of the linearization without a change in magnetic topology. The detachment process is best described as a two part process: First, ions detach by a breakdown of the magnetic moment when the quantity |v/fcLB| becomes of order unity. Second, the turbulent electric field enhances electron transport up to a factor of 4±1 above collisional diffusion; electron cross-field velocities approximate that of the ions and depart on more centralized field lines. Electrons are believed to detach by breakdown of magnetic moment further downstream in the weaker magnetic field.

Le transport anormal de l’energie observé en régime turbulent joue un rôle majeur dans les propriétés de stabilite des plasmas de fusion par confinement magnétique, dans des machines comme ITER. En effet, la turbulence plasma est intimement corrélée au temps de confinement de l’energie, un point clé des recherches en fusion thermonucléaire.<p>Du point de vue théorique, la turbulence plasma est décrite par les équations gyrocinétiques, un ensemble d équations aux dérivées partielles non linéaires couplées. Par suite des très différentes échelles spatiales mises en jeu dans des conditions expérimentales réelles, une simulation numérique directe et complète (DNS) de la turbulence gyrocinétique est totalement hors de portée des plus puissants calculateurs actuels, de sorte que démontrer la faisabilité d’une alternative permettant de réduire l’effort numérique est primordiale. En particulier, les simulations de grandes échelles (”Large-Eddy Simulations” - LES) constituent un candidat pertinent pour permettre une telle r éduction. Les techniques LES ont initialement été développées pour les simulations de fluides turbulents à haut nombre de Reynolds. Dans ces simulations, les plus grandes échelles sont explicitement simulées numériquement, alors que l’influence des plus petites est prise en compte via un modèle implémenté dans le code.<p>Cette thèse présente les premiers développements de techniques LES dans le cadre des équations gyrocinétiques (GyroLES). La modélisation des plus petites échelles est basée sur des bilans d’énergie libre. En effet, l’energie libre joue un rôle important dans la théorie gyrocinétique car elle en est un invariant non lin éaire bien connu. Il est démontré que sa dynamique partage de nombreuses propriétés avec le transfert d’energie dans la turbulence fluide. En particulier, il est montré l’existence d’une cascade d énergie libre, fortement locale et dirigée des grandes échelles vers les petites, dans le plan perpendiculaire â celui du champ magnétique ambiant.<p>La technique GyroLES est aujourd’hui implantée dans le code GENE et a été testée avec succès pour les instabilités de gradient de température ionique (ITG), connues pour jouer un rôle crucial dans la micro-turbulence gyrocinétique. A l’aide des GyroLES, le spectre du flux de chaleur obtenu dans des simulations à très hautes résolutions est correctement reproduit, et ce avec un gain d’un facteur 20 en termes de coût numérique. Pour ces raisons, les simulations gyrocinétiques GyroLES sont potentiellement un excellent candidat pour réduire l’effort numérique des codes gyrocinétiques actuels. <p>/ Anomalous transport due to plasma micro-turbulence is known to play an important role in confinement properties of magnetically confined fusion plasma devices such as ITER. Indeed, plasma turbulence is strongly connected to the energy confinement time, a key issue in thermonuclear fusion research. Plasma turbulence is described by the gyrokinetic equations, a set of nonlinear partial differential equations. Due to the various scales characterizing the turbulent fluctuations in realistic experimental conditions, Direct Numerical Simulations (DNS) of gyrokinetic turbulence remain close to the computational limit of current supercomputers, so that any alternative is welcome to decrease the numerical effort. In particular, Large-Eddy Simulations (LES) are a good candidate for such a decrease. LES techniques have been devised for simulating turbulent fluids at high Reynolds number. In these simulations, the large scales are computed explicitly while the influence of the smallest scales is modeled.<p>In this thesis, we present for the first time the development of the LES for gyrokinetics (GyroLES). The modeling of the smallest scales is based on free energy diagnostics. Indeed, free energy plays an important role in gyrokinetic theory, since it is known to be a nonlinear invariant. It is shown that its dynamics share many properties with the energy transfer in fluid turbulence. In particular, one finds a (strongly) local, forward (from large to small scales) cascade of free energy in the plane perpendicular to the background magnetic field.<p>The GyroLES technique is implemented in the gyrokinetic code Gene and successfully tested for the ion temperature gradient instability (ITG), since ITG is suspected to play a crucial role in gyrokinetic micro-turbulence. Employing GyroLES, the heat flux spectra obtained from highly resolved direct numerical simulations are recovered. It is shown that the gain of GyroLES runs is 20 in terms of computational time. For this reason, Gyrokinetic Large Eddy Simulations can be considered a serious candidate to reduce the numerical cost of gyrokinetic simulations. / Doctorat en Sciences / info:eu-repo/semantics/nonPublished

Turbulence is ubiquitous in the flows of fluids and plasmas. This thesis is devoted to studies of the statistical properties of turbulence in the three-dimensional (3D) Hall magnetohydrodynamic (Hall-MHD) equations, the two-dimensional (2D) MHD equations, the one-dimensional (1D) hyperviscous Burgers equation, and the 3D Navier-Stokes equations. Chapter 1 contains a brief introduction to statistically homogeneous and isotropic turbulence. This is followed by an over-view of the equations we study in the subsequent chapters, the motivation for the studies and a summary of problems we investigate in chapters 2-6. In Chapter 2 we present our study of Hall-MHD turbulence [1]. We show that a shell-model version of the 3D Hall-MHD equations provides a natural theoretical model for investigating the multiscaling behaviors of velocity and magnetic structure functions. We carry out extensive numerical studies of this shell model, obtain the scaling exponents for its structure functions, in both the low-k and high-k power-law ranges of 3D Hall-MHD, and find that the extended-self-similarity procedure is helpful in extracting the multiscaling nature of structure functions in the high-k regime, which otherwise appears to display simple scaling. Our results shed light on intriguing solar-wind measurements. In Chapter 3 we present our study of the inverse-cascade regime in two-dimensional magnetohydrodynamic turbulence [2]. We present a detailed direct numerical simulation (DNS) of statistically steady, homogeneous, isotropic, two-dimensional magnetohydrodynamic (2D MHD) turbulence. Our study concentrates on the inverse cascade of the magnetic vector potential. We examine the dependence of the statistical properties of such turbulence on dissipation and friction coefficients. We extend earlier work significantly by calculating fluid and magnetic spectra, probability distribution functions (PDFs) of the velocity, magnetic, vorticity, current, stream-function, and magnetic-vector-potential fields and their increments. We quantify the deviations of these PDFs from Gaussian ones by computing their flatnesses and hyperflatnesses. We also present PDFs of the Okubo-Weiss parameter, which distinguishes between vortical and extensional flow regions, and its magnetic analog. We show that the hyperflatnesses of PDFs of the increments of the stream-function and the magnetic vector potential exhibit significant scale dependence and we examine the implication of this for the multiscaling of structure functions. We compare our results with those of earlier studies. In Chapter 4 we compare the statistical properties of 2D MHD turbulence for two different energy injection scales. We present systematic DNSs of statistically steady 2D MHD turbulence. Our two DNSs are distinguished by kinj, the wave number at which we inject energy into the system. In our first DNS (run R1), kinj = 2 and, in the second (run R2) kinj = 250. We show that various statistical properties of the turbulent states in the runs R1 and R2 are strikingly different The nature of energy spectrum, probability distribution functions, and topological structures are compared for the two runs R1 and R2 are found to be strikingly different. In Chapter 5 we study the hyperviscous Burgers equation for very high α, order of hyperviscosity [3]. We show, by using direct numerical simulations and theory, how, by increasing α in equations of hydrodynamics, there is a transition from a dissipative to a conservative system. This remarkable result, already conjectured for the asymptotic case α →∞ [U. Frisch et al., Phys. Rev. Lett. 101, 144501 (2008)], is now shown to be true for any large, but finite, value of α greater than a crossover value α crossover. We thus provide a self-consistent picture of how dissipative systems, under certain conditions, start behaving like conservative systems, and hence elucidate the subtle connection between equilibrium statistical mechanics and out-of-equilibrium turbulent flows. In Chapter 6 we show how to use asymptotic-extrapolation and Richardson extrapolation methods to extract the exponents ξ p that characterize the dependence of the order-p moments of the velocity gradients on the Reynolds number Re. To use these extrapolation methods we must have high-precision data for such moments. We obtain these high-precision data by carrying out the most extensive, quadruple precision, pseudospectral DNSs of the Navier-Stokes equation.