Direct and inverse energy cascades in a forced rotating turbulence experiment
Abstract
We present experimental evidence for a double cascade of kinetic energy in a statistically stationary rotating turbulence experiment. Turbulence is generated by a set of vertical flaps which continuously injects velocity fluctuations towards the center of a rotating water tank. The energy transfers are evaluated from twopoint thirdorder threecomponent velocity structure functions, which we measure using stereoscopic particle image velocimetry in the rotating frame. Without global rotation, the energy is transferred from large to small scales, as in classical threedimensional turbulence. For nonzero rotation rates, the horizontal kinetic energy presents a double cascade: a direct cascade at small horizontal scales and an inverse cascade at large horizontal scales. By contrast, the vertical kinetic energy is always transferred from large to small horizontal scales, a behavior reminiscent of the dynamics of a passive scalar in twodimensional turbulence. At the largest rotation rate the flow is nearly twodimensional, and a pure inverse energy cascade is found for the horizontal energy. To describe the scalebyscale energy budget, we consider a generalization of the KármánHowarthMonin equation to inhomogeneous turbulent flows, in which the energy input is explicitly described as the advection of turbulent energy from the flaps through the surface of the control volume where the measurements are performed.
I Introduction
Global rotation is a key ingredient of many geophysical and astrophysical flows. Through the action of the Coriolis force, rotating turbulence tends to approach twodimensionality, i.e. invariance along the rotation axis (hereafter denoted as the vertical axis by convention).CambonBook (); DavidsonBook2013 (); Godeferd2014 () Energetic 2D and 3D flow features therefore coexist in rotating turbulence, and the question of the direction of the energy cascade between spatial scales naturally arises: In 3D, energy is transferred from large to small scalesCambonBook (); DavidsonBook2013 (); FrischBook () whereas it is transferred from small to large scales in 2D, as first proposed by Kraichnan.Kraichnan1967 (); Lindborg1999 (); Tabeling2002 () In rotating turbulence, energy transfers depend on the Rossby number , which compares the rotation period to the turbulent turnover time. In the limit of small , the fluid motions evolving on a time scale much slower than the rotation period are 2D3C (twodimensional, threecomponent), a result known as the TaylorProudman theorem, while the faster motions of frequency up to are in the form of 3D inertial waves.GreenspanBook () In this limit, 3D energy transfers occur through resonant and quasiresonant triadic interactions of inertial waves,Cambon1989 (); Waleffe1992 (); Waleffe1993 (); Bordes2012 (); Smith1999 () which drive energy in a direct cascade, with a net transfer towards slow, smallscale, nearly 2D modes.Waleffe1993 (); Galtier2003 (); Cambon2004 () Exactly resonant triads cannot however drive energy from 3D modes to the exactly 2D mode. In the limit of vanishing Rossby number, only those exact resonances are efficient, so the 2D3C mode is autonomous:Bourouiba2008 () It follows a purely 2D dynamics unaffected by rotation, with an inverse cascade of horizontal energy and a passivescalar mixing of the vertical velocity.FrischBook (); Tabeling2002 () This decoupling implies that, if energy is supplied to the 3D modes only, the 2D mode should not be excited and no inverse cascade should be observed.
In contrast with this asymptotic limit, most experiments and numerical simulations correspond to moderate Rossby numbers. They exhibit the emergence of largescale columnar structures, which suggests a net transfer from the 3D “wave” modes to the 2D3C “vortex” mode.Hopfinger1982 (); Bartello1994 (); Bourouiba2007 (); Mininni2009 (); Mininni2010 (); Pouquet2010 (); Moisy2011 (); Bourouiba2012 (); Pouquet2013 () For such moderate Rossby numbers, nearresonant triadic interactions, which are increasingly important as is increased, allow for nonvanishing energy transfers between 3D and 2D modes,Waleffe1993 (); Chen2005 (); Smith2005 (); Bourouiba2007 (); Bourouiba2012 () thus providing a mechanism for the emergence of inverse energy transfers: once energy is transferred to the 2D vortex mode, local 2D interactions are expected to build an upscale energy cascade. Even for a purely 3D forcing, the vortex mode grows as a result of nearresonant triads involving one 2D mode and two largeverticalscale and smallhorizontalscale 3D modes:Bourouiba2012 () this vortex mode then triggers inverse energy transfers between purely 2D modes. This intermediate Rossby number regime is of first practical interest: the Rossby number of most laboratory experiments and geophysical/astrophysical flows is indeed of the order of . In these situations, a natural question is to what extent direct and inverse cascades may coexist, and what sets their relative amplitudes as the Rossby number is varied.
Inverse energy cascade in rotating turbulence has been mostly investigated numerically, in the simplified configuration of a body force acting at an intermediate wave number in a periodic box.Hossain1994 (); Smith1996 (); Yeung1998 (); Smith1999 (); Chen2005 (); Mininni2009 (); Sen2012 (); Pouquet2013 (); Deusebio2014 () In this setup the inverse cascade is manifested through a growth of the energy spectrum, and hence an inverse spectral transfer, at wave numbers (with the wave number component normal to the rotation axis). As for 2D turbulence, the kinetic energy increases during this transient regime, until energetic domain size structures are formedGallet2013 () or additional largescale dissipation comes into play. Although much weaker, an inverse transfer of energy is also found in numerical simulations of decaying rotating turbulence.Bourouiba2007 (); Thiele2009 (); Teitelbaum2011 () Overall, these simulations indicate that, in addition to the Rossby number, the nature of the forcing, in particular its dimensionality (2D vs. 3D), componentality (2C vs. 3C) and helicity content, play key roles for the existence and intensity of the inverse cascade.Smith2005 (); Sen2012 () In addition, since shallow domains resemble 2D systems, which enhances the inverse cascade, another key parameter in this problem is the vertical confinement: the critical Rossby number under which the inverse cascade appears increases as the ratio of the box height to the forcing scale gets smaller.Smith1996 (); Smith1999 (); Deusebio2014 ()
Although these numerical simulations have provided valuable insight about the conditions under which an inverse cascade takes place in rotating turbulence, the most common assumptions of homogeneity and narrowband spectral forcing are of limited practical interest. More general forcing functions are considered in the simulations of Bourouiba et al.,Bourouiba2012 () with energy input either in a large range of vertical scales and a single horizontal scale, or viceversa. In most flows encountered in the laboratory and in geophysical/astrophysical contexts, energy injection in a given control volume is broadband and results from the spatial gradients of turbulent energy. As a consequence, the wellseparated inverse and direct cascades obtained in numerical simulations with a separating wave number fixed at the forcing wave number are not relevant to describe real flows with boundary forcing. Furthermore, flows of geophysical relevance can often be considered to be in statistically steady state. Such stationary states are easily achieved in laboratory experiments, whereas they generally correspond to prohibitively long integration times for numerical simulation. This provides another justification for considering the problem of the energy cascade directions of rotating turbulence experimentally.
We therefore built an experiment aimed at studying such stationary rotating turbulence. Designing a rotating turbulence experiment which unambiguously exhibits an inverse cascade is however difficult for several reasons. First, in a statistically steady turbulence experiment, an inverse cascade can be identified only from measurements of energy transfers, i.e. from thirdorder velocity correlations. These measurements require very large data sets from advanced imagebased diagnostic such as stereoscopic particle image velocimetry.Lamriben2011 () Second, it is possible to separate the scalebyscale energy fluxes from the spatial transport of energy only under the assumption of weak inhomogeneity of the flow, which is difficult to satisfy with boundary forced experiments.
Because of these difficulties, experimental evidence of inverse cascade in rotating turbulence is scarce. Indirect evidence was first provided by Baroud et al.Baroud2002 () in forced turbulence and later by Morize et al.Morize2005 () in decaying turbulence. In both experiments it is reflected in a change of sign of the thirdorder moment of the longitudinal velocity increments in the plane normal to the rotation axis, (where is the velocity increment projected along the horizontal separation ). Simple relations between and the energy flux exist only either in the 3D3C isotropic case or in the 2D2C isotropic case, but not in the general axisymmetric case, so the change of sign of cannot be unambiguously related to inverse energy transfers in these experiments. More recently, evidence of inverse energy transfers has been reported by Yarom et al.,Yarom2013 () from the transient evolution of the energy spectrum in a forced rotating turbulence experiment. However, because of the unstationnary and inhomogeneous nature of their experiment, it is delicate to distinguish the scalebyscale energy transfers at a given spatial location from the spatial energy transport from the turbulence production device to the measurement area. In all these experiments the aspect ratio is of order unity, so the 2D features of turbulence are essentially due to rotation and not confinement. The extreme case of rotating shallow water experiments is indeed known to produce a purely 2D dynamics with an inverse energy cascade even at small rotation rate (see, e.g., Afanasyev and CraigAfanasyev2013 ()). The integral scales measurements of van Bokhoven et al.,Bokhoven2009 () in which both the fluid height and the rotation rate are varied, also confirm the combined roles of these two parameters in the generation of largescale quasi2D vortices.
In this paper we investigate the interplay between direct and inverse energy cascades in a statistically stationary rotating turbulence experiment from direct measurements of scalebyscale energy transfers. Turbulence is generated by a set of vertical flaps which continuously inject velocity fluctuations towards the center of a rotating water tank. The flaps are vertically invariant, but instabilities in their vicinity induce 3D turbulent fluctuations, so the forcing injects energy both in the 2D and 3D modes. We compute the energy transfers from the divergence of the twopoint thirdorder velocity structure functions extracted from stereoscopic particle image velocimetry measurements in the rotating frame. We observe the emergence of a double cascade of energy, direct at small scales and inverse at large scales, the extension and magnitude of the inverse cascade increasing with global rotation. This overall behavior of the total kinetic energy is the superposition of different behaviors for the horizontal and vertical velocities: for rapid global rotation, the horizontal energy exhibits an inverse cascade, whereas the vertical energy follows a direct cascade. The inverse cascade of horizontal energy is found only at large scale for moderate rotation rate, but gradually spreads down to the smallest scales as the rotation rate is increased. These findings are compatible with a 2D3C dynamics at large rotation rate, with the horizontal velocity following a 2D dynamics and the vertical velocity behaving as a passive scalar.
The energy transfers in homogeneous (but not necessarily isotropic) turbulence can be described in the physical space using the KármánHowarthMonin (KHM) equation.Monin1975 (); FrischBook (); Lamriben2011 () This approach holds for homogeneous decaying turbulence and for stationary turbulence forced by a homogeneous body force. However, it breaks down in boundaryforced experiments, in which inhomogeneities induce a transport of kinetic energy from the forcing region to the region where measurements are performed. Extended versions of the KHM equation including the effects of inhomogeneities have been proposed and proved useful to describe the energy budget in simple configurations, e.g. in windtunnel experiments.Lindborg1999b (); Hill2002 (); Antonia2006 (); Danaila2012a (); Danaila2012b () Here we make use of the inhomogeneous generalization of the KHM equation proposed by Hill.Hill2002 () The measurement of the different terms of this equation in the case of the largest rotation rate, which is closer to the asymptotic 2D3C state, allows us to clarify the effect of the inhomogeneous forcing in this experiment.
Ii Experimental Setup
The experimental setup is similar to the one described in Gallet et al.,Gallet2014 () and only the features specific to the present experiments are described in detail here. The setup consists of a glass tank of cm square base and cm height, filled with cm of water and mounted on a precision rotating platform of 2 m diameter (see Fig. 1(a)). We have carried out experiments at five rotation rates in the range to rad s ( to rpm), together with a reference experiment without rotation (). The rotation rate is constant to better than relative fluctuations. In the central region of the tank, we use a glass lid to avoid the paraboloidal deformation of the free surface and to allow for visualization from above. This lid is cm above the bottom of the tank.
A statistically stationary turbulent flow is produced by a set of ten vortex dipole generators. They are arranged in 5 blocks of 2 generators located around a hexagonal arena of cm width, each of them being oriented towards the center of the arena (Fig. 1(a)). One side of the hexagon is left open to illuminate the center of the arena with a horizontal laser sheet. This forcing device was initially designed to generate turbulence in stratified fluids, and is described in detail in Refs. Billant2000, ; Augier2014, . Each generator consists of a pair of vertical flaps, 60 cm high and cm long, each flap rotating about one of its vertical edges (Fig. 1(b)). Thanks to DC motors and a system of gears and cams, the pairs of flaps are driven in a periodic motion of s duration and 9 amplitude: the two flaps being initially parallel, they first rotate with an angular velocity rad s during s until their tips almost touch each other. They remain motionless during s, before reopening during 1.7 s until they reach the initial parallel configuration again. They finally remain motionless during the last 1.7 s of the cycle. The motions of the two adjacent pairs of flaps of a given block are in phase, but an arbitrary phase shift is set between the five blocks. The rotation of the platform is set long before the start of this forcing device, at least 1 h, in order for transient spinup recirculations to be damped. Once solidbody rotation is reached, we start the generators, and a statistically stationary state is reached in the center of the tank after a few minutes.
The Reynolds number based on the flap length and flap angular velocity is . The flow generated by the closing of the flaps consists of an initially verticallyinvariant vortex dipole (Fig. 1(a)) which quickly becomes unstable and produces smallscale 3D turbulent fluctuations. This turbulent burst selfadvects towards the center of the arena because of the persistent largescale vortex dipole component. The Rossby number based on the flap angular velocity is low, (see Tab. 1), indicating that the flow generated by the flap motion is influenced by rotation right from the generators (except for the nonrotating experiment). Turbulence in the center of the flow can be also characterized locally by the turbulent Reynolds and Rossby numbers based on the r.m.s. velocity and the horizontal integral scale defined in Eq. (4) (see next Section): and (values are given in Tab. 1).
(rpm)  0  2  4  8  12  16 

0.22  0.11  0.055  0.037  0.028  
0.48  0.79  0.89  0.95  0.94  0.97  
0.17  0.15  0.06  0.05  0.04  0.10  
(mm)  24  43  45  44  42  38 
140  230  350  420  400  330  
0.30  0.20  0.13  0.087  0.068 
The three components of the velocity field are measured in a horizontal and a vertical plane in the rotating frame (Fig. 1(a)) using a stereoscopic particle image velocimetry (PIV) system.DaVis (); Pivmat () The two regions of interest are centered with respect to the arena of generators. It is a square of about cm in a vertical plane along the diagonal of the base of the tank, and a square of cm in a horizontal plane at middepth. The flow is seeded with m tracer particles and illuminated by a laser sheet generated by a double mJ Nd:YAG pulsed laser mounted on the rotating platform (Fig. 1(a)). The illuminated flow section is imaged with two doublebuffer cameras aiming at the laser sheet under different incidence angles. Images are taken from above through the glass lid for the measurements in the horizontal plane, and from two adjacent vertical sides of the tank for the measurements in the vertical plane.
Each acquisition consists in 3 600 quadruplets of images (one pair per camera) recorded at 0.35 Hz with a ms time lag between the two images of a given pair. The 3 velocity components are computed in the twodimensional measurement plane using stereoscopic reconstruction. The crosscorrelations are based on interrogation windows of pixels with overlap. The resulting 2D3C velocity fields are sampled on a grid of (resp. ) vectors with a spatial resolution of mm (resp. mm) in the horizontal (resp. vertical) plane.
Iii Local homogeneity and axisymmetry
In this experiment, kinetic energy is injected by the generators located around the region of interest, so an inward transport of energy takes place from the generators to the center of the arena. An important feature of turbulence in this configuration is the presence or not of a mean flow induced by the generators: this indicates whether the transport of energy can be mainly attributed to a reproducible flow or to turbulent fluctuations. This can be addressed by performing a Reynolds decomposition of the velocity field,
(1) 
where is the timeaveraged velocity field and its turbulent part. From this decomposition, we can compute the turbulent and total kinetic energies, and respectively, where is a spatial average over the horizontal region of interest. The turbulence rate is about in the nonrotating experiment but rapidly increases up to 97% as the rotation rate is increased (see Tab. 1), indicating that, under rotation, the turbulent structures are essentially selfadvected from the generators towards the center of the arena. In the following we therefore focus on the turbulent component which dominates the flow in the rotating case.
Although turbulence is necessarily inhomogeneous in this configuration, with more energy near the generators than at the center of the flow, we may expect a reasonable local homogeneity in the measurement area, because of its small size (square of about 13 cm side) compared to the distance to the generators (33 cm from the center of the arena to the tip of the flaps). Before investigating the scalebyscale energy distribution and energy transfers from spatially averaged twopoint statistics, it is therefore important to quantify the degree of homogeneity of the flow. Since most of the energy is turbulent, we can quantify the level of homogeneity by the spatial standard deviation of the turbulent kinetic energy
(2) 
with the local timeaveraged turbulent energy (such that ). This ratio is given in Tab. 1. It is smaller than for rpm, indicating a reasonable degree of homogeneity in the region of interest.
A last singlepoint quantity of interest to characterize the turbulence field in this configuration is the spatially averaged velocity correlation tensor, (the trace of this tensor is twice the turbulence kinetic energy). For axisymmetric turbulence with respect to , one has (i.e. turbulence is isotropic in the horizontal plane), with zero nondiagonal components. For 3D isotropic turbulence the three diagonal components are equal (i.e., ). This tensor therefore characterizes the componentality of turbulence, i.e. the isotropy with respect to the velocity components. It must not be confused with the dimensionality of turbulence, which characterizes the dependence of the twopoint velocity statistics with respect to orientation of the separation vector joining the two points (investigated in next section).
In Fig. 2 we see that turbulence is nearly axisymmetric, with to within 3% in the rotating case and 10% in the nonrotating case, and with the three nondiagonal components less than 10% of the diagonal components. As expected, turbulence is never isotropic, even in the case , for which . This anisotropy originates from the verticallyinvariant forcing by the flaps, which induces significantly weaker vertical velocity fluctuations than horizontal ones. As the rotation rate increases, remains roughly constant, whereas and first increase with and saturates beyond 8 rpm (). At large , the vertical kinetic energy represents about 10% of the total energy.
Iv Scalebyscale energy distribution and transfers
We now focus on the scalebyscale energy distribution and energy transfers. For this we must use twopoint quantities: let us consider two points A and B in the turbulent flow at positions and . We define the midpoint position and the separation vector . Using cylindrical coordinates the separation writes (,,), with and . In homogeneous turbulence, all statistical averages are functions of the separation vector only. However, inhomogeneity plays a key role in boundary forced experiments, and we thus consider the inhomogeneous framework in which ensemble averages remain functions of both and .
The centered velocity increment for separation , midpoint and time is
To perform an energy budget in this inhomogeneous context, we first need to define a control region for the midpoint : we consider only points A and B for which lies in square of side 40 mm centered in the PIV field. The relatively small size of the square allows for the separation to be as large as mm with the two points A and B still lying in the PIV field. The statistical averages are defined as an average over time and over all the positions of inside . In the following, the spatial average over will simply be denoted by .
Defining the statistics from centered or noncentered increments would be equivalent for homogeneous turbulence, and most of the remainder of this section can be understood in this framework. However, the use of centered increments plays a key role in Sec. V where we discuss the scalebyscale energy budget: there is a balance at every scale between viscous dissipation, nonlinear transfers between different scales (flux in space), and advection of kinetic energy at scale through the boundaries of the domain (flux in space). The latter term is the source term of the energy budget, which vanishes if turbulence is assumed to be homogeneous.
iv.1 Energy distribution
We characterize the distribution of the turbulent energy among spatial scales by the anisotropic secondorder structure function, defined as the variance of the centered velocity incrementsfn1 ()
(3) 
The angular average of this quantity, , where is the usual spherical coordinate system, can be interpreted as the energy contained in eddies of size or less, provided that is larger than the dissipative scale.DavidsonBook2004 (); DavidsonBook2013 () For isotropic turbulence, therefore directly measures the cumulative energy from to . For anisotropic turbulence, contains in addition key information on the anisotropic distribution of energy among eddies of characteristic horizontal and vertical scale given by and respectively. For isotropic turbulence, the isosurfaces of are spherical, while for axisymmetric turbulence about the vertical they are invariants with respect to rotations around the axis. Twodimensional turbulence would give exactly cylindrical iso (invariant by translation along ), which is a special case of axisymmetric turbulence.
Figure 3 shows the maps of the normalized energy distributions in the horizontal and vertical planes for and 16 rpm, with taken at and mm the maximum separation.fn2 () In Figs. 3(a,b), the isocontours of are nearly circular in the horizontal plane, both without and with rotation, indicating the good level of twopoint axisymmetry of turbulence. For the largest horizontal scales considered here ( mm), reaches for , indicating that nearly all turbulent energy is contained in the range of scales of interest, whereas it reaches only at rpm, indicating that structures larger than the maximum available scale still contain energy. This is a first indication of the emergence of large horizontal structures in the presence of rotation. This effect can be further quantified by the horizontal integral scale,
(4) 
with , and the twopoint correlation function. The conventional definition is such that , but using here a finite truncation at , chosen such that , is necessary because does not reach at the maximum available scale mm. Values of are given in Tab. 1. In the absence of rotation, mm, which corresponds to the characteristic size of the turbulent fluctuations generated by the flaps. As is increased, grows by nearly a factor 2 compared to the nonrotating case, confirming the generation of largescale structures.
We now turn to the energy distribution in the vertical plane (Figs. 3c,d). is anisotropic both with and without rotation, with a trend towards vertical elongation of the contour lines at large scales. This anisotropy is however weak at and affects preferentially the large scales: this is a direct consequence of the vertical invariance of the forcing device, which creates a nearly 2D flow at large scale carrying smallscale 3D fluctuations. The anisotropy is much more pronounced in the presence of rotation and persists down to the smallest scales, indicating a trend towards quasi2D turbulence. This scale dependence of the anisotropy can be quantified by the ratio
(5) 
with and . This ratio is zero for 3D isotropic turbulence and 1 for 2D turbulence. The plot of in Fig. 4 shows a growth of anisotropy with at all rotation rates. This growth is weak for (with increasing from 0.01 to 0.1), indicating that the 2D nature of the forcing has a weak influence at these scales in the center of the tank. The anisotropy is much stronger when : first grows rapidly from to mm before saturating. For the largest available rotation rate, rpm, turbulence is nearly 2D for mm, with , but remains significantly 3D at smaller scales.
iv.2 Energy transfers
We now consider the scalebyscale energy transfers, defined from thirdorder moments of velocity increments. We start from the KármánHorwarthMonin equation for timedependent homogeneous (but not necessarily isotropic) turbulenceMonin1975 (); FrischBook ()
(6) 
where , and
(7) 
is the energy transfer term in scalespace (with the divergence with respect to the vector separation ), and the instantaneous energy dissipation rate. Here the brackets represent spatial and ensemble average. Similarly to the angular average of , which represents the cumulative energy from scale to , the angular average can be interpreted as the energy flux from scales smaller than towards scales larger than . For isotropic turbulence, the sign of therefore gives the direction of the energy cascade, forward if and inverse if . In the inhomogeneous case, additional terms corresponding to advection of energy between different regions of the turbulent flow appear in Eq. (6). In the absence of body forces, which are not relevant in our experiment, advection of energy from outside the control domain is the only source term to sustain stationary turbulence: we will come back to this point in Sec. V.
In the following we focus on stationary turbulence, and we take in Eq. (6). The ensemble average can be therefore replaced by a temporal average, which we denote as . For axisymmetric turbulence, it is convenient to decompose the flux (7) into its perpendicular (horizontal) and parallel (vertical) contributions,
(8)  
with and . We focus in the following on pure horizontal separations by setting , and we perform an azimuthal average over to improve the statistics. Both contributions from Eq. (8) are then functions of the horizontal separation only.
For strictly 2D turbulence, vertical invariance implies . The vertical flux cannot be measured here, because we cannot access the vertical derivative from measurements in the horizontal plane. In principle one could use the data in the vertical plane, but we found significant departure from axisymmetry for thirdorder quantities (although secondorder quantities are found nearly axisymmetric, as shown in Fig. 3a). This lack of axisymmetry can be circumvented by performing an azimuthal average with respect to , which is possible in the horizontal plane only.
For the nonrotating experiment, Fig. 5 shows that the horizontal flux is negative at all scales, as expected for a direct energy cascade from large to small scales, with as in the viscous range. The observed decrease in for beyond 15 mm might be due to the vertical invariance of the forcing device: the large scales are slightly 2D (see Fig. 4), which may enhance the inverse energy transfers and reduce the direct ones. However, as discussed in section V, boundarydriven flows display some inhomogeneity at large scale, which challenges the interpretation of in terms of pure scalebyscale energy transfers. Therefore, we do not elaborate more on the behavior of at scales larger than mm (see details in section V).
For increasing rotation rates, strongly decreases at intermediate scales in absolute value, and eventually a change of sign is observed at large scales beyond 4 rpm (). This indicates the onset of an inverse energy cascade, which spreads towards smaller scales as increases. Remarkably, the double cascade persists even at the largest rotation rate ( rpm, ), with the coexistence of an inverse flux () at large scales and a direct flux () at small scales. This implies that, on average, energy must be supplied at intermediate scales (of the order of 25 mm for rpm): we return to this point in Sec. V.
iv.3 Horizontal transfers of horizontal and vertical energy
The energy flux contains both the horizontal flux of horizontal energy, , and the horizontal flux of vertical energy, . To get further insight into the double cascade observed in Fig. 5, we decompose as follows
(9)  
(10) 
These two contributions are shown in Fig. 6. Interestingly, we observe that remains negative at all rotation rates, indicating that vertical energy is always transferred from large to small horizontal scales, whereas becomes positive as the rotation rate is increased, a signature of the onset of an inverse cascade for the horizontal energy. In the nonrotating case, this negative flux is compatible with the classical direct cascade framework of 3D turbulence. By contrast, in the presence of rotation, for the scales at which the inverse cascade of is observed, the direct cascade of is reminiscent of the behavior of a passive scalar advected by a twodimensional flow: the stretching and folding of the vertical velocity by the horizontal flow produce small scales through filamentation, inducing a direct horizontal cascade of vertical velocity.Yaglom1949 (); Monin1975 (); Tabeling2002 () We provide in the next section further assessment of this picture.
Figure 6 also indicates that the horizontal flux of vertical energy is a significant contribution to for all rotation rates. For the low rotation rate rpm (), although a significant inverse cascade already takes place at large scale for the horizontal energy (), it is hidden by a stronger direct cascade of vertical energy (). This results in an overall direct cascade of total energy (). For larger rotation rates, the inverse cascade of horizontal energy becomes dominant, eventually leading to at all scales for rpm. The crossover scale separating the direct and inverse cascades of horizontal energy rapidly decreases as increases, going from mm for rpm to zero for rpm (and then over the whole range of scales).
V Scalebyscale energy budget
v.1 Inhomogeneous KármánHowarthMonin equation
To provide a physical interpretation for the sign of the scalebyscale energy flux , we must describe carefully the energy input in the experiment, and in particular its scale dependence. In numerical simulations of homogeneous stationary turbulence, this source term usually originates from a random body force acting on a narrow range of scales. By contrast, here the fluid motion is driven by moving solid boundaries, so the energy injection in a given control volume away from the forcing must originate from the transport of energy through the surface delimiting the control volume. Since a nontrivial stationary state cannot be described by the homogeneous nonforced KHM equation (6), which contains no source term, we must consider explicitly the effects of the inhomogeneities in the following.
We consider here the inhomogeneous generalization of the KHM equation proposed by Hill.Hill2002 () We briefly recall the derivation of this equation, with the addition of the Coriolis force. Let us start from the incompressible NavierStokes (NS) equation in the rotating frame
(11) 
with the pressure modified by the centrifugal force and normalized by the fluid density. Taking the difference between points and and taking the scalar product with yield
(12) 
with . Quantities with denote the average between the two points: , , and , with the local energy dissipation rate. The last term of the equation involves the velocity correlation tensor , with . Importantly, all the terms in Eq. (12) are functions of , and the nonlinear term splits into a scaletoscale transfer term (divergence w.r.t. separation ) and a transport term (divergence w.r.t. midpoint ).
We consider both the spatial average over a control volume and the ensemble average of Eq. (12). Using the divergence theorem to express the inhomogeneous terms as a flux through the closed surface delimiting the control volume , we obtain
(13) 
where the flux term writes
(14) 
The unit vector is directed outward of the control volume by convention. In the scalebyscale budget (13), the energy input (or output) at a given scale is ensured by the term , which originates from the inhomogeneities in the pressure and velocity statistics. For homogeneous turbulence, one has and , so Eq. (13) becomes the usual KHM equation (6).
The flux term (14) contains three contributions:

is the flux of cumulative energy through the surface due to advection by the velocity . It is positive when is directed into the control volume. Note that this term takes a simple form in the classical configuration of a windtunnel: the transport velocity is essentially replaced by the uniform mean velocity , and the energy flux per unit surface becomes . The inward flux of through the upstream face of the control volume is larger than the outward flux through the downstream face, hence a net flux of kinetic energy into the control volume, which is dissipated at the same rate by viscosity. By contrast, the timeaveraged velocity is negligible in the present experiment when , and energy input in the control volume proceeds through advection of turbulent kinetic energy by the turbulent velocity itself.

The term originates from the work of the pressure force through the boundary of the control volume. This term cannot be measured experimentally. However, under the assumption of local axisymmetry and homogeneity at small scale, this term is expected to be much smaller than the advection term for scales much smaller than the characteristic scale of inhomogeneity (see Appendix A).
We stress the fact that the present inhomogeneous KHM equation has a clear interpretation for scales smaller than the characteristic scale of inhomogeneity, for which the contribution from the pressure can be neglected. This is because the twopoint velocitypressure correlation can be written either as a divergence with respect to or with respect to , so the interpretation of as the only scalebyscale transfer term becomes questionable when the velocitypressure correlation is significant, i.e. for large (inhomogeneous) scales. In the following we show that interesting modifications of the energy transfers by global rotation occur in the range of scales for which the velocitypressure correlation is indeed negligible.
v.2 Scalebyscale budget for 2D3C flows
For rapid global rotation, the flow becomes weakly dependent on the vertical coordinate far from the horizontal boundaries. In the following we assume a purely 2D3C velocity field in the bulk of the flow, and write separate equations for the evolution of the horizontal and vertical energies. In this 2D limit, the same analysis as in Sec. V.1, but performed here on the horizontal components of Eq. (11) only, yields
(15) 
with
(16) 
where the subscript in , , and indicates that only the horizontal velocity components are considered, and the indicates that the quantity is a midpoint average. This equation is an inhomogeneous version of the KHM equation for the horizontal flow only. Similarly, from the vertical component of the NavierStokes equation, one can also compute the budget for the vertical energy,
(17) 
with
(18) 
where the subscript refers to the vertical component of the velocity. This equation is an inhomogeneous generalization of Yaglom’s equation for a passive scalar field:Yaglom1949 (); Monin1975 () the tracer enters and leaves the control volume through advection by the horizontal velocity . Inside the control volume, nonlinearities transfer the vertical energy between different scales through stretching and folding by the horizontal field, and viscosity damps the strong gradients created by these processes. We stress the fact that Eq. (17) does not involve pressure: all the terms in this equation can be accessed through stereoPIV measurements in a horizontal plane.
v.3 Experimental assessment of the horizontal and vertical kinetic energy budget
In the following we focus on the highest rotation rate, rpm (), for which we expect the turbulent flow to reach a quasi2D3C state, so that we can apply the vs. decomposition of the inhomogeneous KHM equation derived above. Note that in this case the ensembleaveraged flow is negligible (see Tab. 1), so we simply consider .
We consider for the control domain a centered square in the square PIV field, with a maximum separation mm. One can think about the corresponding control volume as a parallelepiped of arbitrary vertical length, with zero fluxes through the top and bottom boundaries. Because turbulence is stationary in the experiment, we replace the ensemble average by a temporal one, which we denote .
v.3.1 Horizontal kinetic energy budget
Let us first consider the energy budget for the horizontal kinetic energy. In statisticallysteady state, Eq. (15) writes
(19) 
with
(20) 
(21) 
The horizontal flux of horizontal energy is defined in Eq. (9), and is the dissipation of the horizontal velocity by the horizontal shearing (with summation over ), such that . In Eq. (19), contains the viscous contribution from inhomogeneities, which we neglect in the following (it is at least 2 orders of magnitude smaller than the other terms), and the vertical transport, which we cannot compute from 2D3C measurements. Here still denotes the spatial average over . Under the assumption of weak inhomogeneity, the velocitypressure flux is expected to be small compared to the transport (see Appendix A), and the remaining terms in Eq. (19) can be readily computed from the 2D3C PIV measurements.
Figure 7(a) shows the three measurable terms of the r.h.s. of Eq. (19) together with their sum, . These terms are averaged over the azimuthal angle . We observe a good agreement between and the sum for scales smaller than mm, to within 20%. The approximation of locally homogeneous and axisymmetric turbulence, with a negligible velocitypressure flux , thus seems valid at small scales. For larger scales we observe a significant departure between and , up to a factor of 2 at mm: for such large scales the turbulent flow cannot be considered as locally homogeneous anymore: the nonmeasured pressure term cannot be neglected, so the interpretation of as a scalebyscale energy transfer becomes incorrect.
Focusing on scales smaller than mm, the advection of horizontal kinetic energy is the only source term in the energy budget. This source term is maximum at large scales, as expected, but it remains significant over the whole range of scales, suggesting a broadband energy injection in this system. If is interpreted as a cumulative energy for eddies of horizontal scale or less, the scalebyscale energy density has the form , so the corresponding forcing density is . Since we observe at intermediate scales, the forcing density scales as , and it remains significant over the range of scales considered here. This broadband energy injection is an important feature of boundaryforced inhomogeneous turbulence, and is in strong contrast with the narrowband forcing often considered in numerical simulations of homogeneous rotating turbulence.
Although the inverse cascade is evident from this horizontal energy budget, it must be noted that its magnitude remains moderate: is never the dominant contribution to the budget, even at the crossover between the viscous diffusion and the forcing. A well developed inverse cascade over a wide range of quasihomogeneous scales at much larger Reynolds number would be characterized by . Here, the amount of energy transferred to large scales remains at best of the same order as the smallscale viscous dissipation .
v.3.2 Vertical kinetic energy budget
We now consider the vertical kinetic energy budget using Eq. (17) in statisticallysteady state:
(22) 
where
(23) 
(24) 
The horizontal flux of vertical energy is defined in Eq. (9), and is the dissipation of the vertical velocity by the horizontal shearing, with . The viscous contribution from inhomogeneities is once again discarded for simplicity. A key feature of Eq. (22) is that it does not involve pressure: all the terms can therefore be readily measured from stereoPIV measurements and the equation should be satisfied for scales at which quasitwodimensionality holds.
The terms of Eq. (22) are shown in Fig. 7(b), together with the sum of the righthand side terms . Once again, there is a good overall agreement between and for mm, clearly indicating a direct cascade of vertical kinetic energy. The picture here is simpler than for the horizontal energy budget: the vertical energy is advected into the control domain (), it is transferred by the nonlinearities to smaller scales (), and it is dissipated at small scales by viscosity (). We can assume that the instabilities of the vortex dipoles generated near the flaps and away from the measurement domain is the source of vertical energy at large scale. Note that for very large Reynolds numbers, an inertial range of scales with constant flux of vertical energy is expected, characterized by . Although this is not achieved at the moderate Reynolds number of this experiment, we can nonetheless define a significant range of scales, centered around 1020 mm, for which is close to to within 20%.
v.3.3 Largescale energy dissipation
For the largest rotation rate, the scenario of a double cascade, towards small scales for the vertical energy and towards large scales for the horizontal energy, is well established from Fig. 7(a,b), at least for scales mm. At larger scales, the departure between and may originate both from a departure from twodimensionality or from the nonmeasured velocitypressure correlation. Although less pronounced, there is also a discrepancy between and at large scale, which indicates a departure from a pure 2D3C state and possibly an influence of the horizontal top and bottom boundaries.
We are therefore left with the following question: what sink of energy absorbs the inverse horizontal energy flux at large scales? A first candidate is Ekman friction on the horizontal boundaries. Assuming that the flow is 2D3C in the central region of the tank, we consider the laminar scaling for the Ekman layer thickness, , and we deduce a typical energy dissipation , where is the water depth and the characteristic horizontal velocity. For rpm, we obtain m s, which turns out to be of the order of the other terms of Eq. (19) (see Fig. 7). A significant fraction of the input horizontal kinetic energy could therefore be dissipated at large scale through Ekman friction. However, this order of magnitude strongly relies on the boundary layers being laminar, which seems questionable in the present experiment.
As an alternate explanation for the energy sink at large scales, we may invoke a feedback of the largescale flow on the forcing device: the largescale flow resulting from the inverse cascade induces largescale pressure forces that transfer some kinetic energy back to the flaps, reducing the overall energy input in the system. In this scenario, the flaps inject energy at intermediate scales and receive energy from the largescale flow through the work of the pressure forces. In the framework of Eq. (19), the corresponding sink of energy is taken into account by the spatial flux and pressure terms: energy at large scales is transferred outside of the control volume, towards the flaps. Unfortunately, testing this scenario would require a precise measurement of the energy input by the flaps and is beyond the scope of the present study.
Vi Conclusion
In this paper we provide experimental evidence of a double energy cascade, direct at small scales and inverse at large scales, in a forced rotating turbulence experiment. Since turbulence is statistically steady, the inverse cascade does not manifest through a temporal growth of kinetic energy, but it is characterized by a change of sign of the scalebyscale energy flux. As the rotation rate is increased, the inverse cascade becomes more pronounced and spreads down to the smallest scales. As compared to previous experimental observations of an inverse cascade in rotating turbulence,Baroud2002 (); Morize2005 (); Yarom2013 () here we provide for the first time a direct scalebyscale measurement of the energy transfers in the horizontal plane. This allows us to distinguish between the horizontal transfers of vertical and horizontal kinetic energy. At the largest rotation rate, this double cascade of the total energy can be described as the superposition of an inverse cascade of horizontal energy and a direct cascade of vertical energy. This is consistent with the 2D3C dynamics expected in the limit of small Rossby numbers, with the vertical velocity behaving as a passive scalar transported by the horizontal flow.
Contrary to numerical simulations, in which energy is usually supplied by a homogeneous body force acting on a prescribed narrow range of wave numbers, in most experiments and in many natural flows, energy is supplied at the boundaries. For a control domain away from those boundaries, energy is advected from the boundaries into the domain. This spatial flux of energy, which is strongly related to the inhomogeneities of the turbulent statistics, results in an effective broadband energy injection term. In order to interpret the energy transfers in such an experiment, it is therefore necessary to separate the contributions from the spatial transport and from the scalebyscale transfers. We have performed this analysis using the inhomogeneous generalization of the KHM equation proposed by Hill,Hill2002 () and we have measured directly the energy transport term for scales at which the velocitypressure correlations can be neglected (quasihomogeneous approximation). Because of this effective broadband forcing, the inversion scale, which separates the direct and inverse cascades, is not directly prescribed by the geometry of the forcing device and decreases with the imposed rotation rate. Modelling this inversion scale as a function of the Rossby number and forcing geometry remains an open issue of first interest for flows of oceanic and atmospheric relevance, such as convectivelydriven rotating flows.
Acknowledgements.
We acknowledge P. Augier, P. Billant and J.M. Chomaz for kindly providing the flap apparatus, and A. Aubertin, L. Auffray, C. Borget and R. Pidoux for their experimental help. This work is supported by the ANR grant no. 2011BS0400601 “ONLITUR”. The rotating platform “Gyroflow” was funded by the “Triangle de la Physique”.Appendix A Locally homogeneous and axisymmetric turbulence
In this appendix we show how one can neglect the velocitypressure correlations in the inhomogeneous KHM equations (13)(15) under the assumptions of local homogeneity and axisymmetry.
Let us first consider that the turbulence is invariant to a reflection with respect to a horizontal plane (nonhelical turbulence). Let us then denote as the typical scale of the inhomogeneity of the turbulence statistics and focus on small separations , or equivalently . We decompose the velocity and pressure fields into and , where and contain the large scales of the flow and and are the smallscale fluctuations. One can think about and as coarsegrained versions of the velocity and pressure fields on a scale that is smaller than but larger than the range of scales of interest. Then and are the smallscale remainders.
Local axisymmetry and homogeneity consists in assuming that the statistics of and are axisymmetric and homogeneous for separations . Let us evaluate the different terms in the integrand of (14) under this assumption.
The pressure term decomposes into
(25) 
Because and evolve on a spatial scale ,
(26) 
Similarly, . Hence , and . We deal with the term using the assumption of local homogeneity and axisymmetry: under a rotation of angle around a vertical axis passing through the midpoint , followed by a reflection with respect to the horizontal plane containing , becomes , becomes , becomes and becomes . Hence becomes , so that on statistical average this quantity vanishes: . As a consequence, the velocitypressure correlation term is of order in the weakly inhomogeneous limit. We therefore expect to be negligible compared to the transport term , which is of order .
The viscous term in the integrand of (14) is of order
(27) 
It is therefore negligible compared to the advective term of the integrand, because the Reynolds number based on is large. Note that because , derivatives with respect to are much larger than derivatives with respect to . In the ensemble average of Eq. (12), the “inhomogeneous” viscous contribution can be estimated as , hence it is negligible compared to the “homogeneous” viscous contribution . For weakly inhomogeneous turbulence, one can therefore keep the latter while neglecting the former.
To conclude, the advective term is the dominant term in the integrand of Eq. (14). For weakly inhomogeneous turbulence, one can therefore retain only this term of the integrand, which ensures the injection of kinetic energy into the control volume.
References
 (1) P. Sagaut and C. Cambon, Homogeneous turbulence Dynamics (Cambridge University Press, Cambridge, 2008).
 (2) P.A. Davidson, Turbulence in rotating, stratified and electrically conducting fluids (Cambridge University Press, Cambridge, 2013).
 (3) F.S. Godeferd and F. Moisy, “Structure and dynamics of rotating turbulence: a review of recent experimental and numerical results,” Appl. Mech. Rev. (2015), doi:10.1115/1.4029006.
 (4) U. Frisch, Turbulence  The Legacy of A. N. Kolmogorov (Cambridge University Press, Cambridge, 1995).
 (5) R. Kraichnan, “Inertial ranges in twodimensional turbulence,” Phys. Fluids 10, 1417 (1967).
 (6) E. Lindborg, “Can the atmospheric kinetic energy spectrum be explained by twodimensional turbulence?,” J. Fluid Mech. 388, 259 (1999).
 (7) P. Tabeling, “Twodimensional turbulence: a physicist approach,” Physics Reports 362, 1 (2002).
 (8) H. Greenspan The theory of rotating fluids (Cambridge University Press, 1968).
 (9) C. Cambon and L. Jacquin, “Spectral approach to nonisotropic turbulence subjected to rotation,” J. Fluid Mech. 202, 295 (1989).
 (10) F. Waleffe, “Inertial transfers in the helical decomposition,” Phys. Fluids A 5, 677 (1993).
 (11) L.M. Smith and F. Waleffe, “Transfer of energy to twodimensional large scales in forced, rotating threedimensional turbulence,” Phys. Fluids 11, 1608 (1999).
 (12) F. Waleffe, “The nature of triad interactions in homogeneous turbulence,” Phys. Fluids A 4, 350 (1992).
 (13) G. Bordes, F. Moisy, T. Dauxois, P.P. Cortet, “Experimental evidence of a triadic resonance of plane inertial waves in a rotating fluid,” Phys. Fluids 24, 014105 (2012).
 (14) S. Galtier, “Weak inertialwave turbulence theory,” Phys. Rev. E 68, 015301 (2003).
 (15) C. Cambon, R. Rubinstein and F. S. Godeferd, “Advances in wave turbulence: rapidly rotating flows,” New J. Phys. 6, 73 (2004).
 (16) L. Bourouiba, “Model of a truncated fast rotating flow at infinite Reynolds number,” Phys. Fluids 20, 075112 (2008).
 (17) E.J. Hopfinger, F.K. Browand, and Y. Gagne, “Turbulence and waves in a rotating tank,” J. Fluid Mech. 125, 505 (1982).
 (18) P. Bartello, O. Métais, and M. Lesieur, “Coherent structures in rotating threedimensional turbulence,” J. Fluid Mech. 273, 1 (1994).
 (19) L. Bourouiba and P. Bartello, “The intermediate Rossby number range and twodimensional threedimensional transfers in rotating decaying homogeneous turbulence,” J. Fluid Mech. 587, 139 (2007).
 (20) L. Bourouiba, D.N. Straube, and M.L. Waite, “Nonlocal energy transfers in rotating turbulence at intermediate Rossby number,” J. Fluid Mech. 690, 129 (2012).
 (21) F. Moisy, C. Morize, M. Rabaud, and J. Sommeria, “Decay laws, anisotropy and cycloneanticyclone asymmetry in decaying rotating turbulence,” J. Fluid Mech. 666, 5 (2011).
 (22) P.D. Mininni, A. Alexakis, A. Pouquet, “Scale interactions and scaling laws in rotating flows at moderate Rossby numbers and large Reynolds numbers,” Phys. Fluids 21, 015108 (2009).
 (23) P.D. Mininni, and A. Pouquet, “Rotating helical turbulence. II. Intermittency, scale invariance, and structures,” Phys. Fluids 22, 035106 (2010).
 (24) A. Pouquet, P.D. Mininni “The interplay between helicity and rotation in turbulence: implications for scaling laws and smallscale dynamics,” Phil. Trans. R. Soc. A 368, 1635 (2010). revue sur scaling, mais sans cascade inverse
 (25) A. Pouquet, A. Sen, D. Rosenberg, P.D. Mininni, and J. Baerenzung, “Inverse cascades in turbulence and the case of rotating flows,” Phys. Scripta T155, 014032 (2013).
 (26) Q. Chen, S. Chen, G.L. Eyink, and D.D. Holm, “Resonant interactions in rotating homogeneous threedimensional turbulence,” J. Fluid Mech. 542, 139 (2005).
 (27) L. M. Smith and Y. Lee, “On near resonances and symmetry breaking in forced rotating flows at moderate Rossby number, ” J. Fluid Mech. 535, 111 (2005).
 (28) M. Hossain, “Reduction in the dimensionality of turbulence due to a strong rotation,” Phys. Fluids 6, 1077 (1994).
 (29) A. Sen, P.D. Mininni, D. Rosenberg, and A. Pouquet, “Anisotropy and nonuniversality in scaling laws of the largescale energy spectrum in rotating turbulence,” Phys. Rev. E 86, 036319 (2012).
 (30) L.M. Smith, J.R. Chasnov, and F. Waleffe, “Crossover from Two to ThreeDimensional Turbulence,” Phys. Rev. Lett. 77, 2467 (1996).
 (31) E. Deusebio, G. Boffetta, E. Lindborg, and S. Musacchio, “Dimensional transition in rotating turbulence,” Phys. Rev. E 90, 023005 (2014).
 (32) P.K. Yeung and Ye Zhou, “Numerical study of rotating turbulence with external forcing,” Phys. Fluids 10, 2895 (1998).
 (33) B. Gallet and W.R. Young, “A twodimensional vortex condensate at high Reynolds number,” J. Fluid Mech. 715, 359 (2013).
 (34) M. Thiele and W.C. Muller, “Structure and decay of rotating homogeneous turbulence,” J. Fluid Mech. 637, 425 (2009).
 (35) T. Teitelbaum and P.D. Mininni, “The decay of turbulence in rotating flows,” Phys. Fluids 23, 3592325 (2011).
 (36) C. Lamriben, P.P. Cortet, and F. Moisy, “Direct measurements of anisotropic energy transfers in a rotating turbulence experiment,” Phys. Rev. Lett. 107, 024503 (2011).
 (37) C. N. Baroud, B. B. Plapp, Z.S. She, and H. L. Swinney, “Anomalous SelfSimilarity in a Turbulent Rapidly Rotating Fluid,” Phys. Rev. Lett. 88, 114501 (2002).
 (38) C. Morize, F. Moisy, and M. Rabaud, “Decaying gridgenerated turbulence in a rotating tank,” Phys. Fluids 17, 095105 (2005).
 (39) E. Yarom, Y. Vardi, and E. Sharon, “Experimental quantification of inverse energy cascade in deep rotating turbulence,” Phys. Fluids 25, 085105 (2013).
 (40) Y. D. Afanasyev and J. D. C. Craig, “Rotating shallow water turbulence: Experiments with altimetry,” Phys. Fluids 25, 106603 (2013).
 (41) L. J. A. van Bokhoven, H. J. H. Clercx, G. J. F. van Heijst and R. R. Trieling, “Experiments on rapidly rotating turbulent flows,” Phys. Fluids 21, 096601 (2009).
 (42) A.S. Monin and A.M. Yaglom, Statistical Fluid Mechanics, vol. 2 (MIT Press, Cambridge, 1975).
 (43) E. Lindborg, “Correction to the fourfifths law due to variations of the dissipation,” Phys. Fluids 11, 510 (1999).
 (44) R.J. Hill, “Exact secondorder structure function relationships,” J. Fluid Mech. 468, 317 (2002).
 (45) R.A. Antonia, P. Burattini, “Approach to the 4/5 law in homogeneous isotropic turbulence,” J. Fluid Mech. 550, 175 (2006).
 (46) L. Danaila, R.A. Antonia and P. Burattini, “Comparison between kinetic energy and passive scalar energy transfer in locally homogeneous isotropic turbulence,” Physica D: Nonlinear Phenomena 241, 224 (2012).
 (47) L. Danaila, J.F. Krawczynski, F. Thiesset and B. Renou, “Yaglomlike equation in axisymmetric anisotropic turbulence,” Physica D: Nonlinear Phenomena 241, 216 (2012)
 (48) B. Gallet, A. Campagne, P.P. Cortet, and F. Moisy, “Scaledependent cycloneanticyclone asymmetry in a forced rotating turbulence experiment,” Phys. Fluids 26, 035108 (2014).
 (49) P. Billant and J.M. Chomaz, “Experimental evidence for a new instability of a vertical columnar vortex pair in a strongly stratified fluid,” J. Fluid Mech. 418, 167 (2000).
 (50) P. Augier, P. Billant, M.E. Negretti, and J.M. Chomaz, “Experimental study of stratified turbulence forced with columnar dipoles,” Phys. Fluids. 26, 046603 (2014).
 (51) DaVis, LaVision GmbH, AnnaVandenhoeckRing 19, Goettingen, Germany.
 (52) F. Moisy, PIVMat toolbox for Matlab, http://www.fast.upsud.fr/pivmat
 (53) P.A. Davidson, Turbulence: An Introduction for Scientists and Engineers (Oxford University Press, Oxford, 2004).
 (54) A.M. Yaglom, “On the local structure of a temperature field in a turbulent flow,” Dokl. Akad Nauk SSSR 69(6),743 (1949).
 (55) This anisotropic second order structure function (3) is directly related to the 3D turbulent energy spectrum : , with the turbulent kinetic energy and the wavevector.
 (56) The normalization by is chosen so that when an exact decorrelation of and is observed beyond . Note that normalizing by