No Arabic abstract
Two-fluid Braginskii codes have simulated open-field line turbulence for over a decade, and only recently has it become possible to study these systems with continuum gyrokinetic codes. This work presents a first-of-its-kind comparison between fluid and (long-wavelength) gyrokinetic models in open field-lines, using the GDB and Gkeyll codes to simulate interchange turbulence in the Helimak device at the University of Texas (T. N. Bernard, et. al., Phys. of Plasmas 26, 042301 (2019)). Partial agreement is attained in a number of diagnostic channels when the GDB sources and sheath boundary conditions (BCs) are selected carefully, especially the heat-flux BCs which can drastically alter the temperature. The radial profile of the fluctuation levels is qualitatively similar and quantitatively comparable on the low-field side, although statistics such as moments of the probability density function and the high-frequency spectrum show greater differences. This comparison indicates areas for future improvement in both simulations, such as sheath BCs, as well as improvements in GDB like particle conservation and spatially varying thermal conductivity, in order to achieve better fluid-gyrokinetic agreement and increase fidelity when simulating experiments.
The properties of the boundary plasma in a tokamak are now recognized to play a key role in determining the achievable fusion power and the lifetimes of plasma-facing components. Accurate quantitative modeling and improved qualitative understanding of the boundary plasma ultimately require five-dimensional gyrokinetic turbulence simulations, which have been successful in predicting turbulence and transport in the core. The additional challenges of boundary-plasma simulation necessitate the development of new gyrokinetic codes or major modifications to existing core gyrokinetic codes. In this thesis, we develop the first gyrokinetic continuum code capable of simulating plasma turbulence on open magnetic field lines, which is a key feature of a tokamak scrape-off layer. In contrast to prior attempts at this problem, we use an energy-conserving discontinuous Galerkin discretization in space. To model the interaction between the plasma and the wall, we design conducting-sheath boundary conditions that permit local currents into and out of the wall. We start by designing spatially one-dimensional kinetic models of parallel SOL dynamics and solve these systems using novel continuum algorithms. By generalizing these algorithms to higher dimensions and adding a model for collisions, we present results from the first gyrokinetic continuum simulations of turbulence on two types of open-field-line systems. The first simulation features uniform and straight field lines, such as found in some linear plasma devices. The second simulation is of a hypothetical model we developed of the NSTX scrape-off layer featuring helical field lines. These developments comprise a major step towards a gyrokinetic continuum code for quantitative predictions of turbulence and transport in the boundary plasma of magnetic fusion devices.
Curvature-driven turbulence in a helical open-field-line plasma is investigated using electrostatic five-dimensional gyrokinetic continuum simulations in an all-bad-curvature helical-slab geometry. Parameters for a National Spherical Torus Experiment scrape-off-layer plasma are used in the model. The formation and convective radial transport of plasma blobs is observed, and it is shown that the radial particle-transport levels are several times higher than diffusive Bohm-transport estimates. By reducing the strength of the poloidal magnetic field, the profile of the heat flux to the divertor plate is observed to broaden.
A key uncertainty in the design and development of magnetic confinement fusion energy reactors is predicting edge plasma turbulence. An essential step in overcoming this uncertainty is the validation in accuracy of reduced turbulent transport models. Drift-reduced Braginskii two-fluid theory is one such set of reduced equations that has for decades simulated boundary plasmas in experiment, but significant questions exist regarding its predictive ability. To this end, using a novel physics-informed deep learning framework, we demonstrate the first ever direct quantitative comparisons of turbulent field fluctuations between electrostatic two-fluid theory and electromagnetic gyrokinetic modelling with good overall agreement found in magnetized helical plasmas at low normalized pressure. This framework is readily adaptable to experimental and astrophysical environments, and presents a new technique for the numerical validation and discovery of reduced global plasma turbulence models.
3D2V continuum gyrokinetic simulations of electrostatic plasma turbulence in a straight, open-field-line geometry have been performed using the full-$f$ discontinuous-Galerkin code Gkeyll. These simulations include the basic elements of a fusion-device scrape-off layer: localized sources to model plasma outflow from the core, cross-field turbulent transport, parallel flow along magnetic field lines, and parallel losses at the limiter or divertor with sheath model boundary conditions. The set of sheath boundary conditions used in the model allows currents to flow through the walls. In addition to details of the numerical approach, results from numerical simulations of turbulence in the Large Plasma Device (LAPD), a linear device featuring straight magnetic field lines, are presented.
Differential rotation is known to suppress linear instabilities in fusion plasmas. However, even in the absence of growing eigenmodes, subcritical fluctuations that grow transiently can lead to sustained turbulence. Here transient growth of electrostatic fluctuations driven by the parallel velocity gradient (PVG) and the ion temperature gradient (ITG) in the presence of a perpendicular ExB velocity shear is considered. The maximally simplified case of zero magnetic shear is treated in the framework of a local shearing box. There are no linearly growing eigenmodes, so all excitations are transient. The maximal amplification factor of initial perturbations and the corresponding wavenumbers are calculated as functions of q/epsilon (=safety factor/aspect ratio), temperature gradient and velocity shear. Analytical results are corroborated and supplemented by linear gyrokinetic numerical tests. For sufficiently low values of q/epsilon (<7 in our model), regimes with fully suppressed ion-scale turbulence are possible. For cases when turbulence is not suppressed, an elementary heuristic theory of subcritical PVG turbulence leading to a scaling of the associated ion heat flux with q, epsilon, velocity shear and temperature gradient is proposed; it is argued that the transport is much less stiff than in the ITG regime.