Search

While heat transport by baroclinic turbulence in oceans and planetary atmospheres is well described by a two-layer model, the relative depth of the two layers varies greatly depending on the situation of interest, making it an important parameter governing the transport properties of the system. Focusing on the low-drag turbulent regime, we extend the vortex-gas scaling theory to address the case of arbitrary layer depths. To wit, we map the arbitrary-layer-depth system onto an equivalent equal-depth system with rescaled parameters, establishing the asymptotic validity of the mapping for weak bottom drag. This approach leads to quantitative predictions for the turbulent transport by two-layer baroclinic turbulence with arbitrary layer depths, without additional free parameters. We validate these predictions using an extended suite of numerical simulations with either linear or quadratic bottom drag.

We report numerical simulations of surface gravity waves forced at small scale and the subsequent inverse cascade of wave action. We combine the spectral approach to simulating weakly nonlinear waves with the capabilities of modern Graphics Processing Units to reach unprecedented scale separation between the forcing and domain scales. The resulting broad inertial range allows for an unambiguous confirmation of the theoretical prediction for the spectrum in the inverse cascade regime, both in terms of spectral index and dependence of the spectral level on the action flux., Comment: 5 pages, 3 figures

We consider two-dimensional flows above topography, revisiting the selective decay (or minimum-enstrophy) hypothesis of Bretherton and Haidvogel. We derive a 'condensed branch' of solutions to the variational problem where a domain-scale condensate coexists with a flow at the (smaller) scale of the topography. The condensate arises through a supercritical bifurcation as the conserved energy of the initial condition exceeds a threshold value, a prediction that we quantitatively validate using Direct Numerical Simulations (DNS). We then consider the forced-dissipative case, showing how weak forcing and dissipation select a single dissipative state out of the continuum of solutions to the energy-conserving system predicted by selective decay. As the forcing strength increases, the condensate arises through a supercritical bifurcation for topographic-scale forcing and through a subcritical bifurcation for domain-scale forcing, both predictions being quantitatively validated by DNS. This method provides a way of determining the equilibrated state of forced-dissipative flows based on variational approaches to the associated energy-conserving system, such as the statistical mechanics of 2D flows or selective decay.

We experimentally and numerically characterize rapidly rotating radiatively driven thermal convection, beyond the sole heat transport measurements reported in Bouillaut et al. (2021). Based on a suite of direct numerical simulations (DNS) and additional processing of the experimental data collected by Bouillaut et al. (2021), we report the simultaneous validation of the scaling predictions of the `geostrophic turbulence' regime -- the diffusivity-free or `ultimate' regime of rapidly rotating convection -- for the heat transport, the temperature fluctuations, the flow speed and the flow structure. Radiatively driven convection thus appears as a versatile setup for the laboratory observation of the diffusivity-free regimes of various convective flows of geophysical and/or astrophysical interest.

We revisit the vortex-gas scaling theory for heat transport by baroclinic turbulence based on the empirical observation that the vortex core radius departs from the Rossby deformation radius for very low bottom drag coefficient. We derive a scaling prediction for the vortex-core radius. For linear bottom drag this scaling dependence for the vortex-core radius does not affect the vortex-gas predictions for the eddy diffusivity and mixing-length, which remain identical to those in Gallet and Ferrari (Proc. Nat. Acad. Sci. USA, 117, 2020). By contrast, for quadratic drag the scaling dependence of the core radius induces new scaling-laws for the eddy diffusivity and mixing length when the quadratic-drag coefficient becomes asymptotically low. We validate the modified scaling predictions through numerical simulations of the two-layer model with very low quadratic-drag coefficient.

Ocean mesoscale eddies enhance meridional buoyancy transport, notably in the Antarctic Circumpolar Current where they contribute to setting the deep stratification of the neighboring ocean basins. The much-needed parameterization of this buoyancy transport in global climate models requires a theory for the overall flux, but also for its vertical structure inside the fluid column. Based on the quasi-geostrophic dynamics of an idealized patch of ocean hosting an arbitrary vertically sheared zonal flow, we provide a quantitative prediction for the vertical structure of the buoyancy flux without adjustable parameters. The prediction agrees quantitatively with meridional flux profiles obtained through numerical simulations of an idealized patch of ocean with realistic parameter values. This work empowers modelers with an explicit and physically based expression for the vertical profile of buoyancy transport by ocean baroclinic turbulence, as opposed to the common practice of using arbitrary prescriptions for the depth-dependence of the transport coefficients., Comment: Geophysical Research Letters, in press

The transport induced by ocean mesoscale eddies remains unresolved in most state-of-the-art climate models and needs to be parameterized instead. The natural scale separation between the forcing and the emergent turbulent flow calls for a diffusive parameterization, where the eddy-induced fluxes are related to the large-scale gradients by a diffusion tensor. The standard parameterization scheme in climate modeling consists in adopting the Gent-McWilliams/Redi (GM/R) form for the diffusion tensor, initially put forward based on physical intuition and educated guesses before being put on firm analytical footing using thickness-weighted average (TWA). In the present contribution we provide a direct derivation of this diffusion tensor from the quasi-geostrophic (QG) dynamics of a horizontally homogeneous three-dimensional patch of ocean hosting a large-scale vertically-sheared zonal flow on the beta plane. While less general than the TWA approach, the present QG framework leads to rigorous constraints on the diffusion tensor. First, there is no diapycnal diffusivity arising in the QG GM/R tensor for low viscosity and small-scale diffusivities. The diffusion tensor then involves only two vertically dependent coefficients, namely the GM transport coefficient $K_{GM}(z)$ and the Redi diffusivity $K_R(z)$. Secondly, as already identified by previous authors the vertical structures of the two coefficients are related by the so-called Taylor-Bretherton relation. Finally, while the two coefficients generically differ in the interior of the water column, we show that they are equal to one another near the surface and near the bottom of the domain for low-enough dissipative coefficients. We illustrate these findings by numerically simulating the QG dynamics of a horizontally homogeneous patch of ocean hosting a vertically sheared zonal current resembling the Antarctic Circumpolar Current.

Three-dimensional convection driven by internal heat sources and sinks (CISS) leads to experimental and numerical scaling-laws compatible with a mixing-length - or `ultimate' - scaling regime $Nu \sim \sqrt{Ra}$. However, asymptotic analytic solutions and idealized 2D simulations have shown that laminar flow solutions can transport heat even more efficiently, with $Nu \sim Ra$. The turbulent nature of the flow thus has a profound impact on its transport properties. In the present contribution we give this statement a precise mathematical sense. We show that the Nusselt number maximized over all solutions is bounded from above by const.$\times Ra$, before restricting attention to 'fully turbulent branches of solutions', defined as families of solutions characterized by a finite nonzero limit of the dissipation coefficient at large driving amplitude. Maximization of $Nu$ over such branches of solutions yields the better upper-bound $Nu \lesssim \sqrt{Ra}$. We then provide 3D numerical and experimental data of CISS compatible with a finite limiting value of the dissipation coefficient at large driving amplitude. It thus seems that CISS achieves the maximal heat transport scaling over fully turbulent solutions.

The dynamo driven by the magnetorotational instability (MRI) is believed to play an important role in the dynamics of accretion discs and may also explain the origin of the extreme magnetic fields present in magnetars. Its saturation level is an important open question known to be particularly sensitive to the diffusive processes through the magnetic Prandtl number Pm (the ratio of viscosity to resistivity). Despite its relevance to proto-neutron stars and neutron star merger remnants, the numerically challenging regime of high Pm is still largely unknown. Using zero-net flux shearing box simulations in the incompressible approximation, we studied MRI-driven dynamos at unprecedentedly high values of Pm reaching 256. The simulations show that the stress and turbulent energies are proportional to Pm up to moderately high values ($\mathrm{Pm} \sim 50$). At higher Pm, they transition to a new regime consistent with a plateau independent of Pm for $\rm Pm \gtrsim 100$. This trend is independent of the Reynolds number, which may suggest an asymptotic regime where the energy injection and dissipation are independent of the diffusive processes. Interestingly, large values of Pm not only lead to intense small-scale magnetic fields but also to a more efficient dynamo at the largest scales of the box., Comment: 8 pages, 5 figures, accepted for publication in MNRAS

The competition between turbulent convection and global rotation in planetary and stellar interiors governs the transport of heat and tracers, as well as magnetic-field generation. These objects operate in dynamical regimes ranging from weakly rotating convection to the `geostrophic turbulence' regime of rapidly rotating convection. However, the latter regime has remained elusive in the laboratory, despite a worldwide effort to design ever-taller rotating convection cells over the last decade. Building on a recent experimental approach where convection is driven radiatively, we report heat transport measurements in quantitative agreement with this scaling regime, the experimental scaling-law being validated against direct numerical simulations (DNS) of the idealized setup. The scaling exponent from both experiments and DNS agrees well with the geostrophic turbulence prediction. The prefactor of the scaling-law is greater than the one diagnosed in previous idealized numerical studies, pointing to an unexpected sensitivity of the heat transport efficiency to the precise distribution of heat sources and sinks, which greatly varies from planets to stars.

Despite the increasing sophistication of numerical models of hot Jupiter atmospheres, the large time-scale separation required in simulating the wide range in electrical conductivity between the dayside and nightside has made it difficult to run fully consistent magnetohydrodynamic (MHD) models. This has led to many studies that resort to drag parametrizations of MHD. In this study, we revisit the question of the Lorentz force as an effective drag by running a series of direct numerical simulations of a weakly rotating, poorly conducting flow in the presence of a misaligned, strong background magnetic field. We find that the drag parametrization fails once the time-scale associated with the Lorentz force becomes shorter than the dynamical time-scale in the system, beyond which the effective drag coefficient remains roughly constant, despite orders-of-magnitude variation in the Lorentz (magnetic) time-scale. We offer an improvement to the drag parametrization by considering the relevant asymptotic limit of low conductivity and strong background magnetic field, known as the quasi-static MHD approximation of the Lorentz force. This approximation removes the fast time-scale associated with magnetic diffusion, but retains a more complex version of the Lorentz force, which could be utilized in future numerical models of hot Jupiter atmospheric circulation., Comment: Main text 9 pages, total 11 pages, 4 figures, 1 table. v2: corrections and changes made based on referee reports. Published in The Astrophysical Journal

We numerically and theoretically investigate the Boussinesq Eady model, where a rapidly rotating density-stratified layer of fluid is subject to a meridional temperature gradient in thermal wind balance with a uniform vertically sheared zonal flow. Through a suite of numerical simulations, we show that the transport properties of the resulting turbulent flow are governed by quasi-geostrophic (QG) dynamics in the rapidly rotating strongly stratified regime. The 'vortex gas' scaling predictions put forward in the context of the two-layer QG model carry over to this fully 3D system: the functional dependence of the meridional flux on the control parameters is the same, the two ajustable parameters entering the theory taking slightly different values. In line with the QG prediction, the meridional buoyancy flux is depth-independent. The vertical buoyancy flux is such that turbulence transports buoyancy along isopycnals, except in narrow layers near the the top and bottom boundaries, the thickness of which decreases as the diffusivities go to zero. The emergent (re)stratification is set by a simple balance between the vertical buoyancy flux and diffusion along the vertical direction. Overall, this study demonstrates how the vortex-gas scaling theory can be adapted to quantitatively predict the magnitude and vertical structure of the meridional and vertical buoyancy fluxes, and of the emergent stratification, without additional fitting parameters., Comment: Accepted version

We address magnetic-field generation by dynamo action in systems with inhomogeneous electrical conductivity and magnetic permeability. More specifically, we first show that the Taylor-Couette kinematic dynamo undergoes a drastic reduction of its stability threshold when a (zero-mean) modulation of the fluid's electrical conductivity or magnetic permeability is introduced. These results are obtained outside the mean-field regime, for which this effect was initially proposed. Beyond this illustrative example, we extend a duality argument put forward by Favier \& Proctor {(2013)} to show that swapping the distributions of conductivity and permeability and changing ${ u}\to -{ u}$ leaves the dynamo threshold unchanged. This allows one to make connections between {\it a priori} unrelated dynamo studies. Finally, we discuss the possibility of observing such an effect both in laboratory and astrophysical settings.

The meridional temperature profile of the upper layers of planetary atmospheres is set through a balance between differential radiative heating by a nearby star, or by intrinsic heat fluxes emanating from the deep interior, and the redistribution of that heat across latitudes by turbulent flows. These flows spontaneously arise through baroclinic instability of the meridional temperature gradients maintained by the forcing. When planetary curvature is neglected, this turbulence takes the form of coherent vortices that mix the meridional temperature profiles. However, the curvature of the planet favors the emergence of Rossby waves and zonal jets that restrict the meridional wandering of the fluid columns, thereby reducing the mixing efficiency across latitudes. A similar situation arises in the ocean, where the baroclinic instability of zonal currents leads to enhanced meridional heat transport by a turbulent flow consisting of vortices and zonal jets. A recent scaling theory for the turbulent heat transport by vortices is extended to include the impact of planetary curvature, in the framework of the two-layer quasigeostrophic beta-plane model. This leads to a quantitative parameterization providing the meridional temperature profile in terms of the externally imposed heat flux in an idealized model of planetary atmospheres and oceans. In addition, it provides a quantitative prediction for the emergent criticality, i.e., the degree of instability in a canonical model of planetary atmosphere or ocean., Comment: Accepted for publication in AGU Advances

Astrophysical plasmas are often subject to both rotation and large-scale background magnetic fields. Individually, each is known to two-dimensionalize the flow in the perpendicular plane. In realistic flows, both of these effects are simultaneously present and, importantly, need not be aligned. In this work, we numerically investigate three-dimensional forced magnetohydrodynamic (MHD) turbulence subject to the competing effects of global rotation and a perpendicular background magnetic field. We focus on the case of a strong background field and find that increasing the rotation rate from zero produces significant changes in the structure of the turbulent flow. Starting with a two-dimensional inverse cascade at zero rotation, the flow first transitions to a forward cascade of kinetic energy, then to a shear-layer dominated regime, and finally to a second shear-layer regime where the kinetic energy flux is strongly suppressed and the energy transfer is mediated by the induced magnetic field. We show that the first two transitions occur at distinct values of the Rossby number, and the third occurs at a distinct value of the Lehnert number. More generally, our results demonstrate that, when considering the simultaneous limits of strong rotation and a strong guide field, the order in which those limits are taken matters in the misaligned case., Comment: 17 pages, 7 figures

The absorption of light or radiation drives turbulent convection inside stars, supernovae, frozen lakes and the Earth's mantle. In these contexts, the goal of laboratory and numerical studies is to determine the relation between the internal temperature gradients and the heat flux transported by the turbulent flow. This is the constitutive law of turbulent convection, to be input into large-scale models of such natural flows. However, in contrast with the radiative heating of natural flows, laboratory experiments have focused on convection driven by heating and cooling plates: the heat transport is then severely restricted by boundary layers near the plates, which prevents the realization of the mixing-length scaling-law used in evolution models of geophysical and astrophysical flows. There is therefore an important discrepancy between the scaling-laws measured in laboratory experiments and those used eg. in stellar evolution models. Here we provide experimental and numerical evidence that radiatively driven convection spontaneously achieves the mixing-length scaling regime, also known as the `ultimate' regime of thermal convection. This constitutes the first clear observation of this regime of turbulent convection. Our study therefore bridges the gap between models of natural flows and laboratory experiments. It opens an experimental avenue for a priori determinations of the constitutive laws to be implemented into models of geophysical and astrophysical flows, as opposed to empirical fits of these constitutive laws to the scarce observational data.

We report the quantitative experimental observation of the weak inertial-wave turbulence regime of rotating turbulence. We produce a statistically steady homogeneous turbulent flow that consists of nonlinearly interacting inertial waves, using rough top and bottom boundaries to prevent the emergence of a geostrophic flow. As the forcing amplitude increases, the temporal spectrum evolves from a discrete set of peaks to a continuous spectrum. Maps of the bicoherence of the velocity field confirm such a gradual transition between discrete wave interactions at weak forcing amplitude, and the regime described by weak turbulence theory (WTT) for stronger forcing. In the former regime, the bicoherence maps display a near-zero background level, together with sharp localized peaks associated with discrete resonances. By contrast, in the latter regime the bicoherence is a smooth function that takes values of the order of the Rossby number, in line with the infinite-domain and random-phase assumptions of WTT. The spatial spectra then display a power-law behavior, both the spectral exponent and the spectral level being accurately predicted by WTT at high Reynolds number and low Rossby number., Comment: Accepted for publication in Physical Review Letters

The mean state of the atmosphere and ocean is set through a balance between external forcing (winds, radiation, heat and freshwater fluxes) and the emergent turbulence, which transfers energy to dissipative structures. The forcing gives rise to jets in the atmosphere and currents in the ocean, which spontaneously develop turbulent eddies through the baroclinic instability. A critical step in the development of a theory of climate is to properly include the resulting eddy-induced turbulent transport of properties like heat, moisture, and carbon. The baroclinic instability generates flow structures at the Rossby deformation radius, a length scale of order 1000 km in the atmosphere and 100 km in the ocean, smaller than the planetary scale and the typical extent of ocean basins respectively. There is therefore a separation of scales between the large-scale temperature gradient and the smaller eddies that advect it randomly, inducing effective diffusion. Numerical solutions of the two-layer quasi-geostrophic model, the standard model for studies of eddy motions in the atmosphere and ocean, show that such scale separation remains in the strongly nonlinear turbulent regime, provided there is sufficient bottom drag. We compute the scaling-laws governing the eddy-driven transport associated with baroclinic turbulence. First, we provide a theoretical underpinning for empirical scaling-laws reported in previous studies, for different formulations of the bottom drag law. Second, these scaling-laws are shown to provide an important first step toward an accurate local closure to predict the impact of baroclinic turbulence in setting the large-scale temperature profiles in the atmosphere and ocean.

Turbulent flows driven by a vertically invariant body force were proven to become exactly two-dimensional above a critical rotation rate, using upper bound theory. This transition in dimensionality of a turbulent flow has key consequences for the energy dissipation rate. However, its location in parameter space is not provided by the bounding procedure. To determine this precise threshold between exactly 2D and partially 3D flows, we perform a linear stability analysis over a fully turbulent 2D base state. This requires integrating numerically a quasi-2D set of equations over thousands of turnover times, to accurately average the growth rate of the 3D perturbations over the statistics of the turbulent 2Dbase flow. We leverage the capabilities of modern GPUs to achieve this task, which allows us to investigate the parameter space up to Re = 10^5. At Reynolds numbers typical of 3D DNS and laboratory experiments, Re in [10^2, 5x10^3], the turbulent 2D flow becomes unstable to 3D motion through a centrifugal-type instability. However, at even higher Reynolds number another instability takes over. A candidate mechanism for the latter instability is the parametric excitation of inertial waves by the modulated 2D flow, a phenomenon that we illustrate with an oscillatory 2D Kolmogorov flow.

We report on a numerical study of turbulent convection driven by a combination of internal heat sources and sinks. Motivated by a recent experimental realisation (Lepot et al. 2018), we focus on the situation where the cooling is uniform, while the internal heating is localised near the bottom boundary, over approximately one tenth of the domain height. We obtain scaling laws $Nu \sim Ra^\gamma Pr^\chi$ for the heat transfer as measured by the Nusselt number $Nu$ expressed as a function of the Rayleigh number $Ra$ and the Prandtl number $Pr$. After confirming the experimental value $\gamma\approx 1/2$ for the dependence on $Ra$, we identify several regimes of dependence on $Pr$. For a stress-free bottom surface and within a range as broad as $Pr \in [0.003, 10]$, we observe the exponent $\chi\approx 1/2$, in agreement with Spiegel's mixing length theory. For a no-slip bottom surface we observe a transition from $\chi\approx 1/2$ for $Pr \leq 0.04$ to $\chi\approx 1/6$ for $Pr \geq 0.04$, in agreement with scaling predictions by Bouillaut et al. The latter scaling regime stems from heat accumulation in the stagnant layer adjacent to a no-slip bottom boundary, which we characterise by comparing the local contributions of diffusive and convective thermal fluxes., Comment: 11 pages, 3 figures

We report on laboratory experiments of wave-driven rotating turbulence. A set of wavemakers produces inertial-wave beams that interact nonlinearly in the central region of a water tank mounted on a rotating platform. The forcing thus injects energy into inertial waves only. For moderate forcing amplitude, part of the energy of the forced inertial waves is transferred to subharmonic waves, through a standard triadic resonance instability. This first step is broadly in line with the theory of weak turbulence. Surprisingly however, stronger forcing does not lead to an inertial-wave turbulence regime. Instead, most of the kinetic energy condenses into a vertically invariant geostrophic flow, even though the latter is unforced. We show that resonant quartets of inertial waves can trigger an instability -- the ``quartetic instability'' -- that leads to such spontaneous emergence of geostrophy. In the present experiment, this instability sets in as a secondary instability of the classical triadic instability., Comment: Accepted for publication in Physical Review Letters

We explore the near-resonant interaction of inertial waves with geostrophic modes in rotating fluids via numerical and theoretical analysis. When a single inertial wave is imposed, we find that some geostrophic modes are unstable above a threshold value of the Rossby number $kRo$ based on the wavenumber and wave amplitude. We show this instability to be caused by triadic interaction involving two inertial waves and a geostrophic mode such that the sum of their eigen frequencies is non-zero. We derive theoretical scalings for the growth rate of this near-resonant instability. The growth rate scaled by the global rotation rate is proportional to $(kRo)^2$ at low $kRo$ and transitions to a $kRo$ scaling for larger $kRo$. These scalings are in excellent agreement with direct numerical simulations. This instability could explain recent experimental observations of geostrophic instability driven by waves., Comment: Accepted for publication in Journal of Fluid Mechanics Rapids. 13 pages, 4 figures

We report on an instability arising when surface gravity waves propagate in a rotating frame. The Stokes drift associated to the uniform wave field, together with global rotation, drives a mean flow in the form of a horizontally invariant Ekman-Stokes spiral. We show that the latter can be subject to an instability that triggers the appearance of an additional horizontally-structured cellular flow. We determine the instability threshold numerically, in terms of the Rossby number Ro associated to the Stokes drift of the waves and the Ekman number E. We confirm the numerical results through asymptotic expansions at both large and low Ekman number. At large E the instability reduces to that of a standard Ekman spiral driven by the wave-induced surface stress instead of a wind stress, while at low E the Stokes-drift profile crucially determines the shape of the unstable mode. In both limits the instability threshold asymptotes to an Ekman-number-independent critical Rossby number, which in both cases also corresponds to a critical Reynolds number associated to the Lagrangian base-flow velocity profile. Parameter values typical of ocean swell fall into the low-E unstable regime: the corresponding "anti-Stokes" flows are unstable, with possible consequences for particle dispersion and mixing.

We present analytical examples of fluid dynamos that saturate through the action of the Coriolis and inertial terms of the Navier-Stokes equation. The flow is driven by a body force and is subject to global rotation and uniform sweeping velocity. The model can be studied down to arbitrarily low viscosity and naturally leads to the strong-field scaling regime for the magnetic energy produced above threshold: the magnetic energy is proportional to the global rotation rate and independent of the viscosity. Depending on the relative orientations of global rotation and large-scale sweeping, the dynamo bifurcation is either supercritical or subcritical. In the supercritical case, the magnetic energy follows the scaling-law for supercritical strong-field dynamos predicted on dimensional grounds by Petrelis & Fauve (2001). In the subcritical case, the system jumps to a finite-amplitude dynamo branch. The magnetic energy obeys a magneto-geostrophic scaling-law (Roberts & Soward 1972), with a turbulent Elsasser number of the order of unity, where the magnetic diffusivity of the standard Elsasser number appears to be replaced by an eddy diffusivity. In the absence of global rotation, the dynamo bifurcation is subcritical and the saturated magnetic energy obeys the equipartition scaling regime. We consider both the vicinity of the dynamo threshold and the limit of large distance from threshold to put these various scaling behaviors on firm analytical ground.

While the saturated magnetic energy is independent of viscosity in dynamo experiments, it remains viscosity-dependent in state-of-the-art 3D direct numerical simulations (DNS). Extrapolating such viscous scaling-laws to realistic parameter values leads to an underestimation of the magnetic energy by several orders of magnitude. The origin of this discrepancy is that fully 3D DNS cannot reach low enough values of the magnetic Prandtl number $Pm$. To bypass this limitation and investigate dynamo saturation at very low $Pm$, we focus on the vicinity of the dynamo threshold in a rapidly rotating flow: the velocity field then depends on two spatial coordinates only, while the magnetic field consists of a single Fourier mode in the third direction. We perform numerical simulations of the resulting set of reduced equations for $Pm$ down to $2\cdot 10^{-5}$. This parameter regime is currently out of reach to fully 3D DNS. We show that the magnetic energy transitions from a high-$Pm$ viscous scaling regime to a low-$Pm$ turbulent scaling regime, the latter being independent of viscosity. The transition to the turbulent saturation regime occurs at a low value of the magnetic Prandtl number, $Pm \simeq 10^{-3}$, which explains why it has been overlooked by numerical studies so far.

What is the turbulent drag force experienced by an object moving in a rotating fluid? This open and fundamental question can be addressed by measuring the torque needed to drive an impeller at constant angular velocity $\omega$ in a water tank mounted on a platform rotating at a rate $\Omega$. We report a dramatic reduction in drag as $\Omega$ increases, down to values as low as $12$\% of the non-rotating drag. At small Rossby number $Ro = \omega/\Omega$, the decrease in drag coefficient $K$ follows the approximate scaling law $K \sim Ro$, which is predicted in the framework of nonlinear inertial wave interactions and weak-turbulence theory. However, stereoscopic particle image velocimetry measurements indicate that this drag reduction rather originates from a weakening of the turbulence intensity in line with the two-dimensionalization of the large-scale flow., Comment: To appear in Journal of Fluid Mechanics Rapids

We consider the flow of a Newtonian fluid in a three-dimensional domain, rotating about a vertical axis and driven by a vertically invariant horizontal body-force. This system admits vertically invariant solutions that satisfy the 2D Navier-Stokes equation. At high Reynolds number and without global rotation, such solutions are usually unstable to three-dimensional perturbations. By contrast, for strong enough global rotation, we prove rigorously that the 2D (and possibly turbulent) solutions are stable to vertically dependent perturbations: the flow becomes 2D in the long-time limit. These results shed some light on several fundamental questions of rotating turbulence: for arbitrary Reynolds number and small enough Rossby number, the system is attracted towards purely 2D flow solutions, which display no energy dissipation anomaly and no cyclone-anticyclone asymmetry. Finally, these results challenge the applicability of wave turbulence theory to describe stationary rotating turbulence in bounded domains., Comment: To be published in Journal of Fluid Mechanics

This work addresses the question of the stability of stratified, spatially periodic shear flows at low P\'eclet number but high Reynolds number. This little-studied limit is motivated by astrophysical systems, where the Prandtl number is often very small. Furthermore, it can be studied using a reduced set of "low-P\'eclet-number equations" proposed by Lignieres [Astronomy & Astrophysics, 348, 933-939, 1999]. Through a linear stability analysis, we first determine the conditions for instability to infinitesimal perturbations. We formally extend Squire's theorem to the low-P\'eclet-number equations, which shows that the first unstable mode is always two-dimensional. We then perform an energy stability analysis of the low-P\'eclet-number equations and prove that for a given value of the Reynolds number, above a critical strength of the stratification, any smooth periodic shear flow is stable to perturbations of arbitrary amplitude. In that parameter regime, the flow can only be laminar and turbulent mixing does not take place. Finding that the conditions for linear and energy stability are different, we thus identify a region in parameter space where finite-amplitude instabilities could exist. Using direct numerical simulations, we indeed find that the system is subject to such finite-amplitude instabilities. We determine numerically how far into the linearly stable region of parameter space turbulence can be sustained., Comment: To be published in Physics of Fluids

We investigate the behavior of flows, including turbulent flows, driven by a horizontal body-force and subject to a vertical magnetic field, with the following question in mind: for very strong applied magnetic field, is the flow mostly two-dimensional, with remaining weak three-dimensional fluctuations, or does it become exactly 2D, with no dependence along the vertical? We first focus on the quasi-static approximation, i.e. the asymptotic limit of vanishing magnetic Reynolds number Rm << 1: we prove that the flow becomes exactly 2D asymptotically in time, regardless of the initial condition and provided the interaction parameter N is larger than a threshold value. We call this property "absolute two-dimensionalization": the attractor of the system is necessarily a (possibly turbulent) 2D flow. We then consider the full-magnetohydrodynamic equations and we prove that, for low enough Rm and large enough N, the flow becomes exactly two-dimensional in the long-time limit provided the initial vertically-dependent perturbations are infinitesimal. We call this phenomenon "linear two-dimensionalization": the (possibly turbulent) 2D flow is an attractor of the dynamics, but it is not necessarily the only attractor of the system. Some 3D attractors may also exist and be attained for strong enough initial 3D perturbations. These results shed some light on the existence of a dissipation anomaly for magnetohydrodynamic flows subject to a strong external magnetic field., Comment: Journal of Fluid Mechanics, in press

In field theory, particles are waves or excitations that propagate on the fundamental state. In experiments or cosmological models one typically wants to compute the out-of-equilibrium evolution of a given initial distribution of such waves. Wave Turbulence deals with out-of-equilibrium ensembles of weakly nonlinear waves, and is therefore well-suited to address this problem. As an example, we consider the complex Klein-Gordon equation with a Mexican-hat potential. This simple equation displays two kinds of excitations around the fundamental state: massive particles and massless Goldstone bosons. The former are waves with a nonzero frequency for vanishing wavenumber, whereas the latter obey an acoustic dispersion relation. Using wave turbulence theory, we derive wave kinetic equations that govern the coupled evolution of the spectra of massive and massless waves. We first consider the thermodynamic solutions to these equations and study the wave condensation transition, which is the classical equivalent of Bose-Einstein condensation. We then focus on nonlocal interactions in wavenumber space: we study the decay of an ensemble massive particles into massless ones. Under rather general conditions, these massless particles accumulate at low wavenumber. We study the dynamics of waves coexisting with such a strong condensate, and we compute rigorously a nonlocal Kolmogorov-Zakharov solution, where particles are transferred non-locally to the condensate, while energy cascades towards large wave numbers through local interactions. This nonlocal cascading state constitute the intermediate asymptotics between the initial distribution of waves and the thermodynamic state reached in the long-time limit.

We present a spatio-temporal analysis of a statistically stationary rotating turbulence experiment, aiming to extract a signature of inertial waves, and to determine the scales and frequencies at which they can be detected. The analysis uses two-point spatial correlations of the temporal Fourier transform of velocity fields obtained from time-resolved stereoscopic particle image velocimetry measurements in the rotating frame. We quantify the degree of anisotropy of turbulence as a function of frequency and spatial scale. We show that this space-time-dependent anisotropy is well described by the dispersion relation of linear inertial waves at large scale, while smaller scales are dominated by the sweeping of the waves by fluid motion at larger scales. This sweeping effect is mostly due to the low-frequency quasi-two-dimensional component of the turbulent flow, a prominent feature of our experiment which is not accounted for by wave turbulence theory. These results question the relevance of this theory for rotating turbulence at the moderate Rossby numbers accessible in laboratory experiments, which are relevant to most geophysical and astrophysical flows.

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 two-point third-order three-component 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 three-dimensional 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 two-dimensional turbulence. At the largest rotation rate the flow is nearly two-dimensional, and a pure inverse energy cascade is found for the horizontal energy. To describe the scale-by-scale energy budget, we consider a generalization of the K\'arm\'an-Howarth-Monin 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., Comment: to appear in Physics of Fluids

Using recordings of swell from pitch-and-roll buoys, we have reproduced the classic observations of long-range surface wave propagation originally made by Munk et al. (1963) using a triangular array of bottom pressure measurements. In the modern data, the direction of the incoming swell fluctuates by about $\pm 10^\circ$ on a time scale of one hour. But if the incoming direction is averaged over the duration of an event then, in contrast with the observations by Munk et al. (1963), the sources inferred by great-circle backtracking are most often in good agreement with the location of large storms on weather maps of the Southern Ocean. However there are a few puzzling failures of great-circle backtracking e.g., in one case, the direct great-circle route is blocked by the Tuamoto Islands and the inferred source falls on New Zealand. Mirages like this occur more frequently in the bottom-pressure observations of Munk et al. (1963), where several inferred sources fell on the Antarctic continent. Using spherical ray tracing we investigate the hypothesis that the refraction of waves by surface currents produces the mirages. With reconstructions of surface currents inferred from satellite altimetry, we show that mesoscale vorticity significantly deflects swell away from great-circle propagation so that the source and receiver are connected by a bundle of many rays, none of which precisely follow a great circle. The $\pm 10^\circ$ directional fluctuations at the receiver result from the arrival of wave packets that have travelled along the different rays within this multipath. The occasional failure of great-circle backtracking, and the associated mirages, probably results from partial topographic obstruction of the multipath, which biases the directional average at the receiver., Comment: Journal of Marine Research, in press

We characterize the statistical and geometrical properties of the cyclone-anticyclone asymmetry in a statistically-steady forced rotating turbulence experiment. Turbulence is generated by a set of vertical flaps which continuously inject velocity fluctuations towards the center of a tank mounted on a rotating platform. We first characterize the cyclone-anticyclone asymmetry from conventional single-point vorticity statistics. We propose a phenomenological model to explain the emergence of the asymmetry in the experiment, from which we predict scaling laws for the root-mean-square velocity in good agreement with the experimental data. We further quantify the cyclone-anticyclone asymmetry using a set of third-order two-point velocity correlations. We focus on the correlations which are nonzero only if the cyclone-anticyclone symmetry is broken. They offer two advantages over single-point vorticity statistics: first, they are defined from velocity measurements only, so an accurate resolution of the Kolmogorov scale is not required; second, they provide information on the scale-dependence of the cyclone-anticyclone asymmetry. We compute these correlation functions analytically for a random distribution of independent identical vortices. These model correlations describe well the experimental ones, indicating that the cyclone-anticyclone asymmetry is dominated by the large-scale long-lived cyclones., Comment: Accepted for publication in Physics of Fluids

We present a study of several systems in which a large scale field is generated over a turbulent background. These large scale fields usually break a symmetry of the forcing by selecting a direction. Under certain conditions, the large scale field displays reversals so that the symmetry of the forcing is recovered statistically. We present examples of such dynamics in the context of the dynamo instability, of two dimensional turbulent Kolmogorov flows and of turbulent Rayleigh-B\'enard convection. In these systems reversals occur respectively for the dynamo magnetic field, for the large scale circulation generated by a periodic forcing in space and for the large scale roll generated by turbulent thermal convection. We compare the mechanisms involved and show that their properties depend on some symmetries of the system and on the way they are broken.

In a mixture of two species of grains of equal size but different mass, placed in a vertically vibrated shallow box, there is spontaneous segregation. Once the system is at least partly segregated and clusters of the heavy particles have formed, there are sudden peaks of the horizontal kinetic energy of the heavy particles, that is otherwise small. Together with the energy peaks the clusters rapidly expand and the segregation is partially lost. The process repeats once segregation has taken place again. Depending on the experimental or numerical parameters, the energy bursts can occur either randomly or with some regularity in time. An explanation for these events is provided based on the existence of a fixed point for an isolated particle bouncing with only vertical motion. The horizontal energy peaks occur when the energy stored in the vertical motion is partly transferred into horizontal energy through a chain reaction of collisions between heavy particles. A necessary condition for the observed regularity of the events is that chain reactions involve most of the heavy particles., Comment: 4 pages, 3 figures, Physical Review Letters (accepted)

We consider the effect of radial fluid injection and suction on Taylor-Couette flow. Injection at the outer cylinder and suction at the inner cylinder generally results in a linearly unstable steady spiralling flow, even for cylindrical shears that are linearly stable in the absence of a radial flux. We study nonlinear aspects of the unstable motions with the energy stability method. Our results, though specialized, may have implications for drag reduction by suction, accretion in astrophysical disks, and perhaps even in the flow in the earth's polar vortex., Comment: 34 pages, 9 figures

We report an experimental investigation on the influence of an external magnetic field on forced 3D turbulence of liquid gallium in a closed vessel. We observe an exponential damping of the turbulent velocity fluctuations as a function of the interaction parameter N (ratio of Lorentz force over inertial terms of the Navier-Stokes equation). The flow structures develop some anisotropy but do not become bidimensional. From a dynamical viewpoint, the damping first occurs homogeneously over the whole spectrum of frequencies. For larger values of N, a very strong additional damping occurs at the highest frequencies. However, the injected mechanical power remains independent of the applied magnetic field. The simultaneous measurement of induced magnetic field and electrical potential differences shows a very weak correlation between magnetic field and velocity fluctuations. The observed reduction of the fluctuations is in agreement with a previously proposed mechanism for the saturation of turbulent dynamos and with the order of magnitude of the Von Karman Sodium dynamo magnetic field., Comment: 34 pages, 20 figures

We show that hemispherical dynamos can result from weak equatorial symmetry breaking of the flow in the interior of planets and stars. Using a model of spherical dynamo, we observe that the interaction between a dipolar and a quadrupolar mode can localize the magnetic field in only one hemisphere when the equatorial symmetry is broken. This process is shown to be related to the one that is responsible for reversals of the magnetic field. These seemingly very different behaviors are thus understood in a unified framework.

The magnetic field of planets or stars is generated by the motion of a conducting fluid through a dynamo instability. The saturation of the magnetic field occurs through the reaction of the Lorentz force on the flow. In relation to this phenomenon, we study the effect of a magnetic field on a turbulent flow of liquid Gallium. The measurement of electric potential differences provides a signal related to the local velocity fluctuations. We observe a reduction of velocity fluctuations at all frequencies in the spectrum when the magnetic field is increased., Comment: accepted for Physical Review E