Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userGyrokinetic continuum simulations of plasma turbulence in the Texas Helimak
58
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The first gyrokinetic simulations of plasma turbulence in the Texas Helimak device, a simple magnetized torus, are presented. The device has features similar to the scrape-off layer region of tokamaks, such as bad-curvature-driven instabilities and sheath boundary conditions on the end plates, which are included in these simulations. Comparisons between simulations and measurements from the experiment show similarities, including equilibrium profiles and fluctuation amplitudes that approach experimental values, but also some important quantitative differences. Both experimental and simulation results exhibit turbulence statistics that are characteristic of blob transport.
rate research
Read More
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.
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.
We analyze plasma heating in weakly collisional kinetic Alfven wave (KAW) turbulence using high resolution gyrokinetic simulations spanning the range of scales between the ion and the electron gyroradii. Real space structures that have a higher than average heating rate are shown not to be confined to current sheets. This novel result is at odds with previous studies, which use the electromagnetic work in the local electron fluid frame, i.e. $mathbf{J} !cdot! (mathbf{E} + mathbf{v}_etimesmathbf{B})$, as a proxy for turbulent dissipation to argue that heating follows the intermittent spatial structure of the electric current. Furthermore, we show that electrons are dominated by parallel heating while the ions prefer the perpendicular heating route. We comment on the implications of the results presented here.
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.
We are developing a new continuum gyrokinetic code, Gkeyll, for use in edge plasma simulations, and here present initial simulations of turbulence on open field lines with model sheath boundary conditions. The code implements an energy conserving discontinuous Galerkin scheme, applicable to a general class of Hamiltonian equations. Several applications to test problems have been done, including a calculation of the parallel heat-flux on divertor plates resulting from an ELM crash in JET, for a 1x/1v SOL scenario explored previously, where the ELM is modeled as a time-dependent intense upstream source. Here we present initial simulations of turbulence on open field lines in the LAPD linear plasma device. We have also done simulations in a helical open-field-line geometry. While various simplifications have been made at present, this still includes some of the key physics of SOL turbulence, such as bad-curvature drive for instabilities and rapid parallel losses with sheath boundary conditions. This is useful for demonstrating the overall feasibility of this approach and for initial physics studies of SOL turbulence. We developed a novel version of DG that uses Maxwellian-weighted basis functions while still preserving exact particle and energy conservation. The Maxwellian-weighted DG method achieves the same error with 4 times less computational cost in 1v, or 16 times lower cost in the 2 velocity dimensions of gyrokinetics (assuming memory bandwidth is the limiting factor).
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
