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

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

BOUT++ : Recent and current developments

BOUT++ is a 3D nonlinear finite-difference plasma simulation code, capable of solving quite general systems of PDEs, but targeted particularly on studies of the edge region of tokamak plasmas. BOUT++ is publicly available, and has been adopted by a growing number of researchers worldwide. Here we present improvements which have been made to the code since its original release, both in terms of structure and its capabilities. Some recent applications of these methods are reviewed, and areas of active development are discussed. We also present algorithms and tools which have been developed to enable creation of inputs from analytic expressions and experimental data, and for processing and visualisation of output results. This includes a new tool Hypnotoad for the creation of meshes from experimental equilibria. Algorithms have been implemented in BOUT++ to solve a range of linear algebraic problems encountered in the simulation of reduced MHD and gyro-fluid models: A preconditioning scheme is presented which enables the plasma potential to be calculated efficiently using iterative methods supplied by the PETSc library, without invoking the Boussinesq approximation. Scaling studies are also performed of a linear solver used as part of physics-based preconditioning to accelerate the convergence of implicit time-integration schemes

Hermes : global plasma edge fluid turbulence simulations

The transport of heat and particles in the relatively collisional edge regions of magnetically confined plasmas is a scientifically challenging and technologically important problem. Understanding and predicting this transport requires the self-consistent evolution of plasma fluctuations, global profiles and flows, but the numerical tools capable of doing this in realistic (diverted) geometry are only now being developed. Here a 5-field reduced 2-fluid plasma model for the study of instabilities and turbulence in magnetised plasmas is presented, built on the BOUT++ framework. This cold ion model allows the evolution of global profiles, electric fields and flows on transport timescales, with flux-driven cross-field transport determined self-consistently by electromagnetic turbulence. Developments in the model formulation and numerical implementation are described, and simulations are performed in poloidally limited and diverted tokamak configurations

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

Influence of plasma turbulence on microwave propagation

It is not fully understood how electromagnetic waves propagate through plasma density fluctuations when the size of the fluctuations is comparable with the wavelength of the incident radiation. In this paper, the perturbing effect of a turbulent plasma density layer on a traversing microwave beam is simulated with full-wave simulations. The deterioration of the microwave beam is calculated as a function of the characteristic turbulence structure size, the turbulence amplitude, the depth of the interaction zone and the size of the waist of the incident beam. The maximum scattering is observed for a structure size on the order of half the vacuum wavelength. The scattering and beam broadening was found to increase linearly with the depth of the turbulence layer and quadratically with the fluctuation strength. Consequences for experiments and 3D effects are considered.Comment: 16 pages, 13 figures. This is an author-created, un-copyedited version of an article submitted for publication in Plasma Physics and Controlled Fusion. IoP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from i

Helium exhaust experiments on JET with Type I ELMs in H-mode and with type III ELMs in ITB discharges

An analysis of helium exhaust experiments on JET in the MkII-GB divertor configuration is presented. Helium is pumped by applying an argon frost layer on the divertor cryo pump. Measurement of the helium retention time, tau(He)(*),, is performed in two ways: by the introduction of helium in gas puffs and measurement of the subsequent decay time constant of the helium content, tau(He)(d*); and by helium beam injection and measurement of the helium replacement time, tau(He)(r*). In ELMy H-mode, with plasma configuration optimized for pumping, tau(He)(d*) approximate to 7.2 x tau(E)(th) is achieved, where tau(E)(th) is the thermal energy replacement time. For quasi-steady internal transport barrier (ITB) discharges, the achieved tau(He)(r*) approximate to 4.1 x tau(E)(th) is significantly lower. The achieved helium recycling coefficient, confirmed by an independent measurement to be R-eff approximate to 0.91, is the same in both scenarios. None of the discharges are dominated by core confinement. The difference in tau(He)(*)/tau(E)(th) is instead due to the confinement properties of the edge plasma, which is characterized by Type I ELMs for the H-mode discharges studied, and Type III ELMs for the quasi-steady ITB discharges. This difference is quantified by an independent measurement of the ratio of the helium replacement time with a helium edge source to the energy confinement time