No Arabic abstract
We present a method for studying the evolution of plasma turbulence by tracking dispersion relations in the energy spectrum in the wavenumber-frequency domain. We apply hybrid plasma simulations in a simplified two-dimensional geometry to demonstrate our method and its applicability to plasma turbulence in the ion kinetic regime. We identify four dispersion relations: ion-Bernstein waves, oblique whistler waves, oblique Alfven/ion-cyclotron waves, and a zero-frequency mode. The energy partition and frequency broadening are evaluated for these modes. The method allows us to determine the evolution of decaying plasma turbulence in our restricted geometry and shows that it cascades along the dispersion relations during the early phase with an increasing broadening around the dispersion relations.
We explore the multi-faceted important features of turbulence (e.g., anisotropy, dispersion, diffusion) in the three-dimensional (3D) wavenumber domain ($k_parallel$, $k_{perp,1}$, $k_{perp,2}$), by employing the k-filtering technique to the high-quality measurements of fields and particles from the MMS multi-spacecraft constellation. We compute the 3D power spectral densities (PSDs) of magnetic and electric fluctuations (marked as $rm{PSD}(delta mathbf{B}(mathbf{k}))$ and $rm{PSD}(delta mathbf{E}_{langlemathbf{v}_mathrm{i}rangle}(mathbf{k}))$), both of which show a prominent spectral anisotropy in the sub-ion range. We give the first 3D image of the bifurcation between power spectra of the electric and magnetic fluctuations, by calculating the ratio between $rm{PSD}(delta mathbf{E}_{ langlemathbf{v}_mathrm{i}rangle}(mathbf{k}))$ and $rm{PSD}(delta mathbf{B}(mathbf{k}))$, the distribution of which is related to the non-linear dispersion relation. We also compute the ratio between electric spectra in different reference frames defined by the ion bulk velocity, that is $mathrm{PSD}(delta{mathbf{E}_{mathrm{local} mathbf{v}_mathrm{i}}})/mathrm{PSD}(delta{mathbf{E}_{ langlemathbf{v}_mathrm{i}rangle}})$, to visualize the turbulence ion diffusion region (T-IDR) in wavenumber space. The T-IDR has an anisotropy and a preferential direction of wavevectors, which is generally consistent with the plasma wave theory prediction based on the dominance of kinetic Alfven waves (KAW). This work manifests the worth of the k-filtering technique in diagnosing turbulence comprehensively, especially when the electric field is involved.
We present the first study of the formation and dissipation of current sheets at electron scales in a wave-driven, weakly collisional, 3D kinetic turbulence simulation. We investigate the relative importance of dissipation associated with collisionless damping via resonant wave-particle interactions versus dissipation in small-scale current sheets in weakly collisional plasma turbulence. Current sheets form self-consistently from the wave-driven turbulence, and their filling fraction is well correlated to the electron heating rate. However, the weakly collisional nature of the simulation necessarily implies that the current sheets are not significantly dissipated via Ohmic dissipation. Rather, collisionless damping via the Landau resonance with the electrons is sufficient to account for the measured heating as a function of scale in the simulation, without the need for significant Ohmic dissipation. This finding suggests the possibility that the dissipation of the current sheets is governed by resonant wave-particle interactions and that the locations of current sheets correspond spatially to regions of enhanced heating.
A scaling theory of long-wavelength electrostatic turbulence in a magnetised, weakly collisional plasma (e.g., ITG turbulence) is proposed, with account taken both of the nonlinear advection of the perturbed particle distribution by fluctuating ExB flows and of its phase mixing, which is caused by the streaming of the particles along the mean magnetic field and, in a linear problem, would lead to Landau damping. It is found that it is possible to construct a consistent theory in which very little free energy leaks into high velocity moments of the distribution function, rendering the turbulent cascade in the energetically relevant part of the wave-number space essentially fluid-like. The velocity-space spectra of free energy expressed in terms of Hermite-moment orders are steep power laws and so the free-energy content of the phase space does not diverge at infinitesimal collisionality (while it does for a linear problem); collisional heating due to long-wavelength perturbations vanishes in this limit (also in contrast with the linear problem, in which it occurs at the finite rate equal to the Landau-damping rate). The ability of the free energy to stay in the low velocity moments of the distribution function is facilitated by the anti-phase-mixing effect, whose presence in the nonlinear system is due to the stochastic version of the plasma echo (the advecting velocity couples the phase-mixing and anti-phase-mixing perturbations). The partitioning of the wave-number space between the (energetically dominant) region where this is the case and the region where linear phase mixing wins its competition with nonlinear advection is governed by the critical balance between linear and nonlinear timescales (which for high Hermite moments splits into two thresholds, one demarcating the wave-number region where phase mixing predominates, the other where plasma echo does).
Kinetic-range turbulence in magnetized plasmas and, in particular, in the context of solar-wind turbulence has been extensively investigated over the past decades via numerical simulations. Among others, one of the widely adopted reduced plasma model is the so-called hybrid-kinetic model, where the ions are fully kinetic and the electrons are treated as a neutralizing (inertial or massless) fluid. Within the same model, different numerical methods and/or approaches to turbulence development have been employed. In the present work, we present a comparison between two-dimensional hybrid-kinetic simulations of plasma turbulence obtained with two complementary approaches spanning about two decades in wavenumber - from MHD inertial range to scales well below the ion gyroradius - with a state-of-the-art accuracy. One approach employs hybrid particle-in-cell (HPIC) simulations of freely-decaying Alfvenic turbulence, whereas the other consists of Eulerian hybrid Vlasov-Maxwell (HVM) simulations of turbulence continuously driven with partially-compressible large-scale fluctuations. Despite the completely different initialization and injection/drive at large scales, the same properties of turbulent fluctuations at $k_perprho_igtrsim1$ are observed. The system indeed self-consistently reprocesses the turbulent fluctuations while they are cascading towards smaller and smaller scales, in a way which actually depends on the plasma beta parameter. Small-scale turbulence has been found to be mainly populated by kinetic Alfven wave (KAW) fluctuations for $betageq1$, whereas KAW fluctuations are only sub-dominant for low-$beta$.
Understanding the nature of the turbulent fluctuations below the ion gyroradius in solar-wind turbulence is a great challenge. Recent studies have been mostly in favor of kinetic Alfven wave (KAW) type of fluctuations, but other kinds of fluctuations with characteristics typical of magnetosonic, whistler and ion Bernstein modes, could also play a role depending on the plasma parameters. Here we investigate the properties of the sub-proton-scale cascade with high-resolution hybrid-kinetic simulations of freely-decaying turbulence in 3D3V phase space, including electron inertia effects. Two proton plasma beta are explored: the intermediate $beta_p=1$ and low $beta_p=0.2$ regimes, both typically observed in solar wind and corona. The magnetic energy spectum exhibits $k_perp^{-8/3}$ and $k_|^{-7/2}$ power laws at $beta_p=1$, while they are slightly steeper at $beta_p=0.2$. Nevertheless, both regimes develop a spectral anisotropy consistent with $k_|sim k_perp^{2/3}$ at $k_perprho_p>1$, and pronounced small-scale intermittency. In this context, we find that the kinetic-scale cascade is dominated by KAW-like fluctuations at $beta_p=1$, whereas the low-$beta$ case presents a more complex scenario suggesting the simultaneous presence of different types of fluctuations. In both regimes, however, a non-negligible role of ion Bernstein type of fluctuations at the smallest scales seems to emerge.