Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA reanalysis of a strong-flow gyrokinetic formalism
96
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We reanalyse an arbitrary-wavelength gyrokinetic formalism [A. M. Dimits, Phys. Plasmas $bf17$, 055901 (2010)], which orders only the vorticity to be small and allows strong, time-varying flows on medium and long wavelengths. We obtain a simpler gyrocentre Lagrangian up to second order. In addition, the gyrokinetic Poisson equation, derived either via variation of the system Lagrangian or explicit density calculation, is consistent with that of the weak-flow gyrokinetic formalism [T. S. Hahm, Phys. Fluids $bf31$, 2670 (1988)] at all wavelengths in the weak flow limit. The reanalysed formalism has been numerically implemented as a particle-in-cell code. An iterative scheme is described which allows for numerical solution of this system of equations, given the implicit dependence of the Euler-Lagrange equations on the time derivative of the potential.
rate research
Read More
In order to predict and analyze turbulent transport in tokamaks, it is important to model transport that arises from microinstabilities. For this task, quasilinear codes have been developed that seek to calculate particle, angular momentum, and heat fluxes both quickly and accurately. In this tutorial, we present a derivation of one such code known as QuaLiKiz, a quasilinear gyrokinetic transport code. The goal of this derivation is to provide a self-contained and complete description of the underlying physics and mathematics of QuaLiKiz from first principles. This work serves both as a comprehensive overview of QuaLiKiz specifically as well as an illustration for deriving quasilinear models in general.
A gyrokinetic model is presented that can properly describe strong flows, large and small amplitude electromagnetic fluctuations occurring on scale lengths ranging from the electron Larmor radius to the equilibrium perpendicular pressure gradient scale length, and large deviations from thermal equilibrium. The formulation of the gyrokinetic model is based on a second order description of the single charged particle dynamics, derived from Lie perturbation theory, where the fast particle gyromotion is decoupled from the slow drifts, assuming that the ratio of the ion sound Larmor radius to the perpendicular equilibrium pressure scale length is small. The collective behavior of the plasma is obtained by a gyrokinetic Boltzmann equation that describes the evolution of the gyroaveraged distribution function and includes a non-linear gyrokinetic Dougherty collision operator. The gyrokinetic model is then developed into a set of coupled fluid equations referred to as the gyrokinetic moment hierarchy. To obtain this hierarchy, the gyroaveraged distribution function is expanded onto a velocity-space Hermite-Laguerre polynomial basis and the gyrokinetic equation is projected onto the same basis, obtaining the spatial and temporal evolution of the Hermite-Laguerre expansion coefficients. The Hermite-Laguerre projection is performed accurately at arbitrary perpendicular wavenumber values. Finally, the self-consistent evolution of the electromagnetic fields is described by a set of gyrokinetic Maxwells equations derived from a variational principle, with the velocity integrals of the gyroaveraged distribution function explicitly evaluated.
Previous limiter-biasing experiments on the Texas Helimak, a simple magnetized torus, have been inconclusive on the effect of flow shear on turbulence levels. To investigate this, the first gyrokinetic simulations of limiter biasing in the Helimak using the plasma physics code Gkeyll have been carried out, and results are presented here. For the scenarios considered, turbulence is mostly driven by the interchange instability, which depends on gradients of equilibrium density profiles. An analysis of both experimental and simulation data demonstrates that shear rates are mostly less than than local linear growth rates, and not all requirements for shear stabilization are met. Rather, the mostly vertical shear flow has an important effect on bulk transport and experimental equilibrium density profiles, and changes to the gradients correspond to changes in turbulence levels.
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.
Magnetic confinement fusion reactors suffer severely from heat and particle losses through turbulent transport, which has inspired the construction of ever larger and more expensive reactors. Numerical simulations are vital to their design and operation, but particle collisions are too infrequent for fluid descriptions to be valid. Instead, strongly magnetised fusion plasmas are described by the gyrokinetic equations, a nonlinear integro-differential system for evolving the particle distribution functions in a five-dimensional position and velocity space, and the consequent electromagnetic field. Due to the high dimensionality, simulations of small reactor sections require hundreds of thousands of CPU hours on High Performance Computing platforms. We develop a Hankel-Hermite spectral representation for velocity space that exploits structural features of the gyrokinetic system. The representation exactly conserves discrete free energy in the absence of explicit dissipation, while our Hermite hypercollision operator captures Landau damping with few variables. Calculation of the electromagnetic fields becomes purely local, eliminating inter-processor communication in, and vastly accelerating, searches for linear instabilities. We implement these ideas in SpectroGK, an efficient parallel code. Turbulent fusion plasmas may dissipate free energy through linear phase mixing to fine scales in velocity space, as in Landau damping, or through a nonlinear cascade to fine scales in physical space, as in hydrodynamic turbulence. Using SpectroGK to study saturated electrostatic drift-kinetic turbulence, we find that the nonlinear cascade suppresses linear phase mixing at energetically-dominant scales, so the turbulence is fluid-like. We use this observation to derive Fourier-Hermite spectra for the electrostatic potential and distribution function, and confirm these spectra with simulations.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
