This comment clarifies the relation of the research in a recently published article [Phys. Plasmas 14, 042503 (2007)] to other prior publications addressing the inclusion of electromagnetic and drift-kinetic electron physics in gyrokinetic simulation, raises a concern related to the inclusion of kinetic electrons in a system with magnetic shear, and discusses alternatives in the face of an important limitation on the general applicability of the algorithm described therein.
Nonlinear gyrokinetic simulations have been conducted to investigate turbulent transport in tokamak plasmas with rotational shear. At sufficiently large flow shears, linear instabilities are suppressed, but transiently growing modes drive subcritical turbulence whose amplitude increases with flow shear. This leads to a local minimum in the heat flux, indicating an optimal E x B shear value for plasma confinement. Local maxima in the momentum fluxes are also observed, allowing for the possibility of bifurcations in the E x B shear. The sensitive dependence of heat flux on temperature gradient is relaxed for large flow shear values, with the critical temperature gradient increasing at lower flow shear values. The turbulent Prandtl number is found to be largely independent of temperature and flow gradients, with a value close to unity.
Global electromagnetic gyrokinetic simulations show the existence of near threshold conditions for both a high-$n$ kinetic ballooning mode (KBM) and an intermediate-$n$ kinetic version of peeling-ballooning mode (KPBM) in the edge pedestal of two DIII-D H-mode discharges. When the magnetic shear is reduced in a narrow region of steep pressure gradient, the KPBM is significantly stabilized, while the KBM is weakly destabilized and hence becomes the most-unstable mode. Collisions decrease the KBMs critical $beta$ and increase the growth rate.
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 report on first evidence of wave activity during neutral beam heating in KSTAR plasmas: 40 kHz magnetic fluctuations with a toroidal mode number of n=1. Our analysis suggests this a beta-induced Alfven eigenmode resonant with the q=1 surface. A kinetic analysis, when coupled with electron temperature measurements from electron cyclotron emission and ion/electron temperature ratios from crystallography, enables calculation of the frequency evolution, which is in agreement with observations. Complementary detailed MHD modelling of the magnetic configuration and wave modes supports the BAE mode conclusion, by locating an n=1 mode separated from the continuum in the core region. Finally, we have computed the threshold to marginal stability for a range of ion temperature profiles. These suggest the BAE can be driven unstable by energetic ions when the ion temperature radial gradient is sufficiently large. Our findings suggest that mode existence could be used as a form of inference for temperature profile consistency in the radial interval of the mode, thereby extending the tools of MHD spectroscopy.
The Hall term has often been neglected in MHD codes as it is difficult to compute. Nevertheless setting it aside for numerical reasons led to ignoring it altogether. This is especially problematic when dealing with tokamak physics as the Hall term cannot be neglected as this paper shows.
Bruce I. Cohen
,Andris M. Dimits
,
.
(2007)
.
"Comment on Electromagnetic Global Gyrokinetic Simulation of Shear Alfven Wave Dynamics in Tokamak Plasmas [Phys. Plasmas 14, 042503 (2007)]"
.
Bruce Cohen I
هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا