Spectral energy scaling in precessing turbulence.

We study precessing turbulence, which appears in several geophysical and astrophysical systems, by direct numerical simulations of homogeneous turbulence where precessional instability is triggered due to the imposed background flow. We show that the time development of kinetic energy K occurs in two main phases associated with different flow topologies: (i) an exponential growth characterizing three-dimensional turbulence dynamics and (ii) nonlinear saturation during which K remains almost time independent, the flow becoming quasi-two-dimensional. The latter stage, wherein the development of K remains insensitive to the initial state, shares an important common feature with other quasi-two-dimensional rotating flows such as rotating Rayleigh-Bénard convection, or the large atmospheric scales: in the plane k_{∥}=0, i.e., the plane associated to an infinite wavelength in the direction parallel to the principal rotation axis, the kinetic energy spectrum scales as k_{⊥}^{-3}. We show that this power law is observed for wave numbers ranging between the Zeman "precessional" and "rotational" scales, k_{S}^{-1} and k_{Ω}^{-1}, respectively, at which the associated background shear or inertial timescales are equal to the eddy turnover time. In addition, an inverse cascade develops for (k_{⊥},k)

Energy partition, scale by scale, in magnetic Archimedes Coriolis weak wave turbulence.

Magnetic Archimedes Coriolis (MAC) waves are omnipresent in several geophysical and astrophysical flows such as the solar tachocline. In the present study, we use linear spectral theory (LST) and investigate the energy partition, scale by scale, in MAC weak wave turbulence for a Boussinesq fluid. At the scale k^{-1}, the maximal frequencies of magnetic (Alfvén) waves, gravity (Archimedes) waves, and inertial (Coriolis) waves are, respectively, V_{A}k,N, and f. By using the induction potential scalar, which is a Lagrangian invariant for a diffusionless Boussinesq fluid [Salhi et al., Phys. Rev. E 85, 026301 (2012)PLEEE81539-375510.1103/PhysRevE.85.026301], we derive a dispersion relation for the three-dimensional MAC waves, generalizing previous ones including that of f-plane MHD "shallow water" waves [Schecter et al., Astrophys. J. 551, L185 (2001)AJLEEY0004-637X10.1086/320027]. A solution for the Fourier amplitude of perturbation fields (velocity, magnetic field, and density) is derived analytically considering a diffusive fluid for which both the magnetic and thermal Prandtl numbers are one. The radial spectrum of kinetic, S_{κ}(k,t), magnetic, S_{m}(k,t), and potential, S_{p}(k,t), energies is determined considering initial isotropic conditions. For magnetic Coriolis (MC) weak wave turbulence, it is shown that, at large scales such that V_{A}k/f≪1, the Alfvén ratio S_{κ}(k,t)/S_{m}(k,t) behaves like k^{-2} if the rotation axis is aligned with the magnetic field, in agreement with previous direct numerical simulations [Favier et al., Geophys. Astrophys. Fluid Dyn. (2012)] and like k^{-1} if the rotation axis is perpendicular to the magnetic field. At small scales, such that V_{A}k/f≫1, there is an equipartition of energy between magnetic and kinetic components. For magnetic Archimedes weak wave turbulence, it is demonstrated that, at large scales, such that (V_{A}k/N≪1), there is an equipartition of energy between magnetic and potential components, while at small scales (V_{A}k/N≫1), the ratio S_{p}(k,t)/S_{κ}(k,t) behaves like k^{-1} and S_{κ}(k,t)/S_{m}(k,t)=1. Also, for MAC weak wave turbulence, it is shown that, at small scales (V_{A}k/sqrt[N^{2}+f^{2}]≫1), the ratio S_{p}(k,t)/S_{κ}(t) behaves like k^{-1} and S_{κ}(k,t)/S_{m}(k,t)=1.

Cross-helicity in rotating homogeneous shear-stratified turbulence.

We consider homogeneous shear-stratified turbulence in a rotating frame, that exhibits complex nonlinear dynamics. Since the analysis of relative orientation between coupled fluctuating fields helps us to understand turbulence dynamics, we focus on the alignment properties of both the velocity and gravity fields with the potential vorticity gradient. With the help of statistical mechanics, we define a vector field which plays a role in the analogous so-called cross-helicity in magnetohydrodynamics. High-resolution direct numerical simulations of developed homogeneous baroclinic turbulence are performed, and a detailed analysis of probability density functions for cross-helicity is provided. A net preference for positive cross-helicity is shown to be related to a new alignment mechanism. We argue that the analysis of cross-helicity is crucial for understanding the dynamics of buoyancy driven flows.

Magnetized stratified rotating shear waves.

We present a spectral linear analysis in terms of advected Fourier modes to describe the behavior of a fluid submitted to four constraints: shear (with rate S), rotation (with angular velocity Ω), stratification, and magnetic field within the linear spectral theory or the shearing box model in astrophysics. As a consequence of the fact that the base flow must be a solution of the Euler-Boussinesq equations, only radial and/or vertical density gradients can be taken into account. Ertel's theorem no longer is valid to show the conservation of potential vorticity, in the presence of the Lorentz force, but a similar theorem can be applied to a potential magnetic induction: The scalar product of the density gradient by the magnetic field is a Lagrangian invariant for an inviscid and nondiffusive fluid. The linear system with a minimal number of solenoidal components, two for both velocity and magnetic disturbance fields, is eventually expressed as a four-component inhomogeneous linear differential system in which the buoyancy scalar is a combination of solenoidal components (variables) and the (constant) potential magnetic induction. We study the stability of such a system for both an infinite streamwise wavelength (k(1) = 0, axisymmetric disturbances) and a finite one (k(1) ≠ 0, nonaxisymmetric disturbances). In the former case (k(1) = 0), we recover and extend previous results characterizing the magnetorotational instability (MRI) for combined effects of radial and vertical magnetic fields and combined effects of radial and vertical density gradients. We derive an expression for the MRI growth rate in terms of the stratification strength, which indicates that purely radial stratification can inhibit the MRI instability, while purely vertical stratification cannot completely suppress the MRI instability. In the case of nonaxisymmetric disturbances (k(1) ≠ 0), we only consider the effect of vertical stratification, and we use Levinson's theorem to demonstrate the stability of the solution at infinite vertical wavelength (k(3) = 0): There is an oscillatory behavior for τ > 1+|K(2)/k(1)|, where τ = St is a dimensionless time and K(2) is the radial component of the wave vector at τ = 0. The model is suitable to describe instabilities leading to turbulence by the bypass mechanism that can be relevant for the analysis of magnetized stratified Keplerian disks with a purely azimuthal field. For initial isotropic conditions, the time evolution of the spectral density of total energy (kinetic + magnetic + potential) is considered. At k(3) = 0, the vertical motion is purely oscillatory, and the sum of the vertical (kinetic + magnetic) energy plus the potential energy does not evolve with time and remains equal to its initial value. The horizontal motion can induce a rapid transient growth provided K(2)/k(1)>>1. This rapid growth is due to the aperiodic velocity vortex mode that behaves like K(h)/k(h) where k(h)(τ)=[k(1)(2) + (K(2) - k(1)τ)(2)](1/2) and K(h) =k(h)(0). After the leading phase (τ > K(2)/k(1)>>1), the horizontal magnetic energy and the horizontal kinetic energy exhibit a similar (oscillatory) behavior yielding a high level of total energy. The contribution to energies coming from the modes k(1) = 0 and k(3) = 0 is addressed by investigating the one-dimensional spectra for an initial Gaussian dense spectrum. For a magnetized Keplerian disk with a purely vertical field, it is found that an important contribution to magnetic and kinetic energies comes from the region near k(1) = 0. The limit at k(1) = 0 of the streamwise one-dimensional spectra of energies, or equivalently, the streamwise two-dimensional (2D) energy, is then computed. The comparison of the ratios of these 2D quantities with their three-dimensional counterparts provided by previous direct numerical simulations shows a quantitative agreement.