The effect of small deviations from a Maxwellian equilibrium on turbulent momentum transport in tokamak plasmas is considered. These non-Maxwellian features, arising from diamagnetic effects, introduce a strong dependence of the radial flux of co-cur
rent toroidal angular momentum on collisionality: As the plasma goes from nearly collisionless to weakly collisional, the flux reverses direction from radially inward to outward. This indicates a collisionality-dependent transition from peaked to hollow rotation profiles, consistent with experimental observations of intrinsic rotation.
Scaling laws for the transport and heating of trace heavy ions in low-frequency, magnetized plasma turbulence are derived and compared with direct numerical simulations. The predicted dependences of turbulent fluxes and heating on ion charge and mass
number are found to agree with numerical results for both stationary and differentially rotating plasmas. Heavy ion momentum transport is found to increase with mass, and heavy ions are found to be preferentially heated, implying a mass-dependent ion temperature for very weakly collisional plasmas and for partially-ionized heavy ions in strongly rotating plasmas.
Scaling laws for ion temperature gradient driven turbulence in magnetized toroidal plasmas are derived and compared with direct numerical simulations. Predicted dependences of turbulence fluctuation amplitudes, spatial scales, and resulting heat flux
es on temperature gradient and magnetic field line pitch are found to agree with numerical results in both the driving and inertial ranges. Evidence is provided to support the critical balance conjecture that parallel streaming and nonlinear perpendicular decorrelation times are comparable at all spatial scales, leading to a scaling relationship between parallel and perpendicular spatial scales. This indicates that even strongly magnetized plasma turbulence is intrinsically three-dimensional.
Nonlinear gyrokinetic simulations have been conducted to investigate turbulent transport in tokamak plasmas with rotational shear. At sufficiently large flow shears, linear instabilities are suppressed, but transiently growing modes drive subcritical
turbulence whose amplitude increases with flow shear. This leads to a local minimum in the heat flux, indicating an optimal E x B shear value for plasma confinement. Local maxima in the momentum fluxes are also observed, allowing for the possibility of bifurcations in the E x B shear. The sensitive dependence of heat flux on temperature gradient is relaxed for large flow shear values, with the critical temperature gradient increasing at lower flow shear values. The turbulent Prandtl number is found to be largely independent of temperature and flow gradients, with a value close to unity.
Direct coupling between a transport solver and local, nonlinear gyrokinetic calculations using the multiscale gyrokinetic code TRINITY [M. Barnes, Ph.D. thesis, arxiv:0901.2868] is described. The coupling of the microscopic and macroscopic physics is
done within the framework of multiscale gyrokinetic theory, of which we present the assumptions and key results. An assumption of scale separation in space and time allows for the simulation of turbulence in small regions of the space-time grid, which are embedded in a coarse grid on which the transport equations are implicitly evolved. This leads to a reduction in computational expense of several orders of magnitude, making first-principles simulations of the full fusion device volume over the confinement time feasible on current computing resources. Numerical results from TRINITY simulations are presented and compared with experimental data from JET and ASDEX Upgrade plasmas.
Many plasmas of interest to the astrophysical and fusion communities are weakly collisional. In such plasmas, small scales can develop in the distribution of particle velocities, potentially affecting observable quantities such as turbulent fluxes. C
onsequently, it is necessary to monitor velocity space resolution in gyrokinetic simulations. In this paper, we present a set of computationally efficient diagnostics for measuring velocity space resolution in gyrokinetic simulations and apply them to a range of plasma physics phenomena using the continuum gyrokinetic code GS2. For the cases considered here, it is found that the use of a collisionality at or below experimental values allows for the resolution of plasma dynamics with relatively few velocity space grid points. Additionally, we describe implementation of an adaptive collision frequency which can be used to improve velocity space resolution in the collisionless regime, where results are expected to be independent of collision frequency.
To faithfully simulate ITER and other modern fusion devices, one must resolve electron and ion fluctuation scales in a five-dimensional phase space and time. Simultaneously, one must account for the interaction of this turbulence with the slow evolut
ion of the large-scale plasma profiles. Because of the enormous range of scales involved and the high dimensionality of the problem, resolved first-principles global simulations are very challenging using conventional (brute force) techniques. In this thesis, the problem of resolving turbulence is addressed by developing velocity space resolution diagnostics and an adaptive collisionality that allow for the confident simulation of velocity space dynamics using the approximate minimal necessary dissipation. With regard to the wide range of scales, a new approach has been developed in which turbulence calculations from multiple gyrokinetic flux tube simulations are coupled together using transport equations to obtain self-consistent, steady-state background profiles and corresponding turbulent fluxes and heating. This approach is embodied in a new code, Trinity, which is capable of evolving equilibrium profiles for multiple species, including electromagnetic effects and realistic magnetic geometry, at a fraction of the cost of conventional global simulations. Furthermore, an advanced model physical collision operator for gyrokinetics has been derived and implemented, allowing for the study of collisional turbulent heating, which has not been extensively studied. To demonstrate the utility of the coupled flux tube approach, preliminary results from Trinity simulations of the core of an ITER plasma are presented.
A set of key properties for an ideal dissipation scheme in gyrokinetic simulations is proposed, and implementation of a model collision operator satisfying these properties is described. This operator is based on the exact linearized test-particle co
llision operator, with approximations to the field-particle terms that preserve conservation laws and an H-Theorem. It includes energy diffusion, pitch-angle scattering, and finite Larmor radius effects corresponding to classical (real-space) diffusion. The numerical implementation in the continuum gyrokinetic code GS2 is fully implicit and guarantees exact satisfaction of conservation properties. Numerical results are presented showing that the correct physics is captured over the entire range of collisionalities, from the collisionless to the strongly collisional regimes, without recourse to artificial dissipation.
