Magnetized turbulence is ubiquitous in many astrophysical and terrestrial plasmas but no universal theory exists. Even the detailed energy dynamics in magnetohydrodynamic (MHD) turbulence are still not well understood. We present a suite of subsonic, super-Alfvénic, high plasma beta MHD turbulence simulations that only vary in their dynamical range, i.e., in their separation between the large-scale forcing and dissipation scales, and their dissipation mechanism (implicit large eddy simulation, ILES, and direct numerical simulation (DNS)). Using an energy transfer analysis framework we calculate the effective numerical viscosities and resistivities, and demonstrate that all ILES calculations of MHD turbulence are resolved and correspond to an equivalent visco-resistive MHD turbulence calculation. Increasing the number of grid points used in an ILES corresponds to lowering the dissipation coefficients, i.e., larger (kinetic and magnetic) Reynolds numbers for a constant forcing scale. Independently, we use this same framework to demonstrate that—contrary to hydrodynamic turbulence—the cross-scale energy fluxes are not constant in MHD turbulence. This applies both to different mediators (such as cascade processes or magnetic tension) for a given dynamical range as well as to a dependence on the dynamical range itself, which determines the physical properties of the flow. We do not observe anymore »
Despite decades of study of high-temperature weakly collisional plasmas, a complete understanding of how energy is transferred between particles and fields in turbulent plasmas remains elusive. Two major questions in this regard are how fluid-scale energy transfer rates, associated with turbulence, connect with kinetic-scale dissipation, and what controls the fraction of dissipation on different charged species. Although the rate of cascade has long been recognized as a limiting factor in the heating rate at kinetic scales, there has not been direct evidence correlating the heating rate with MHD-scale cascade rates. Using kinetic simulations and in situ spacecraft data, we show that the fluid-scale energy flux indeed accounts for the total energy dissipated at kinetic scales. A phenomenology, based on disruption of proton gyromotion by fluctuating electric fields that are produced in turbulence at proton scales, argues that the proton versus electron heating is controlled by the ratio of the nonlinear timescale to the proton cyclotron time and by the plasma beta. The proposed scalings are supported by the simulations and observations.
- Publication Date:
- NSF-PAR ID:
- 10386251
- Journal Name:
- The Astrophysical Journal
- Volume:
- 941
- Issue:
- 2
- Page Range or eLocation-ID:
- Article No. 137
- ISSN:
- 0004-637X
- Publisher:
- DOI PREFIX: 10.3847
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract -
ABSTRACT The physical foundations of the dissipation of energy and the associated heating in weakly collisional plasmas are poorly understood. Here, we compare and contrast several measures that have been used to characterize energy dissipation and kinetic-scale conversion in plasmas by means of a suite of kinetic numerical simulations describing both magnetic reconnection and decaying plasma turbulence. We adopt three different numerical codes that can also include interparticle collisions: the fully kinetic particle-in-cell vpic, the fully kinetic continuum Gkeyll, and the Eulerian Hybrid Vlasov–Maxwell (HVM) code. We differentiate between (i) four energy-based parameters, whose definition is related to energy transfer in a fluid description of a plasma, and (ii) four distribution function-based parameters, requiring knowledge of the particle velocity distribution function. There is an overall agreement between the dissipation measures obtained in the PIC and continuum reconnection simulations, with slight differences due to the presence/absence of secondary islands in the two simulations. There are also many qualitative similarities between the signatures in the reconnection simulations and the self-consistent current sheets that form in turbulence, although the latter exhibits significant variations compared to the reconnection results. All the parameters confirm that dissipation occurs close to regions of intense magnetic stresses, thusmore »
-
Abstract Magnetic fields grow quickly, even at early cosmological times, suggesting the action of a small-scale dynamo (SSD) in the interstellar medium (ISM) of galaxies. Many studies have focused on idealized, isotropic, homogeneous, turbulent driving of the SSD. Here we analyze more realistic simulations of supernova-driven turbulence to understand how it drives an SSD. We find that SSD growth rates are intermittently variable as a result of the evolving multiphase ISM structure. Rapid growth in the magnetic field typically occurs in hot gas, with the highest overall growth rates occurring when the fractional volume of hot gas is large. SSD growth rates correlate most strongly with vorticity and fluid Reynolds number, which also both correlate strongly with gas temperature. Rotational energy exceeds irrotational energy in all phases, but particularly in the hot phase while SSD growth is most rapid. Supernova rate does not significantly affect the ISM average kinetic energy density. Rather, higher temperatures associated with high supernova rates tend to increase SSD growth rates. SSD saturates with total magnetic energy density around 5% of equipartition to kinetic energy density, increasing slightly with magnetic Prandtl number. While magnetic energy density in the hot gas can exceed that of the othermore »
-
Abstract Subsonic, compressive turbulence transfers energy to cosmic rays (CRs), a process known as nonresonant reacceleration. It is often invoked to explain the observed ratios of primary to secondary CRs at ∼GeV energies, assuming wholly diffusive CR transport. However, such estimates ignore the impact of CR self-confinement and streaming. We study these issues in stirring box magnetohydrodynamic (MHD) simulations using Athena++, with field-aligned diffusive and streaming CR transport. For diffusion only, we find CR reacceleration rates in good agreement with analytic predictions. When streaming is included, reacceleration rates depend on plasma
β . Due to streaming-modified phase shifts between CR and gas variables, they are slower than canonical reacceleration rates in low-β environments like the interstellar medium but remain unchanged in high-β environments like the intracluster medium. We also quantify the streaming energy-loss rate in our simulations. For sub-Alfvénic turbulence, it is resolution dependent (hence unconverged in large-scale simulations) and heavily suppressed compared to the isotropic loss ratev A· ∇P CR/P CR∼v A/L 0, due to misalignment between the mean field and isotropic CR gradients. Unlike acceleration efficiencies, CR losses are almost independent of magnetic field strength overβ ∼ 1–100 and are, therefore, not the primary factor behind lower acceleration rates when streaming is included. While this paper is primarilymore » -
ABSTRACT Models for cosmic ray (CR) dynamics fundamentally depend on the rate of CR scattering from magnetic fluctuations. In the ISM, for CRs with energies ∼MeV-TeV, these fluctuations are usually attributed either to ‘extrinsic turbulence’ (ET) – a cascade from larger scales – or ‘self-confinement’ (SC) – self-generated fluctuations from CR streaming. Using simple analytic arguments and detailed ‘live’ numerical CR transport calculations in galaxy simulations, we show that both of these, in standard form, cannot explain even basic qualitative features of observed CR spectra. For ET, any spectrum that obeys critical balance or features realistic anisotropy, or any spectrum that accounts for finite damping below the dissipation scale, predicts qualitatively incorrect spectral shapes and scalings of B/C and other species. Even if somehow one ignored both anisotropy and damping, observationally required scattering rates disagree with ET predictions by orders of magnitude. For SC, the dependence of driving on CR energy density means that it is nearly impossible to recover observed CR spectral shapes and scalings, and again there is an orders-of-magnitude normalization problem. But more severely, SC solutions with super-Alfvénic streaming are unstable. In live simulations, they revert to either arbitrarily rapid CR escape with zero secondary production, ormore »
