Multiscale gyrokinetic turbulence simulations with the real ion-to-electron mass ratio and β value are realized for the first time, where the β value is given by the ratio of plasma pressure to ...magnetic pressure and characterizes electromagnetic effects on microinstabilities. Numerical analysis at both the electron scale and the ion scale is used to reveal the mechanism of their cross-scale interactions. Even with the real-mass scale separation, ion-scale turbulence eliminates electron-scale streamers and dominates heat transport, not only of ions but also of electrons. Suppression of electron-scale turbulence by ion-scale eddies, rather than by long-wavelength zonal flows, is also demonstrated by means of direct measurement of nonlinear mode-to-mode coupling. When the ion-scale modes are stabilized by finite-β effects, the contribution of the electron-scale dynamics to the turbulent transport becomes non-negligible and turns out to enhance ion-scale turbulent transport. Damping of the ion-scale zonal flows by electron-scale turbulence is responsible for the enhancement of ion-scale transport.
It is found that, in the studies of heavy ion transport with gyrokinetic simulations, the ion parallel drift frequency can reverse sign in velocity space when the amplitude variation of the ...electrostatic potential fluctuation is strong along the magnetic field line. As a result, the particle transport related to the parallel dynamics is strongly enhanced. It is noted that, while parallel gradient of the fluctuation amplitude can be instigated by a large magnetic shear or safety factor in a tokamak, the generic mechanism is independent of its cause, which suggests broader applications to kinetic plasma problems. Some relevant topics are briefly addressed in the end.
•A contour dynamics method is applied to the Vlasov-Poisson system.•An efficient implementation of the periodic boundary condition is proposed.•The new method is benchmarked for the linear and ...nonlinear Landau damping.•Energy and particle conservation is examined for the benchmark test.•Particle trapping process is well captured by the contour dynamics method.
We revisit the contour dynamics (CD) simulation method which is applicable to large deformation of distribution function in the Vlasov-Poisson plasma with the periodic boundary, where contours of distribution function are traced without using spatial grids. Novelty of this study lies in application of CD to the one-dimensional Vlasov-Poisson plasma with the periodic boundary condition. A major difficulty in application of the periodic boundary is how to deal with contours when they cross the boundaries. It has been overcome by virtue of a periodic Green's function, which effectively introduces the periodic boundary condition without cutting nor reallocating the contours. The simulation results are confirmed by comparing with an analytical solution for the piece-wise constant distribution function in the linear regime and a linear analysis of the Landau damping. Also, particle trapping by Langmuir wave is successfully reproduced in the nonlinear regime.
A linear eigenmode analysis of the magnetosphere‐ionosphere coupling shows that the feedback instability leading to spontaneous formation of auroral arc structures remains unstable even in a case ...with strong vertical shear of horizontal ion flows induced by ion‐neutral collisions in the E layer. For low‐order Alfvén harmonics, the linear frequency and growth rate obtained by means of a height‐resolved model of the ionosphere are comparable to those resulted from the height‐integrated one, where fine vertical structures of the ionospheric density and magnetic field perturbations are well resolved in the former. The present result confirms validity of the height‐integrated ionosphere model in the long parallel wavelength limit, as assumed in previous studies on the feedback instability.
Plain Language Summary
Spontaneous formation of auroral arcs can be explained in terms of a plasma instability in the magnetosphere‐ionosphere coupling. The present study confirms validity and robustness of the basic instability theory even in cases with inhomogeneous profiles of the ionospheric parameters such as the ion‐neutral collision.
Key Points
Feedback instability in the magnetosphere‐ionosphere coupling is unstable even in a case with strong vertical flow shear in the E layer
The linear eigenmode analysis captures entire profiles of the electromagnetic fields and the density perturbations along a field line
A conventional magnetosphere‐ionosphere coupling model with a height‐integrated ionosphere remains valid in a low‐frequency regime
Gyrokinetic turbulence simulations are applied for the first time to the cross-scale interactions of microtearing modes (MTMs) and electron-temperature-gradient (ETG) modes. The investigation of the ...fluctuation response in a multiscale simulation including both types of instabilities indicates that MTMs are suppressed by ETG turbulence. A detailed analysis of nonlinear mode coupling reveals that radially localized current-sheet structures of MTMs are strongly distorted by fine-scale E×B flows of ETG turbulence. Consequently, electron heat transport caused by the magnetic flutter of MTMs is significantly reduced and ETG turbulence dominates electron heat transport.
Abstract
This article discusses triad transfer analysis via quadratic nonlinearity. To avoid fictitious interactions, symmetrization of the triad transfer is reviewed, including arbitrary orthogonal ...decomposition and coarse graining. The directional representation of the symmetrized triad transfer is proposed by minimizing the number of edges in a network graph of triad interactions with keeping the energy consistency. The directional representation simplifies visualization and allows us to reduce the energy transfer into a one-to-one relation, while eliminating fictitious interactions in non-symmetrized triad transfer functions. Energy transfer processes among plasma turbulent fluctuations that decompose by the singular value decomposition are analyzed as an application. A network graph visualization clearly demonstrates the importance of symmetrization and the consistency between the symmetrized triad transfer and its directional representation.
Plasma β dependence of electromagnetic turbulent transport is investigated by means of gyrokinetic simulations with self-consistent change of the equilibrium magnetic field. It is found that energy ...transport due to ion-temperature-gradient (ITG) driven turbulence does not decrease with increasing β; that is, the ion energy diffusivity does not decrease, and the electron energy diffusivity increases with β. This is because magnetic fluctuations are significantly influenced by the background magnetic field structure change with β by the Pfirsch-Schluter current. The magnetic field change weakens the suppression effect of magnetic perturbations on the growth of the ITG mode, and it also suppresses nonlinear zonal flow production. The influence of the magnetic field change is significant as the global magnetic shear increases.