Search

We perform numerical simulations of decaying rotating stratified turbulence and show, in the Boussinesq framework, that helicity (velocity-vorticity correlation), as observed in supercell storms and hurricanes, is spontaneously created due to an interplay between buoyancy and rotation common to large-scale atmospheric and oceanic flows. Helicity emerges from the joint action of eddies and of inertia-gravity waves (with inertia and gravity with respective associated frequencies f and N), and it occurs when the waves are sufficiently strong. For N/f < 3 the amount of helicity produced is correctly predicted by a quasilinear balance equation. Outside this regime, and up to the highest Reynolds number obtained in this study, namely Re ≈ 10000, helicity production is found to be persistent for N/f as large as ≈ 17, and for ReFr² and ReRo², respectively, as large as ≈ 100 and ≈ 24000. [ABSTRACT FROM AUTHOR]

The decay rate of isotropic and homogeneous turbulence is known to be affected by the large-scale spectrum of the initial perturbations, associated with at least two canonical self-preserving solutions of the von Kármán-Howarth equation: the so-called Batchelor and Saffman spectra. The effect of long-range correlations in the decay of anisotropic flows is less clear, and recently it has been proposed that the decay rate of rotating turbulence may be independent of the large-scale spectrum of the initial perturbations. We analyze numerical simulations of freely decaying rotating turbulence with initial energy spectra ~k4 (Batchelor turbulence) and ~k² (Saffman turbulence) and show that, while a self-similar decay can not be identified for the total energy, the decay is indeed affected by long-range correlations. The decay of two- and three-dimensional modes follows distinct power laws in each case, which are consistent with predictions derived from the anisotropic von Kármán-Howarth equation, and with conservation of anisotropic integral quantities by the flow evolution [ABSTRACT FROM AUTHOR]

Rapidly rotating turbulent flow is characterized by the emergence of columnar structures that are representative of quasi-two-dimensional behavior of the flow. It is known that when energy is injected into the fluid at an intermediate scale Lf, it cascades towards smaller as well as larger scales. In this paper we analyze the flow in the inverse cascade range, at a small but fixed Rossby number, Rof ≈ 0.05. Several numerical simulations with helical and nonhelical forcing functions are considered in periodic boxes with unit aspect ratio. In order to resolve the inverse cascade range with reasonably large Reynolds number, the analysis is based on large eddy simulations which include the effect of helicity on eddy viscosity and eddy noise. Thus, we model the small scales and resolve explicitly the large scales. We show that the large-scale energy spectrum has at least two solutions: one that is consistent with Kolmogorov-Kraichnan-Batchelor-Leith phenomenology for the inverse cascade of energy in two-dimensional (2D) turbulence with a ˜K-5/3 scaling, and the other that corresponds to a steeper ˜k-3 spectrum in which the three-dimensional (3D) modes release a substantial fraction of their energy per unit time to the 2D modes. The spectrum that emerges depends on the anisotropy of the forcing function, the former solution prevailing for forcings in which more energy is injected into the 2D modes while the latter prevails for isotropic forcing. In the case of anisotropic forcing, whence the energy goes from the 2D to the 3D modes at low wave numbers, large-scale shear is created, resulting in a time scale τsh, associated with shear, thereby producing a ˜ k-1 spectrum for the total energy awith the horizontal energy of the 2D modes still following ˜K-5/3 scaling. [ABSTRACT FROM AUTHOR]

We derive statistical equilibrium solutions of the truncated inviscid surface quasigeostrophic (SQG) equations, and verify the validity of these solutions at late times in numerical simulations. The results indicate the pseudoenstrophy thermalizes while the pseudoenergy can condense at the gravest modes, in agreement with previous indications of a direct cascade of pseudoenstrophy and an inverse cascade of pseudoenergy in forced-dissipative SQG systems. At early times, the truncated inviscid SQG simulations show a behavior reminiscent of forced-dissipative SQG turbulence, and we identify spectral scaling laws for the pseudoenergy and pseudoenstrophy spectra. More importantly, a comparison between viscous and inviscid simulations allows us to identify free-decay laws for the pseudoenstrophy in SQG turbulence at very large Reynolds number. [ABSTRACT FROM AUTHOR]

We study the stability of shear flows in a fully ionized plasma. Kelvin-Helmholtz is a well-known macroscopic and ideal shear-driven instability. In sufficiently low-density plasmas, also the microscopic Hall magnetoshear instability can take place. We performed three-dimensional simulations of the Hall-magnetohydrodynamic equations where these two instabilities are present, and carried out a comparative study. We find that when the shear flow is so intense that its vorticity surpasses the ion-cyclotron frequency of the plasma, the Hall magnetoshear instability is not only non-negligible, but it actually displays growth rates larger than those of the Kelvin-Helmholtz instability. [ABSTRACT FROM AUTHOR]

Understanding the in situ amplification of large-scale magnetic fields in turbulent astrophysical rotators has been a core subject of dynamo theory. When turbulent velocities are helical, large-scale dynamos that substantially amplify fields on scales that exceed the turbulent forcing scale arise, but the minimum sufficient fractional kinetic helicity fh,c has not been previously well quantified. Using direct numerical simulations for a simple helical dynamo, we show that fh,c decreases as the ratio of forcing to large-scale wave numbers kF/kmin increases. From the condition that a large-scale helical dynamo must overcome the back reaction from any nonhelical field on the large scales, we develop a theory that can explain the simulations. For kF/kmin &sce; 8 we find fh,c ≲ 3%, implying that very small helicity fractions strongly influence magnetic spectra for even moderate-scale separation. [ABSTRACT FROM AUTHOR]

Scalar transport by turbulent flows is best described in terms of Lagrangian parcel motions. Here we measure the Eulerian distance travel along Lagrangian trajectories in a simple point vortex flow to determine the probabilistic impulse response function for scalar transport in the absence of molecular diffusion. As expected, the mean squared Eulerian displacement scales ballistically at very short times and diffusively for very long times, with the displacement distribution at any given time approximating that of a random walk. However, significant deviations in the displacement distributions from Rayleigh are found. The probability of long distance transport is reduced over inertial range time scales due to spatial and temporal intermittency. This can be modeled as a series of trapping events with durations uniformly distributed below the Eulerian integral time scale. The probability of long distance transport is, on the other hand, enhanced beyond that of the random walk for both times shorter than the Lagrangian integral time and times longer than the Eulerian integral time. The very short-time enhancement reflects the underlying Lagrangian velocity distribution, while that at very long times results from the spatial and temporal variation of the flow at the largest scales. The probabilistic impulse response function, and with it the expectation value of the scalar concentration at any point in space and time, can be modeled using only the evolution of the lowest spatial wave number modes (the mean and the lowest harmonic) and an eddy based constrained random walk that captures the essential velocity phase relations associated with advection by vortex motions. Preliminary examination of Lagrangian tracers in three-dimensional homogeneous isotropic turbulence suggests that transport in that setting can be similarly modeled. [ABSTRACT FROM AUTHOR]

We present numerical evidence of a critical-like transition in an out-of-equilibrium mean-field description of a quantum system. By numerically solving the Gross-Pitaevskii equation we show that quantum turbulence displays an abrupt change between three-dimensional (3D) and two-dimensional (2D) behavior. The transition is observed both in quasi-2D flows in cubic domains (controlled by the amplitude of a 3D perturbation to the flow), as well as in flows in thin domains (controlled by the domain aspect ratio) in a configuration that mimics systems realized in laboratory experiments. In one regime the system displays a transfer of the energy towards smaller scales, while in the other the system displays a transfer of the energy towards larger scales and a coherent self-organization of the quantized vortices. [ABSTRACT FROM AUTHOR]

Understanding the in situ amplification of large-scale magnetic fields in turbulent astrophysical rotators has been a core subject of dynamo theory. When turbulent velocities are helical, large-scale dynamos that substantially amplify fields on scales that exceed the turbulent forcing scale arise, but the minimum sufficient fractional kinetic helicity f(h,C) has not been previously well quantified. Using direct numerical simulations for a simple helical dynamo, we show that f(h,C) decreases as the ratio of forcing to large-scale wave numbers k(F)/k(min) increases. From the condition that a large-scale helical dynamo must overcome the back reaction from any nonhelical field on the large scales, we develop a theory that can explain the simulations. For k(F)/k(min)≥8 we find f(h,C)≲3%, implying that very small helicity fractions strongly influence magnetic spectra for even moderate-scale separation.

Magnetic field generation by dynamo action is often studied within the theoretical framework of magnetohydrodynamics (MHD). However, for sufficiently diffuse media, the Hall effect may become non-negligible. We present results from three-dimensional simulations of the Hall-MHD equations subjected to random nonhelical forcing. We study the role of the Hall effect in the dynamo efficiency for different values of the Hall parameter. For small values of the Hall parameter, the small-scale dynamo is more efficient, displaying faster growth and saturating at larger amplitudes of the magnetic field. For larger values of the Hall parameter, saturation of the magnetic field is reached at smaller amplitudes than in the MHD case. We also study energy transfer rates among spatial scales and show that the Hall effect produces a reduction of the direct energy cascade at scales larger than the Hall scale, therefore leading to smaller energy dissipation rates. Finally, we present results stemming from simulations at large magnetic Prandtl numbers, which is the relevant regime in the hot and diffuse interstellar medium. In the range of magnetic Prandtl numbers considered, the Hall effect moves the peak of the magnetic energy spectrum as well as other relevant magnetic length scales toward the Hall scale.

We demonstrate that, for the case of quasiequipartition between the velocity and the magnetic field, the Lagrangian-averaged magnetohydrodynamics (LAMHD) alpha model reproduces well both the large-scale and the small-scale properties of turbulent flows; in particular, it displays no increased (superfilter) bottleneck effect with its ensuing enhanced energy spectrum at the onset of the subfilter scales. This is in contrast to the case of the neutral fluid in which the Lagrangian-averaged Navier-Stokes alpha model is somewhat limited in its applications because of the formation of spatial regions with no internal degrees of freedom and subsequent contamination of superfilter-scale spectral properties. We argue that, as the Lorentz force breaks the conservation of circulation and enables spectrally nonlocal energy transfer (associated with Alfvén waves), it is responsible for the absence of a viscous bottleneck in magnetohydrodynamics (MHD), as compared to the fluid case. As LAMHD preserves Alfvén waves and the circulation properties of MHD, there is also no (superfilter) bottleneck found in LAMHD, making this method capable of large reductions in required numerical degrees of freedom; specifically, we find a reduction factor of approximately 200 when compared to a direct numerical simulation on a large grid of 1536;{3} points at the same Reynolds number.

We compute solutions of the Lagrangian-averaged Navier-Stokes alpha - (LANS alpha ) model for significantly higher Reynolds numbers (up to Re approximately 8300 ) than have previously been accomplished. This allows sufficient separation of scales to observe a Navier-Stokes inertial range followed by a second inertial range specific to the LANS alpha model. Both fully helical and nonhelical flows are examined, up to Reynolds numbers of approximately 1300. Analysis of the third-order structure function scaling supports the predicted l3 scaling; it corresponds to a k-1 scaling of the energy spectrum for scales smaller than alpha. The energy spectrum itself shows a different scaling, which goes as k1. This latter spectrum is consistent with the absence of stretching in the subfilter scales due to the Taylor frozen-in hypothesis employed as a closure in the derivation of the LANS alpha model. These two scalings are conjectured to coexist in different spatial portions of the flow. The l3 [E(k) approximately k-1] scaling is subdominant to k1 in the energy spectrum, but the l3 scaling is responsible for the direct energy cascade, as no cascade can result from motions with no internal degrees of freedom. We demonstrate verification of the prediction for the size of the LANS alpha attractor resulting from this scaling. From this, we give a methodology either for arriving at grid-independent solutions for the LANS alpha model, or for obtaining a formulation of the large eddy simulation optimal in the context of the alpha models. The fully converged grid-independent LANS alpha model may not be the best approximation to a direct numerical simulation of the Navier-Stokes equations, since the minimum error is a balance between truncation errors and the approximation error due to using the LANS alpha instead of the primitive equations. Furthermore, the small-scale behavior of the LANS alpha model contributes to a reduction of flux at constant energy, leading to a shallower energy spectrum for large alpha. These small-scale features, however, do not preclude the LANS alpha model from reproducing correctly the intermittency properties of the high-Reynolds-number flow.

We present direct numerical simulations of dynamo action in a forced Roberts flow. The behavior of the dynamo is followed as the mechanical Reynolds number is increased, starting from the laminar case until a turbulent regime is reached. The critical magnetic Reynolds for dynamo action is found, and in the turbulent flow it is observed to be nearly independent on the magnetic Prandtl number in the range from approximately 0.3 to approximately 0.1. Also the dependence of this threshold with the amount of mechanical helicity in the flow is studied. For the different regimes found, the configuration of the magnetic and velocity fields in the saturated steady state are discussed.

We investigate the transfer of energy from large scales to small scales in fully developed forced three-dimensional magnetohydrodynamics (MHD) turbulence by analyzing the results of direct numerical simulations in the absence of an externally imposed uniform magnetic field. Our results show that the transfer of kinetic energy from large scales to kinetic energy at smaller scales and the transfer of magnetic energy from large scales to magnetic energy at smaller scales are local, as is also found in the case of neutral fluids and in a way that is compatible with the Kolmogorov theory of turbulence. However, the transfer of energy from the velocity field to the magnetic field is a highly nonlocal process in Fourier space. Energy from the velocity field at large scales can be transferred directly into small-scale magnetic fields without the participation of intermediate scales. Some implications of our results to MHD turbulence modeling are also discussed.

We present direct numerical simulations and Lagrangian averaged (also known as alpha model) simulations of forced and free decaying magnetohydrodynamic turbulence in two dimensions. The statistics of sign cancellations of the current at small scales is studied using both the cancellation exponent and the fractal dimension of the structures. The alpha model is found to have the same scaling behavior between positive and negative contributions as the direct numerical simulations. The alpha model is also able to reproduce the time evolution of these quantities in free decaying turbulence. At large Reynolds numbers, an independence of the cancellation exponent with the Reynolds numbers is observed.

We present direct numerical simulations and alpha -model simulations of four familiar three-dimensional magnetohydrodynamic (MHD) turbulence effects: selective decay, dynamic alignment, inverse cascade of magnetic helicity, and the helical dynamo effect. The MHD alpha model is shown to capture the long-wavelength spectra in all these problems, allowing for a significant reduction of computer time and memory at the same kinetic and magnetic Reynolds numbers. In the helical dynamo, not only does the alpha model correctly reproduce the growth rate of magnetic energy during the kinematic regime, it also captures the nonlinear saturation level and the late generation of a large scale magnetic field by the helical turbulence.