No Arabic abstract
Direct numerical simulations are employed to reveal three distinctly different flow regions in rotating spherical Rayleigh-Benard convection. In the low-latitude region $mathrm{I}$ vertical (parallel to the axis of rotation) convective columns are generated between the hot inner and the cold outer sphere. The mid-latitude region $mathrm{II}$ is dominated by vertically aligned convective columns formed between the Northern and Southern hemispheres of the outer sphere. The diffusion-free scaling, which indicates bulk-dominated convection, originates from this mid-latitude region. In the equator region $mathrm{III}$ the vortices are affected by the outer spherical boundary and are much shorter than in region $mathrm{II}$. Thermally driven turbulence with background rotation in spherical Rayleigh-Benard convection is found to be characterized by three distinctly different flow regions. The diffusion-free scaling, which indicates the heat transfer is bulk-dominated, originates from the mid-latitude region in which vertically aligned vortices are stretched between the Northern and Southern hemispheres of the outer sphere. These results show that the flow physics in rotating convection are qualitatively different in planar and spherical geometries. This finding underlines that it is crucial to study convection in spherical geometries to better understand geophysical and astrophysical flow phenomena.
We study numerically the melting of a horizontal layer of a pure solid above a convecting layer of its fluid rotating about the vertical axis. In the rotating regime studied here, with Rayleigh numbers of order $10^7$, convection takes the form of columnar vortices, the number and size of which depend upon the Ekman and Prandtl numbers, as well as the geometry -- periodic or confined. As the Ekman and Rayleigh numbers vary, the number and average area of vortices vary in inverse proportion, becoming thinner and more numerous with decreasing Ekman number. The vortices transport heat to the phase boundary thereby controlling its morphology, characterized by the number and size of the voids formed in the solid, and the overall melt rate, which increases when the lower boundary is governed by a no-slip rather than a stress-free velocity boundary condition. Moreover, the number and size of voids formed are relatively insensitive to the Stefan number, here inversely proportional to the latent heat of fusion. For small values of the Stefan number, the convection in the fluid reaches a slowly evolving geostrophic state wherein columnar vortices transport nearly all the heat from the lower boundary to melt the solid at an approximately constant rate. In this quasi-steady state, we find that the Nusselt number, characterizing the heat flux, co-varies with the interfacial roughness, for all the flow parameters and Stefan numbers considered here. This confluence of processes should influence the treatment of moving boundary problems, particularly those in astrophysical and geophysical problems where rotational effects are important.
The effect of rotation on the boundary layers (BLs) in a Rayleigh-Benard (RB) system at a relatively low Rayleigh number, i.e. $Ra = 4times10^7$, is studied for different Pr by direct numerical simulations and the results are compared with laminar BL theory. In this regime we find a smooth onset of the heat transfer enhancement as function of increasing rotation rate. We study this regime in detail and introduce a model based on the Grossmann-Lohse theory to describe the heat transfer enhancement as function of the rotation rate for this relatively low Ra number regime and weak background rotation $Rogtrsim 1$. The smooth onset of heat transfer enhancement observed here is in contrast to the sharp onset observed at larger $Ra gtrsim 10^8$ by Stevens {it{et al.}} [Phys. Rev. Lett. {bf{103}}, 024503, 2009], although only a small shift in the Ra-Ro-Pr phase space is involved.
For rapidly rotating turbulent Rayleigh--Benard convection in a slender cylindrical cell, experiments and direct numerical simulations reveal a boundary zonal flow (BZF) that replaces the classical large-scale circulation. The BZF is located near the vertical side wall and enables enhanced heat transport there. Although the azimuthal velocity of the BZF is cyclonic (in the rotating frame), the temperature is an anticyclonic traveling wave of mode one whose signature is a bimodal temperature distribution near the radial boundary. The BZF width is found to scale like $Ra^{1/4}Ek^{2/3}$ where the Ekman number $Ek$ decreases with increasing rotation rate.
Using direct numerical simulations, we study rotating Rayleigh-Benard convection in a cylindrical cell for a broad range of Rayleigh, Ekman, and Prandtl numbers from the onset of wall modes to the geostrophic regime, an extremely important one in geophysical and astrophysical contexts. We connect linear wall-mode states that occur prior to the onset of bulk convection with the boundary zonal flow that coexists with turbulent bulk convection in the geostrophic regime through the continuity of length and time scales and of convective heat transport. We quantitatively collapse drift frequency, boundary length, and heat transport data from numerous sources over many orders of magnitude in Rayleigh and Ekman numbers. Elucidating the heat transport contributions of wall modes and of the boundary zonal flow are critical for characterizing the properties of the geostrophic regime of rotating convection in finite, physical containers and is crucial for connecting the geostrophic regime of laboratory convection with geophysical and astrophysical systems.
Double-diffusive convection driven by both thermal and compositional buoyancy in a rotating spherical shell can exhibit a rather large number of behaviours often distinct from that of the single diffusive system. In order to understand how the differences in thermal and compositional molecular diffusivities determine the dynamics of thermo-compositional convection we investigate numerically the linear onset of convective instability in a double-diffusive setup. We construct an alternative equivalent formulation of the non-dimensional equations where the linearised double-diffusive problem is described by an effective Rayleigh number, $text{Ra}$, measuring the amplitude of the combined buoyancy driving, and a second parameter, $alpha$, measuring the mixing of the thermal and compositional contributions. This formulation is useful in that it allows for the analysis of several limiting cases and reveals dynamical similarities in the parameters space which are not obvious otherwise. We analyse the structure of the critical curves in this $text{Ra}-alpha$ space, explaining asymptotic behaviours in $alpha$, transitions between inertial and diffusive regimes, and transitions between large scale (fast drift) and small scale (slow drift) convection. We perform this analysis for a variety of diffusivities, rotation rates and shell aspect ratios showing where and when new modes of convection take place.