More Like this

1. Shearing-box simulations of MRI-driven turbulence in weakly collisional accretion discs

   https://doi.org/10.1093/mnras/stz1111  

   Kempski, Philipp ; Quataert, Eliot ; Squire, Jonathan ; Kunz, Matthew W ( July 2019 , Monthly Notices of the Royal Astronomical Society)

   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. Kelvin–Helmholtz Instability “Tube” and “Knot” Dynamics. Part II: KHI T&K Dynamics in a Multiscale Gravity Wave Direct Numerical Simulation

   https://doi.org/10.1175/JAS-D-22-0193.1  

   Fritts, David C. ; Wang, Ling ( October 2023 , Journal of the Atmospheric Sciences)

   A companion paper by Fritts et al. reviews evidence for Kelvin–Helmholtz instability (KHI) “tube” and “knot” (T&K) dynamics that appear to be widespread throughout the atmosphere. Here we describe the results of an idealized direct numerical simulation of multiscale gravity wave dynamics that reveals multiple larger- and smaller-scale KHI T&K events. The results enable assessments of the environments in which these dynamics arise and their competition with concurrent gravity wave breaking in driving turbulence and energy dissipation. A larger-scale event is diagnosed in detail and reveals diverse and intense T&K dynamics driving more intense turbulence than occurs due to gravity wave breaking in the same environment. Smaller-scale events reveal that KHI T&K dynamics readily extend to weaker, smaller-scale, and increasingly viscous shear flows. Our results suggest that KHI T&K dynamics should be widespread, perhaps ubiquitous, wherever superposed gravity waves induce intensifying shear layers, because such layers are virtually always present. A second companion paper demonstrates that KHI T&K dynamics exhibit elevated turbulence generation and energy dissipation rates extending to smaller Reynolds numbers for relevant KHI scales wherever they arise. These dynamics are suggested to be significant sources of turbulence and mixing throughout the atmosphere that are currently ignored or underrepresented in turbulence parameterizations in regional and global models.

   Atmospheric observations reveal that Kelvin–Helmholtz instabilities (KHI) often exhibit complex interactions described as “tube” and “knot” (T&K) dynamics in the presence of larger-scale gravity waves (GWs). These dynamics may prove to make significant contributions to energy dissipation and mixing that are not presently accounted for in large-scale modeling and weather prediction. We explore here the occurrence of KHI T&K dynamics in an idealized model that describes their behavior and character arising at larger and smaller scales due to superposed, large-amplitude GWs. The results reveal that KHI T&K dynamics arise at larger and smaller scales, and that their turbulence intensities can be comparable to those of the GWs.

    

   more » « less
3. Asymptotic behaviour of rotating convection-driven dynamos in the plane layer geometry

   https://doi.org/10.1017/jfm.2022.848  

   Yan, Ming ; Calkins, Michael A. ( November 2022 , Journal of Fluid Mechanics)

   Dynamos driven by rotating convection in the plane layer geometry are investigated numerically for a range of Ekman number ( $E$ ), magnetic Prandtl number ( $Pm$ ) and Rayleigh number ( $Ra$ ). The primary purpose of the investigation is to compare results of the simulations with previously developed asymptotic theory that is applicable in the limit of rapid rotation. We find that all of the simulations are in the quasi-geostrophic regime in which the Coriolis and pressure gradient forces are approximately balanced at leading order, whereas all other forces, including the Lorentz force, act as perturbations. Agreement between simulation output and asymptotic scalings for the energetics, flow speeds, magnetic field amplitude and length scales is found. The transition from large-scale dynamos to small-scale dynamos is well described by the magnetic Reynolds number based on the small convective length scale, $\widetilde {Rm}$ , with large-scale dynamos preferred when $\widetilde {Rm} \lesssim O(1)$ . The magnitude of the large-scale magnetic field is observed to saturate and become approximately constant with increasing Rayleigh number. Energy spectra show that all length scales present in the flow field and the small-scale magnetic field are consistent with a scaling of $E^{1/3}$ , even in the turbulent regime. For a fixed value of $E$ , we find that the viscous dissipation length scale is approximately constant over a broad range of $Ra$ ; the ohmic dissipation length scale is approximately constant within the large-scale dynamo regime, but transitions to a $\widetilde {Rm}^{-1/2}$ scaling in the small-scale dynamo regime. 

   more » « less
4. Extending a Physics-informed Machine-learning Network for Superresolution Studies of Rayleigh–Bénard Convection

   https://doi.org/10.3847/1538-4357/ad1c55  

   Salim, Diane M. ; Burkhart, Blakesley ; Sondak, David ( March 2024 , The Astrophysical Journal)

   Advancing our understanding of astrophysical turbulence is bottlenecked by the limited resolution of numerical simulations that may not fully sample scales in the inertial range. Machine-learning (ML) techniques have demonstrated promise in upscaling resolution in both image analysis and numerical simulations (i.e., superresolution). Here we employ and further develop a physics-constrained convolutional neural network ML model called “MeshFreeFlowNet” (MFFN) for superresolution studies of turbulent systems. The model is trained on both the simulation images and the evaluated partial differential equations (PDEs), making it sensitive to the underlying physics of a particular fluid system. We develop a framework for 2D turbulent Rayleigh–Bénard convection generated with theDedaluscode by modifying the MFFN architecture to include the full set of simulation PDEs and the boundary conditions. Our training set includes fully developed turbulence sampling Rayleigh numbers (Ra) ofRa= 10⁶–10¹⁰. We evaluate the success of the learned simulations by comparing the power spectra of the directDedalussimulation to the predicted model output and compare both ground-truth and predicted power spectral inertial range scalings to theoretical predictions. We find that the updated network performs well at allRastudied here in recovering large-scale information, including the inertial range slopes. The superresolution prediction is overly dissipative at smaller scales than that of the inertial range in all cases, but the smaller scales are better recovered in more turbulent than laminar regimes. This is likely because more turbulent systems have a rich variety of structures at many length scales compared to laminar flows.

    

   more » « less
5. Compressible two-dimensional turbulence: cascade reversal and sensitivity to imposed magnetic field

   https://doi.org/10.1088/1367-2630/ad0172  

   Fouxon, Itzhak ; Kritsuk, Alexei G. ; Mond, Michael ( November 2023 , New Journal of Physics)

   We study the impact of compressibility on two-dimensional turbulent flows, such as those modeling astrophysical disks. We demonstrate that the direction of cascade undergoes continuous transition as the Mach numberMaincreases, from inverse atMa = 0, to direct at$Ma=\mathrm{\infty }$. Thus, at$Ma\sim 1$comparable amounts of energy flow from the pumping scale to large and small scales, in accord with previous data. For supersonic turbulence with$Ma\gg 1$the cascade is direct, as in three dimensions, which results in multifractal density field. For that regime ($Ma\gg 1$) we derive a Kolmogorov-type law for potential forcing and obtain an explicit expression for the third order correlation tensor of the velocity. We further show that all third order structure functions are zero up to first order in the inertial range scales, which is in sharp contrast with incompressible turbulence where the third order structure function, that describes the energy flux associated with the energy cascade is non-zero. The properties of compressible turbulence have significant implications on the amplification of magnetic fields in conducting fluids. We thus demonstrate that imposing external magnetic field on compressible flows of conducting fluids allows to manipulate the flow producing possibly large changes even at small Mach numbers. Thus Zeldovich’s antidynamo theorem, by which atMa = 0 the magnetic field is zero in the steady state, must be used with caution. Real flows have finiteMaand, however small it is, for large enough values ofI, the magnetic flux through the disk, the magnetic field changes the flow appreciably, or rearranges it completely. This renders the limitMa → 0 singular for non-zero values ofI. Of particular interest is the effect of the density multifractality, at$Ma\gg 1$which is relevant for astrophysical disks. We demonstrate that in that regime, in the presence of non-zeroIthe magnetic field energy is enhanced by a large factor as compared to its estimates based on the mean field. Finally, based on the insights described above, we propose a novel two-dimensional Burgers’ turbulence, whose three-dimensional counterpart is used for studies of the large-scale structure of the Universe, as a model for supersonic two-dimensional magnetohydrodynamic flows.

    

   more » « less