Search

Mesoscale convection covers an intermediate scale range between small-scale turbulence and the global organization of the convection flow. It is often characterized by an order of the convection patterns despite very high Rayleigh numbers and strong turbulent fluctuations. In this review, we discuss several aspects of mesoscale convection, which have been investigated by three-dimensional direct numerical simulations. The numerical studies are performed in a characteristic configuration of a plane layer that is heated from below and cooled from above or subject to constant heat flux at the top and bottom boundaries. We discuss the role of the thermal and mechanical boundary conditions for structure formation and study the impact of the domain shape as well as the Prandtl number. With respect to the latter, we focus on low values that arise in astrophysical convection and are partly not anymore accessible in laboratory experiments with liquid metals. Beside these experiments in the Boussinesq approximation, we report studies of non-Boussinesq mesoscale convection. This is done by investigating effects of compressibility and temperature dependence of material properties. The kinetic energy dissipation rate turns out to remain a central quantity for the turbulent mixing in compressible convection. Their different components, statistics, relation to the turbulent viscosity, and the multifractal properties are discussed., Comment: 26 pages, 16 figures (some in reduced quality)

The impact of variable material properties, such as temperature-dependent thermal conductivity and dynamical viscosity, on the dynamics of a fully compressible turbulent convection flow beyond the anelastic limit are studied in the present work by two series of three-dimensional direct numerical simulations in a layer of aspect ratio 4 with periodic boundary conditions in both horizontal directions. One simulation series is for a weakly stratified adiabatic background, one for a strongly stratified one. The Rayleigh number is $10^5$ and the Prandtl number is 0.7 throughout this study. The temperature dependence of material parameters is imposed as a power law with an exponent $\beta$. It generates a superadiabaticity $\varepsilon(z)$ that varies across the convection layer. Central statistical quantities of the flow, such as the mean superadiabatic temperature, temperature and density fluctuations, or turbulent Mach numbers are compared in the form of horizontal plane-time averaged profiles. It is found that the additional material parameter dependence causes systematic quantitative changes of all these quantities, but no qualitative ones. A growing temperature power law exponent $\beta$ also enhances the turbulent momentum transfer in the weak stratification case by 40\%, it reduces the turbulent heat transfer by up to 50\% in the strong stratification case.

Buoyancy-driven turbulent convection leads to a fully compressible flow with a strong top-down asymmetry of first- and second-order statistics when the adiabatic equilibrium profiles of temperature, density and pressure decay very strongly across the convection layer. The growth of this asymmetry and the formation of an increasingly thicker stabilized sublayer with a negative mean convective heat flux at the top of the convection zone is reported here by a series highly resolved three-dimensional direct numerical simulations beyond the Oberbeck-Boussinesq and anelastic limits for dimensionless dissipation numbers $0.1 \le D\le 0.8$ at fixed Rayleigh number $Ra=10^6$ and superadiabaticity. The highly stratified compressible convection regime appears for $D > D_{\rm crit}\approx 0.65$, when density fluctuations collapse to those of pressure; it is characterized by an up to nearly 50\% reduced global turbulent heat transfer and a sparse network of focused thin thermal plumes falling through the top sublayer deep into the bulk.

Fully compressible turbulent convection beyond the Oberbeck-Boussinesq limit and anelastic regime is studied in three-dimensional numerical simulations. Superadiabaticity $\epsilon$ and dissipation number $D$, which measures the strength of stratification of adiabatic equilibria, cause two limits of compressible convection -- nearly top-down-symmetric, strongly superadiabatic and highly top-down-asymmetric, strongly stratified convection. Highest turbulent Mach numbers $M_t$ follow for a symmetric blend of these two limits which we term the fully compressible case. Particularly, the strongly stratified convection case leads to a fluctuation-reduced top layer in the convection zone, a strongly reduced global heat transfer, and differing boundary layer dynamics between top and bottom. We detect this asymmetry for growing dissipation number $D$ also in the phase plane which is spanned by the turbulent Mach number $M_t$ and the dilatation parameter $\delta$ which relates the dilatational velocity fluctuations to the solenoidal ones. A detailed analysis of the different transport currents in the fully compressible energy budget relates the low-$D$ convection cases to the standard definition of the dimensionless Nusselt number in the Oberbeck-Boussinesq limit.