The RF-Transpond code couples a fluid plasma transport solver with a frequency domain cold plasma RF wave solver in a 1D domain parallel to a strong background magnetic field. A ponderomotive force ...term proportional to parallel gradients in the electric field strength is included in the transport model in order to describe ponderomotive effects in the scrape-off layer (SOL) of fusion plasmas. The transport and wave codes are verified independently and a coupled case corresponding to experimental parameters from the LArge Plasma Device (LAPD) is presented. The density perturbation ratio Rn, calculated to describe ponderomotive force driven modifications, is up to 20% for the simulation inputs used.
Program Title: rf-transpond
CPC Library link to program files:https://doi.org/10.17632/xn3y2yx9wj.1
Developer's repository link:https://github.com/rhealbarnett/rf-transpond.git
Licensing provisions: MIT
Programming language: Matlab
Nature of problem: Self consistent coupled model describing ponderomotive force driven density modification in the near field of RF antennas.
Solution method: The density and velocity solutions are calculated from the continuity and momentum transport equations, solved using a finite difference time domain method, which include a ponderomotive force term that depends on the radio frequency electric field. The frequency domain electric field solution is calculated from the cold plasma wave equation, solved using a finite difference frequency domain method, where the cold plasma dielectric tensor is a function of the density. The electric field and density couple the two models, providing self consistency.
In-depth experimental characterisation of spontaneous shear flow patterning into a so-called E×B staircase-named after its planetary analogue-is shown in magnetised plasma turbulence, using ...ultrafast-sweeping reflectometry in the Tore Supra tokamak. Staircase signatures are found in a large variety of L-mode plasma conditions. Sensitivity to the dominant source of free energy is highlighted for the first time. A connection between staircase shear layer permeability and deviation from gyro-Bohm confinement scaling is strongly suggested, opening new routes to understanding confinement in drift-wave turbulence.
Gyrokinetic tokamak plasmas can exhibit intrinsic toroidal rotation driven by the residual stress. While most studies have attributed the residual stress to the parallel-momentum flux from the ...turbulent E×B motion, the parallel-momentum flux from the drift-orbit motion (denoted Π_{∥}^{D}) and the E×B-momentum flux from the E×B motion (denoted Π_{E×B}) are often neglected. Here, we use the global total-f gyrokinetic code XGC to study the residual stress in the core and the edge of a DIII-D H-mode plasma. Numerical results show that both Π_{∥}^{D} and Π_{E×B} make up a significant portion of the residual stress. In particular, Π_{∥}^{D} in the core is higher than the collisional neoclassical level in the presence of turbulence, while in the edge it represents an outflux of countercurrent momentum even without turbulence. Using a recently developed "orbit-flux" formulation, we show that the higher-than-neoclassical-level Π_{∥}^{D} in the core is driven by turbulence, while the outflux of countercurrent momentum from the edge is mainly due to collisional ion orbit loss. These results suggest that Π_{∥}^{D} and Π_{E×B} can be important for the study of intrinsic toroidal rotation.Gyrokinetic tokamak plasmas can exhibit intrinsic toroidal rotation driven by the residual stress. While most studies have attributed the residual stress to the parallel-momentum flux from the turbulent E×B motion, the parallel-momentum flux from the drift-orbit motion (denoted Π_{∥}^{D}) and the E×B-momentum flux from the E×B motion (denoted Π_{E×B}) are often neglected. Here, we use the global total-f gyrokinetic code XGC to study the residual stress in the core and the edge of a DIII-D H-mode plasma. Numerical results show that both Π_{∥}^{D} and Π_{E×B} make up a significant portion of the residual stress. In particular, Π_{∥}^{D} in the core is higher than the collisional neoclassical level in the presence of turbulence, while in the edge it represents an outflux of countercurrent momentum even without turbulence. Using a recently developed "orbit-flux" formulation, we show that the higher-than-neoclassical-level Π_{∥}^{D} in the core is driven by turbulence, while the outflux of countercurrent momentum from the edge is mainly due to collisional ion orbit loss. These results suggest that Π_{∥}^{D} and Π_{E×B} can be important for the study of intrinsic toroidal rotation.
This work presents a code to simulate collisions and laser-plasma interactions of multiple-species plasmas by solving 1D3V Vlasov-Fokker-Planck-Maxwell equations. The distribution functions of all ...species are expanded into spherical harmonic series of arbitrary order. Each electron or ion species possesses its distinct and suitable velocity mesh and the number of harmonics. Collisions between different ion species are considered and treated with the Fokker-Planck collision term. The code employs the total variation diminishing (TVD) scheme to avoid numerical instabilities in the region of steep density, rarefied plasma-vacuum interface, and non-periodic boundary conditions, which makes it possible to simulate cases that are close to actual physical scenes. The number of particles is globally conserved for each species, the energy for each species is conserved for collisions between the same species, and total energy is conserved for collisions between different species.
Abstract
Doubly peaked density distribution is expected not only to affect the plasma-wetted area at divertor plates, but also to correlate with the upstream density profile and hence characteristics ...of magnetohydrodynamic activities in tokamak plasmas (Wang
et al
2020
Phys. Rev. Lett.
124
195002). Clarifying its origin is important to understand the compatibility between power/particle exhausts in divertor and high-performance core plasmas required by present-day and future tokamak devices. In this paper, we analyze the double-peak density profile appearing in the modeling during the physics design phase of the new lower tungsten divertor for EAST by using a comprehensive 2D SOLPS-ITER code package, including full drifts and currents, with a concentration on an unfavorable magnetic field (ion
B
× ∇
B
drift is directed away from the primary
X
-point). The results indicate that
E
×
B
drift induced by the plasma potential gradient near the target, which is closely related to the divertor state, plays essential roles in the formation of a double-peak profile at the target: (1) large enough radial
E
p
×
B
drift produces a broadened high-density region; (2) strong poloidal
E
r
×
B
drift drives a significant particle sink and creates a valley on the high-density profile. Thus, the simulation results can explain why this kind of doubly peaked density profile is usually observed at the high-recycling divertor regime. In addition, features of the double-peak ion saturation current distribution measured in preliminary experiments testing the new lower tungsten divertor are qualitatively consistent with the simulations.
Neoclassical transport processes are important to the understanding of plasma confinement physics in doubly periodic magnetized toroidal plasmas, especially, after the impact of the momentum ...confinement on the particle and energy confinement is recognized. Real doubly periodic tori in general are non-axisymmetric, with symmetric tori as a special case. An eight-moment approach to transport theory with plasma density N, plasma pressure p, mass flow velocity V and heat flow q as independent variables is adopted. Transport processes are dictated by the solutions of the momentum and heat flux balance equations. For toroidal plasma confinement devices, the first order (in the gyro-radius ordering) plasma flows are on the magnetic surface to guarantee good plasma confinement and are thus two-dimensional. Two linearly independent components of the momentum equation are required to determine the flows completely. Once this two-dimensional flow is relaxed, i.e. the momentum equation reaches a steady state, plasmas become ambipolar, and all the transport fluxes are determined through the flux-force relation. The flux-force relation is derived both from the kinetic definitions for the transport fluxes and from the manipulation of the momentum and heat flux balance equations to illustrate the nature of the transport fluxes by examining their corresponding driven forces and their roles in the momentum and heat flux balance equations. Steady-state plasma flows are determined by the components of the stress and heat stress tensors in the momentum and heat flux balance equations. This approach emphasizes the pivotal role of the momentum equation in the transport processes and is particularly useful in modelling plasma flows in experiments. The methodology for neoclassical transport theory is applied to fluctuation-driven transport fluxes in the quasilinear theory to unify these two theories. Experimental observations in tokamaks and stellarators for the physics discussed are presented.
Edge localized modes (ELMs) have a detrimental effect on the plasma facing components and pose one of the most serious obstacles for steady-state operation in a future fusion device. For future ...fusion machines, the control or even full suppression of ELMs is mandatory. In the past years, extensive effort has been directed to the development of operational regimes that maintain the high confinement and good performance of the H-mode, while at the same time ELMs are suppressed or mitigated. Several natural ELM-free and small-ELM regimes, such as the quiescent H-mode, the improved energy confinement mode, the type-II and the grassy ELM-regime, have been obtained in various tokamaks. The state-of-the-art and recent advances of these ELM-free and small-ELM scenarios are reviewed, and the access and sustainment as well as their applicability to ITER are discussed.
A bottom-up Einstein-Maxwell-dilaton holographic model is used to compute, for the first time, the behavior of several transport coefficients of the hot and baryon-rich strongly coupled quark-gluon ...plasma at the critical point and also across the first-order phase transition line in the phase diagram. The observables under study are the shear and bulk viscosities, the baryon and thermal conductivities, the baryon diffusion, the jet quenching parameter $\hat{q}$, as well as the heavy-quark drag force and the Langevin diffusion coefficients. These calculations provide a phenomenologically promising estimate for these coefficients, given that our model quantitatively reproduces lattice QCD thermodynamics results, both at zero and finite baryon density, besides naturally incorporating the nearly perfect fluidity of the quark-gluon plasma. Here we find that the diffusion of baryon charge, and also the shear and bulk viscosities, are suppressed with increasing baryon density, indicating that the medium becomes even closer to perfect fluidity at large densities. On the other hand, the jet quenching parameter and the heavy-quark momentum diffusion are enhanced with increasing density. The observables display a discontinuity gap when crossing the first-order phase transition line, while developing an infinite slope at the critical point. The transition temperatures associated with different transport coefficients differ in the crossover region but are found to converge at the critical point.
We observe plasma heating due to collisional diffusion across a separatrix when a magnesium ion column in a Penning-Malmberg trap is cyclically pushed back and forth across a partial trapping ...barrier. The barrier is an externally applied axisymmetric “squeeze” potential, which creates a velocity separatrix between trapped and passing particles. Weak ion-ion collisions then cause separatrix crossings, leading to irreversible heating. Furthermore, the heating rate scales as the square root of the oscillation rate times the collision frequency and thus can be dominant for low-collisionality plasmas. The particle velocity distribution function is measured with coherent laser induced fluorescence and shows passing and trapped particles having an out-of-phase response to the forced plasma oscillations.