The coupling state between ions and neutrals in the interstellar medium plays a key role in the dynamics of magnetohydrodynamic (MHD) turbulence, but is challenging to study numerically. In this work, we investigate the damping of MHD turbulence in a partially ionized medium using 3D two-fluid (ions + neutrals) simulations generated with the athenak code. Specifically, we examine the velocity, density, and magnetic field statistics of the two-fluid MHD turbulence in different regimes of neutral-ion coupling. Our results demonstrate that when ions and neutrals are strongly coupled, the velocity statistics resemble those of single-fluid MHD turbulence. Both the velocity structures and kinetic energy spectra of ions and neutrals are similar, while their density structures can be significantly different. With an excess of small-scale sharp density fluctuations in ions, the density spectrum in ions is shallower than that of neutrals. When ions and neutrals are weakly coupled, the turbulence in ions is more severely damped due to the ion-neutral collisional friction than that in neutrals, resulting in a steep kinetic energy spectrum and density spectrum in ions compared to the Kolmogorov spectrum. We also find that the magnetic energy spectrum basically follows the shape of the kinetic energy spectrum of ions, irrespective of the coupling regime. In addition, we find large density fluctuations in ions and neutrals and thus spatially inhomogeneous ionization fractions. As a result, the neutral-ion decoupling and damping of MHD turbulence take place over a range of length-scales.
- Award ID(s):
- 2019828
- NSF-PAR ID:
- 10332448
- Date Published:
- Journal Name:
- Physics of Plasmas
- Volume:
- 29
- Issue:
- 5
- ISSN:
- 1070-664X
- Page Range / eLocation ID:
- 052501
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
In this paper we present an experimental study of edge turbulence in the Large Plasma Device at UCLA. We utilize a scan of discharge power and prefill pressure (neutral density) to show experimentally that turbulent density fluctuations decrease with decreasing density gradient, as predicted for resistive drift-wave turbulence (RDWT). As expected for RDWT, we observe that the cross-phase between the density and potential fluctuations is close to 0. Moreover, the addition of an electron temperature gradient leads to a reduction in the amplitude of the density fluctuations, as expected for RDWT. However, counter to theoretical expectations, we find that the potential fluctuations do not follow the same trends as the density fluctuations for changes either in density gradients or the addition of a temperature gradient. The disconnect between the density and potential fluctuations is connected to changes in the parallel flows as a result of differences in the prefill pressure, i.e. neutral density. Further analysis of the density and potential fluctuation spectra show that the electron temperature gradient reduces the low frequency fluctuations up to $10 \,{\rm kHz}$ and the introduction of a temperature gradient leads to an unexpected ${\sim }{\rm \pi}$ shift of the density–potential cross-phase at ${\sim }10\,{\rm kHz}$ , while maintaining the typical resistive drift-wave cross-phase at lower frequencies. These experiments partly confirm existing knowledge on resistive drift-wave turbulence, but also introduce new observations that indicate a need for dedicated nonlinear three-dimensional turbulence simulations that include neutrals.more » « less
-
more » « less
Abstract The matter in an accretion disk must lose angular momentum when moving radially inwards but how this works has long been a mystery. By calculating the trajectories of individual colliding neutrals, ions, and electrons in a weakly ionized 2D plasma containing gravitational and magnetic fields, we numerically simulate accretion disk dynamics at the particle level. As predicted by Lagrangian mechanics, the fundamental conserved global quantity is the total canonical angular momentum, not the ordinary angular momentum. When the Kepler angular velocity and the magnetic field have opposite polarity, collisions between neutrals and charged particles cause: (i) ions to move radially inwards, (ii) electrons to move radially outwards, (iii) neutrals to lose ordinary angular momentum, and (iv) charged particles to gain canonical angular momentum. Neutrals thus spiral inward due to their decrease of ordinary angular momentum while the accumulation of ions at small radius and accumulation of electrons at large radius produces a radially outward electric field. In 3D, this radial electric field would drive an out-of-plane poloidal current that produces the magnetic forces that drive bidirectional astrophysical jets. Because this neutral angular momentum loss depends only on neutrals colliding with charged particles, it should be ubiquitous. Quantitative scaling of the model using plausible disk density, temperature, and magnetic field strength gives an accretion rate of 3 × 10−8solar mass per year, which is in good agreement with observed accretion rates.
-
more » « less
Abstract Subauroral Polarization Streams (SAPS) are associated with closure of region 2 field‐aligned current (R2 FAC) through the low conductivity region. Although SAPS have often been studied from a magnetosphere‐ionosphere coupling perspective, recent observations suggest strong interaction also exists between SAPS and the thermosphere. Our study focuses on thermospheric wind driving and its impact on SAPS and R2 FAC during the 17 March 2013 geomagnetic storm using both observations and the physics‐based Rice Convection Model‐Coupled Thermosphere, Ionosphere, Plasmasphere, electrodynamics (RCM‐CTIPe) model that self‐consistently couples the magnetosphere‐ionosphere‐thermosphere system. Defense Meteorological Satellite Program (DMSP)‐18 and Gravity Field and Steady‐State Ocean Circulation Explorer (GOCE) satellite observations show that, as the storm progresses, sunward ion flows intensify and expand equatorward and are accompanied by strengthening of subauroral neutral winds with some delay. Our model successfully reproduces time evolution and overall structure of the sunward ion drift and neutral wind. A force term analysis is performed to investigate the momentum transfer to the neutrals from the ions. Contrary to previous studies showing that Coriolis force is the main driver of neutrals during storm time, we find that the ion drag is the largest force driving westward neutral wind in the SAPS region where the ion density is low in the trough region. Furthermore, simulations with and without the neutral wind dynamo effect are compared to quantify the effect of the neutral to plasma flow. The comparison shows that the self‐consistent active ionosphere thermosphere coupling increases the R2 FAC and the westward ion drift equatorward of the SAPS region by 20% and 40% by the flywheel effect, respectively.
-
null (Ed.)ABSTRACT In the cold neutral medium, high out-of-equilibrium temperatures are created by intermittent dissipation processes, including shocks, viscous heating, and ambipolar diffusion. The high-temperature excursions are thought to explain the enhanced abundance of CH+ observed along diffuse molecular sightlines. Intermittent high temperatures should also have an impact on H2 line luminosities. We carry out simulations of magnetohydrodynamic (MHD) turbulence in molecular clouds including heating and cooling, and post-process them to study H2 line emission and hot-gas chemistry, particularly the formation of CH+. We explore multiple magnetic field strengths and equations of state. We use a new H2 cooling function for $n_{\text{H}}\le 10^5\, {\text{cm}}^{-3}$, $T\le 5000\, {\text{K}}$, and variable H2 fraction. We make two important simplifying assumptions: (i) the H2/H fraction is fixed everywhere and (ii) we exclude from our analysis regions where the ion–neutral drift velocity is calculated to be greater than 5 km s−1. Our models produce H2 emission lines in accord with many observations, although extra excitation mechanisms are required in some clouds. For realistic root-mean-square (rms) magnetic field strengths (≈10 μG) and velocity dispersions, we reproduce observed CH+ abundances. These findings contrast with those of Valdivia et al. (2017) Comparison of predicted dust polarization with observations by Planck suggests that the mean field is ≳5 µG, so that the turbulence is sub-Alfvénic. We recommend future work treating ions and neutrals as separate fluids to more accurately capture the effects of ambipolar diffusion on CH+ abundance.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1063/5.0087131
