Search

Nonlinear phenomena and turbulence are central to our understanding and modeling the dynamics of fluids and plasmas, and yet they still resist analytical resolutions in many instances. However, progress has been made recently, displaying a richness of phenomena which was somewhat unexpected a few years back, such as the double constant-flux cascades of a same invariant to both the large and to the small scales, or the presence of non-Gaussian wings in the large-scale fields, for fluids and plasmas. Here, I will concentrate on the direct measurement of the magnitude of dissipation and an evaluation of intermittency in a turbulent plasma using exact laws stemming from invariance principles and involving cross-correlation tensors with both the velocity and the magnetic fields. I will illustrate these points through scaling laws, together with data analysis from existing experiments, observations and numerical simulations. Finally, I will also briefly explore the possible implications for validity and use of several modeling strategies. To appear, Plasma Physics and Controlled Fusion, 2023.

Space plasmas are three-dimensional dynamic entities. Except under very special circumstances, their structure in space and their behavior in time are not related in any simple way. Therefore, single spacecraft in situ measurements cannot unambiguously unravel the full space-time structure of the heliospheric plasmas of interest in the inner heliosphere, in the Geospace environment, or the outer heliosphere. This shortcoming leaves numerous central questions incompletely answered. Deficiencies remain in at least two important subjects, Space Weather and fundamental plasma turbulence theory, due to a lack of a more complete understanding of the space-time structure of dynamic plasmas. Only with multispacecraft measurements over suitable spans of spatial separation and temporal duration can these ambiguities be resolved. We note that these characterizations apply to turbulence across a wide range of scales, and also equally well to shocks, flux ropes, magnetic clouds, current sheets, stream interactions, etc. In the following, we will describe the basic requirements for resolving space-time structure in general, using turbulence' as both an example and a principal target or study. Several types of missions are suggested to resolve space-time structure throughout the Heliosphere., Comment: White Paper submitted to: Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033. arXiv admin note: substantial text overlap with arXiv:1903.06890

We observe the emergence of strong vertical drafts in direct numerical simulations of the Boussinesq equations in a range of parameters of geophysical interest. These structures, which appear intermittently in space and time, generate turbulence and enhance kinetic and potential energy dissipation, providing a possible explanation for the observed variability of the local energy dissipation in the bulk of oceanic flows, and the modulation of its probability distribution function. We show how, due to the extreme drafts, in runs with Froude numbers observable in geophysical scenarios, roughly 10% of the domain flow can account for up to 50% of the global volume dissipation, reminiscent of estimates based on oceanic models., Comment: 16 pages, 6 figures

Non-Gaussian statistics of large-scale fields are routinely observed in data from atmospheric and oceanic campaigns and global models. Recent direct numerical simulations (DNSs) showed that large-scale intermittency in stably stratified flows is due to the emergence of sporadic, extreme events in the form of bursts in the vertical velocity and the temperature. This phenomenon results from the interplay between waves and turbulent motions, affecting mixing. We provide evidence of the enhancement of the classical small-scale (or internal) intermittency due to the emergence of large-scale drafts, connecting large- and small-scale bursts. To this aim we analyze a large set of DNSs of the stably stratified Boussinesq equations over a wide range of values of the Froude number ($Fr\approx 0.01-1$). The variation of the buoyancy field kurtosis with $Fr$ is similar to (though with smaller values than) the kurtosis of the vertical velocity, both showing a non-monotonic trend. We present a mechanism for the generation of extreme vertical drafts and vorticity enhancements which follows from the exact equations for field gradients., Comment: 7 pages, 4 figures

Helicity, a measure of the breakage of reflectional symmetry representing the topology of turbulent flows, contributes in a crucial way to their dynamics and to their fundamental statistical properties. We review several of their main features, both new and old, as the discovery of bi-directional cascades or the role of helical vortices in the enhancement of large-scale magnetic fields in the dynamo problem. The dynamical contribution in magnetohydrodynamic (MHD) of the cross-correlation between velocity and induction is discussed as well. We consider next how turbulent transport is affected by helical constraints, in particular in the context of magnetic reconnection and fusion plasmas {under one- and two-fluid approximations}. Central issues on how to construct turbulence models for non-reflectionally symmetric helical flows are reviewed, including in the presence of shear, and we finally briefly mention the possible role of helicity in the development of strongly localized quasi-singular structures {at small scale)., Comment: 14 pages, 0 figure 0 table

In turbulence, for neutral or conducting fluids, a large ratio of scales is excited because of the possible occurrence of inverse cascades to large, global scales together with direct cascades to small, dissipative scales, as observed in the atmosphere and oceans, or in the solar environment. In this context, using direct numerical simulations with forcing, we analyze scale dynamics in the presence of magnetic fields with a generalized Ohm's law including a Hall current. The ion inertial length epsilon_H serves as the control parameter at fixed Reynolds number. Both the magnetic and generalized helicity -- invariants in the ideal case -- grow linearly with time, as expected from classical arguments. The cross-correlation between the velocity and magnetic field grows as well, more so in relative terms for a stronger Hall current. We find that the helical growth rates vary exponentially with epsilon_H, provided the ion inertial scale resides within the inverse cascade range. These exponential variations are recovered phenomenologically using simple scaling arguments. They are directly linked to the wavenumber power-law dependence of generalized and magnetic helicity, k^(-2), in their inverse ranges. This illustrates and confirms the important role of the interplay between large and small scales in the dynamics of turbulent flows.

We study the spatio-temporal behavior of the Els\"asser variables describing magnetic and velocity field fluctuations, using direct numerical simulations of three-dimensional magnetohydrodynamic turbulence. We consider cases with relatively small, intermediate, and large values of a mean background magnetic field, and with null, small, and high cross-helicity (correlations between the velocity and the magnetic field). Wavenumber-dependent time correlation functions are computed for the different simulations. From these correlation functions, the decorrelation time is computed and compared with different theoretical characteristic times: the local non-linear time, the random-sweeping time, and the Alfv\'enic time. It is found that decorrelation times are dominated by sweeping effects for low values of the mean magnetic field and for low values of the cross-helicity, while for large values of the background field or of the cross-helicity and for wave vectors sufficiently aligned with the guide field, decorrelation times are controlled by Alfv\'enic effects. Finally, we observe counter-propagation of Alfv\'enic fluctuations due to reflections produced by inhomogeneities in the total magnetic field. This effect becomes more prominent in flows with large cross-helicity, strongly modifying the propagation of waves in turbulent magnetohydrodynamic flows., Comment: 16 pages, 40 figures

Geophysical fluid flows are predominantly turbulent and often strongly affected by the Earth's rotation, as well as by stable density stratification. Using direct numerical simulations of forced Boussinesq equations, we study the influence of these effects on the motion of fluid particles, focusing on cases where the frequencies associated with rotation and stratification (RaS), $N$ and $f$ respectively, are held at a fixed ratio $N/f=5$. As the intensity of RaS increases, a sharp transition is observed between a regime dominated by eddies to a regime dominated by waves, which can also be seemingly described by simply comparing the time scale $1/N$ and $\tau_\eta$ (the Kolmogorov time scale). We perform a detailed study of Lagrangian statistics of acceleration, velocity and related quantities in the two regimes. The flow anisotropy induces a clear difference between particle motion in the horizontal and vertical directions. In the regime $N\tau_\eta <1$, acceleration statistics in both horizontal and vertical directions, exhibit well known characteristics of isotropic turbulence. In contrast for $N\tau_\eta >1$, they are directly influenced by imposed RaS. The Lagrangian velocity statistics exhibit visible anisotropy for all runs; nevertheless the degree of anisotropy becomes very strong in the regime $N\tau_\eta >1$. We find that in the regime $N\tau_\eta <1$, rotation enhances the mean displacement of particles in horizontal planes at short times, but inhibits them at longer times. This inhibition of horizontal displacement becomes stronger for $N\tau_\eta >1$, with no clear diffusive behavior. Displacements in the vertical direction are always inhibited. The inhibition becomes extremely strong when $N\tau_\eta >1$, with the particles almost being trapped horizontally., Comment: 24 pages, 10 figures

Analyzing a large data base of high-resolution three-dimensional direct numerical simulations of decaying rotating stratified flows, we show that anomalous mixing and dissipation, marked anisotropy, and strong intermittency are all observed simultaneously in an intermediate regime of parameters in which both waves and eddies interact nonlinearly. A critical behavior governed by the stratification occurs at Richardson numbers of order unity, and with the flow close to being in a state of instability. This confirms the central dynamical role, in rotating stratified turbulence, of large-scale intermittency, which occurs in the vertical velocity and temperature fluctuations, as an adjustment mechanism of the energy transfer in the presence of strong waves.

This paper briefly reviews a number of fundamental measurements that need to be made in order to characterize turbulence in space plasmas such as the solar wind. It has long been known that many of these quantities require simultaneous multipoint measurements to attain a proper characterization that would reveal the fundamental physics of plasma turbulence. The solar wind is an ideal plasma for such an investigation, and it now appears to be technologically feasible to carry out such an investigation, following the pioneering Cluster and MMS missions. Quantities that need to be measured using multipoint measurements include the two-point, two-time second correlation function of velocity, magnetic field and density, and higher order statistical objects such as third and fourth order structure functions. Some details of these requirements are given here, with a eye towards achieving closure on fundamental questions regarding the cascade rate, spectral anisotropy, characteristic coherent structures, intermittency, and dissipation mechanisms that describe plasma turbuelence, as well as its variability with plasma parameters in the solar wind. The motivation for this discussion is the current planning for a proposed Helioswarm mission that would be designed to make these measurements,leading to breakthrough understanding of the physics of space and astrophysical turbulence., Comment: White paper submitted to the PLASMA 2020 Decadal Survey Committee

An existing hybrid MPI-OpenMP scheme is augmented with a CUDA-based fine grain parallelization approach for multidimensional distributed Fourier transforms, in a well-characterized pseudospectral fluid turbulence code. Basics of the hybrid scheme are reviewed, and heuristics provided to show a potential benefit of the CUDA implementation. The method draws heavily on the CUDA runtime library to handle memory management, and on the cuFFT library for computing local FFTs. The manner in which the interfaces are constructed to these libraries, and ISO bindings utilized to facilitate platform portability, are discussed. CUDA streams are implemented to overlap data transfer with cuFFT computation. Testing with a baseline solver demonstrates significant aggregate speed-up over the hybrid MPI-OpenMP solver by offloading to GPUs on an NVLink-based test system. While the batch streamed approach provides little benefit with NVLink, we see a performance gain of 30% when tuned for the optimal number of streams on a PCIe-based system. It is found that strong GPU scaling is ideal, or slightly better than ideal, in all cases. In addition to speed-up measurements for the fiducial solver, we also consider several other solvers with different numbers of transform operations and find that aggregate speed-ups are nearly constant for all solvers., Comment: 17 pages, 8 figures; submitted to Parallel Computing

We briefly review helicity dynamics, inverse and bi-directional cascades in fluid and magnetohydrodynamic (MHD) turbulence, with an emphasis on the latter. The energy of a turbulent system, an invariant in the non-dissipative case, is transferred to small scales through nonlinear mode coupling. Fifty years ago, it was realized that, for a two-dimensional fluid, energy cascades instead to larger scales, and so does magnetic helicity in three-dimensional MHD. However, evidence obtained recently indicates that in fact, for a range of governing parameters, there are systems for which their ideal invariants can be transferred, with constant fluxes, to both the large scales and the small scales, as for MHD or rotating stratified flows, in the latter case including with quasi-geostrophic forcing. Such bi-directional, split, cascades directly affect the rate at which mixing and dissipation occur in these flows in which nonlinear eddies interact with fast waves with anisotropic dispersion laws, due for example to imposed rotation, stratification or uniform magnetic fields. The directions of cascades can be obtained in some cases through the use of phenomenological arguments, one of which we derive here following classical lines in the case of the inverse magnetic helicity cascade in electron MHD. With more highly-resolved data sets stemming from large laboratory experiments, high-performance computing and in-situ satellite observations, machine-learning tools are bringing novel perspectives to turbulence research, e.g. in helping devise new explicit sub-grid scale parameterizations, which may lead to enhanced physical insight, including in the future in the case of these new bi-directional cascades.

We investigate the large-scale intermittency of vertical velocity and temperature, and the mixing properties of stably stratified turbulent flows using both Lagrangian and Eulerian fields from direct numerical simulations, in a parameter space relevant for the atmosphere and the oceans. Over a range of Froude numbers of geophysical interest ($\approx 0.05-0.3$) we observe very large fluctuations of the vertical components of the velocity and the potential temperature, localized in space and time, with a sharp transition leading to non-Gaussian wings of the probability distribution functions. This behavior is captured by a simple model representing the competition between gravity waves on a fast time-scale and nonlinear steepening on a slower time-scale. The existence of a resonant regime characterized by enhanced large-scale intermittency, as understood within the framework of the proposed model, is then linked to the emergence of structures in the velocity and potential temperature fields, localized overturning and mixing. Finally, in the same regime we observe a linear scaling of the mixing efficiency with the Froude number and an increase of its value of roughly one order of magnitude.

We present a model for the scaling of mixing in weakly rotating stratified flows characterized by their Rossby, Froude and Reynolds numbers Ro, Fr, Re. It is based on quasi-equipartition between kinetic and potential modes, sub-dominant vertical velocity and lessening of the energy transfer to small scales as measured by the ratio rE of kinetic energy dissipation to its dimensional expression. We determine their domains of validity for a numerical study of the unforced Boussinesq equations mostly on grids of 10243 points, with Ro/Fr> 2.5 and with 1600< Re<1.9x104; the Prandtl number is one, initial conditions are either isotropic and at large scale for the velocity, and zero for the temperature {\theta}, or in geostrophic balance. Three regimes in Fr are observed: dominant waves, eddy-wave interactions and strong turbulence. A wave-turbulence balance for the transfer time leads to rE growing linearly with Fr in the intermediate regime, with a saturation at ~0.3 or more, depending on initial conditions for larger Froude numbers. The Ellison scale is also found to scale linearly with Fr, and the flux Richardson number Rf transitions for roughly the same parameter values as well. Putting together the 3 relationships of the model allows for the prediction of mixing efficiency scaling as Fr-2~RB-1 in the low and intermediate regimes, whereas for higher Fr, it scales as RB-1/2, as already observed: as turbulence strengthens, rE~1, the velocity is isotropic and smaller buoyancy fluxes altogether correspond to a decoupling of velocity and temperature fluctuations, the latter becoming passive.

In this paper, we present an overview of concepts and data concerning inverse cascades of excitation towards scales larger than the forcing scale in a variety of contexts, from two-dimensional fluids and wave turbulence, to geophysical flows in the presence of rotation and stratification. We briefly discuss the role of anisotropy in the occurrence and properties of such cascades. We then show that the cascade of some invariant, for example the total energy, may be transferred through nonlinear interactions to both the small scales and the large scales, with in each case a constant flux. This is in contrast to the classical picture, and we illustrate such a dual cascade in the context of atmospheric and oceanic observations, direct numerical simulations and modeling. We also show that this dual cascade of total energy can in fact be decomposed in some cases into separate cascades of the kinetic and potential energies, provided the Froude and Rossby numbers are small enough. In all cases, the potential energy flux remains small, of the order of 10% or less relative to the kinetic energy flux. Finally, we demonstrate that, in the small-scale inertial range, approximate equipartition between potential and kinetic modes is obtained, leading to an energy ratio close to one, with strong departure at large scales due to the dominant kinetic energy inverse cascade and piling-up at the lowest spatial frequency, and at small scales due to unbalanced dissipation processes, even though the Prandtl number is equal to one.

The large-scale structures in the ocean and the atmosphere are in geostrophic balance, and a conduit must be found to channel the energy to the small scales where it can be dissipated. In turbulence this takes the form of an energy cascade, whereas one possible mechanism in a balanced flow at large scales is through the formation of fronts, a common occurrence in geophysical dynamics. We show in this paper that an iconic configuration in laboratory and numerical experiments for the study of turbulence, that of the Taylor-Green or von K\'arm\'an swirling flow, can be suitably adapted to the case of fluids with large aspect ratios, leading to the creation of an imposed large-scale vertical shear. To this effect we use direct numerical simulations of the Boussinesq equations without net rotation and with no small-scale modeling, and with this idealized Taylor-Green set-up. Various grid spacings are used, up to $2048^2\times 256$ spatial points. The grids are always isotropic, with box aspect ratios of either $1:4$ or $1:8$. We find that when shear and stratification are comparable, the imposed shear layer resulting from the forcing leads to the formation of multiple fronts and filaments which destabilize and further evolve into a turbulent flow in the bulk, with a sizable amount of dissipation and mixing, and with a cycle of front creation, instability, and development of turbulence. The results depend on the vertical length scales for shear and for stratification, with stronger large-scale gradients being generated when the two length scales are comparable., Comment: 19 pages, 15 figures, several simulations added in this new version

Kraichnan seminal ideas on inverse cascades yielded new tools to study common phenomena in geophysical turbulent flows. In the atmosphere and the oceans, rotation and stratification result in a flow that can be approximated as two-dimensional at very large scales, but which requires considering three-dimensional effects to fully describe turbulent transport processes and non-linear phenomena. Motions can thus be classified into two classes: fast modes consisting of inertia-gravity waves, and slow quasi-geostrophic modes for which the Coriolis force and horizontal pressure gradients are close to balance. In this paper we review previous results on the strength of the inverse cascade in rotating and stratified flows, and then present new results on the effect of varying the strength of rotation and stratification (measured by the ratio $N/f$ of the Brunt-V\"ais\"ala frequency to the Coriolis frequency) on the amplitude of the waves and on the flow quasi-geostrophic behavior. We show that the inverse cascade is more efficient in the range of $N/f$ for which resonant triads do not exist, $1/2 \le N/f \le 2$. We then use the spatio-temporal spectrum, and characterization of the flow temporal and spatial scales, to show that in this range slow modes dominate the dynamics, while the strength of the waves (and their relevance in the flow dynamics) is weaker., Comment: 17 pages, 13 figures, new version has all changes in the paper published in the Two Dimensional Turbulence Focus Issue of Physics of Fluids

The effects of changing the orientation of the rotation axis on homogeneous turbulence is considered. We perform direct numerical simulations on a periodic box of $1024^3$ grid points, where the orientation of the rotation axis is changed (a) at a fixed time instant (b) regularly at time intervals commensurate with the rotation time scale. The former is characterized by a dominant inverse energy cascade whereas in the latter, the inverse cascade is stymied due to the recurrent changes in the rotation axis resulting in a strong forward energy transfer and large scale structures that resemble those of isotropic turbulence., Comment: 7 pages, 8 figures, The European Physical Journal E (EPJ E)

The interplay between waves and eddies in stably stratified rotating flows is investigated by means of world-class direct numerical simulations using up to $3072^3$ grid points. Strikingly, we find that the shift from vortex to wave dominated dynamics occurs at a wavenumber $k_R$ which does not depend on Reynolds number, suggesting that partition of energy between wave and vortical modes is not sensitive to the development of turbulence at the smaller scales. We also show that $k_R$ is comparable to the wavenumber at which exchanges between kinetic and potential modes stabilize at close to equipartition, emphasizing the role of potential energy, as conjectured in the atmosphere and the oceans. Moreover, $k_R$ varies as the inverse of the Froude number as explained by the scaling prediction proposed, consistent with recent observations and modeling of the Mesosphere-Lower Thermosphere and of the ocean.

We study the partition of energy between waves and vortices in stratified turbulence, with or without rotation, for a variety of parameters, focusing on the behavior of the waves and vortices in the inverse cascade of energy towards the large scales. To this end, we use direct numerical simulations in a cubic box at a Reynolds number Re=1000, with the ratio between the Brunt-V\"ais\"al\"a frequency N and the inertial frequency f varying from 1/4 to 20, together with a purely stratified run. The Froude number, measuring the strength of the stratification, varies within the range 0.02 < Fr < 0.32. We find that the inverse cascade is dominated by the slow quasi-geostrophic modes. Their energy spectra and fluxes exhibit characteristics of an inverse cascade, even though their energy is not conserved. Surprisingly, the slow vortices still dominate when the ratio N/f increases, also in the stratified case, although less and less so. However, when N/f increases, the inverse cascade of the slow modes becomes weaker and weaker, and it vanishes in the purely stratified case. We discuss how the disappearance of the inverse cascade of energy with increasing N/f can be interpreted in terms of the waves and vortices, and identify three major effects that can explain this transition based on inviscid invariants arguments.

We build on recent developments in the study of fluid turbulence [Gibbon \textit{et al.} Nonlinearity 27, 2605 (2014)] to define suitably scaled, order-$m$ moments, $D_m^{\pm}$, of $\omega^\pm= \omega \pm j$, where $\omega$ and $j$ are, respectively, the vorticity and current density in three-dimensional magnetohydrodynamics (MHD). We show by mathematical analysis, for unit magnetic Prandtl number $P_M$, how these moments can be used to identify three possible regimes for solutions of the MHD equations; these regimes are specified by inequalities for $D_m^{\pm}$ and $D_1^{\pm}$. We then compare our mathematical results with those from our direct numerical simulations (DNSs) and thus demonstrate that 3D MHD turbulence is like its fluid-turbulence counterpart insofar as all solutions, which we have investigated, remain in \textit{only one of these regimes}; this regime has depleted nonlinearity. We examine the implications of our results for the exponents $q^{\pm}$ that characterize the power-law dependences of the energy spectra $\mathcal{E}^{\pm}(k)$ on the wave number $k$, in the inertial range of scales. We also comment on (a) the generalization of our results to the case $P_M \neq 1$ and (b) the relation between $D_m^{\pm}$ and the order-$m$ moments of gradients of hydrodynamic fields, which are used in characterizing intermittency in turbulent flows., Comment: 14 pages, 3 figures

We live in an age in which high-performance computing is transforming the way we do science. Previously intractable problems are now becoming accessible by means of increasingly realistic numerical simulations. One of the most enduring and most challenging of these problems is turbulence. Yet, despite these advances, the extreme parameter regimes encountered in space physics and astrophysics (as in atmospheric and oceanic physics) still preclude direct numerical simulation. Numerical models must take a Large Eddy Simulation (LES) approach, explicitly computing only a fraction of the active dynamical scales. The success of such an approach hinges on how well the model can represent the subgrid-scales (SGS) that are not explicitly resolved. In addition to the parameter regime, heliophysical and astrophysical applications must also face an equally daunting challenge: magnetism. The presence of magnetic fields in a turbulent, electrically conducting fluid flow can dramatically alter the coupling between large and small scales, with potentially profound implications for LES/SGS modeling. In this review article, we summarize the state of the art in LES modeling of turbulent magnetohydrodynamic (MHD) flows. After discussing the nature of MHD turbulence and the small-scale processes that give rise to energy dissipation, plasma heating, and magnetic reconnection, we consider how these processes may best be captured within an LES/SGS framework. We then consider several specific applications in heliophysics and astrophysics, assessing triumphs, challenges, and future directions., Comment: 51 pages, 6 figures (Figs 2, 3, and 4 color), accepted to Space Science Reviews (in press). The paper is a product of a workshop on "LES Modeling in MHD Turbulence" held in Boulder, CO in May, 2013, sponsored by the Geophysical Turbulence Program at the National Center for Atmospheric Research

We perform two high resolution direct numerical simulations of stratified turbulence for Reynolds number equal to Re~25000 and Froude number respectively of Fr~0.1 and Fr~0.03. The flows are forced at large scale and discretized on an isotropic grid of 2048^3 points. Stratification makes the flow anisotropic and introduces two extra characteristic scales with respect to homogeneous isotropic turbulence: the buoyancy scale, L_B, and the Ozmidov scale, l_{oz}. The former is related to the number of layers that the flow develops in the direction of gravity, the latter is regarded as the scale at which isotropy is recovered. The values of L_B and l_{oz} depend on the Froude number and their absolute and relative size affect the repartition of energy among Fourier modes. By contrasting the behavior of the two simulated flows we identify some surprising similarities: after an initial transient the two flows evolve towards comparable values of the kinetic and potential enstrophy, and energy dissipation rate. Further similarities emerge at large scales: the same ratio between potential and total energy (~0.1) is spontaneously selected by the flows, and slow modes grow monotonically in both regimes causing a slow increase of the total energy in time. The axisymmetric total energy spectrum shows a wide variety of spectral slopes as a function of the angle between the imposed stratification and the wave vector. One-dimensional energy spectra computed in the direction parallel to gravity are flat from the forcing up to buoyancy scale. At intermediate scales a ~ k^{-3} parallel spectrum develops for the Fr ~ 0.03 run, whereas for weaker stratification, the saturation spectrum does not have enough scales to develop and instead one observes a power law compatible with Kolmogorov scaling. Finally, the spectrum of helicity is flat until L_B, as observed in the nocturnal planetary boundary layer.

Employing a simple ideal magnetohydrodynamic model in spherical geometry,we show that the presence of either rotation or finite magnetic helicity is sufficient to induce dynamical reversals of the magnetic dipole moment. The statistical character of the model is similar to that of terrestrial magnetic field reversals, with the similarity being stronger when rotation is present.The connection between long time correlations, $1/f$ noise, and statistics of reversals is supported, consistent with earlier suggestions., Comment: accepted in Physical Review E

We report results on rotating stratified turbulence in the absence of forcing, with large-scale isotropic initial conditions, using direct numerical simulations computed on grids of up to 4096^3 points. The Reynolds and Froude numbers are respectively equal to Re=5.4 x 10^4 and Fr=0.0242. The ratio of the Brunt-V\"ais\"al\"a to the inertial wave frequency, N/f, is taken to be equal to 4.95, a choice appropriate to model the dynamics of the southern abyssal ocean at mid latitudes. This gives a global buoyancy Reynolds number R_B=ReFr^2=32, a value sufficient for some isotropy to be recovered in the small scales beyond the Ozmidov scale, but still moderate enough that the intermediate scales where waves are prevalent are well resolved. We concentrate on the large-scale dynamics, for which we find a spectrum compatible with the Bolgiano-Obukhov scaling, and confirm that the Froude number based on a typical vertical length scale is of order unity, with strong gradients in the vertical. Two characteristic scales emerge from this computation, and are identified from sharp variations in the spectral distribution of either total energy or helicity. A spectral break is also observed at a scale at which the partition of energy between the kinetic and potential modes changes abruptly, and beyond which a Kolmogorov-like spectrum recovers. Large slanted layers are ubiquitous in the flow in the velocity and temperature fields, with local overturning events indicated by small Richardson numbers, and a small large-scale enhancement of energy directly attributable to the effect of rotation is also observed., Comment: 19 pages, 9 figures (including compound figures)

We present results from direct numerical simulations of the Boussinesq equations in the presence of rotation and/or stratification, both in the vertical direction. The runs are forced isotropically and randomly at small scales and have spatial resolutions of up to $1024^3$ grid points and Reynolds numbers of $\approx 1000$. We first show that solutions with negative energy flux and inverse cascades develop in rotating turbulence, whether or not stratification is present. However, the purely stratified case is characterized instead by an early-time, highly anisotropic transfer to large scales with almost zero net isotropic energy flux. This is consistent with previous studies that observed the development of vertically sheared horizontal winds, although only at substantially later times. However, and unlike previous works, when sufficient scale separation is allowed between the forcing scale and the domain size, the total energy displays a perpendicular (horizontal) spectrum with power law behavior compatible with $\sim k_\perp^{-5/3}$, including in the absence of rotation. In this latter purely stratified case, such a spectrum is the result of a direct cascade of the energy contained in the large-scale horizontal wind, as is evidenced by a strong positive flux of energy in the parallel direction at all scales including the largest resolved scales.

Most of the turbulent flows appearing in nature (e.g. geophysical and astrophysical flows) are subjected to strong rotation and stratification. These effects break the symmetries of classical, homogenous isotropic turbulence. In doing so, they introduce a natural decomposition of phase space in terms of wave modes and potential vorticity modes. The appearance of a new time scale associated to the propagation of waves, in addition to the eddy turnover time, increases the complexity of the energy transfers between the various scales; nonlinearly interacting waves may dominate at some scales while balanced motion may prevail at others. In the end, it is difficult to predict \emph{a priori} if the energy cascades downscale as in homogeneous isotropic turbulence, upscale as expected from balanced dynamics, or follows yet another phenomenology. In this paper, we suggest a theoretical approach based on equilibrium statistical mechanics for the ideal system, inspired from the restricted partition function formalism introduced in metastability studies. In this framework, we show analytically that in the presence of rotation, when the dynamics is restricted to the slow modes, the equilibrium energy spectrum features an infrared divergence characteristic of an inverse cascade regime, whereas this is not the case for purely stratified flows.

The ocean and the atmosphere, and hence the climate, are governed at large scale by interactions between pressure gradient, Coriolis and buoyancy forces. This leads to a quasi-geostrophic balance in which, in a two-dimensional-like fashion, the energy injected by solar radiation, winds or tides goes to large scales in what is known as an inverse cascade. Yet, except for Ekman friction, energy dissipation and turbulent mixing occur at small scale implying the formation of such scales associated with breaking of geostrophic dynamics through wave-eddy interactions \cite{ledwell_00, vanneste_13} or frontogenesis \cite{hoskins_72, mcwilliams_10}, in opposition to the inverse cascade. Can it be both at the same time? We exemplify here this dual behavior of energy with the help of three-dimensional direct numerical simulations of rotating stratified Boussinesq turbulence. We show that efficient small-scale mixing and large-scale coherence develop simultaneously in such geophysical and astrophysical flows, both with constant flux as required by theoretical arguments, thereby clearly resolving the aforementioned contradiction., Comment: 5 pages, 2 figures

There is a clear distinction between simple laminar and complex turbulent fluids. But in some cases, as for the nocturnal planetary boundary layer, a stable and well-ordered flow can develop intense and sporadic bursts of turbulent activity which disappear slowly in time. This phenomenon is ill-understood and poorly modeled; and yet, it is central to our understanding of weather and climate dynamics. We present here a simple model which shows that in stably stratified turbulence, the stronger bursts can occur when the flow is expected to be more stable. The bursts are generated by a rapid non-linear amplification of energy stored in waves, and are associated with energetic interchanges between vertical velocity and temperature (or density) fluctuations. Direct numerical simulations on grids of 2048^3 points confirm this somewhat paradoxical result of measurably stronger events for more stable flows, displayed not only in the temperature and vertical velocity derivatives, but also in the amplitude of the fields themselves.

We present results from an ensemble of 50 runs of two-dimensional hydrodynamic turbulence with spatial resolution of 2048^2 grid points, and from an ensemble of 10 runs with 4096^2 grid points. All runs in each ensemble have random initial conditions with same initial integral scale, energy, enstrophy, and Reynolds number. When both ensemble- and time-averaged, inverse energy cascade behavior is observed, even in the absence of external mechanical forcing: the energy spectrum at scales larger than the characteristic scale of the flow follows a k^(-5/3) law, with negative flux, together with a k^(-3) law at smaller scales, and a positive flux of enstrophy. The source of energy for this behavior comes from the modal energy around the energy containing scale at t=0. The results shed some light into connections between decaying and forced turbulence, and into recent controversies in experimental studies of two-dimensional and magnetohydrodynamic turbulent flows., Comment: 7 pages, 6 figures

We examine the scaling laws of MHD turbulence for three different types of forcing functions and imposing at all times the four-fold symmetries of the Taylor-Green (TG) vortex generalized to MHD; no uniform magnetic field is present and the magnetic Prandtl number is equal to unity. We also include a forcing in the induction equation, and we take the three configurations studied in the decaying case in [E. Lee et al. Phys. Rev.E {\bf 81}, 016318 (2010)]. To that effect, we employ direct numerical simulations up to an equivalent resolution of $2048^3$ grid points. We find that, similarly to the case when the forcing is absent, different spectral indices for the total energy spectrum emerge, corresponding to either a Kolmogorov law, an Iroshnikov-Kraichnan law that arises from the interactions of turbulent eddies and Alfv\'en waves, or to weak turbulence when the large-scale magnetic field is strong. We also examine the inertial range dynamics in terms of the ratios of kinetic to magnetic energy, and of the turn-over time to the Alfv\'en time, and analyze the temporal variations of these quasi-equilibria.

A numerical study of decaying stably-stratified flows is performed. Relatively high stratification and moderate Reynolds numbers are considered, and a particular emphasis is placed on the role of helicity (velocity-vorticity correlations). The problem is tackled by integrating the Boussinesq equations in a periodic cubical domain using different initial conditions: a non-helical Taylor-Green (TG) flow, a fully helical Beltrami (ABC) flow, and random flows with a tunable helicity. We show that for stratified ABC flows helicity undergoes a substantially slower decay than for unstratified ABC flows. This fact is likely associated to the combined effect of stratification and large scale coherent structures. Indeed, when the latter are missing, as in random flows, helicity is rapidly destroyed by the onset of gravitational waves. A type of large-scale dissipative "cyclostrophic" balance can be invoked to explain this behavior. When helicity survives in the system it strongly affects the temporal energy decay and the energy distribution among Fourier modes. We discover in fact that the decay rate of energy for stratified helical flows is much slower than for stratified non-helical flows and can be considered with a phenomenological model in a way similar to what is done for unstratified rotating flows. We also show that helicity, when strong, has a measurable effect on the Fourier spectra, in particular at scales larger than the buoyancy scale for which it displays a rather flat scaling associated with vertical shear.

We perform numerical simulations of decaying rotating stratified turbulence and show, in the Boussinesq framework, that helicity (velocity-vorticity correlation), as observed in super-cell 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 quasi-linear balance equation. Outside this regime, and up to the highest Reynolds number obtained in this study, namely $Re\approx 10000$, helicity production is found to be persistent for $N/f$ as large as $\approx 17$, and for $ReFr^2$ and $ReRo^2 $ respectively as large as $\approx 100$ and $\approx 24000$., Comment: 10 pages, 5 figures

A review of some of the issues that have arisen over the years concerning the energy distribution among scales for magnetohydrodynamics (MHD) turbulence is given here. A variety of tools are employed to that effect, and a central role is played by taking into consideration the ideal (non-dissipative) invariants, namely the total energy, the magnetic helicity and the cross-correlations between the velocity and the magnetic field (concentrating on the three-dimensional case). These concepts, based mostly on theory, models and direct numerical simulations, are briefly put in the context of observations, in particular the solar wind, and some of the remaining open questions are delineated as well. New results on ideal MHD dynamics in three dimensions on equivalent grids of up to $6144^3$ points using the Taylor-Green flow generalized to MHD are also mentioned., Comment: 17 pages, 4 figures, review

We investigate the ideal and incompressible magnetohydrodynamic (MHD) equations in three space dimensions for the development of potentially singular structures. The methodology consists in implementing the four-fold symmetries of the Taylor-Green vortex generalized to MHD, leading to substantial computer time and memory savings at a given resolution; we also use a re-gridding method that allows for lower-resolution runs at early times, with no loss of spectral accuracy. One magnetic configuration is examined at an equivalent resolution of $6144^3$ points, and three different configurations on grids of $4096^3$ points. At the highest resolution, two different current and vorticity sheet systems are found to collide, producing two successive accelerations in the development of small scales. At the latest time, a convergence of magnetic field lines to the location of maximum current is probably leading locally to a strong bending and directional variability of such lines. A novel analytical method, based on sharp analysis inequalities, is used to assess the validity of the finite-time singularity scenario. This method allows one to rule out spurious singularities by evaluating the rate at which the logarithmic decrement of the analyticity-strip method goes to zero. The result is that the finite-time singularity scenario cannot be ruled out, and the singularity time could be somewhere between $t=2.33$ and $t=2.70.$ More robust conclusions will require higher resolution runs and grid-point interpolation measurements of maximum current and vorticity., Comment: 18 pages, 13 figures, 2 tables; submitted to Physical Review E

We investigate long-time properties of three-dimensional MHD turbulence in the absence of forcing and examine in particular the role played by the quadratic invariants of the system and by the symmetries of the initial configurations. We observe that, when sufficient accuracy is used, initial conditions with a high degree of symmetries, as in the absence of helicity, do not travel through parameter space over time whereas by perturbing these solutions either explicitly or implicitly using for example single precision for long times, the flows depart from their original behavior and can become either strongly helical, or have a strong alignment between the velocity and the magnetic field. When the symmetries are broken, the flows evolve towards different end states, as predicted by statistical arguments for non-dissipative systems with the addition of an energy minimization principle, as already analyzed in \cite{stribling_90} for random initial conditions using a moderate number of Fourier modes. Furthermore, the alignment properties of these flows, between velocity, vorticity, magnetic potential, induction and current, correspond to the dominance of two main regimes, one helically dominated and one in quasi-equipartition of kinetic and magnetic energy. We also contrast the scaling of the ratio of magnetic energy to kinetic energy as a function of wavenumber to the ratio of eddy turn-over time to Alfv\'en time as a function of wavenumber. We find that the former ratio is constant with an approximate equipartition for scales smaller than the largest scale of the flow whereas the ratio of time scales increases with increasing wavenumber., Comment: 14 pages, 6 figures

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 $L_f$, it cascades towards smaller as well as larger scales. In this paper we analyze the flow in the \textit{inverse cascade} range at a small but fixed Rossby number, {$\mathcal{R}o_f \approx 0.05$}. Several {numerical simulations with} helical and non-helical 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 {$\sim k_{\perp}^{-5/3}$} scaling, and the other that corresponds to a steeper {$\sim k_{\perp}^{-3}$} spectrum in which the three-dimensional (3D) modes release a substantial fraction of their energy per unit time to 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 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 wavenumbers, large-scale shear is created resulting in another time scale $\tau_{sh}$, associated with shear, {thereby producing} a $\sim k^{-1}$ spectrum for the {total energy} with the 2D modes still following a {$\sim k_{\perp}^{-5/3}$} scaling.

We first summarize briefly several properties concerning the dynamics of two-dimensional (2D) turbulence, with an emphasis on the inverse cascade of energy to the largest accessible scale of the system. In order to study a similar phenomenon in three-dimensional (3D) turbulence undergoing strong solid-body rotation, we test a previously developed Large Eddy Simulation (LES) model against a high-resolution direct numerical simulation of rotating turbulence on a grid of $3072^3$ points. We then describe new numerical results on the inverse energy cascade in rotating flows using this LES model and contrast the case of 2D versus 3D forcing, as well as non-helical forcing (i.e., with weak overall alignment between velocity and vorticity) versus the fully helical Beltrami case, both for deterministic and random forcing. The different scaling of the inverse energy cascade can be attributed to the dimensionality of the forcing, with, in general, either a $k_{\perp}^{-3}$ or a $k_{\perp}^{-5/3}$ energy spectrum of slow modes at large scales, perpendicular referring to the direction of rotation. We finally invoke the role of shear in the case of a strongly anisotropic deterministic forcing, using the so-called ABC flow., Comment: 10 pages, 3 figures

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 backreaction from any non-helical 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., Comment: 5 pages, 4 figures

This study shows the connection between three previously observed but seemingly unrelated phenomena in hydrodynamic (HD) and magnetohydrodynamic (MHD) turbulent flows, involving the emergence of fluctuations occurring on very long time scales: the low-frequency 1/f noise in the power frequency spectrum, the delayed ergodicity of complex valued amplitude fluctuations in wavenumber space, and the spontaneous flippings or reversals of large scale fields. Direct numerical simulations of ideal MHD and HD are employed in three space dimensions, at low resolution, for long periods of time, and with high accuracy to study several cases: Different geometries, presence of rotation and/or a uniform magnetic field, and different values of the associated conserved global quantities. It is conjectured that the origin of all these long-time phenomena is rooted in the interaction of the longest wavelength fluctuations available to the system with fluctuations at much smaller scales. The strength of this non-local interaction is controlled either by the existence of conserved global quantities with a back-transfer in Fourier space, or by the presence of a slow manifold in the dynamics., Comment: 17 pages, 21 figures, accepted for publication in Physical Review E

We present numerical evidence of how three-dimensionalization occurs at small scale in rotating turbulence with Beltrami (ABC) forcing, creating helical flow. The Zeman scale $\ell_{\Omega}$ at which the inertial and eddy turn-over times are equal is more than one order of magnitude larger than the dissipation scale, with the relevant domains (large-scale inverse cascade of energy, dual regime in the direct cascade of energy $E$ and helicity $H$, and dissipation) each moderately resolved. These results stem from the analysis of a large direct numerical simulation on a grid of $3072^3$ points, with Rossby and Reynolds numbers respectively equal to 0.07 and $2.7\times 10^4$. At scales smaller than the forcing, a helical wave-modulated inertial law for the energy and helicity spectra is followed beyond $\ell_{\Omega}$ by Kolmogorov spectra for $E$ and $H$. Looking at the two-dimensional slow manifold, we also show that the helicity spectrum breaks down at $\ell_{\Omega}$, a clear sign of recovery of three-dimensionality in the small scales., Comment: 13 pages, 6 figures

We examine three--dimensional turbulent flows in the presence of solid-body rotation and helical forcing in the framework of stochastic Schramm-L\"owner evolution curves (SLE). The data stems from a run on a grid of $1536^3$ points, with Reynolds and Rossby numbers of respectively 5100 and 0.06. We average the parallel component of the vorticity in the direction parallel to that of rotation, and examine the resulting $<\omega_\textrm{z}>_\textrm{z}$ field for scaling properties of its zero-value contours. We find for the first time for three-dimensional fluid turbulence evidence of nodal curves being conformal invariant, belonging to a SLE class with associated Brownian diffusivity $\kappa=3.6\pm 0.1$. SLE behavior is related to the self-similarity of the direct cascade of energy to small scales in this flow, and to the partial bi-dimensionalization of the flow because of rotation. We recover the value of $\kappa$ with a heuristic argument and show that this value is consistent with several non-trivial SLE predictions., Comment: 4 pages, 3 figures, submitted to PRL

With the help of a model of magnetohydrodynamic (MHD) turbulence tested previously, we explore high Reynolds number regimes up to equivalent resolutions of 6000^3 grid points in the absence of forcing and with no imposed uniform magnetic field. For the given initial condition chosen here, with equal kinetic and magnetic energy, the flow ends up being dominated by the magnetic field, and the dynamics leads to an isotropic Iroshnikov-Kraichnan energy spectrum. However, the locally anisotropic magnetic field fluctuations perpendicular to the local mean field follow a Kolmogorov law. We find that the ratio of the eddy turnover time to the Alfven time increases with wavenumber, contrary to the so-called critical balance hypothesis. Residual energy and helicity spectra are also considered; the role played by the conservation of magnetic helicity is studied, and scaling laws are found for the magnetic helicity and residual helicity spectra. We put these results in the context of the dynamics of a globally isotropic MHD flow which is locally anisotropic because of the influence of the strong large-scale magnetic field, leading to a partial equilibration between kinetic and magnetic modes for the energy and the helicity., Comment: 11 pages, 7 figures, submitted to Phys. Rev. E