Gyrokinetic Simulation of Global Turbulent Transport Properties in Tokamak Experiments

Gyrokinetic Simulation of Global Turbulent Transport Properties in Tokamak Experiments
Author:
Publisher:
Total Pages: 92505
Release: 2006
Genre:
ISBN:

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.

Comprehensive Gyrokinetic Simulation of Tokamak Turbulence at Finite Relative Gyroradius

Comprehensive Gyrokinetic Simulation of Tokamak Turbulence at Finite Relative Gyroradius
Author:
Publisher:
Total Pages: 8
Release: 2002
Genre:
ISBN:

OAK B202 COMPREHENSIVE GYROKINETIC SIMULATION OF TOKAMAK TURBULENCE AT FINITE RELATIVE GYRORADIUS. A continuum global gyrokinetic code GYRO has been developed to comprehensively simulate turbulent transport in actual experimental profiles and allow 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, and all in real tokamak geometry. Most importantly the code operates at finite relative gyroradius ([rho]*) so as to treat the profile shear stabilization effects which break gyroBohm scaling. The code operates in a cyclic flux tube limit which allows only gyroBohm scaling and a noncyclic radial annulus with physical profile variation. The later requires an adaptive source to maintain equilibrium profiles. Simple ITG simulations demonstrate the broken gyroBohm scaling depends on the actual rotational velocity shear rates competing with mode growth rates, direct comprehensive simulations of the DIII-D [rho]*-scaled L-mode experiments are presented as a quantitative test of gyrokinetics and the paradigm.

ADVANCES IN COMPREHENSIVE GYROKINETIC SIMULATIONS OF TRANSPORT IN TOKAMAKS.

ADVANCES IN COMPREHENSIVE GYROKINETIC SIMULATIONS OF TRANSPORT IN TOKAMAKS.
Author: R. E. WALTZ
Publisher:
Total Pages:
Release: 2004
Genre:
ISBN:

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.

Advances in Comprehensive Gyrokinetic Simulations of Transport in Tokamaks

Advances in Comprehensive Gyrokinetic Simulations of Transport in Tokamaks
Author:
Publisher:
Total Pages:
Release: 2004
Genre:
ISBN:

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.

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

Magnetic Fluctuations in Gyrokinetic Simulations of Tokamak Scrape-Off Layer Turbulence
Author: Noah Roth Mandell
Publisher:
Total Pages: 0
Release: 2021
Genre:
ISBN:

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.

Kinetic Simulation of Edge Instability in Fusion Plasmas

Kinetic Simulation of Edge Instability in Fusion Plasmas
Author: Daniel Patrick Fulton
Publisher:
Total Pages: 103
Release: 2015
Genre:
ISBN: 9781321995824

In this work, gyrokinetic simulations in edge plasmas of both tokamaks and field reversed configurations (FRC) have been carried out using the Gyrokinetic Toroidal Code (GTC) and A New Code (ANC) has been formulated for cross-separatrix FRC simulation. In the tokamak edge, turbulent transport in the pedestal of an H-mode DIII-D plasma is studied via simulations of electrostatic driftwaves. Annulus geometry is used and simulations focus on two radial locations corresponding to the pedestal top with mild pressure gradient and steep pressure gradient. A reactive trapped electron instability with typical ballooning mode structure is excited in the pedestal top. At the steep gradient, the electrostatic instability exhibits unusual mode structure, peaking at poloidal angles theta=+- pi/2. Simulations find this unusual mode structure is due to steep pressure gradients in the pedestal but not due to the particular DIII-D magnetic geometry. Realistic DIII-D geometry has a stabilizing effect compared to a simple circular tokamak geometry. Driftwave instability in FRC is studied for the first time using gyrokinetic simulation. GTC is upgraded to treat realistic equilibrium calculated by an MHD equilibrium code. Electrostatic local simulations in outer closed flux surfaces find ion-scale modes are stable due to the large ion gyroradius and that electron drift-interchange modes are excited by electron temperature gradient and bad magnetic curvature. In the scrape-off layer (SOL) ion-scale modes are excited by density gradient and bad curvature. Collisions have weak effects on instabilities both in the core and SOL. Simulation results are consistent with density fluctuation measurements in the C-2 experiment using Doppler backscattering (DBS). The critical density gradients measured by the DBS qualitatively agree with the linear instability threshold calculated by GTC simulations. One outstanding critical issue in the FRC is the interplay between turbulence in the FRC core and SOL regions. While the magnetic flux coordinates used by GTC provide a number of computational advantages, they present unique challenges at the magnetic field separatrix. To address this limitation, a new code, capable of coupled core-SOL simulations, is formulated, implemented, and successfully verified.

Advances in Quasilinear Gyrokinetic Modeling of Turbulent Transport

Advances in Quasilinear Gyrokinetic Modeling of Turbulent Transport
Author: Cole Darin Stephens
Publisher:
Total Pages: 312
Release: 2021
Genre:
ISBN:

The quest to harness fusion energy requires the successful modeling of plasma turbulence and transport in magnetic confinement devices. For such modeling, the requisite length and time scales span many orders of magnitude. Integrated modeling approaches are constructed to account for the wide range of physics involved in turbulent transport by coupling separate physical models together. The primary physical models used in this work are kinetic and designed to simulate microturbulence on the smallest scales associated with turbulent transport. However, high precision nonlinear kinetic simulations often cannot be easily coupled to integrated modeling suites due to the extreme computational costs that would be involved. Model reduction which drastically reduces the computational complexity of the problem is therefore necessary. One must of course ensure that the reduced model does not severely diminish the accuracy of the calculation; the model reduction itself must be founded on more exact computational approaches as well as fundamental theoretical principles. One of the most successful approaches in model reduction is quasilinear gyrokinetics. There are two fundamental assumptions for the quasilinear model examined in this work. First, the three adiabatic invariants (the magnetic moment, the longitudinal invariant, and the poloidal flux) must be appropriately conserved and their associated single charged particle motions (the gyromotion, the bounce-transit motion, and the toroidal drift motion) must be characterized accurately. Second, the quasilinear approximation must hold such that the coherent linear response is adequate enough to compute the quasilinear fluxes without full calculation of the nonlinear physics. The particular model used, QuaLiKiz, has been proven successful in reproducing local gyrokinetic fluxes in the tokamak core while remaining computationally tractable. There are three primary goals of this dissertation project. The first is to examine the fundamental physics underlying gyrokinetic and reduced model approaches at the single charged particle scale. To achieve this goal, we examine the assumption of magnetic moment invariance in a wide variety of electromagnetic fields. We successfully identify the dimensionless parameters that determine magnetic moment conservation in each scenario and then proceed to quantify the degree to which magnetic moment conservation is broken. In doing so, we confirm that the magnetic moment is sufficiently conserved for a wide range of regimes relevant to tokamak plasmas. In addition, we derive new analytic formulas for quantities associated with bounce-transit motion in circular tokamak fields. We compare these new, more exact calculations to approximations commonly used in reduced models (including QuaLiKiz) and determine the conditions such that the approximations break down. We then also confirm that the approximations are valid in the tokamak core for conventional, large aspect ratio devices. The second goal of this dissertation project is to rederive and compile the model equations for QuaLiKiz from first principles. Over the years of QuaLiKiz's development, there has never been a complete manuscript that sketches the derivation of QuaLiKiz from start to finish. The lack of such a document makes it difficult to extend the physics of QuaLiKiz to new parameter regimes of interest. Various possible extensions such as including electromagnetic effects or more realistic tokamak geometries require the adjustment of several different assumptions that would affect the derivation in key ways. As such, correct implementations of new physics would require an existing derivation as a reference point lest the implementation be handled in an incoherent fashion. In addition, a step-by-step outline of how each assumption of QuaLiKiz affects the derivation can be helpful in determining which assumptions can be relaxed for a more accurate model. The successful completion of this derivation, included in this dissertation, will be immensely useful for future QuaLiKiz improvement and validation. With the derivation in hand, we proceed to the third goal of this project: improving the collisional model of QuaLiKiz. Collisions play an essential role in characterizing the transport associated with trapped electron modes. It has become evident in recent studies that the collisional model in QuaLiKiz requires improvement; in integrated modeling, the imprecise treatment of collisional trapped electron modes leads to incorrect density profile predictions near the tokamak core for highly collisional regimes. We revisit the collision model implemented in QuaLiKiz and use the more exact gyrokinetic code GENE (Gyrokinetic Electromagnetic Numerical Experiment) to make improvements to QuaLiKiz's collision operator. We then use the new version of QuaLiKiz in integrated modeling to compare density profiles predicted by the old and new collision operators. We confirm that the new collision operator leads to density profiles that more accurately match the experimental profiles.