Measures of intermittency in driven supersonic flows

Turbulent stagnation in a Z-pinch plasma

The ion kinetic energy in a stagnating plasma was previously determined by Kroupp et al. [Phys. Rev. Lett. 107, 105001 (2011)PRLTAO0031-900710.1103/PhysRevLett.107.105001] from Doppler-dominated line shapes augmented by measurements of plasma properties and assuming a uniform-plasma model. Notably, the energy was found to be dominantly stored in hydrodynamic flow. Here we advance a new description of this stagnation as supersonically turbulent. Such turbulence implies a nonuniform density distribution. We demonstrate how to reanalyze the spectroscopic data consistent with the turbulent picture and show that this leads to better concordance of the overconstrained spectroscopic measurements, while also substantially lowering the inferred mean density.

Variations of characteristic time scales in rotating stratified turbulence using a large parametric numerical study

We study rotating stratified turbulence (RST) making use of numerical data stemming from a large parametric study varying the Reynolds, Froude and Rossby numbers, Re, Fr and Ro in a broad range of values. The computations are performed using periodic boundary conditions on grids of 1024(3) points, with no modeling of the small scales, no forcing and with large-scale random initial conditions for the velocity field only, and there are altogether 65 runs analyzed in this paper. The buoyancy Reynolds number defined as R(B) = ReFr2 varies from negligible values to ≈ 10(5), approaching atmospheric or oceanic regimes. This preliminary analysis deals with the variation of characteristic time scales of RST with dimensionless parameters, focusing on the role played by the partition of energy between the kinetic and potential modes, as a key ingredient for modeling the dynamics of such flows. We find that neither rotation nor the ratio of the Brunt-Väisälä frequency to the inertial frequency seem to play a major role in the absence of forcing in the global dynamics of the small-scale kinetic and potential modes. Specifically, in these computations, mostly in regimes of wave turbulence, characteristic times based on the ratio of energy to dissipation of the velocity and temperature fluctuations, T(V) and T(P), vary substantially with parameters. Their ratio γ=T(V)/T(P) follows roughly a bell-shaped curve in terms of Richardson number Ri. It reaches a plateau - on which time scales become comparable, γ≈0.6 - when the turbulence has significantly strengthened, leading to numerous destabilization events together with a tendency towards an isotropization of the flow.

Small-scale behavior of Hall magnetohydrodynamic turbulence

Decaying Hall magnetohydrodynamic (HMHD) turbulence is studied using three-dimensional (3D) direct numerical simulations with grids up to 768(3) points and two different types of initial conditions. Results are compared to analogous magnetohydrodynamic (MHD) runs and both Laplacian and Laplacian-squared dissipative operators are examined. At scales below the ion inertial length, the ratio of magnetic to kinetic energy as a function of wave number transitions to a magnetically dominated state. The transition in behavior is associated with the advection term in the momentum equation becoming subdominant to dissipation. Examination of autocorrelation functions reveals that, while current and vorticity structures are similarly sized in MHD, HMHD current structures are narrower and vorticity structures are wider. The electric field autocorrelation function is significantly narrower in HMHD than in MHD and is similar to the HMHD current autocorrelation function at small separations. HMHD current structures are found to be significantly more intense than in MHD and appear to have an enhanced association with strong alignment between the current and magnetic field, which may be important in collisionless plasmas where field-aligned currents can be unstable. When hyperdiffusivity is used, a longer region consistent with a k(-7/3) scaling is present for right-polarized fluctuations when compared to Laplacian dissipation runs.

Scaling and statistics in three-dimensional compressible turbulence

The scaling and statistical properties of three-dimensional compressible turbulence are studied using high-resolution numerical simulations and a heuristic model. The two-point statistics of the solenoidal component of the velocity field are found to be not significantly different from those of incompressible turbulence, while the scaling exponents of the velocity structure function for the compressive component become saturated at high orders. Both the simulated flow and the heuristic model reveal the presence of a power-law tail in the probability density function of negative velocity divergence (high compression regime). The power-law exponent is different from that in Burgers turbulence, and this difference is shown to have a major contribution from the pressure effect, which is absent in the Burgers turbulence.

Exact relation for correlation functions in compressible isothermal turbulence

Compressible isothermal turbulence is analyzed under the assumption of homogeneity and in the asymptotic limit of a high Reynolds number. An exact relation is derived for some two-point correlation functions which reveals a fundamental difference with the incompressible case. The main difference resides in the presence of a new type of term which acts on the inertial range similarly as a source or a sink for the mean energy transfer rate. When isotropy is assumed, compressible turbulence may be described by the relation -2/3ε(eff)r = F(r)(r), where F(r) is the radial component of the two-point correlation functions and ε(eff) is an effective mean total energy injection rate. By dimensional arguments, we predict that a spectrum in k(-5/3) may still be preserved at small scales if the density-weighted fluid velocity ρ(1/3)u is used.

Passive scalar structures in supersonic turbulence

We conduct a systematic numerical study of passive scalar structures in supersonic turbulent flows. We find that the degree of intermittency in the scalar structures increases only slightly as the flow changes from transonic to highly supersonic, while the velocity structures become significantly more intermittent. This difference is due to the absence of shocklike discontinuities in the scalar field. The structure functions of the scalar field are well described by the intermittency model of She and Lévêque [Phys. Rev. Lett 72, 336 (1994)], and the most intense scalar structures are found to be sheetlike at all Mach numbers.

Intermittency and universality in fully developed inviscid and weakly compressible turbulent flows

We perform high-resolution numerical simulations of homogenous and isotropic compressible turbulence, with an average 3D Mach number close to 0.3. We study the statistical properties of intermittency for velocity, density, and entropy. For the velocity field, which is the only quantity that can be compared to the isotropic incompressible case, we find no statistical differences in its behavior in the inertial range due either to the slight compressibility or to the different dissipative mechanism. For the density field, we find evidence of "frontlike" structures, although no shocks are produced by the simulation.

HDR VolVis: high dynamic range volume visualization

In this paper, we present an interactive high dynamic range volume visualization framework (HDR VolVis) for visualizing volumetric data with both high spatial and intensity resolutions. Volumes with high dynamic range values require high precision computing during the rendering process to preserve data precision. Furthermore, it is desirable to render high resolution volumes with low opacity values to reveal detailed internal structures, which also requires high precision compositing. High precision rendering will result in a high precision intermediate image (also known as high dynamic range image). Simply rounding up pixel values to regular display scales will result in loss of computed details. Our method performs high precision compositing followed by dynamic tone mapping to preserve details on regular display devices. Rendering high precision volume data requires corresponding resolution in the transfer function. To assist the users in designing a high resolution transfer function on a limited resolution display device, we propose a novel transfer function specification interface with nonlinear magnification of the density range and logarithmic scaling of the color/ opacity range. By leveraging modern commodity graphics hardware, multiresolution rendering techniques and out-of-core acceleration, our system can effectively produce an interactive visualization of large volume data, such as 2,048(3).

Inertial range scaling in numerical turbulence with hyperviscosity

Numerical turbulence with hyperviscosity is studied and compared with direct simulations using ordinary viscosity and data from wind tunnel experiments. It is shown that the inertial range scaling is similar in all three cases. Furthermore, the bottleneck effect is approximately equally broad (about one order of magnitude) in these cases and only its height is increased in the hyperviscous case-presumably as a consequence of the steeper decent of the spectrum in the hyperviscous subrange. The mean normalized dissipation rate is found to be in agreement with both wind tunnel experiments and direct simulations. The structure function exponents agree with the She-Leveque model. Decaying turbulence with hyperviscosity still gives the usual t(-1.25) decay law for the kinetic energy, and also the bottleneck effect is still present and about equally strong.

Simulations of nonhelical hydromagnetic turbulence

Nonhelical hydromagnetic forced turbulence is investigated using large scale simulations on up to 256 processors and 1024(3) mesh points. The magnetic Prandtl number is varied between 1/8 and 30, although in most cases it is unity. When the magnetic Reynolds number is based on the inverse forcing wave number, the critical value for dynamo action is shown to be around 35 for magnetic Prandtl number of unity. For small magnetic Prandtl numbers we find the critical magnetic Reynolds number to increase with decreasing magnetic Prandtl number. The Kazantsev k(3/2) spectrum for magnetic energy is confirmed for the kinematic regime, i.e., when nonlinear effects are still unimportant and when the magnetic Prandtl number is unity. In the nonlinear regime, the energy budget converges for large Reynolds numbers (around 1000) such that for our parameters about 70% is in kinetic energy and about 30% is in magnetic energy. The energy dissipation rates are converged to 30% viscous dissipation and 70% resistive dissipation. Second-order structure functions of the Elsasser variables give evidence for a k(-5/3) spectrum. Nevertheless, the three-dimensional spectrum is close to k(-3/2), but we argue that this is due to the bottleneck effect. The bottleneck effect is shown to be equally strong both for magnetic and nonmagnetic turbulence, but it is far weaker in one-dimensional spectra that are normally studied in laboratory turbulence. Structure function exponents for other orders are well described by the She-Leveque formula, but the velocity field is significantly less intermittent and the magnetic field is more intermittent than the Elsasser variables.