Impurity and Trace Tritium Transport in Tokamak Edge Turbulence

The turbulent transport of impurity or minority species, as for example Tritium, is investigated in drift-Alfv\'en edge turbulence. The full effects of perpendicular and parallel convection are kept for the impurity species. The impurity density develops a granular structure with steep gradients and locally exceeds its initial values due to the compressibility of the flow. An approximate decomposition of the impurity flux into a diffusive part and an effective convective part (characterized by a pinch velocity) is performed and a net inward pinch effect is recovered. The pinch velocity is explained in terms of Turbulent Equipartition and is found to vary poloidally. The results show that impurity transport modeling needs to be two-dimensional, considering besides the radial direction also the strong poloidal variation in the transport coefficients.Comment: 12 Pages, 5 Figure

Cold pulse and rotation reversals with turbulence spreading and residual stress

Transport modeling based on inclusion of turbulence spreading and residual stresses shows internal rotation reversals and polarity reversal of cold pulses, with a clear indication of nonlocal transport effects due to fast spreading in the turbulence intensity field. The effects of turbulence spreading and residual stress are calculated from the gradient of the turbulence intensity. In the model presented in this paper, the flux is carried by the turbulence intensity field, which in itself is subject to radial transport effects. The pulse polarity inversion and the rotation profile reversal positions are close to the radial location of the stable/unstable transition. Both effects have no direct explanation within the framework of classical transport modeling, where the fluxes are related directly to the linear growth rates, the turbulence intensity profile is not considered and the corresponding residual stress is absent. Our simulations are in qualitative agreement with measurements from ohmically heated plasmas. Rotation reversal at a finite radius is found in situations not displaying saturated confinement, which we identify as situations where the plasma is nearly everywhere unstable. As an additional and new effect, the model predicts a perturbation of the velocity profile following a cold pulse from the edge. This allows direct experimental confirmation of both the existence of residual stress caused by turbulence intensity profiles and fundamental ideas of transport modeling presented here. Published by AIP Publishing

Collisional transport across the magnetic field in drift-fluid models

Drift ordered fluid models are widely applied in studies of low-frequency turbulence in the edge and scrape-off layer regions of magnetically confined plasmas. Here, we show how collisional transport across the magnetic field is self-consistently incorporated into drift-fluid models without altering the drift-fluid energy integral. We demonstrate that the inclusion of collisional transport in drift-fluid models gives rise to diffusion of particle density, momentum and pressures in drift-fluid turbulence models and thereby obviate the customary use of artificial diffusion in turbulence simulations. We further derive a computationally efficient, two-dimensional model which can be time integrated for several turbulence de-correlation times using only limited computational resources. The model describes interchange turbulence in a two-dimensional plane perpendicular to the magnetic field located at the outboard midplane of a tokamak. The model domain has two regions modeling open and closed field lines. The model employs a computational expedient model for collisional transport. Numerical simulations show good agreement between the full and the simplified model for collisional transport

Simulation of transition dynamics to high confinement in fusion plasmas

The transition dynamics from the low (L) to the high (H) confinement mode in magnetically confined plasmas is investigated using a first-principles four-field fluid model. Numerical results are in close agreement with measurements from the Experimental Advanced Superconducting Tokamak - EAST. Particularly, the slow transition with an intermediate dithering phase is well reproduced by the numerical solutions. Additionally, the model reproduces the experimentally determined L-H transition power threshold scaling that the ion power threshold increases with increasing particle density. The results hold promise for developing predictive models of the transition, essential for understanding and optimizing future fusion power reactors

Anomalous diffusion, clustering, and pinch of impurities in plasma edge turbulence

The turbulent transport of impurity particles in plasma edge turbulence is investigated. The impurities are modeled as a passive fluid advected by the electric and polarization drifts, while the ambient plasma turbulence is modeled using the two-dimensional Hasegawa--Wakatani paradigm for resistive drift-wave turbulence. The features of the turbulent transport of impurities are investigated by numerical simulations using a novel code that applies semi-Lagrangian pseudospectral schemes. The diffusive character of the turbulent transport of ideal impurities is demonstrated by relative-diffusion analysis of the evolution of impurity puffs. Additional effects appear for inertial impurities as a consequence of compressibility. First, the density of inertial impurities is found to correlate with the vorticity of the electric drift velocity, that is, impurities cluster in vortices of a precise orientation determined by the charge of the impurity particles. Second, a radial pinch scaling linearly with the mass--charge ratio of the impurities is discovered. Theoretical explanation for these observations is obtained by analysis of the model equations.Comment: This article has been submitted to Physics of Plasmas. After it is published, it will be found at http://pop.aip.org/pop

Shear Flow Generation and Energetics in Electromagnetic Turbulence

Zonal flows are recognised to play a crucial role for magnetised plasma confinement. The genesis of these flows out of turbulent fluctuations is therefore of significant interest. We investigate the relative importance of zonal flow generation mechanisms via the Reynolds stress, Maxwell stress, and geodesic acoustic mode (GAM) transfer in drift-Alfv\'en turbulence. By means of numerical computations we quantify the energy transfer into zonal flows owing to each of these effects. The importance of the three driving ingredients in electrostatic and electromagnetic turbulence for conditions relevant to the edge of fusion devices is revealed for a broad range of parameters. The Reynolds stress is found to provide a flow drive, while the electromagnetic Maxwell stress is in the cases considered a sink for the flow energy. In the limit of high plasma beta, where electromagnetic effects and Alfv\'en dynamics are important, the Maxwell stress is found to cancel the Reynolds stress to a high degree. The geodesic oscillations, related to equilibrium pressure profile modifications due to poloidally asymmetric transport, can act as both sinks as drive terms, depending on the parameter regime. For high beta cases the GAMs are the main drive of the flow. This is also reflected in the frequency dependence of the flow, showing a distinct peak at the GAM frequency in that regime.Comment: 16 pages, 12 Figure

Comparison of Edge and Internal Transport Barriers in Drift Wave Predictive Simulations

We have simulated the formation of an internal transport barrier on JET including a self-consistent treatment of ion and electron temperatures and poloidal and toroidal momentum. Similar simulations of edge transport barriers, including the L-H transition have also been made. However, here only polodal momentum and the temperatures were simulated. The internal barrier included an anomalous spinup of poloidal momentum similar to that in the experiment. Also the edge barrier was accompanied by a spinup of poloidal momentum. The experimental density (with no barrier) was used and kept fixed for the internal barrier. For the edge barrier the edge density was varied and it turned out that a lower edge density gave a stronger barrier. Electromagnetic and nonlocal effects were important for both types of barriers