skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Reynolds number dependence of Lagrangian dispersion in direct numerical simulations of anisotropic magnetohydrodynamic turbulence
Large-scale magnetic fields thread through the electrically conducting matter of the interplanetary and interstellar medium, stellar interiors and other astrophysical plasmas, producing anisotropic flows with regions of high-Reynolds-number turbulence. It is common to encounter turbulent flows structured by a magnetic field with a strength approximately equal to the root-mean-square magnetic fluctuations. In this work, direct numerical simulations of anisotropic magnetohydrodynamic (MHD) turbulence influenced by such a magnetic field are conducted for a series of cases that have identical resolution, and increasing grid sizes up to $2048^3$ . The result is a series of closely comparable simulations at Reynolds numbers ranging from 1400 up to 21 000. We investigate the influence of the Reynolds number from the Lagrangian viewpoint by tracking fluid particles and calculating single-particle and two-particle statistics. The influence of Alfvénic fluctuations and the fundamental anisotropy on the MHD turbulence in these statistics is discussed. Single-particle diffusion curves exhibit mildly superdiffusive behaviours that differ in the direction aligned with the magnetic field and the direction perpendicular to it. Competing alignment processes affect the dispersion of particle pairs, in particular at the beginning of the inertial subrange of time scales. Scalings for relative dispersion, which become clearer in the inertial subrange for a larger Reynolds number, can be observed that are steeper than indicated by the Richardson prediction.  more » « less
Award ID(s):
1907876
PAR ID:
10422989
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
944
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    ABSTRACT We present a systematic shearing-box investigation of magnetorotational instability (MRI)-driven turbulence in a weakly collisional plasma by including the effects of an anisotropic pressure stress, i.e. anisotropic (Braginskii) viscosity. We constrain the pressure anisotropy (Δp) to lie within the stability bounds that would be otherwise imposed by kinetic microinstabilities. We explore a broad region of parameter space by considering different Reynolds numbers and magnetic-field configurations, including net vertical flux, net toroidal-vertical flux, and zero net flux. Remarkably, we find that the level of turbulence and angular-momentum transport are not greatly affected by large anisotropic viscosities: the Maxwell and Reynolds stresses do not differ much from the MHD result. Angular-momentum transport in Braginskii MHD still depends strongly on isotropic dissipation, e.g. the isotropic magnetic Prandtl number, even when the anisotropic viscosity is orders of magnitude larger than the isotropic diffusivities. Braginskii viscosity nevertheless changes the flow structure, rearranging the turbulence to largely counter the parallel rate of strain from the background shear. We also show that the volume-averaged pressure anisotropy and anisotropic viscous transport decrease with increasing isotropic Reynolds number (Re); e.g. in simulations with net vertical field, the ratio of anisotropic to Maxwell stress (αA/αM) decreases from ∼0.5 to ∼0.1 as we move from Re ∼ 103 to Re ∼ 104, while 〈4$$\pi$$Δp/B2〉 → 0. Anisotropic transport may thus become negligible at high Re. Anisotropic viscosity nevertheless becomes the dominant source of heating at large Re, accounting for $${\gtrsim } 50 {{\ \rm per\ cent}}$$ of the plasma heating. We conclude by briefly discussing the implications of our results for radiatively inefficient accretion flows on to black holes. 
    more » « less
  2. 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
  3. The inertial subrange of turbulent scales is commonly reflected by a power law signature in ensemble statistics such as the energy spectrum and structure functions – both in theory and from observations. Despite promising findings on the topic of fractal geometries in turbulence, there is no accepted image for the physical flow features corresponding to this statistical signature in the inertial subrange. The present study uses boundary layer turbulence measurements to evaluate the self-similar geometric properties of velocity isosurfaces and investigate their influence on statistics for the velocity signal. The fractal dimension of streamwise velocity isosurfaces, indicating statistical self-similarity in the size of ‘wrinkles’ along each isosurface, is shown to be constant only within the inertial subrange of scales. For the transition between the inertial subrange and production range, it is inferred that the largest wrinkles become increasingly confined by the overall size of large-scale coherent velocity regions such as uniform momentum zones. The self-similarity of isosurfaces yields power-law trends in subsequent one-dimensional statistics. For instance, the theoretical 2/3 power-law exponent for the structure function can be recovered by considering the collective behaviour of numerous isosurface level sets. The results suggest that the physical presence of inertial subrange eddies is manifested in the self-similar wrinkles of isosurfaces. 
    more » « less
  4. The magnetohydrodynamic (MHD) equations, as a collisional fluid model that remains in local thermodynamic equilibrium (LTE), have long been used to describe turbulence in myriad space and astrophysical plasmas. Yet, the vast majority of these plasmas, from the solar wind to the intracluster medium (ICM) of galaxy clusters, are only weakly collisional at best, meaning that significant deviations from LTE are not only possible but common. Recent studies have demonstrated that the kinetic physics inherent to this weakly collisional regime can fundamentally transform the evolution of such plasmas across a wide range of scales. Here, we explore the consequences of pressure anisotropy and Larmor-scale instabilities for collisionless,$$\beta \gg 1$$, turbulence, focusing on the role of a self-organizational effect known as ‘magneto-immutability’. We describe this self-organization analytically through a high-$$\beta$$, reduced ordering of the Chew–Goldberger–Low-MHD (CGL-MHD) equations, finding that it is a robust inertial-range effect that dynamically suppresses magnetic-field-strength fluctuations, anisotropic-pressure stresses and dissipation due to heat fluxes. As a result, the turbulent cascade of Alfvénic fluctuations continues below the putative viscous scale to form a robust, nearly conservative, MHD-like inertial range. These findings are confirmed numerically via Landau-fluid CGL-MHD turbulence simulations that employ a collisional closure to mimic the effects of microinstabilities. We find that microinstabilities occupy a small ($${\sim }5\,\%$$) volume-filling fraction of the plasma, even when the pressure anisotropy is driven strongly towards its instability thresholds. We discuss these results in the context of recent predictions for ion-vs-electron heating in low-luminosity accretion flows and observations implying suppressed viscosity in ICM turbulence. 
    more » « less
  5. We report an investigation into random-jet-stirred homogeneous turbulence generated in a vertical octagonal prism-shaped tank where there are jet arrays on four of the eight vertical faces. We show that the turbulence is homogeneous at all scales in the central region of the tank that span multiple integral scales in all directions. The jet forcing from four sides in the horizontal direction guarantees isotropy in horizontal planes but leads to more energy in the horizontal fluctuations com- pared with the vertical fluctuations. This anisotropy between the horizontal and vertical fluctuations decreases at smaller scales, so that the inertial and dissipation range statistics show isotropic behavior. Using four jet arrays allows us to achieve higher turbulence intensity and Reynolds number with a shorter jet merging distance compared to previous facilities with two-facing arrays. By changing the array-to-array distance, the parameters of the algorithm that drives random-jet stirring, and attachments to the exits of each jet, we show that we are able to vary the turbulence scales and Reynolds number. We provide scaling relations for the turbulent fluctuation velocity, integral scale, and dissipation rate, and we show how these scales of motion are primarily determined by the properties of individual jets and the diffusion of their momentum with distance from the nozzles. Finally, we examine the signatures of individual jets in the turbulent velocity spectra and report the conditions under which individual jet flows, not fully mixed with the background turbulence, produce a spectral peak and the corresponding frequency associated with the jet forcing timescale. 
    more » « less