This content will become publicly available on May 1, 2023

Kinetic modeling of neutral transport for a continuum gyrokinetic code
We present the first-of-its-kind coupling of a continuum full- f gyrokinetic turbulence model with a 6D continuum model for kinetic neutrals, carried out using the Gkeyll code. Our objective is to improve the first-principle understanding of the role of neutrals in plasma fueling, detachment, and their interaction with edge plasma profiles and turbulence statistics. Our model includes only atomic hydrogen and incorporates electron-impact ionization, charge exchange, and wall recycling. These features have been successfully verified with analytical predictions and benchmarked with the DEGAS2 Monte Carlo neutral code. We carry out simulations for a scrape-off layer (SOL) with simplified geometry and National Spherical Torus Experiment parameters. We compare these results to a baseline simulation without neutrals and find that neutral interactions reduce the normalized density fluctuation levels and associated skewness and kurtosis, while increasing auto-correlation times. A flatter density profile is also observed, similar to the SOL density shoulder formation in experimental scenarios with high fueling.
Authors:
; ; ; ; ; ; ; ;
Award ID(s):
Publication Date:
NSF-PAR ID:
10332448
Journal Name:
Physics of Plasmas
Volume:
29
Issue:
5
Page Range or eLocation-ID:
052501
ISSN:
1070-664X
1. 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\,{\rmmore » 2. 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 ofmore » 3. 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 thatmore »