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.

Inspection of available data on the decay exponent for the kinetic energy of homogeneous and isotropic turbulence (HIT) shows that it varies by as much as 100\%. Measurements and simulations often show no correspondence with theoretical arguments, which are themselves varied. This situation is unsatisfactory given that HIT is a building block of turbulence theory and modeling. We take recourse to a large base of direct numerical simulations and study decaying HIT for a variety of initial conditions. We show that the Kolmogorov decay exponent and the Birkhoff-Saffman decay are both readily observed, albeit approximately, for long periods of time if the initial conditions are appropriately arranged. We also present, for both cases, other turbulent statistics such as the velocity derivative skewness, energy spectra and dissipation, and show that the decay and growth laws are approximately as expected theoretically, though the wavenumber spectrum near the origin begins to change relatively quickly, suggesting that the invariants do not strictly exist. We comment briefly on why the decay exponent has varied so widely in past experiments and simulations., Comment: 14 pages, 11 figures

Mixing of passive scalars in compressible turbulence does not obey the same classical Reynolds number scaling as its incompressible counterpart. We first show from a large database of direct numerical simulations that even the solenoidal part of the velocity field fails to follow the classical incompressible scaling when the forcing includes a substantial dilatational component. Though the dilatational effects on the flow remain significant, our main results are that both the solenoidal energy spectrum and the passive scalar spectrum scale assume incompressible forms, and that the scalar gradient aligns with the most compressive eigenvalue of the solenoidal part, provided that only the solenoidal components are used for scaling in a consistent manner. Minor modifications to this statement are also pointed out., Comment: 6 pages and 7 figures

Turbulent flows, ubiquitous in nature and engineering, comprise fluctuations over a wide range of spatial and temporal scales. While flows with fluctuations in thermodynamic variables are much more common, much less is known about these flows than their incompressible counterparts in which thermodynamics is decoupled from hydrodynamics. A critical element in the study of the latter has been the concept of universal scaling laws which provides fundamental as well as practical information about the spatio-temporal behavior of these complex systems. Part of this success is due to the dependence on a single non-dimensional parameter, that is the Reynolds number. Universality in compressible flows, on the other hand, have proven to be more elusive as no unifying set of parameters were found to yield universal scaling laws. This severely limits our understanding of these flows and the successful development of theoretically sound models. Using results in specific asymptotic limits of the governing equations we show that universal scaling is indeed observed when the set of governing parameters is expanded to include internally generated dilatational scales which are the result of the driving mechanisms that produce the turbulence. The analysis demonstrate why previous scaling laws fail in the general case, and opens up new venues to identify physical processes of interest and create turbulence models needed for simulations of turbulent flows at realistic conditions. We support our results with a new massive database of highly-resolved direct numerical simulations along with all data available in the literature. In search of universal features, we suggest classes which bundle the evolution of flows in the new parameter space. An ultimate asymptotic regime predicted independently by renormalization group theories and statistical mechanics is also assssed with available data.