Study of mirror effect on scrape-off layer-divertor plasma based on a generalized fluid model incorporating ion temperature anisotropy

2018 ◽  
Vol 58 (6-8) ◽  
pp. 556-562 ◽  
Author(s):  
S. Togo ◽  
T. Takizuka ◽  
D. Reiser ◽  
K. Hoshino ◽  
K. Ibano ◽  
...  
Author(s):  
Menglong Zhao ◽  
Tom Rognlien ◽  
Aaro Einari Jaervinen ◽  
Ilon Joseph

Abstract Understanding tokamak exhaust-power heat loads on divertor plates depends critically on having a realistic model of the scrape-off layer (SOL) plasma. The Braginskii fluid model is often solved to understand the SOL plasma behavior. This model is based on the collisional limit for transport along the magnetic field B⃗. The ions and electron gyrofrequencies are assumed to be much larger than the Coulomb collision frequencies, which are nonetheless, sufficiently large to yield common parallel and perpendicular temperatures for each species, i.e., the temperatures are assumed to be isotropic. In certain circumstances such as encountered for the tokamak H-mode, the ion temperature can be quite anisotropic. In this work, the anisotropy effects are implemented in the 2D transport code UEDGE. Various geometries (1D slab, 2D slab and a toroidal tokamak geometry) are used to study the 2D structure of ion temperature anisotropy and its effects on plasma transport in detail. Results show that the effects of ion temperature anisotropy on the plasma parallel transport are substantial near the magnetic X-point, which leads to different steady state density profiles in the divertor regions. The extra mirror force introduced by ion temperature anisotropy can be one of the main forces contributing to the plasma flow in the SOL.


2014 ◽  
Vol 81 (1) ◽  
Author(s):  
E. E. Scime ◽  
P. A. Keiter ◽  
M. M. Balkey ◽  
J. L. Kline ◽  
X. Sun ◽  
...  

The West Virginia University Hot hELIcon eXperiment (HELIX) provides variable density and ion temperature plasmas, with controllable levels of thermal anisotropy, for space relevant laboratory experiments in the Large Experiment on Instabilities and Anisotropy (LEIA) as well as fundamental studies of helicon source physics in HELIX. Through auxiliary ion heating, the ion temperature anisotropy (T⊥/T∥) is variable from 1 to 20 for parallel plasma beta (β = 8πnkTi∥/B2) values that span the range of 0.0001 to 0.01 in LEIA. The ion velocity distribution function is measured throughout the discharge volume in steady-state and pulsed plasmas with laser induced fluorescence (LIF). The wavelengths of very short wavelength electrostatic fluctuations are measured with a coherent microwave scattering system. Operating at low neutral pressures triggers spontaneous formation of a current-free electric double layer. Ion acceleration through the double layer is detected through LIF. LIF-based velocity space tomography of the accelerated beam provides a two-dimensional mapping of the bulk and beam ion distribution functions. The driving frequency for the m = 1 helical antenna is continuously variable from 8.5 to 16 MHz and frequency dependent variations of the RF coupling to the plasma allow the spontaneously appearing double layers to be turned on and off without modifying the plasma collisionality or magnetic field geometry. Single and multi-species plasmas are created with argon, helium, nitrogen, krypton, and xenon. The noble gas plasmas have steep neutral density gradients, with ionization levels reaching 100% in the core of the plasma source. The large plasma density in the source enables the study of Aflvén waves in the HELIX device.


Author(s):  
Debing Zhang ◽  
Limin Yu ◽  
Erbing Xue ◽  
Xianmei Zhang ◽  
Haijun Ren

Abstract In the nowadays and future fusion devices such as ITER and CFETR, as the use of various heating schemes, the parallel and perpendicular temperature of plasmas can be different; this temperature anisotropy may have significant effects on the turbulence. In this work, the anomalous transport driven by the ion temperature gradient instability is investigated in an anisotropic deuterium-tritium (D-T) plasma. The anisotropic factor $\alpha$, defined as the ratio of perpendicular temperature to parallel temperature, is introduced to describe the temperature anisotropy in the equilibrium distribution function of D. The linear dispersion relation in local kinetic limit is derived, and then numerically evaluated to study the dependence of mode frequency on the anisotropic factor $\alpha$ and the proportion for T particle $\vareT$ by choosing three sets of typical parameters, denoted as the cyclone base case (CBC), ITER and CFETR cases. Based on the linear results, the mixing length model approximation is adopted to analyze the quasi-linear particle and energy fluxes for D and T. It is found that choosing small $\alpha$ and large $\vareT$ is beneficial for the confinement of particle and energy for D and T. This work may be helpful for the estimation of turbulent transport level in the ITER and CFETR devices.


2019 ◽  
Vol 85 (6) ◽  
Author(s):  
Alessandro Geraldini ◽  
F. I. Parra ◽  
F. Militello

The magnetic presheath is a boundary layer occurring when magnetized plasma is in contact with a wall and the angle $\unicode[STIX]{x1D6FC}$ between the wall and the magnetic field $\boldsymbol{B}$ is oblique. Here, we consider the fusion-relevant case of a shallow-angle, $\unicode[STIX]{x1D6FC}\ll 1$ , electron-repelling sheath, with the electron density given by a Boltzmann distribution, valid for $\unicode[STIX]{x1D6FC}/\sqrt{\unicode[STIX]{x1D70F}+1}\gg \sqrt{m_{\text{e}}/m_{\text{i}}}$ , where $m_{\text{e}}$ is the electron mass, $m_{\text{i}}$ is the ion mass, $\unicode[STIX]{x1D70F}=T_{\text{i}}/ZT_{\text{e}}$ , $T_{\text{e}}$ is the electron temperature, $T_{\text{i}}$ is the ion temperature and $Z$ is the ionic charge state. The thickness of the magnetic presheath is of the order of a few ion sound Larmor radii $\unicode[STIX]{x1D70C}_{\text{s}}=\sqrt{m_{\text{i}}(ZT_{\text{e}}+T_{\text{i}})}/ZeB$ , where e is the proton charge and $B=|\boldsymbol{B}|$ is the magnitude of the magnetic field. We study the dependence on $\unicode[STIX]{x1D70F}$ of the electrostatic potential and ion distribution function in the magnetic presheath by using a set of prescribed ion distribution functions at the magnetic presheath entrance, parameterized by $\unicode[STIX]{x1D70F}$ . The kinetic model is shown to be asymptotically equivalent to Chodura’s fluid model at small ion temperature, $\unicode[STIX]{x1D70F}\ll 1$ , for $|\text{ln}\,\unicode[STIX]{x1D6FC}|>3|\text{ln}\,\unicode[STIX]{x1D70F}|\gg 1$ . In this limit, despite the fact that fluid equations give a reasonable approximation to the potential, ion gyro-orbits acquire a spatial extent that occupies a large portion of the magnetic presheath. At large ion temperature, $\unicode[STIX]{x1D70F}\gg 1$ , relevant because $T_{\text{i}}$ is measured to be a few times larger than $T_{\text{e}}$ near divertor targets of fusion devices, ions reach the Debye sheath entrance (and subsequently the wall) at a shallow angle whose size is given by $\sqrt{\unicode[STIX]{x1D6FC}}$ or $1/\sqrt{\unicode[STIX]{x1D70F}}$ , depending on which is largest.


2008 ◽  
Vol 15 (6) ◽  
pp. 957-964 ◽  
Author(s):  
S. M. Finnegan ◽  
M. E. Koepke ◽  
D. J. Knudsen

Abstract. The stationary inertial Alfvén (StIA) wave (Knudsen, 1996) was predicted for cold, collisionless plasma. The model was generalized (Finnegan et al., 2008) to include nonzero values of electron and ion collisional resistivity and thermal pressure. Here, the two-fluid model is further generalized to include anisotropic thermal pressure. A bounded range of values of parallel electron drift velocity is found that excludes periodic stationary Alfvén wave solutions. This exclusion region depends on the value of the local Alfvén speed VA, plasma beta perpendicular to the magnetic field β⊥ and electron temperature anisotropy.


2000 ◽  
Vol 64 (2) ◽  
pp. 131-153 ◽  
Author(s):  
R. N. FRANKLIN ◽  
J. SNELL

This paper reports the results of computations to obtain the spatial distributions of the charged particles in a bounded active plasma dominated by negative ions. Using the fluid model with a constant collision frequency for electrons, positive ions and negative ions the cases of both detachment-dominated gases (such as oxygen) and recombination-dominated gases (such as chlorine) are examined. It is concluded that it is valid to use a Boltzmann relation ne = ne0exp(eV/kT) for the electrons of density ne, where the temperature T is approximately the electron temperature Te, and that the density nn of the negative ions at low pressures obeys nn = nn0exp(eV/kTn), where Tn is the negative-ion temperature. However, at high pressure in detachment-dominated gases where the ratio of negative-ion density to electron density is constant and greater than unity, and when the attachment rate is larger than the ionization rate, the negative ions are distributed with the same effective temperature as the electrons. In all other cases there is no simple relationship. Thus to put nn/ne = const, nn = ne0exp(eV/kTe) and nn = nn0exp(eV/kTn) simultaneously is mathematically inconsistent and physically unsound. Accordingly, expressions deduced for ambipolar diffusion coefficients based on these assumptions have no validity. The correct expressions for the situation where nn/ne = const are obtained without invoking a Boltzmann relation for the negative ions.


2020 ◽  
Author(s):  
Neeraj Jain ◽  
Joerg Buechner

<p>Spacecraft observations show the radial dependence of the solar wind temperature to be slower than what is expected from the adiabatic cooling of the solar wind expanding radially outwards from the sun. The most viable process considered to explain the observed slower-than-adiabatic cooling is the heating of the solar wind plasma by dissipation of the turbulent fluctuations. In solar wind which is  a collisionless plasma in turbulent state, macroscopic energy is cascaded down to kinetic scales where kinetic plasma processes can finally dissipate the energy into heat. The kinetic scale plasma processes responsible  for the dissipation of energy are, however, not well understood. A number of observational and simulation studies have shown that the heating is concentrated in and around current sheets self-consistently formed at kinetic scales. The current sheets contain free energy sources for the growth of plasma instabilities which can serve as the mechanism of the collisionless dissipation. A detailed information on the free energy sources contained in these current sheets of plasma turbulence is lacking but essential to understand the role of  plasma instabilities in collisionless dissipation.</p><p>We carry out 2-D hybrid simulations of kinetic plasma turbulence to study in detail free energy sources available in the current sheets formed in the turbulence. We focus on three free energy sources, namely, plasma density gradient, velocity gradients for both ions and electrons and ion temperature anisotropy. Our simulations show formation of current sheets in which electric current parallel to the externally applied magnetic field flows in a thickness of the order of an ion inertial length. Inside a current sheet, electron flow velocity dominates ion flow velocity in the parallel direction resulting in a larger cross-gradient of the former. The perpendicular electron velocity inside a current sheet also has variations sharper than the corresponding ion velocity. Cross gradients in plasma density are weak (under 10 % variation inside current sheets). Ion temperature is anisotropic in current sheets. Thus the current in the sheets is primarily due to electron shear flow. A theoretical model to explain the difference between electron and ion velocities in current sheets is developed. Spacecraft observations of electron shear flow in space plasma turbulence will be pointed out.   </p><p>These results suggest that the current sheets formed in kinetic plasma turbulence are close to the force free equilibrium rather than the often assumed Harris equilibrium.  This demands investigations of the linear stability properties and nonlinear evolution of force free current sheets with temperature anisotropy. Such studies can provide effective dissipation coefficients to be included in macroscopic model of the solar wind evolution.   </p>


Sign in / Sign up

Export Citation Format

Share Document