Gyrokinetic Simulations of Turbulent Impurity Transport in Tokamaks


Book Description

Understanding impurity transport in the core of tokamak plasmas is central to achieving controlled fusion. Indeed impurities are ubiquitous in these devices and their presence in the core are detrimental to plasma confinement (fuel dilution, Bremsstrahlung). Recently, specific attention was given to the convective mechanism related to the gradient of the toroidal rotation to explain experimental flat/hollow impurity profiles in the plasma core. In this thesis, up-to-date modelling tools (NEO for neoclassical transport and GKW for turbulent transport) including the impact of toroidal rotation are used to study both the neoclassical and turbulent contributions to impurity fluxes. A comparison of the experimental and modelled carbon density peaking factor (R/LnC) is performed for a large number of baseline and hybrid H-mode plasmas (increased confinement regimes) with modest to high toroidal rotation from the European tokamak JET. Confrontation of experimental and modelled carbon peaking factor yields two main results. First roto-diffusion is found to have a nonnegligible impact on the carbon peaking factor at high values of the toroidal rotation frequency gradient. Second, there is a tendency to overpredict the experimental R/LnC in the core inner region where the carbon density profiles are hollow. This disagreement between experimental and modelled R/LnC, closely related to the collisionality, is also observed for the momentum transport channel which hints at a common parallel symmetry breaking mechanism lacking in the simulations.




Impurity Transport in Tokamak Plasmas


Book Description

Impurity transport is an issue of utmost importance for tokamaks. Indeed high-Z materials are only partially ionized in the plasma core, so that they can lead to prohibitive radiative losses even at low concentrations, and impact dramatically plasma performance and stability. On-axis accumulation of tungsten has been widely observed in tokamaks.While the very core impurity peaking is generally attributed to neoclassical effects, turbulent transport could well dominate in the gradient region at ITER relevant collisionality. Up to recently, first principles simulations of corresponding fluxes were performed with different dedicated codes, implicitly assuming that both transport channels are separable and therefore additive. The validity of this assumption is questionned. Simulations obtained with the gyrokinetic code GYSELA have shown clear evidences of a neoclassical-turbulence synergy for impurity transport and allowed the identification of a mechanism that underly this synergy.An analytical work allows to compute the level and the structure of the axisymmetric part of the electric potential knowing the turbulence intensity. Two mechanisms are found for the generation of poloidal asymmetries of the electric potential: flow compressibility and the ballooning of the turbulence. A new prediction for the neoclassical impurity flux in presence of large poloidal asymmetries and pressure anisotropies has been derived. A fair agreement has been found between the new theoretical prediction for neoclassical impurity flux and the results of a GYSELA simulation displaying large poloidal asymmetries and pressure anisotropies induced by the presence of turbulence.




Advances in Comprehensive Gyrokinetic Simulations of Transport in Tokamaks


Book Description

A continuum global gyrokinetic code GYRO has been developed to comprehensively simulate core turbulent transport in actual experimental profiles and enable direct quantitative comparisons to the experimental transport flows. GYRO not only treats the now standard ion temperature gradient (ITG) mode turbulence, but also treats trapped and passing electrons with collisions and finite [beta], equilibrium ExB shear stabilization, and all in real tokamak geometry. Most importantly the code operates at finite relative gyroradius ([rho]{sub *}) so as to treat the profile shear stabilization and nonlocal effects which can break gyroBohm scaling. The code operates in either a cyclic flux-tube limit (which allows only gyroBohm scaling) or a globally with physical profile variation. Rohm scaling of DIII-D L-mode has been simulated with power flows matching experiment within error bars on the ion temperature gradient. Mechanisms for broken gyroBohm scaling, neoclassical ion flows embedded in turbulence, turbulent dynamos and profile corrugations, plasma pinches and impurity flow, and simulations at fixed flow rather than fixed gradient are illustrated and discussed.







ADVANCES IN COMPREHENSIVE GYROKINETIC SIMULATIONS OF TRANSPORT IN TOKAMAKS.


Book Description

A continuum global gyrokinetic code GYRO has been developed to comprehensively simulate core turbulent transport in actual experimental profiles and enable direct quantitative comparisons to the experimental transport flows. GYRO not only treats the now standard ion temperature gradient (ITG) mode turbulence, but also treats trapped and passing electrons with collisions and finite {beta}, equilibrium ExB shear stabilization, and all in real tokamak geometry. Most importantly the code operates at finite relative gyroradius ({rho}{sub *}) so as to treat the profile shear stabilization and nonlocal effects which can break gyroBohm scaling. The code operates in either a cyclic flux-tube limit (which allows only gyroBohm scaling) or globally with physical profile variation. Bohm scaling of DIII-D L-mode has been simulated with power flows matching experiment within error bars on the ion temperature gradient. Mechanisms for broken gyroBohm scaling, neoclassical ion flows embedded in turbulence, turbulent dynamos and profile corrugations, are illustrated.










Transport in Gyrokinetic Tokamaks


Book Description

A comprehensive study of transport in full-volume gyrokinetic (gk) simulations of ion temperature gradient driven turbulence in core tokamak plasmas is presented. Though this g̀̀yrokinetic tokamak ̀̀is much simpler than experimental tokamaks, such simplicity is an asset, because a dependable nonlinear transport theory for such systems should be more attainable. Toward this end, we pursue two related lines of inquiry. (1) We study the scalings of gk tokamaks with respect to important system parameters. In contrast to real machines, the scalings of larger gk systems (a/?{sub s} ≳ 64) with minor radius, with current, and with a/?{sub s} are roughly consistent with the approximate theoretical expectations for electrostatic turbulent transport which exist as yet. Smaller systems manifest quite different scalings, which aids in interpreting differing mass-scaling results in other work. (2) With the goal of developing a first-principles theory of gk transport, we use the gk data to infer the underlying transport physics. The data indicate that, of the many modes k present in the simulation, only a modest number (N{sub k} ∼ 10) of k dominate the transport, and for each, only a handful (N{sub p} ∼ 5) of couplings to other modes p appear to be significant, implying that the essential transport physics may be described by a far simpler system than would have been expected on the basis of earlier nonlinear theory alone. Part of this analysis is the inference of the coupling coefficients M{sub kpq} governing the nonlinear mode interactions, whose measurement from tokamak simulation data is presented here for the first time.




Gyrokinetic Simulation of Global Turbulent Transport Properties in Tokamak Experiments


Book Description

A general geometry gyro-kinetic model for particle simulation of plasma turbulence in tokamak experiments is described. It incorporates the comprehensive influence of noncircular cross section, realistic plasma profiles, plasma rotation, neoclassical (equilibrium) electric fields, and Coulomb collisions. An interesting result of global turbulence development in a shaped tokamak plasma is presented with regard to nonlinear turbulence spreading into the linearly stable region. The mutual interaction between turbulence and zonal flows in collisionless plasmas is studied with a focus on identifying possible nonlinear saturation mechanisms for zonal flows. A bursting temporal behavior with a period longer than the geodesic acoustic oscillation period is observed even in a collisionless system. Our simulation results suggest that the zonal flows can drive turbulence. However, this process is too weak to be an effective zonal flow saturation mechanism.




Magnetic Fluctuations in Gyrokinetic Simulations of Tokamak Scrape-Off Layer Turbulence


Book Description

Understanding turbulent transport physics in the tokamak edge and scrape-off layer (SOL) is critical to developing a successful fusion reactor. The dynamics in these regions plays a key role in achieving high fusion performance by determining the edge pedestal that suppresses turbulence in the high-confinement mode (H-mode). Additionally, the survivability of a reactor is set by the heat load to the vessel walls, making it important to understand turbulent spreading of heat as it flows along open magnetic field lines in the SOL. Large-amplitude fluctuations, magnetic X-point geometry, and plasma interactions with material walls make simulating turbulence in the edge/SOL more challenging than in the core region, necessitating specialized gyrokinetic codes. Further, the inclusion of electromagnetic effects in gyrokinetic simulations that can handle the unique challenges of the boundary plasma is critical to the understanding of phenomena such as the pedestal and edge-localized modes, for which electromagnetic dynamics are expected to be important.In this thesis, we develop the first capability to simulate electromagnetic gyrokinetic turbulence on open magnetic field lines. This is an important step towards comprehensive electromagnetic gyrokinetic simulations of the coupled edge/SOL system. By using a continuum full-f approach via an energy-conserving discontinuous Galerkin (DG) discretization scheme that avoids the Ampere cancellation problem, we show that electromagnetic fluctuations can be handled in a robust, stable, and efficient manner in the gyrokinetic module of the Gkeyll code. We then present results which roughly model the scrape-off layer of the National Spherical Torus Experiment (NSTX), and show that electromagnetic effects can affect blob dynamics and transport. We also formulate the gyrokinetic system in field-aligned coordinates for modeling realistic edge and scrape-off layer geometries in experiments. A novel DG algorithm for maintaining positivity of the distribution function while preserving conservation laws is also presented.