No Arabic abstract
In the preparation of Cafe Latte, spectacular layer formation can occur between the expresso shot in a glass of milk and the milk itself. Xue et al. (Nat. Commun., vol. 8, 2017, pp. 1-6) showed that the injection velocity of expresso determines the depth of coffee-milk mixture. After a while when a stable stratification forms in the mixture, the layering process can be modelled as a double diffusive convection system with a stably-stratified coffee-milk mixture cooled from the side. More specifically, we perform (two-dimensional) direct numerical simulations of laterally cooled double diffusive convection for a wide parameter range, where the convective flow is driven by a lateral temperature gradient while stabilized by a vertical concentration gradient. When the thermal driving force dominates over the stabilizing force, the flow behaves like vertical convection in which a large-scale circulation develops. However, with increasing strength of the stabilizing force, a meta-stable layered regime emerges. Initially, several vertically-stacked convection rolls develop, and these well-mixed layers are separated by sharp interfaces with large concentration gradients. The initial thickness of these emerging layers can be estimated by balancing the work exerted by thermal driving and the required potential energy to bring fluid out of its equilibrium position in the stably stratified fluid. In the layered regime, we further observe successive layer merging, and eventually only a single convection roll remains. We elucidate the following merging mechanism: As weakened circulation leads to accumulation of hot fluid adjacent to the hot sidewall, larger buoyancy forces associated with hotter fluid eventually break the layer interface. Then two layers merge into a larger layer, and circulation establishes again within the merged structure.
Three dimensional roll-type double-diffusive convection in a horizontally infinite layer of an uncompressible liquid is considered in the neighborhood of Hopf bifurcation points. A system of amplitude equations for the variations of convective rolls amplitude is derived by multiple-scaled method. An attention is paid to an interaction of convection and horizontal vortex. Different cases of the derived equations are discussed.
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.
Vortices play an unique role in heat and momentum transports in astro- and geo-physics, and it is also the origin of the Earths dynamo. A question existing for a long time is whether the movement of vortices can be predicted or understood based on their historical data. Here we use both the experiments and numerical simulations to demonstrate some generic features of vortex motion and distribution. It can be found that the vortex movement can be described on the framework of Brownian particles where they move ballistically for the time shorter than some critical timescales, and then move diffusively. Traditionally, the inertia of vortex has often been neglected when one accounts for their motion, our results imply that vortices actually have inertial-induced memory such that their short term movement can be predicted. Extending to astro- and geo-physics, the critical timescales of transition are in the order of minutes for vortices in atmosphere and ocean, in which this inertial effect may often be neglected compared to other steering sources. However, the timescales for vortices are considerably larger which range from days to a year. It infers the new concept that not only the external sources alone, for example the solar wind, but also the internal source, which is the vortex inertia, can contribute to the short term Earths magnetic field variation.
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.
The process referred to as semi-convection in astrophysics and double-diffusive convection in the diffusive regime in Earth and planetary sciences, occurs in stellar and planetary interiors in regions which are stable according to the Ledoux criterion but unstable according to the Schwarzschild criterion. In this series of papers, we analyze the results of an extensive suite of 3D numerical simulations of the process, and ultimately propose a new 1D prescription for heat and compositional transport in this regime which can be used in stellar or planetary structure and evolution models. In a preliminary study of the phenomenon, Rosenblum et al. (2011) showed that, after saturation of the primary instability, a system can evolve in one of two possible ways: the induced turbulence either remains homogeneous, with very weak transport properties, or transitions into a thermo-compositional staircase where the transport rate is much larger (albeit still smaller than in standard convection). In this paper, we show that this dichotomous behavior is a robust property of semi-convection across a wide region of parameter space. We propose a simple semi-analytical criterion to determine whether layer formation is expected or not, and at what rate it proceeds, as a function of the background stratification and of the diffusion parameters (viscosity, thermal diffusivity and compositional diffusivity) only. The theoretical criterion matches the outcome of our numerical simulations very adequately in the numerically accessible planetary parameter regime, and can easily be extrapolated to the stellar parameter regime. Subsequent papers will address more specifically the question of quantifying transport in the layered case and in the non-layered case.