Energy dynamics in a simulation of LAPD turbulence

Energy dynamics calculations in a 3D fluid simulation of drift wave turbulence in the linear Large Plasma Device (LAPD) [W. Gekelman et al., Rev. Sci. Inst. 62, 2875 (1991)] illuminate processes that drive and dissipate the turbulence. These calculations reveal that a nonlinear instability dominates the injection of energy into the turbulence by overtaking the linear drift wave instability that dominates when fluctuations about the equilibrium are small. The nonlinear instability drives flute-like ($k_\parallel = 0$) density fluctuations using free energy from the background density gradient. Through nonlinear axial wavenumber transfer to $k_\parallel \ne 0$ fluctuations, the nonlinear instability accesses the adiabatic response, which provides the requisite energy transfer channel from density to potential fluctuations as well as the phase shift that causes instability. The turbulence characteristics in the simulations agree remarkably well with experiment. When the nonlinear instability is artificially removed from the system through suppressing $k_\parallel=0$ modes, the turbulence develops a coherent frequency spectrum which is inconsistent with experimental data

The effects of non-uniform drive on plasma filaments

Wendelstein 7-X core fueling is primarily achieved through pellet injection. The trajectory of plasmoids from an ablating pellet is an ongoing research question, which is complicated by the complex magnetic geometry of W7-X; curvature drive varies significantly toroidally, including a change in the drift drive direction. Here we use the Hermes model in BOUT++ to simulate cold plasma filaments in slab geometries where the magnetic drift drive is non-uniform along the field line. It is shown that if the field-line-averaged curvature drive is non-zero, a filament will propagate coherently in the direction of average drive. It is also shown that a non-uniform drive will provide a non-uniform propagation; an effect which is reduced at higher temperatures due to an increased sound speed along the field line. Finally, simulations with curvature similar to that found in Wendelstein 7-X are performed which indicate a slow radial propagation of plasma filaments despite a highly non-uniform drive profile

Intrinsic suppression of turbulence in linear plasma devices

Plasma turbulence is the dominant transport mechanism for heat and particles in magnetized plasmas in linear devices and tokamaks, so the study of turbulence is important in limiting and controlling this transport. Linear devices provide an axial magnetic field that serves to confine a plasma in cylindrical geometry as it travels along the magnetic field from the source to the strike point. Due to perpendicular transport, the plasma density and temperature have a roughly Gaussian radial profile with gradients that drive instabilities, such as resistive drift-waves and Kelvin-Helmholtz. If unstable, these instabilities cause perturbations to grow resulting in saturated turbulence, increasing the cross-field transport of heat and particles. When the plasma emerges from the source, there is a time, τk, that describes the lifetime of the plasma based on parallel velocity and length of the device. As the plasma moves down the device, it also moves azimuthally according to E × B and diamagnetic velocities. There is a balance point in these parallel and perpendicular times that sets the stabilisation threshold. We simulate plasmas with a variety of parallel lengths and magnetic fields to vary the parallel and perpendicular lifetimes, respectively, and find that there is a clear correlation between the saturated RMS density perturbation level and the balance between these lifetimes. The threshold of marginal stability is seen to exist where τk ≈ 11τ⊥. This is also associated with the product τkγ∗, where γ∗ is the drift-wave linear growth rate, indicating that the instability must exist for roughly 100 times the growth time for the instability to enter the non-linear growth phase. We explore the root of this correlation and the implications for linear device design

Verification of BOUT++ by the method of manufactured solutions

BOUT++ is a software package designed for solving plasma fluid models. It has been used to simulate a wide range of plasma phenomena ranging from linear stability analysis to 3D plasma turbulence and is capable of simulating a wide range of drift-reduced plasma fluid and gyro-fluid models. A verification exercise has been performed as part of a EUROfusion Enabling Research project, to rigorously test the correctness of the algorithms implemented in BOUT++, by testing order-of-accuracy convergence rates using the Method of Manufactured Solutions (MMS). We present tests of individual components including time-integration and advection schemes, non-orthogonal toroidal field-aligned coordinate systems and the shifted metric procedure which is used to handle highly sheared grids. The flux coordinate independent approach to differencing along magnetic field-lines has been implemented in BOUT++ and is here verified using the MMS in a sheared slab configuration. Finally, we show tests of three complete models: 2-field Hasegawa-Wakatani in 2D slab, 3-field reduced magnetohydrodynamics (MHD) in 3D field-aligned toroidal coordinates, and 5-field reduced MHD in slab geometry

A novel flexible field-aligned coordinate system for tokamak edge plasma simulation

Tokamak plasmas are confined by a magnetic field that limits the particle and heat transport perpendicular to the field. Parallel to the field the ionised particles can move freely, so to obtain confinement the field lines are “closed” (ie.form closed surfaces of constant poloidal flux) in the core of a tokamak. Towards, the edge, however, the field lines intersect physical surfaces, leading to interaction between neutral and ionised particles, and the potential melting of the material surface. Simulation of this interaction is important for predicting the performance and lifetime of future tokamak devices such as ITER. Field-aligned coordinates are commonly used in the simulation of tokamak plasmas due to the geometry and magnetic topology of the system. However, these coordinates are limited in the geometry they allow in the poloidal plane due to orthogonality requirements. A novel 3D coordinate system is proposed herein that relaxes this constraint so that any arbitrary, smoothly varying geometry can be matched in the poloidal plane while maintaining a field-aligned coordinate. This system is implemented in BOUT++ and tested for accuracy using the method of manufactured solutions. A MAST edge cross-section is simulated using a fluid plasma model and the results show expected behaviour for density, temperature, and velocity. Finally, simulations of an isolated divertor leg are conducted with and without neutrals to demonstrate the ion-neutral interaction near the divertor plate and the corresponding beneficial decrease in plasma temperature