Search

We consider inertial waves propagating in a fluid contained in a non-axisymmetric three-dimensional rotating cavity. We focus on the particular case of a fluid enclosed inside a truncated cone or frustum, which is the volume that lies between two horizontal parallel planes cutting an upright cone. While this geometry has been studied in the past, we generalise it by breaking its axisymmetry and consider the case of a truncated elliptical cone for which the horizontal sections are elliptical instead of circular. The problem is first tackled using ray tracing where local wave packets are geometrically propagated and reflected within the closed volume without attenuation. We complement these results with a local asymptotic analysis and numerical simulations of the original linear viscous problem. We show that the attractors, well-known in two dimensional or axisymmetric domains, can be trapped in a particular plane in three-dimension provided that the axisymmetry of the domain is broken. Contrary to previous examples of attractors in three-dimensional domains, all rays converge towards the same limit cycle regardless of initial conditions, and it is localised in the bulk of the fluid., Comment: 22 pages, 11 figures

Following our previous work on periodic ray paths (He et al, 2022), we study asymptotically and numerically the structure of internal shear layers for very small Ekman numbers in a three-dimensional (3D) spherical shell and in a two-dimensional (2D) cylindrical annulus when the rays converge towards an attractor. We first show that the asymptotic solution obtained by propagating the self-similar solution generated at the critical latitude on the librating inner core describes the main features of the numerical solution. The internal shear layer structure and the scaling for its width and velocity amplitude in $E^{1/3}$ and $E^{1/12}$ respectively are recovered. The amplitude of the asymptotic solution is shown to decrease to $E^{1/6}$ when it reaches the attractor, as it is also observed numerically. However, some discrepancies are observed close to the particular attractors along which the phase of the wave beam remains constant. Another asymptotic solution close to those attractors is then constructed using the model of Ogilvie (2005). The solution obtained for the velocity has an $O(E^{1/6})$ amplitude, but a different self-similar structure than the critical-latitude solution. It also depends on the Ekman pumping at the contact points of the attractor with the boundaries. We demonstrate that it reproduces correctly the numerical solution. Surprisingly, the solution close to an attractor with phase shift (that is an attractor that touches the axis in 3D or in 2D with a symmetric forcing) is found to be much weaker., Comment: 31 pages, 13 figures

The heat transport by rapidly-rotating Rayleigh-B\'enard convection is of fundamental importance to many geophysical flows. Laboratory measurements are impeded by robust wall modes which develop along vertical walls, significantly perturbing the heat flux. We show that narrow horizontal fins along the vertical walls efficiently suppress wall modes ensuring that their contribution to the global heat flux is negligible compared with bulk convection in the geostrophic regime, thereby paving the way for new experimental studies of geophysically relevant regimes of rotating convection., Comment: 6 pages, 6 figures

We experimentally investigate the behaviour of instantaneous localised releases of heavy particles falling as turbulent clouds in quiescent water, both with and without background rotation. We present the results of 514 systematic experiments for no rotation and for three rotation rates $\Omega = 5,10,20$rpm, and for the size of particles in the range $5\mu m$ to $1mm$, exploring four decades of the Rouse number $\mathcal{R} \in [6\times 10^{-4},4]$ which quantifies the inertia of particles. In the canonical framework of turbulent thermals described by Morton \textit{et al.}, [Proc. R. Soc. A: Math. Phys. Sci. \textbf{234}, 1 (1956)], we compare particle clouds with salt-water thermals to highlight specificities due to the particulate nature of the turbulence forcing. In the absence of rotation, particle clouds initially behave as salty thermals with a modulation of their entrainment capacity, which is optimally enhanced for a finite inertia $\mathcal{R}\simeq 0.3$ due to particulate effects. However this regime of turbulence is limited in time due to the inertial decoupling between turbulent eddies and particles. For the three values of $\Omega$ explored here, the particulate enhancement of entrainment is inhibited. Moreover the cloud's expansion is interrupted when the Coriolis force overcomes its inertia, forcing the cloud to transform into vortical columnar flows which considerably increase the residence time of particles.

Intense east-west winds called zonal jets are observed in the atmospheres of Jupiter and Saturn and extend in their deep interior. We present experimental results from a fully three-dimensional laboratory analog of deep gas giants zonal jets. We use a rapidly rotating deep cylindrical tank, filled with water, and forced by a small-scale hydraulic circulation at the bottom. A topographic $\beta$-effect is naturally present because of the curvature of the free surface. Instantaneous turbulent zonal jets spontaneously emerge from the small-scale forcing, equilibrate at large scale, and can contain up to 70% of the total kinetic energy of the flow once in a quasi-steady state. We show that the spectral properties of the experimental flows are consistent with the theoretical predictions in the zonostrophic turbulence regime, argued to be relevant to gas giants. This constitutes the first fully-experimental validation of the zonostrophic theory in a completely three-dimensional framework. Complementary, quasi-geostrophic (QG) simulations show that this result is not sensitive to the forcing scale. Next, we quantify the potential vorticity (PV) mixing. While PV staircasing should emerge in the asymptotic regime of the gas giants, only a moderate PV mixing occurs because of the strong forcing and dissipation, as confirmed by QG simulations at smaller Ekman number. We quantify the local PV mixing by measuring the equivalent of a Thorpe scale, and confirm that it can be used to estimate the upscale energy transfer rate of the flow, which otherwise needs to be estimated from a much more demanding spectral analysis., Comment: 40 pages, 12 figures, 6 tables

This article presents a numerical analysis of the instability developing in horizontally sheared Poiseuille flow, when stratification extends along the vertical direction. Our study builds up on the previous work that originally detected the linear instability of such configuration, by means of experiments, theoretical analysis and numerical simulations \citep{G21}. We extend hereafter this former investigation beyond linear theory, investigating nonlinear regimes with direct numerical simulations. We notice that the flow loses its vertical homogeneity through a secondary bifurcation, due to harmonic resonances, and further describe this symmetry-breaking mechanism in the vicinity of the instability threshold. When departing away from this limit, we observe a series of bursting events that break down the flow into disordered motions driven by localized shear instabilities. This intermittent dynamics leads to the coexistence of horizontal localized layers of stratified turbulence surrounded by quiescent regions of meandering waves.

We investigate the dynamics of a fluid layer subject to an imposed bottom heat flux and a top monotonically-increasing temperature profile driving horizontal convection. We use direct numerical simulations and consider a large range of flux-based Rayleigh numbers $10^6 \leq Ra_F \leq 10^9$ and imposed top horizontal to bottom vertical heat flux ratios $0 \leq \Lambda \leq 1$. The fluid domain is a closed two-dimensional box with aspect ratio $4\leq \Gamma \leq 16$ and we consider no-slip boundaries and adiabatic side walls. We demonstrate a regime transition from Rayleigh--B\'enard convection (RB) to horizontal convection (HC) at $\Lambda\approx 10^{-2}$, which is independent of $Ra_F$ and $\Gamma$. At small $\Lambda$, the flow is organized in multiple overturning cells with approximately unit aspect ratio, while at large $\Lambda$ a single cell is obtained. The RB-relevant Nusselt number scaling with $Ra_F$ and the HC-relevant Nusselt number scaling with the horizontal Rayleigh number $Ra_L=Ra_F\Lambda\Gamma^4$ are in good agreement with previous results from classical RB convection and HC studies in the limit $\Lambda \ll 10^{-2}$ and $\Lambda \gg 10^{-2}$, respectively. We demonstrate that the system is multi-stable near the transition $\Lambda\approx10^{-2}$, i.e. the exact number of cells not only depends on $\Lambda$ but also on the system's history. Our results suggest that subglacial lakes, which motivated this study, are likely to be dominated by RB convection, unless the slope of the ice-water interface, which controls the horizontal temperature gradient via the pressure-dependence of the freezing point, is greater than unity., Comment: Accepted for publication in Journal of Fluid Mechanics. 30 pages, 15 figures

Instabilities of fluid flows often generate turbulence. Using extensive direct numerical simulations, we study two-dimensional turbulence driven by a wavenumber-localised instability superposed on stochastic forcing, in contrast to previous studies of state-independent forcing. As the contribution of the instability forcing, measured by a parameter $\gamma$, increases, the system undergoes two transitions. For $\gamma$ below a first threshold, a regular large-scale vortex condensate forms. Above this threshold, shielded vortices (SVs) emerge within the condensate. At a second, larger value of $\gamma$, the condensate breaks down, and a gas of weakly interacting vortices with broken symmetry spontaneously emerges, characterised by preponderance of vortices of one sign only and suppressed inverse energy cascade. The latter transition is shown to depend on the damping mechanism. The number density of SVs in the broken symmetry state slowly increases via a random nucleation process. Bistability is observed between the condensate and mixed SV-condensate states. Our findings provide new evidence for a strong dependence of two-dimensional turbulence phenomenology on the forcing.

In marine plankton, many swimming species can perceive their environment with flow sensors. Can they use this flow information to travel faster in turbulence? To address this question, we consider plankters swimming at constant speed, whose goal is to move upwards. We propose a robust analytical behavior that allows plankters to choose a swimming direction according to the local flow gradients. We show numerically that such plankters can "surf" on turbulence and reach net vertical speeds up to twice their swimming speed. This new physics-based model suggests that planktonic organisms can exploit turbulence features for navigation., Comment: Article: 5 pages, 5 figures. Supplemental: 11 pages, 6 figures. For associated open-source code, see https://github.com/C0PEP0D/sheld0n

Recent photometric observations of massive stars show ubiquitous low-frequency "red-noise" variability, which has been interpreted as internal gravity waves (IGWs). Simulations of IGWs generated by convection show smooth surface wave spectra, qualitatively matching the observed red-noise. On the other hand, theoretical calculations by Shiode et al (2013) and Lecoanet et al (2019) predict IGWs should manifest at the surface as regularly-spaced peaks associated with standing g-modes. In this work, we compare these theoretical approaches to simplified 2D numerical simulations. The simulations show g-mode peaks at their surface, and are in good agreement with Lecoanet et al (2019). The amplitude estimates of Shiode et al (2013) did not take into account the finite width of the g-mode peaks; after correcting for this finite width, we find good agreement with simulations. However, simulations need to be run for hundreds of convection turnover times for the peaks to become visible; this is a long time to run a simulation, but a short time in the life of a star. The final spectrum can be predicted by calculating the wave energy flux spectrum in much shorter simulations, and then either applying the theory of Shiode et al (2013) or Lecoanet et al (2019)., Comment: submitted to MNRAS

We consider the problem of an inextensible but flexible fiber advected by a steady chaotic flow, and ask the simple question whether the fiber can spontaneously knot itself. Using a 1D Cosserat model, a simple local viscous drag model and discrete contact forces, we explore the probability of finding knots at any given time when the fiber is interacting with the ABC class of flows. The bending rigidity is shown to have a marginal effect compared to that of increasing the fiber length. Complex knots are formed up to 11 crossings, but some knots are more probable than others. The finite-time Lyapunov exponent of the flow is shown to have a positive effect on the knot probability. Finally, contact forces appear to be crucial since knotted configurations can remain stable for times much longer than the turnover time of the flow, something that is not observed when the fiber can freely cross itself., Comment: 13 pages, 8 figures, accepted for publication in Phys. Rev. E

This paper is associated with a poster winner of a 2019 American Physical Society Division of Fluid Dynamics (DFD) Milton van Dyke Award for work presented at the DFD Gallery of Fluid Motion. The original poster is available online at the Gallery of Fluid Motion, https://doi.org/10.1103/APS.DFD.2019.GFM.P0015.

Jupiter's dynamics shapes its cloud patterns but remains largely unknown below this natural observational barrier. Unraveling the underlying three-dimensional flows is thus a primary goal for NASA's ongoing Juno mission that was launched in 2011. Here, we address the dynamics of large Jovian vortices using laboratory experiments complemented by theoretical and numerical analyses. We determine the generic force balance responsible for their three-dimensional pancake-like shape. From this, we define scaling laws for their horizontal and vertical aspect ratios as a function of the ambient rotation, stratification and zonal wind velocity. For the Great Red Spot in particular, our predicted horizontal dimensions agree well with measurements at the cloud level since the Voyager mission in 1979. We additionally predict the Great Red Spot's thickness, inaccessible to direct observation: it has surprisingly remained constant despite the observed horizontal shrinking. Our results now await comparison with upcoming Juno observations.

Rotating turbulence is commonly known for being dominated by geostrophic vortices that are invariant along the rotation axis and undergo inverse cascade. Yet, it has recently been shown to sustain fully three-dimensional states with a downscale energy cascade. In this letter, we investigate the statistical properties of three-dimensional rotating turbulence by the means of direct numerical simulations in a triply periodic box where geostrophic vortices are specifically damped. The resulting turbulent flow is an inertial wave turbulence that verifies the Zakharov-Kolmogorov spectrum derived analytically by Galtier (Phys. Rev. E, 68, 2003), thus offering numerical proof of the relevance of wave turbulence theory for three-dimensional, anisotropic waves. Lastly, we show that the same forcing leads to either geostrophic or wave turbulence depending on the initial condition. Our results thus bring further evidence for bi-stability in rotating turbulent flows., Comment: 7 pages, 4 figures

Coupled mixed convective and stratified systems are common in natural flows. To study experimentally the associated dynamics, we use a singular property of water: its non-linear equation of state is characterised by a maximum density close to $4^{\circ}$C. By heating the top of a tank at $35^{\circ}$C and cooling the bottom at $0^{\circ}$C, a two-layer configuration spontaneously appears. The convective motion in the bottom layer consists mostly of a large-scale circulation and rising cold plumes. This turbulent flow generates propagating internal gravity waves in the upper stably-stratified layer. PIV measurements are performed and spectral characteristics of the convection and internal gravity waves are presented. An horizontal large-scale reversing flow in the stratified layer is observed which is viscously driven by a third, intermediate layer. This buffer layer is located between the convective and stratified layers and is thermally coupled with the convective one, hence sustaining a strong horizontal shear. Three dimensional direct numerical simulations with geometry and physical parameters close to the experimental ones corroborate our experimental results., Comment: 21 pages. 22 figures

The dynamics, structure and stability of zonal jets in planetary flows are still poorly understood, especially in terms of coupling with the small-scale turbulent flow. Here, we use an experimental approach to address the questions of zonal jets formation and long-term evolution. A strong and uniform topographic $\beta$-effect is obtained inside a water-filled rotating tank thanks to the paraboloidal fluid free upper surface combined with a specifically designed bottom plate. A small-scale turbulent forcing is performed by circulating water through the base of the tank. Time-resolving PIV measurements reveal the self-organization of the flow into multiple zonal jets with strong instantaneous signature. We identify a subcritical bifurcation between two regimes of jets depending on the forcing intensity. In the first regime, the jets are steady, weak in amplitude, and directly forced by the local Reynolds stresses due to our forcing. In the second one, we observe highly energetic and dynamic jets of width larger than the forcing scale. An analytical modeling based on the quasi-geostrophic approximation reveals that this subcritical bifurcation results from the resonance between the directly forced Rossby waves and the background zonal flow.

We develop and analyse the first second-order phase-field model to combine melting and dissolution in multi-component flows. This provides a simple and accurate way to simulate challenging phase-change problems in existing codes. Phase-field models simplify computation by describing separate regions using a smoothed phase field. The phase field eliminates the need for complicated discretisations that track the moving phase boundary. However standard phase-field models are only first-order accurate. They often incur an error proportional to the thickness of the diffuse interface. We eliminate this dominant error by developing a general framework for asymptotic analysis of diffuse-interface methods in arbitrary geometries. With this framework we can consistently unify previous second-order phase-field models of melting and dissolution and the volume-penalty method for fluid-solid interaction. We finally validate second-order convergence of our model in two comprehensive benchmark problems using the open-source spectral code Dedalus., Comment: 24 pages, 6 figures

We report an idealized numerical study of a melting and freezing solid adjacent to a turbulent, buoyancy-affected shear flow, in order to improve our understanding of topography generation by phase changes in the environment. We use the phase-field method to dynamically couple the heat equation for the solid with the Navier-Stokes equations for the fluid. We investigate the evolution of an initially-flat and horizontal solid boundary overlying a pressure-driven turbulent flow. We assume a linear equation of state for the fluid and change the sign of the thermal expansion coefficient, such that the background density stratification is either stable, neutral or unstable. We find that channels aligned with the direction of the mean flow are generated spontaneously by phase changes at the fluid-solid interface. Streamwise vortices in the fluid, the interface topography and the temperature field in the solid influence each other and adjust until a statistical steady state is obtained. The crest-to-trough amplitude of the channels are larger than about 10$\delta_{\nu}$ in all cases, with $\delta_{\nu}$ the viscous length scale, but are much larger and more persistent for an unstable stratification than for a neutral or stable stratification. This happens because a stable stratification makes the cool melt fluid buoyant such that it shields the channel from further melting, whereas an unstable stratification makes the cool melt fluid sink, inducing further melting by rising hot plumes. The statistics of flow velocities and melt rates are investigated, and we find that channels and keels emerging in our simulations do not significantly change the mean drag coefficient., Comment: 35 pages, 15 figures

We show, using direct numerical simulations with experimentally realizable boundary conditions, that wall modes in Rayleigh-B\'enard convection in a rapidly rotating cylinder persist even very far from their linear onset. These nonlinear wall states survive in the presence of turbulence in the bulk and are robust with respect to changes in the shape of the boundary of the container. In this sense these modes behave much like the topologically protected states present in two-dimensional chiral systems even though rotating convection is a three-dimensional nonlinear driven dissipative system. We suggest that the robustness of this nonlinear mode may provide an explanation for the strong zonal flows observed recently in simulations and experiments on rapidly rotating convection at high Rayleigh number., Comment: 14 pages, 7 figures

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

A pure and incompressible material is confined between two plates such that it is heated from below and cooled from above. When its melting temperature is comprised between these two imposed temperatures, an interface separating liquid and solid phases appears. Depending on the initial conditions, freezing or melting occurs until the interface eventually converges towards a stationary state. This evolution is studied numerically in a two-dimensional configuration using a phase-field method coupled with the Navier-Stokes equations. Varying the control parameters of the model, we exhibit two types of equilibria: diffusive and convective. In the latter case, Rayleigh-B\'enard convection in the liquid phase shapes the solid-liquid front, and a macroscopic topography is observed. A simple way of predicting these equilibrium positions is discussed and then compared with the numerical simulations. In some parameter regimes, we show that multiple equilibria can coexist depending on the initial conditions. We also demonstrate that, in this bi-stable regime, transitioning from the diffusive to the convective equilibrium is inherently a nonlinear mechanism involving finite-amplitude perturbations.

Bowman et al. (2019) reported low-frequency photometric variability in 164 O- and B-type stars observed with K2 and TESS. They interpret these motions as internal gravity waves, which could be excited stochastically by convection in the cores of these stars. The detection of internal gravity waves in massive stars would help distinguish between massive stars with convective or radiative cores, determine core size, and would provide important constraints on massive star structure and evolution. In this work, we study the observational signature of internal gravity waves generated by core convection. We calculate the \textit{wave transfer function}, which links the internal gravity wave amplitude at the base of the radiative zone to the surface luminosity variation. This transfer function varies by many orders of magnitude for frequencies $\lesssim 1 \, {\rm d}^{-1}$, and has regularly-spaced peaks near $1 \, {\rm d}^{-1}$ due to standing modes. This is inconsistent with the observed spectra which have smooth ``red noise'' profiles, without the predicted regularly-spaced peaks. The wave transfer function is only meaningful if the waves stay predominately linear. We next show that this is the case: low frequency traveling waves do not break unless their luminosity exceeds the radiative luminosity of the star; and, the observed luminosity fluctuations at high frequencies are so small that standing modes would be stable to nonlinear instability. These simple calculations suggest that the observed low-frequency photometric variability in massive stars is not due to internal gravity waves generated in the core of these stars. We finish with a discussion of (sub)surface convection, which produces low-frequency variability in low-mass stars, very similar to that observed in Bowman et al. (2019) in higher mass stars., Comment: Accepted to ApJL. Supplementary information available at: http://www.princeton.edu/~lecoanet/data/SupplementaryInfo_Lecoanet2019.pdf

Planetary rotation organizes fluid motions into coherent, long-lived swirls, known as large-scale vortices (LSVs), which play an important role in the dynamics and long-term evolution of geophysical and astrophysical fluids. Here, using direct numerical simulations, we show that LSVs in rapidly rotating mixed convective and stably stratified fluids, which approximates the two-layer, turbulent-stratified dynamics of many geophysical and astrophysical fluids, have a generic shape and that their size can be predicted. We show that LSVs emerge in the convection zone from upscale energy transfers and can penetrate into the stratified layer. At the convective-stratified interface, the LSV cores have a positive buoyancy anomaly. Due to the thermal wind constraint, this buoyancy anomaly leads to winds in the stratified layer that decay over a characteristic vertical length scale. Thus LSVs take the shape of a depth-invariant cylinder with a finite-size radius in the turbulent layer and of a penetrating half dome in the stratified layer. Importantly, we demonstrate that when LSVs penetrate all the way through the stratified layer and reach a boundary that is no-slip, they saturate by boundary friction. We provide a prediction for the penetration depth and maximum radius of LSVs as a function of the LSV vorticity, the stratified layer depth, and the stratification. Our results, which apply for cyclonic LSVs, suggest that LSVs in slowly rotating stars and Earth's liquid core are confined to the convective layer, while in Earth's atmosphere and oceans they can penetrate far into the stratified layer., Comment: 15 pages, 12 figures

In this paper, we present an experimental investigation of the turbulent saturation of the flow driven by parametric resonance of inertial waves in a rotating fluid. In our set-up, a half-meter wide ellipsoid filled with water is brought to solid body rotation, and then undergoes sustained harmonic modulation of its rotation rate. This triggers the exponential growth of a pair of inertial waves via a mechanism called the libration-driven elliptical instability. Once the saturation of this instability is reached, we observe a turbulent state for which energy is injected into the resonant inertial waves only. Depending on the amplitude of the rotation rate modulation, two different saturation states are observed. At large forcing amplitudes, the saturation flow mainly consists of a steady, geostrophic anticyclone. Its amplitude vanishes as the forcing amplitude is decreased while remaining above the threshold of the elliptical instability. Below this secondary transition, the saturation flow is a superposition of inertial waves which are in weakly non-linear resonant interaction, a state that could asymptotically lead to inertial wave turbulence. In addition to being a first experimental observation of a wave-dominated saturation in unstable rotating flows, the present study is also an experimental confirmation of the model of Le Reun et al, PRL (2017) who introduced the possibility of these two turbulent regimes. The transition between these two regimes and their relevance to geophysical applications are finally discussed., Comment: Published in Journal of Fluid Mechanics, 30 pages, 19 pages

The possibility of subcritical behaviour in the geostrophic turbulence regime of rapidly rotating thermally driven convection is explored. In this regime a non-local inverse energy transfer may compete with the more traditional and local direct cascade. We show that, even for control parameters for which no inverse cascade has previously been observed, a subcritical transition towards a large-scale vortex state can occur when the system is initialized with a vortex dipole of finite amplitude. This new example of bistability in a turbulent flow, which may not be specific to rotating convection, opens up new avenues for studying energy transfer in strongly anisotropic three-dimensional flows., Comment: 12 pages, 6 figures

We study the evolution of a melting front between the solid and liquid phases of a pure incompressible material where fluid motions are driven by unstable temperature gradients. In a plane layer geometry, this can be seen as classical Rayleigh-B\'enard convection where the upper solid boundary is allowed to melt due to the heat flux brought by the fluid underneath. This free-boundary problem is studied numerically in two dimensions using a phase-field approach, classically used to study the melting and solidification of alloys, which we dynamically couple with the Navier-Stokes equations in the Boussinesq approximation. The advantage of this approach is that it requires only moderate modifications of classical numerical methods. We focus on the case where the solid is initially nearly isothermal, so that the evolution of the topography is related to the inhomogeneous heat flux from thermal convection, and does not depend on the conduction problem in the solid. From a very thin stable layer of fluid, convection cells appears as the depth -- and therefore the effective Rayleigh number of the layer increases. The continuous melting of the solid leads to dynamical transitions between different convection cell sizes and topography amplitudes. The Nusselt number can be larger than its value for a planar upper boundary, due to the feedback of the topography on the flow, which can stabilize large-scale laminar convection cells., Comment: 36 pages, 16 figures

Dynamo action in planetary cores has been extensively studied in the context of convectively-driven flows. We show in this letter that mechanical forcings, namely tides, libration and precession, are also able to kinematically sustain a magnetic field against ohmic diffusion. Previous attempts published in the literature focused on the laminar response or considered idealized spherical configurations. In contrast, we focus here on the developed turbulent regime and we self-consistently solve the magnetohydrodynamic (MHD) equations in an ellipsoidal container. Our results open new avenues of research in dynamo theory where both convection and mechanical forcing can play a role, independently or simultaneously.

We present three-dimensional direct numerical simulations of internal waves excited by turbulent convection in a self-consistent, Boussinesq and Cartesian model of convective--stably-stratified fluids. We demonstrate that in the limit of large Rayleigh number ($Ra\in [4\times 10^7,10^9]$) and large stratification (Brunt-V\"{a}is\"{a}l\"{a} frequencies $f_N \gg f_c$, where $f_c$ is the convective frequency), simulations are in good agreement with a theory that assumes waves are generated by Reynolds stresses due to eddies in the turbulent region (Lecoanet \& Quataert 2013 MNRAS 430 (3) 2363-2376). Specifically, we demonstrate that the wave energy flux spectrum scales like $k_{\perp}^4f^{-13/2}$ for weakly-damped waves (with $k_{\perp}$ and $f$ the waves' horizontal wavenumbers and frequencies), and that the total wave energy flux decays with $z$, the distance from the convective region, like $z^{-13/8}$., Comment: 13 pages, 6 figures

We demonstrate via direct numerical simulations that a periodic, oscillating mean flow spontaneously develops from turbulently generated internal waves. We consider a minimal physical model where the fluid self-organizes in a convective layer adjacent to a stably stratified one. Internal waves are excited by turbulent convective motions, then nonlinearly interact to produce a mean flow reversing on timescales much longer than the waves' period. Our results demonstrate for the first time that the three-scale dynamics due to convection, waves, and mean flow is generic and hence can occur in many astrophysical and geophysical fluids. We discuss efforts to reproduce the mean flow in reduced models, where the turbulence is bypassed. We demonstrate that wave intermittency, resulting from the chaotic nature of convection, plays a key role in the mean-flow dynamics, which thus cannot be captured using only second-order statistics of the turbulent motions.

We investigate the stability of stratified fluid layers undergoing homogeneous and periodic tidal deformation. We first introduce a local model which allows to study velocity and buoyancy fluctuations in a Lagrangian domain periodically stretched and sheared by the tidal base flow. While keeping the key physical ingredients only, such a model is efficient to simulate planetary regimes where tidal amplitudes and dissipation are small. With this model, we prove that tidal flows are able to drive parametric subharmonic resonances of internal waves, in a way reminiscent of the elliptical instability in rotating fluids. The growth rates computed via Direct Numerical Simulations (DNS) are in very good agreement with WKB analysis and Floquet theory. We also investigate the turbulence driven by this instability mechanism. With spatio-temporal analysis, we show that it is a weak internal wave turbulence occurring at small Froude and buoyancy Reynolds numbers. When the gap between the excitation and the Brunt-V\"ais\"al\"a frequencies is increased, the frequency spectrum of this wave turbulence displays a -2 power law reminiscent of the high-frequency branch of the Garett and Munk spectrum (Garrett & Munk 1979) which has been measured in the oceans. In addition, we find that the mixing efficiency is altered compared to what is computed in the context of DNS of stratified turbulence excited at small Froude and large buoyancy Reynolds numbers and is consistent with a superposition of waves., Comment: Accepted for publication in Journal of Fluid Mechanics, 27 pages, 21 figures

We present the stability analysis of a plane Couette flow which is stably stratified in the vertical direction orthogonally to the horizontal shear. Interest in such a flow comes from geophysical and astrophysical applications where background shear and vertical stable stratification commonly coexist. We perform the linear stability analysis of the flow in a domain which is periodic in the stream-wise and vertical directions and confined in the cross-stream direction. The stability diagram is constructed as a function of the Reynolds number Re and the Froude number Fr, which compares the importance of shear and stratification. We find that the flow becomes unstable when shear and stratification are of the same order (i.e. Fr $\sim$ 1) and above a moderate value of the Reynolds number Re$\gtrsim$700. The instability results from a resonance mechanism already known in the context of channel flows, for instance the unstratified plane Couette flow in the shallow water approximation. The result is confirmed by fully non linear direct numerical simulations and to the best of our knowledge, constitutes the first evidence of linear instability in a vertically stratified plane Couette flow. We also report the study of a laboratory flow generated by a transparent belt entrained by two vertical cylinders and immersed in a tank filled with salty water linearly stratified in density. We observe the emergence of a robust spatio-temporal pattern close to the threshold values of F r and Re indicated by linear analysis, and explore the accessible part of the stability diagram. With the support of numerical simulations we conclude that the observed pattern is a signature of the same instability predicted by the linear theory, although slightly modified due to streamwise confinement.

The combination of elliptical deformation of streamlines and vorticity can lead to the destabilisation of any rotating flow via the elliptical instability. Such a mechanism has been invoked as a possible source of turbulence in planetary cores subject to tidal deformations. The saturation of the elliptical instability has been shown to generate turbulence composed of non-linearly interacting waves and strong columnar vortices with varying respective amplitudes, depending on the control parameters and geometry. In this paper, we present a suite of numerical simulations to investigate the saturation and the transition from vortex-dominated to wave-dominated regimes. This is achieved by simulating the growth and saturation of the elliptical instability in an idealised triply periodic domain, adding a frictional damping to the geostrophic component only, to mimic its interaction with boundaries. We reproduce several experimental observations within one idealised local model and complement them by reaching more extreme flow parameters. In particular, a wave-dominated regime that exhibits many signatures of inertial wave turbulence is characterised for the first time. This regime is expected in planetary interiors., Comment: 5 pages, 4 figures, accepted for publication in Phys. Rev. Let

Regions of quiet Sun generally exhibit a complex distribution of small-scale magnetic field structures, which interact with the near-surface turbulent convective motions. Furthermore, it is probable that some of these magnetic fields are generated locally by a convective dynamo mechanism. In addition to the well-known granular and supergranular convective scales, various observations have indicated that there is an intermediate scale of convection, known as mesogranulation, with vertical magnetic flux concentrations accumulating preferentially at mesogranular boundaries. Our aim is to investigate the small-scale dynamo properties of a convective flow that exhibits both granulation and mesogranulation, comparing our findings with solar observations. Adopting an idealised model for a localised region of quiet Sun, we use numerical simulations of compressible magnetohydrodynamics, in a 3D Cartesian domain, to investigate the parametric dependence of this system (focusing particularly upon the effects of varying the aspect ratio and the Reynolds number). In purely hydrodynamic convection, we find that mesogranulation is a robust feature of this system provided that the domain is wide enough to accommodate these large-scale motions. The mesogranular peak in the kinetic energy spectrum is more pronounced in the higher Reynolds number simulations. We investigate the dynamo properties of this system in both the kinematic and the nonlinear regimes and we find that the dynamo is always more efficient in larger domains, when mesogranulation is present. Furthermore, we use a filtering technique in Fourier space to demonstrate that it is indeed the larger scales of motion that are primarily responsible for driving the dynamo. In the nonlinear regime, the magnetic field distribution compares very favourably to observations, both in terms of the spatial distribution and the measured field strengths., Comment: 12 pages, 11 figures, accepted for publication in Astronomy & Astrophysics

Mean-field dynamo theory suggests that turbulent convection in a rotating layer of electrically-conducting fluid produces a significant alpha-effect, which is one of the key ingredients in any mean-field dynamo model. Provided that this alpha-effect operates more efficiently than (turbulent) magnetic diffusion, such a system should be capable of sustaining a large-scale dynamo. However, in the Boussinesq model that was considered by Cattaneo&Hughes (2006) the dynamo produced small-scale, intermittent magnetic fields with no significant large-scale component. In this paper, we consider the compressible analogue of the rotating convective layer that was considered by Cattaneo&Hughes (2006). Varying the horizontal scale of the computational domain, we investigate the dependence of the dynamo upon the rotation rate. Our simulations indicate that these turbulent compressible flows can drive a small-scale dynamo but, even when the layer is rotating very rapidly (with a mid-layer Taylor number of Ta=10^8), we find no evidence for the generation of a significant large-scale component of the magnetic field on a dynamical timescale. Like Cattaneo&Hughes(2006), we measure a negligible (time-averaged) alpha-effect when a uniform horizontal magnetic field is imposed across the computational domain. Although the total horizontal magnetic flux is a conserved quantity in these simulations, the (depth-dependent) horizontally-averaged magnetic field always exhibits strong fluctuations. If these fluctuations are artificially suppressed within the code, we measure a significant mean electromotive force that is comparable to that found in related calculations in which the alpha-effect is measured using the test-field method, even though we observe no large-scale dynamo action., Comment: 25 pages, 10 figures, to appear in J. Fluid Mech

We study dynamo action in a convective layer of electrically-conducting, compressible fluid, rotating about the vertical axis. At the upper and lower bounding surfaces, perfectly-conducting boundary conditions are adopted for the magnetic field. Two different levels of thermal stratification are considered. If the magnetic diffusivity is sufficiently small, the convection acts as a small-scale dynamo. Using a definition for the magnetic Reynolds number $R_M$ that is based upon the horizontal integral scale and the horizontally-averaged velocity at the mid-layer of the domain, we find that rotation tends to reduce the critical value of $R_M$ above which dynamo action is observed. Increasing the level of thermal stratification within the layer does not significantly alter the critical value of $R_M$ in the rotating calculations, but it does lead to a reduction in this critical value in the non-rotating cases. At the highest computationally-accessible values of the magnetic Reynolds number, the saturation levels of the dynamo are similar in all cases, with the mean magnetic energy density somewhere between 4 and 9% of the mean kinetic energy density. To gain further insights into the differences between rotating and non-rotating convection, we quantify the stretching properties of each flow by measuring Lyapunov exponents. Away from the boundaries, the rate of stretching due to the flow is much less dependent upon depth in the rotating cases than it is in the corresponding non-rotating calculations. It is also shown that the effects of rotation significantly reduce the magnetic energy dissipation in the lower part of the layer. We also investigate certain aspects of the saturation mechanism of the dynamo., Comment: 23 pages, 12 figures

This article is focused on the dynamics of a rotating electrically conducting fluid in a turbulent state. As inside the Earth's core or in various industrial processes, a flow is altered by the presence of both background rotation and a large scale magnetic field. In this context, we present a set of 3D direct numerical simulations of incompressible decaying turbulence. We focus on parameters similar to the ones encountered in geophysical and astrophysical flows, so that the Rossby number is small, the interaction parameter is large, but the Elsasser number, defining the ratio between Coriolis and Lorentz forces, is about unity. These simulations allow to quantify the effect of rotation and thus inertial waves on the growth of magnetic fluctuations due to Alfv\'en waves. Rotation prevents the occurrence of equipartition between kinetic and magnetic energies, with a reduction of magnetic energy at decreasing Elsasser number {\Lambda}. It also causes a decrease of energy transfer mediated by cubic correlations. In terms of flow structure, a decrease of {\Lambda} corresponds to an increase in the misalignment of velocity and magnetic field., Comment: 18 pages, 12 figures

We analyze the anisotropy of turbulence in an electrically conducting fluid in the presence of a uniform magnetic field, for low magnetic Reynolds number, using the quasi-static approximation. In the linear limit, the kinetic energy of velocity components normal to the magnetic field decays faster than the kinetic energy of component along the magnetic field [Moffatt, JFM 28, 1967]. However, numerous numerical studies predict a different behaviour, wherein the final state is characterized by dominant horizontal energy. We investigate the corresponding nonlinear phenomenon using Direct Numerical Simulations. The initial temporal evolution of the decaying flow indicates that the turbulence is very similar to the so-called "two-and-a-half-dimensional" flow [Montgomery & Turner, Phys. Fluids 25(2), 1982] and we offer an explanation for the dominance of horizontal kinetic energy., Comment: 17 pages, 8 figures