No Arabic abstract
Monte Carlo methods are often employed to numerically integrate kinetic equations, such as the particle-in-cell method for the plasma kinetic equation, but these methods suffer from the introduction of counting noise to the solution. We report on a cautionary tale of counting noise modifying the nonlinear saturation of kinetic instabilities driven by unstable beams of plasma. We find a saturated magnetic field in under-resolved particle-in-cell simulations due to the sampling error in the current density. The noise-induced magnetic field is anomalous, as the magnetic field damps away in continuum kinetic and increased particle count particle-in-cell simulations. This modification of the saturated state has implications for a broad array of astrophysical phenomena beyond the simple plasma system considered here, and it stresses the care that must be taken when using particle methods for kinetic equations.
A Hamiltonian two-field gyrofluid model for kinetic Alfven waves (KAWs) in a magnetized electron-proton plasma, retaining ion finite-Larmor-radius corrections and parallel magnetic field fluctuations, is used to study the inverse cascades that develop when turbulence is randomly driven at sub-ion scales. In the directions perpendicular to the ambient field, the dynamics of the cascade turns out to be nonlocal and the ratio $chi_f$ of the wave period to the characteristic nonlinear time at the driving scale affect some of its properties. For example, at small values of $chi_f$, parametric decay instability of the modes driven by the forcing can develop, enhancing for a while inverse transfers. The balanced state, obtained at early time when the two counter-propagating waves are equally driven, also becomes unstable at small $chi_f$, leading to an inverse cascade. For $beta_e$ smaller than a few units, the cascade slows down when reaching the low-dispersion spectral range. For higher $beta_e$, the ratio of the KAW to the Alfven frequencies displays a local minimum. At the corresponding transverse wavenumber, a condensate is formed, and the cascade towards larger scales is then inhibited. Depending on the parameters, a parallel inverse cascade can develop, enhancing the elongation of the ion-scale magnetic vortices that generically form.
Magnetic helicity is robustly conserved in systems with large magnetic Reynolds numbers, including most systems of astrophysical interest. This plays a major role in suppressing the kinematic large scale dynamo and driving the large scale dynamo through the magnetic helicity flux. Numerical simulations of astrophysical systems typically lack sufficient resolution to enforce global magnetic helicity over several dynamical times. Errors in the internal distribution of magnetic helicity are equally serious and possibly larger. Here we propose an algorithm for enforcing strict local conservation of magnetic helicity in the Coulomb gauge in numerical simulations.
In the present paper, we investigate the power-law behaviour of the magnetic field spectra in the Earths magnetosheath region using Cluster spacecraft data under solar minimum condition. The power spectral density of the magnetic field data and spectral slopes at various frequencies are analysed. Propagation angle and compressibility are used to test the nature of turbulent fluctuations. The magnetic field spectra have the spectral slopes between -1.5 to 0 down to spatial scales of 20 ion gyroradius and show clear evidence of a transition to steeper spectra for small scales with a second power-law, having slopes between -2.6 to -1.8. At low frequencies, f_sc<0.3f_ci(where f_ci is ion gyro-frequency), propagation angle approximately 90 degrees to the mean magnetic field, B_0, and compressibility shows a broad distribution, 0.1 < R > 0.9. On the other hand at f_sc>10f_ci, the propagation angle exhibits a broad range between 30-90 degree while R has a small variation: 0.2 < R > 0.5. We conjecture that at high frequencies, the perpendicularly propagating Alfven waves could partly explain the statistical analysis of spectra. To support our prediction of kinetic Alfven wave-dominated spectral slope behaviour at high frequency, we also present a theoretical model and simulate the magnetic field turbulence spectra due to the nonlinear evolution of kinetic Alfven waves. The present study also shows the analogy between the observational and simulated spectra.
We use 3D fully kinetic particle-in-cell simulations to study the occurrence of magnetic reconnection in a simulation of decaying turbulence created by anisotropic counter-propagating low-frequency Alfven waves consistent with critical-balance theory. We observe the formation of small-scale current-density structures such as current filaments and current sheets as well as the formation of magnetic flux ropes as part of the turbulent cascade. The large magnetic structures present in the simulation domain retain the initial anisotropy while the small-scale structures produced by the turbulent cascade are less anisotropic. To quantify the occurrence of reconnection in our simulation domain, we develop a new set of indicators based on intensity thresholds to identify reconnection events in which both ions and electrons are heated and accelerated in 3D particle-in-cell simulations. According to the application of these indicators, we identify the occurrence of reconnection events in the simulation domain and analyse one of these events in detail. The event is related to the reconnection of two flux ropes, and the associated ion and electron exhausts exhibit a complex three-dimensional structure. We study the profiles of plasma and magnetic-field fluctuations recorded along artificial-spacecraft trajectories passing near and through the reconnection region. Our results suggest the presence of particle heating and acceleration related to small-scale reconnection events within magnetic flux ropes produced by the anisotropic Alfvenic turbulent cascade in the solar wind. These events are related to current structures of order a few ion inertial lengths in size.
A familiar problem in space and astrophysical plasmas is to understand how dissipation and heating occurs. These effects are often attributed to the cascade of broadband turbulence which transports energy from large scale reservoirs to small scale kinetic degrees of freedom. When collisions are infrequent, local thermodynamic equilibrium is not established. In this case the final stage of energy conversion becomes more complex than in the fluid case, and both pressure-dilatation and pressure strain interactions (Pi-D $equiv -Pi_{ij} D_{ij}$) become relevant and potentially important. Pi-D in plasma turbulence has been studied so far primarily using simulations. The present study provides a statistical analysis of Pi-D in the Earths magnetosheath using the unique measurement capabilities of the Magnetospheric Multiscale (MMS) mission. We find that the statistics of Pi-D in this naturally occurring plasma environment exhibit strong resemblance to previously established fully kinetic simulations results. The conversion of energy is concentrated in space and occurs near intense current sheets, but not within them. This supports recent suggestions that the chain of energy transfer channels involves regional, rather than pointwise, correlations.