Zhu, Hongxuan; Stoltzfus-Dueck, T; Hager, R; Ku, S; Chang, C. S.
Abstract:
Ion orbit loss has been used to model the formation of a strong negative radial electric field Er in the tokamak edge, as well as edge momentum transport and toroidal rotation. To quantitatively measure ion orbit loss, an orbit-flux formulation has been developed and numerically applied to the gyrokinetic particle-in-cell code XGC. We study collisional ion orbit loss in an axisymmetric DIII-D L-mode plasma using gyrokinetic ions and drift-kinetic electrons. Numerical simulations, where the plasma density and temperature profiles are maintained through neutral ionization and heating, show the formation of a quasisteady negative Er in the edge. We have measured a radially outgoing ion gyrocenter flux due to collisional scattering of ions into the loss orbits, which is balanced by the radially incoming ion gyrocenter flux from confined orbits on the collisional time scale. This suggests that collisional ion orbit loss can shift Er in the negative direction compared to that in plasmas without orbit loss. It is also found that collisional ion orbit loss can contribute to a radially outgoing (counter-current) toroidal-angular-momentum flux, which is not balanced by the toroidal-angular-momentum flux carried by ions on the confined orbits. Therefore, the edge toroidal rotation shifts in the co-current direction on the collisional time scale.
Coronal mass ejections (CMEs) are some of the most energetic and violent events in our solar system. The prediction and understanding of CMEs is of particular importance due to the impact that they can have on Earth-based satellite systems, and in extreme cases, ground-based electronics. CMEs often occur when long-lived magnetic flux ropes (MFRs) anchored to the solar surface destabilize and erupt away from the Sun. One potential cause for these eruptions is an ideal magnetohydrodynamic (MHD) instability such as the kink or torus instability. Previous experiments on the Magnetic Reconnection eXperiment (MRX) revealed a class of MFRs that were torus-unstable but kink-stable, which failed to erupt. These “failed-tori” went through a process similar to Taylor relaxation where the toroidal current was redistributed before the eruption ultimately failed. We have investigated this behavior through additional diagnostics that measure the current distribution at the foot points and the energy distribution before and after an event. These measurements indicate that ideal MHD effects are sufficient to explain the energy distribution changes during failed torus events. This excludes Taylor relaxation as a possible mechanism of current redistribution during an event. A new model that only requires non-ideal effects in a thin layer above the electrodes is presented to explain the observed phenomena. This work broadens our understanding of the stability of MFRs and the mechanism behind the failed torus through the improved prediction of the torus instability and through new diagnostics to measure the energy inventory and current profile at the foot points.
Trinczek, Silvia; Parra, Felix I.; Catto, Peter J.; Calvo, Iván; Landreman, Matt
Abstract:
We present a new neoclassical transport model for large aspect ratio tokamaks where the gradient scale lengths are of the size of the ion poloidal gyroradius. Previous work on neoclassical transport across transport barriers assumed large density and potential gradients but a small temperature gradient, or neglected the gradient of the mean parallel flow. Using large aspect ratio and low collisionality expansions, we relax these restrictive assumptions. We define a new set of variables based on conserved quantities, which simplifies the drift kinetic equation whilst keeping strong gradients, and derive equations describing the transport of particles, parallel momentum and energy by ions in the banana regime. The poloidally varying parts of density and electric potential are included. Studying contributions from both passing and trapped particles, we show that the resulting transport is dominated by trapped particles. We find that a non-zero neoclassical particle flux requires parallel momentum input which could be provided through interaction with turbulence or impurities. We derive upper and lower bounds for the energy flux across a transport barrier in both temperature and density and present example profiles and fluxes.
Cheng, C. Z.; Kim, E.-H.; Boardsen, S. A.; Cairns, I. H.; Johnson, J. R.
Abstract:
This includes data for a journal article "O-Mode Instability for Non-thermal Continuum Radiation Generated near the Plasmapause". It includes Matlab data and files shown in the forthcoming Journal article.
The dynamic interplay between the core and the edge plasma has important consequences in the confinement and heating of fusion plasma. The transport of the Scrape-Off-Layer (SOL) plasma imposes boundary conditions on the core plasma, and neutral transport through the SOL influences the core plasma sourcing. In order to better study these effects in a self-consistent, time-dependent fashion with reasonable turn-around time, a reduced model is needed. In this paper we introduce the SOL Box Model, a reduced SOL model that calculates the plasma temperature and density in the SOL given the core-to-edge particle and power fluxes and recycling coefficients. The analytic nature of the Box Model allows one to readily incorporate SOL physics in time-dependent transport solvers for pulse design applications in the control room. Here we demonstrate such a coupling with the core transport solver TRANSP and compare the results with density and temperature measurements, obtained through Thomson scattering and Langmuir probes, of an NSTX discharge. Implications for future interpretive and predictive simulations are discussed.
The data set consists of the figures in a manuscript titled "Self-consistent simulation of compressional Alfvén
eigenmodes excitation and runaway electron diffusion in tokamak disruptions", and plotting script used for figure generation.
