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.

Transport in Gyrokinetic Tokamaks

Transport in Gyrokinetic Tokamaks
Author:
Publisher:
Total Pages: 45
Release: 1995
Genre:
ISBN:

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.

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.

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.

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.