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. 
                        more » 
                        « less   
                    
                            
                            Damping of MHD turbulence in a partially ionized medium
                        
                    
    
            ABSTRACT 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. 
        more » 
        « less   
        
    
                            - Award ID(s):
- 2307840
- PAR ID:
- 10475297
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 527
- Issue:
- 2
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 3945-3961
- Size(s):
- p. 3945-3961
- Sponsoring Org:
- National Science Foundation
More Like this
- 
            
- 
            Abstract We study the properties of sub-Alfvénic magnetohydrodynamic (MHD) turbulence, i.e., turbulence with Alfvén mach numberMA=VL/VA< 1, whereVLis the velocity at the injection scale andVAis the Alfvén velocity. We demonstrate that MHD turbulence can have different properties, depending on whether it is driven by velocity or magnetic fluctuations. If the turbulence is driven by isotropic bulk forces acting upon the fluid, i.e., is velocity driven, in an incompressible conducting fluid we predict that the kinetic energy is times larger than the energy of magnetic fluctuations. This effect arises from the long parallel wavelength tail of the forcing, which excites modes withk∥/k⊥<MA. We also predict that as the MHD turbulent cascade reaches the strong regime, the energy of slow modes exceeds the energy of Alfvén modes by a factor . These effects are absent if the turbulence is driven through magnetic fluctuations at the injection scale. We confirm our predictions with numerical simulations. Since the assumption of magnetic and kinetic energy equipartition is at the core of the Davis–Chandrasekhar–Fermi (DCF) approach to measuring magnetic field strength in sub-Alfvénic turbulence, we conclude that the DCF technique is not universally applicable. In particular, we suggest that the dynamical excitation of long azimuthal wavelength modes in the galactic disk may compromise the use of the DCF technique. We discuss alternative expressions that can be used to obtain magnetic field strength from observations and consider ways of distinguishing the cases of velocity and magnetically driven turbulence using observational data.more » « less
- 
            Abstract In a strongly magnetized, magnetically dominated relativistic plasma, Alfvénic turbulence can extend to scales much smaller than the particle inertial scales. It leads to an energy cascade somewhat analogous to inertial- or kinetic-Alfvén turbulent cascades existing in nonrelativistic space and astrophysical plasmas. Based on phenomenological modeling and particle-in-cell numerical simulations, we propose that the energy spectrum of such relativistic kinetic-scale Alfvénic turbulence is close tok−3or slightly steeper than that due to intermittency corrections or Landau damping. We note the analogy of this spectrum with the Kraichnan spectrum corresponding to the enstrophy cascade in 2D incompressible fluid turbulence. Such turbulence strongly energizes particles in the direction parallel to the background magnetic field, leading to nearly one-dimensional particle momentum distributions. We find that these distributions have universal log-normal statistics.more » « less
- 
            Abstract Magnetic fields and turbulence are fundamental to the evolutions of galaxies, yet their precise measurement and analysis present significant challenges. The recently developed Velocity Gradient Technique (VGT), which capitalizes on the anisotropy inherent in magnetohydrodynamic (MHD) turbulence, represents a new method for mapping magnetic fields in galaxies using spectroscopic observations. Most validations of VGT thus far have relied upon idealized MHD turbulence simulations, however, which lack the more complex dynamics found in galaxies and galaxy mergers. In this study, we scrutinize VGT using an AREPO-based cosmological galaxy merger simulation, testing its effectiveness across pre-merger, merging, and post-merger stages. We examine the underlying assumptions of VGT and probe the statistics of gas density, velocity, and magnetic fields over time. We find that the velocity fluctuations are indeed anisotropic at each stage, being larger in the direction perpendicular to the local magnetic field, as required by VGT. We find additionally that galaxy mergers substantially intensify the velocity and density fluctuations and amplify the magnetic fields at all scales. The observed scaling of the velocity fluctuations shows a steeper trend thanr1/2between 0.6 and 3 kpc and a shallower trend at larger scales. The scaling of the magnetic field and density fluctuations at scales ≲1.0 kpc also predominantly aligns withr1/2. Finally, we compare results from VGT to those derived from polarization-like mock magnetic field measurements, finding consistent and statistically significant global agreement in all cases.more » « less
- 
            Magnetohydrodynamic (MHD) turbulence is a cross-field process relevant to many systems. A prerequisite for understanding these systems is to constrain the role of MHD turbulence, and in particular, the energy exchange between kinetic and magnetic forms. The energetics of strongly magnetized and compressible turbulence has so far resisted attempts to understand them. Numerical simulations reveal that kinetic energy can be orders of magnitude higher than fluctuating magnetic energy. We solved this lack-of-balance puzzle by calculating the energetics of compressible and sub-Alfvénic turbulence based on the dynamics of coherent cylindrical fluid parcels. Using the MHD Lagrangian, we proved analytically that the bulk of the magnetic energy transferred to kinetic energy is the energy that is stored in the coupling between the ordered and fluctuating magnetic field. The analytical relations are in strikingly good agreement with numerical data, up to second-order terms.more » « less
 An official website of the United States government
An official website of the United States government 
				
			 
					 
					
