skip to main content


Title: On the Conservation of Turbulence Energy in Turbulence Transport Models
Abstract Zank et al. developed models describing the transport of low-frequency incompressible and nearly incompressible turbulence in inhomogeneous flows. The formalism was based on expressing the fluctuating variables in terms of the Elsässar variables and then taking “moments” subject to various closure hypotheses. The turbulence transport models are different according to whether the plasma beta regime is large, of order unity, or small. Here, we show explicitly that the three sets of turbulence transport models admit a conservation representation that resembles the well-known WKB transport equation for Alfvén wave energy density after introducing appropriate definitions of the “pressure” associated with the turbulent fluctuations. This includes introducing a distinct turbulent pressure tensor for 3D incompressible turbulence (the large plasma beta limit) and pressure tensors for quasi-2D and slab turbulence (the plasma beta order-unity or small regimes) that generalize the form of the WKB pressure tensor. Various limits of the different turbulent pressure tensors are discussed. However, the analogy between the conservation form of the turbulence transport models and the WKB model is not close for multiple reasons, including that the turbulence models express fully nonlinear physical processes unlike the strictly linear WKB description. The analysis presented here both serves as a check on the validity and correctness of the turbulence transport models and also provides greater transparency of the energy dissipation term and the “turbulent pressure” in our models, which is important for many practical applications.  more » « less
Award ID(s):
1655280
NSF-PAR ID:
10356041
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
928
Issue:
2
ISSN:
0004-637X
Page Range / eLocation ID:
176
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. This work aims at comparing the accuracy and overall performance of a low-Mach CFD solver and a fully-compressible CFD solver for direct numerical simulation (DNS) of nonequilibrium plasma assisted ignition (PAI) using a phenomenological model described in Castela et al. [1]. The phenomenological model describes the impact of nanosecond pulsed plasma discharges by introducing source terms in the reacting flow equations, instead of solving the detailed plasma kinetics at every time step of the discharge. Ultra-fast gas heating and dissociation ofO2 are attributed to the electronic excitation ofN2 and the subsequent quenching to ground state. This process is highly exothermic, and is responsible for dissociation of O2 to form O radicals; both of which promote faster ignition. Another relatively slower process of gas heating associated with vibrational-to-translational relaxation is also accounted for, by solving an additional vibrational energy transport equation. A fully-compressible CFD solver for high Mach (M>0.2) reacting flows, developed by extending the default rhoCentralFoam solver in OpenFOAM, is used to perform DNS of PAI in a 2D domain representing a cross section of a pin-to-pin plasma discharge configuration. The same case is also simulated using a low-Mach, pressure-based CFD solver, built by extending the default reactingFoam solver. The lack of flow or wave dominated transport after the plasma-induced weak shock wave leaves the domain causes inaccurate computation of all the transport variables, with a rather small time step dictated by the CFL condition, with the fully-compressible solver. These issues are not encountered in the low-Mach solver. Finally, the low-Mach solver is used to perform DNS of PAI in lean, premixed, isotropic turbulent mixtures of CH4-air at two different Reynolds numbers of 44 and 395. Local convection of the radicals and vibrational energy from the discharge domain, and straining of the high temperature reaction zones resulted in slower ignition of the case with the higher Re. A cascade effect of temperature reduction in the more turbulent case also resulted in a five - six times smaller value of the vibrational to translational gas heating source term, which further inhibited ignition. Two pulses were sufficient for ignition of the Re = 44 case, whereas three pulses were required for the Re = 395 case; consistent with the results of Ref. [1]. 
    more » « less
  2. Abstract

    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 atMa=. Thus, atMa1comparable amounts of energy flow from the pumping scale to large and small scales, in accord with previous data. For supersonic turbulence withMa1the cascade is direct, as in three dimensions, which results in multifractal density field. For that regime (Ma1) 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, atMa1which 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
  3. ABSTRACT

    The nature of cosmic ray (CR) transport in the Milky Way remains elusive. The predictions of current microphysical CR transport models in magnetohydrodynamic (MHD) turbulence are drastically different from what is observed. These models usually focus on MHD turbulence with a strong guide field and ignore the impact of turbulent intermittency on particle propagation. This motivates our studying the alternative regime of large-amplitude turbulence with δB/B0 ≫ 1, in which intermittent small-scale magnetic field reversals are ubiquitous. We study particle transport in such turbulence by integrating trajectories in stationary snapshots. To quantify spatial diffusion, we use a set-up with continuous particle injection and escape, which we term the turbulent leaky box. We find that particle transport is very different from the strong guide-field case. Low-energy particles are better confined than high-energy particles, despite less efficient pitch-angle isotropization at small energies. In the limit of weak guide field, energy-dependent confinement is driven by the energy-dependent (in)ability to follow reversing magnetic field lines exactly and by the scattering in regions of ‘resonant curvature’, where the field line bends on a scale that is of the order of the local particle gyro-radius. We derive a heuristic model of particle transport in magnetic folds that approximately reproduces the energy dependence of transport found numerically. We speculate that CR propagation in the Galaxy is regulated by the intermittent field reversals highlighted here and discuss the implications of our findings for CR transport in the Milky Way.

     
    more » « less
  4. Abstract Nearly incompressible magnetohydrodynamic (NI MHD) theory for β ∼ 1 (or β ≪ 1) plasma has been developed and applied to the study of solar wind turbulence. The leading-order term in β ∼ 1 or β ≪ 1 plasma describes the majority of 2D turbulence, while the higher-order term describes the minority of slab turbulence. Here, we develop new NI MHD turbulence transport model equations in the high plasma beta regime. The leading-order term in a β ≫ 1 plasma is fully incompressible and admits both structures (flux ropes or magnetic islands) and slab (Alfvén waves) fluctuations. This paper couples the NI MHD turbulence transport equations with three fluid (proton, electron, and pickup ion) equations, and solves the 1D steady-state equations from 1–75 au. The model is tested against 27 yr of Voyager 2 data, and Ulysses and NH SWAP data. The results agree remarkably well, with some scatter, about the theoretical predictions. 
    more » « less
  5. ABSTRACT

    We present the first empirical constraints on the turbulent velocity field of the diffuse circumgalactic medium around four luminous quasi-stellar objects (QSOs) at z ≈ 0.5–1.1. Spatially extended nebulae of ≈50–100 physical kpc in diameter centred on the QSOs are revealed in [O ii] $\lambda \lambda \, 3727,3729$ and/or [O iii] $\lambda \, 5008$ emission lines in integral field spectroscopic observations obtained using Multi-Unit Spectroscopic Explorer on the Very Large Telescope. We measure the second- and third-order velocity structure functions (VSFs) over a range of scales, from ≲5 kpc to ≈20–50 kpc, to quantify the turbulent energy transfer between different scales in these nebulae. While no constraints on the energy injection and dissipation scales can be obtained from the current data, we show that robust constraints on the power-law slope of the VSFs can be determined after accounting for the effects of atmospheric seeing, spatial smoothing, and large-scale bulk flows. Out of the four QSO nebulae studied, one exhibits VSFs in spectacular agreement with the Kolmogorov law, expected for isotropic, homogeneous, and incompressible turbulent flows. The other three fields exhibit a shallower decline in the VSFs from large to small scales. However, with a limited dynamic range in the spatial scales in seeing-limited data, no constraints can be obtained for the VSF slopes of these three nebulae. For the QSO nebula consistent with the Kolmogorov law, we determine a turbulence energy cascade rate of ≈0.2 cm2 s−3. We discuss the implication of the observed VSFs in the context of QSO feeding and feedback in the circumgalactic medium.

     
    more » « less