Gyrokinetic treatment of a grazing angle magnetic field

We develop a gyrokinetic treatment for ions in the magnetic presheath, close to the plasma-wall boundary. We focus on magnetic presheaths with a small magnetic field to wall angle, $\alpha \ll 1$. Characteristic lengths perpendicular to the wall in such a magnetic presheath scale with the typical ion Larmor orbit size, $\rho_{\text{i}}$. The smallest scale length associated with variations parallel to the wall is taken to be across the magnetic field, and ordered $l = \rho_{\text{i}} / \delta$, where $ \delta \ll 1$ is assumed. The scale lengths along the magnetic field line are assumed so long that variations associated with this direction are neglected. These orderings are consistent with what we expect close to the divertor target of a tokamak. We allow for a strong electric field $\vec{E}$ in the direction normal to the electron repelling wall, with strong variation in the same direction. The large change of the electric field over an ion Larmor radius distorts the orbit so that it is not circular. We solve for the lowest order orbits by identifying coordinates, which consist of constants of integration, an adiabatic invariant and a gyrophase, associated with periodic ion motion in the system with $\alpha = \delta = 0$. By using these new coordinates as variables in the limit $\alpha \sim \delta \ll 1$, we obtain a generalized ion gyrokinetic equation. We find another quantity that is conserved to first order and use this to simplify the gyrokinetic equation, solving it in the case of a collisionless magnetic presheath. Assuming a Boltzmann response for the electrons, a form of the quasineutrality equation that exploits the change of variables is derived. The gyrokinetic and quasineutrality equations give the ion distribution function and electrostatic potential in the magnetic presheath if the entrance boundary condition is specified.

Comments: 49 pages including appendices and references

Similar Publications

The anti-Stokes scattering and Stokes scattering in stimulated Brillouin scattering (SBS) cascade have been researched by the Vlasov-Maxwell simulation. In the high-intensity laser-plasmas interaction, the stimulated anti-Stokes Brillouin scattering (SABS) will occur after the second stage SBS rescattering. The mechanism of SABS has been put forward to explain this phenomenon. Read More


In this work, we study the outward propagation of temperature perturbations. For this purpose, we apply an advanced analysis technique, the Transfer Entropy, to ECE measurements performed in ECR heated discharges at the low-shear stellarator TJ-II. We observe that the propagation of these perturbations is not smooth, but is slowed down at specific radial positions, near 'trapping zones' characterized by long time lags with respect to the perturbation origin. Read More


We discuss the dynamics of seeded large amplitude plasma blobs and depletions in an (effective) gravitational field. If the flows are incompressible, the radial center of mass velocity of blobs and depletions is proportional to the square root of their initial cross-field size and amplitude. For compressible flows this scaling holds only for ratios of amplitude to size larger than a critical value. Read More


The nonlinear thin-shell instability (NTSI) may explain some of the turbulent hydrodynamic structures that are observed close to the collision boundary of energetic astrophysical outflows. It develops in nonplanar shells that are bounded on either side by a hydrodynamic shock, provided that the amplitude of the seed oscillations is sufficiently large. The hydrodynamic NTSI has a microscopic counterpart in collisionless plasma. Read More


A scheme with gold cone-capillary is proposed to improve the protons acceleration and involved problems are investigated by using the two-dimensional particle-in-cell simulations. It is demonstrated that the cone-capillary can efficiently guide and collimate the protons to a longer distance and lead to a better beam quality with a dense density $\geq10n_c$, monoenergetic peak energy $E_k \sim 1.51~\mathrm{GeV}$, spatial emittance $\sim0. Read More


A recently proposed technique correlating electric fields and particle velocity distributions is applied to single-point time series extracted from linearly unstable, electrostatic numerical simulations. The form of the correlation, which measures the transfer of phase-space energy density between the electric field and plasma distributions and had previously been applied to damped electrostatic systems, is modified to include the effects of drifting equilibrium distributions of the type that drive counter-streaming and bump-on-tail instabilities. By using single-point time series, the correlation is ideal for diagnosing dynamics in systems where access to integrated quantities, such as energy, is observationally infeasible. Read More


Using two-dimensional hybrid expanding box simulations we study the competition between the continuously driven parallel proton temperature anisotropy and fire hose instabilities in collisionless homogeneous plasmas. For quasi radial ambient magnetic field the expansion drives $T_{\mathrm{p}\|}>T_{\mathrm{p}\perp}$ and the system becomes eventually unstable with respect to the dominant parallel fire hose instability. This instability is generally unable to counteract the induced anisotropization and the system typically becomes unstable with respect to the oblique fire hose instability later on. Read More


From numerical simulations, we show that non-rotating magnetohydrodynamic shear flows are unstable to finite amplitude velocity perturbations and become turbulent, leading to the growth and sustenance of magnetic energy, including large scale fields. This supports the concept that sustained magnetic energy from turbulence is independent of the driving mechanism for large enough magnetic Reynolds numbers. Read More


We present a 2.5-dimensional charge-conservative electromagnetic particle-in-cell (EM-PIC) algorithm optimized for the analysis of vacuum electronic devices (VED) with cylindrical symmetry (axisymmetry). We explore the axisymmetry present in the device geometry, fields, and sources to reduce the dimensionality of the problem from 3D to 2D. Read More


The pilot system development in metre-scale negative laboratory discharges is studied with ns-fast photography. The systems appear as bipolar structures in the vicinity of the negative high-voltage electrode. They appear as a result of a single negative streamer propagation and determine further discharge development. Read More