Search

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 1,400 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\'enic fluctuations and the fundamental anisotropy on the MHD turbulence in these statistics is discussed. Single-particle diffusion curves exhibit mildly superdiffusive behaviors 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 larger Reynolds number, can be observed that are steeper than indicated by the Richardson prediction., Comment: 23 pages, 11 figures

We investigate the utility of the convex hull of many Lagrangian tracers to analyze transport properties of turbulent flows with different anisotropy. In direct numerical simulations of statistically homogeneous and stationary Navier-Stokes turbulence, neutral fluid Boussinesq convection, and MHD Boussinesq convection a comparison with Lagrangian pair dispersion shows that convex hull statistics capture the asymptotic dispersive behavior of a large group of passive tracer particles. Moreover, convex hull analysis provides additional information on the sub-ensemble of tracers that on average disperse most efficiently in the form of extreme value statistics and flow anisotropy via the geometric properties of the convex hulls. We use the convex hull surface geometry to examine the anisotropy that occurs in turbulent convection. Applying extreme value theory, we show that the maximal square extensions of convex hull vertices are well described by a classic extreme value distribution, the Gumbel distribution. During turbulent convection, intermittent convective plumes grow and accelerate the dispersion of Lagrangian tracers. Convex hull analysis yields information that supplements standard Lagrangian analysis of coherent turbulent structures and their influence on the global statistics of the flow., Comment: 18 pages, 10 figures, preprint

Local regions of anomalous particle dispersion, and intermittent events that occur in turbulent flows can greatly influence the global statistical description of the flow. These local behaviors can be identified and analyzed by comparing the growth of neighboring convex hulls of Lagrangian tracer particles. Although in our simulations of homogeneous turbulence the convex hulls generally grow in size, after the Lagrangian particles that define the convex hulls begin to disperse, our analysis reveals short periods when the convex hulls of the Lagrangian particles shrink, evidence that particles are not dispersing simply. Shrinkage can be associated with anisotropic flows, since it occurs most frequently in the presence of a mean magnetic field or thermal convection. We compare dispersion between a wide range of statistically homogeneous and stationary turbulent flows ranging from homogeneous isotropic Navier-Stokes turbulence over different configurations of magnetohydrodynamic turbulence and Boussinesq convection., Comment: The following article has been submitted to Physics of Fluids. After it is published, it will be found at (http://scitation.aip.org/content/aip/journal/pof2). Copyright 2014 J. Pratt, A. Busse, W.-C. Mueller, S.C. Chapman, and N.W. Watkins. This article is distributed under a Creative Commons Attribution 3.0 Unported License

Intermittent large-scale high-shear flows are found to occur frequently and spontaneously in direct numerical simulations of statistically stationary turbulent Boussinesq magnetohydrodynamic (MHD) convection. The energetic steady-state of the system is sustained by convective driving of the velocity field and small-scale dynamo action. The intermittent emergence of flow structures with strong velocity and magnetic shearing generates magnetic energy at an elevated rate over time-scales longer than the characteristic time of the large-scale convective motion. The resilience of magnetic energy amplification suggests that intermittent shear-bursts are a significant driver of dynamo action in turbulent magnetoconvection.

A key prediction of turbulence theories is frame-invariance, and in magnetohydrodynamic (MHD) turbulence, axisymmetry of fluctuations with respect to the background magnetic field. Paradoxically the power in fluctuations in the turbulent solar wind are observed to be ordered with respect to the bulk macroscopic flow as well as the background magnetic field. Here, non- axisymmetry across the inertial and dissipation ranges is quantified using in-situ observations from Cluster. The observed inertial range non- axisymmetry is reproduced by a 'fly through' sampling of a Direct Numerical Simulation of MHD turbulence. Furthermore, 'fly through' sampling of a linear superposition of transverse waves with axisymmetric fluctuations generates the trend in non- axisymmetry with power spectral exponent. The observed non-axisymmetric anisotropy may thus simply arise as a sampling effect related to Taylor's hypothesis and is not related to the plasma dynamics itself., Comment: 4 pages, 2 figures, accepted for publication in PRL

We quantify possible differences between turbulent dynamo action in the Sun and the dynamo action studied in idealized simulations. For this purpose we compare Fourier-space shell-to-shell energy transfer rates of three incrementally more complex dynamo simulations: an incompressible, periodic simulation driven by random flow, a simulation of Boussinesq convection, and a simulation of fully compressible convection that includes physics relevant to the near-surface layers of the Sun. For each of the simulations studied, we find that the dynamo mechanism is universal in the kinematic regime because energy is transferred from the turbulent flow to the magnetic field from wavenumbers in the inertial range of the energy spectrum. The addition of physical effects relevant to the solar near-surface layers, including stratification, compressibility, partial ionization, and radiative energy transport, does not appear to affect the nature of the dynamo mechanism. The role of inertial-range shear stresses in magnetic field amplification is independent from outer-scale circumstances, including forcing and stratification. Although the shell-to-shell energy transfer functions have similar properties to those seen in mean-flow driven dynamos in each simulation studied, the saturated states of these simulations are not universal because the flow at the driving wavenumbers is a significant source of energy for the magnetic field., Comment: 16 pages, 9 figures, accepted for publication in ApJ

Navier-Stokes turbulence subject to solid-body rotation is studied by high-resolution direct numerical simulations (DNS) of freely decaying and stationary flows. Setups characterized by different Rossby numbers are considered. In agreement with experimental results strong rotation is found to lead to anisotropy of the direct nonlinear energy flux, which is attenuated primarily in the direction of the rotation axis. In decaying turbulence the evolution of kinetic energy follows an approximate power law with a distinct dependence of the decay exponent on the rotation frequency. A simple phenomenological relation between exponent and rotation rate reproduces this observation. Stationary turbulence driven by large-scale forcing exhibits $k_\perp^{-2}$-scaling in the rotation-dominated inertial range of the one-dimensional energy spectrum taken perpendicularly to the rotation axis. The self-similar scaling is shown to be the cumulative result of individual spectral contributions which, for low rotation rate, display $k_\perp^{-3}$-scaling near the $k_\parallel=0$ plane. A phenomenology which incorporates the modification of the energy cascade by rotation is proposed. In the observed regime the nonlinear turbulent interactions are strongly influenced by rotation but not suppressed. {Longitudinal two-point velocity structure functions taken perpendicularly to the axis of rotation indicate weak intermittency of the $k_\parallel=0$ (2D) component of the flow while the intermittent scaling of $k_\parallel\neq 0$ (3D) fluctuations is well captured by a modified She-L\'ev\^eque intermittency model which yields the expression $\zeta_p=p/6 + 2(1-(2/3)^{p/2})$ for the structure function scaling exponents., Comment: accepted for publication in Journal of Fluid Mechanics, 2 figure in reduced quality (data reduction)

We present a collection of eight data sets, from state-of-the-art experiments and numerical simulations on turbulent velocity statistics along particle trajectories obtained in different flows with Reynolds numbers in the range $R_\lambda \in [120:740]$. Lagrangian structure functions from all data sets are found to collapse onto each other on a wide range of time lags, revealing a universal statistics, and calling for a unified theoretical description. Parisi-Frisch Multifractal theory, suitable extended to the dissipative scales and to the Lagrangian domain, is found to capture intermittency of velocity statistics over the whole three decades of temporal scales here investigated., Comment: 5 pages, 1 figure; content changed, references updated

We report on a comparison of high-resolution numerical simulations of Lagrangian particles advected by incompressible turbulent hydro- and magnetohydrodynamic (MHD) flows. Numerical simulations were performed with up to $1024^3$ collocation points and 10 million particles in the Navier-Stokes case and $512^3$ collocation points and 1 million particles in the MHD case. In the hydrodynamics case our findings compare with recent experiments from Mordant et al. [1] and Xu et al. [2]. They differ from the simulations of Biferale et al. [3] due to differences of the ranges choosen for evaluating the structure functions. In Navier-Stokes turbulence intermittency is stronger than predicted by a multifractal approach of [3] whereas in MHD turbulence the predictions from the multifractal approach are more intermittent than observed in our simulations. In addition, our simulations reveal that Lagrangian Navier-Stokes turbulence is more intermittent than MHD turbulence, whereas the situation is reversed in the Eulerian case. Those findings can not consistently be described by the multifractal modeling. The crucial point is that the geometry of the dissipative structures have different implications for Lagrangian and Eulerian intermittency. Application of the multifractal approach for the modeling of the acceleration PDFs works well for the Navier-Stokes case but in the MHD case just the tails are well described., Comment: to appear in J. Plasma Phys