We studied flow organization and heat transfer properties in mixed turbulent convection within Poiseuille–Rayleigh–Bénard channels subjected to temporally modulated sinusoidal wall temperatures. Three-dimensional direct numerical simulations were performed for Rayleigh numbers in the range $10^6 \leqslant Ra \leqslant 10^8$, a Prandtl number $Pr = 0.71$ and a bulk Reynolds number $Re_b \approx 5623$. We found that high-frequency wall temperature oscillations had minimal impact on flow structures, while low-frequency oscillations induced adaptive changes, forming stable stratified layers during cooling. Proper orthogonal decomposition (POD) analysis revealed a dominant streamwise unidirectional shear flow mode. Large-scale rolls oriented in the streamwise direction appeared as higher POD modes and were significantly influenced by lower-frequency wall temperature variations. Long-time-averaged statistics showed that the Nusselt number increased with decreasing frequency by up to 96 %, while the friction coefficient varied by less than 15 %. High-frequency modulation predominantly influenced near-wall regions, enhancing convective effects, whereas low frequencies reduced these effects via stable stratified layer formation. Phase-averaged statistics showed that high-frequency modulation resulted in phase-stable streamwise velocity and temperature profiles, while low-frequency modulation caused significant variations due to weakened turbulence. Turbulent kinetic energy (TKE) profiles remained high near the wall during both heating and cooling at high frequency, but decreased during cooling at low frequencies. A TKE budget analysis revealed that during heating, TKE production was dominated by shear near the wall and by buoyancy in the bulk region; while during cooling, the production, distribution and dissipation of TKE were all nearly zero.

This study introduces a novel approach to investigate the Reynolds analogy in complex flow scenarios. It is shown that the total mechanical energy $\mathit {B}$, viz. the sum of kinetic energy and pressure work, and the field $\Gamma =\theta ^2/2$ (where $\theta$ is the transported passive scalar) are governed by two equations that are similar in form, when time-averaged for statistically stationary flows. For fully developed channel flows the integral energy balance links the mean bulk velocity and scalar with the volume averages of the respective dissipation rates, allowing the assessment of the Reynolds analogy in terms of the dissipation fields. This approach is tested on direct numerical simulation data of rough-wall turbulent channel flow at two different roughness Reynolds numbers, namely $k^+=15$ and $k^+=90$. For a unit Prandtl number, the same qualitative behaviour is observed for the mean wall-normal distributions of the budget-equation terms of $B$ and $\Gamma$, the latter being larger than the corresponding terms in the mechanical-energy budget. The Reynolds decomposition of the flow into temporal mean and stochastic parts reveals that roughness primarily affects the mean-flow dissipation. For the $k^+=90$ case, the analysis shows that attached-flow and high-shear regions dominate the integral mean scalar and momentum transfer and exhibit the greatest differences between the mean mechanical and scalar dissipation rates. In contrast, well-mixed regions, sheltered by large roughness elements, contribute similarly and minimally to the integral scalar and momentum transfer.

We analysed the performance of convolutional autoencoders in generating reduced-order representations of the temperature field of two-dimensional Rayleigh–Bénard flows at $\textit{Pr} =1$ and Rayleigh numbers extending from $10^6$ to $10^8$, capturing the range where the flow transitions to turbulence. We present a way of estimating the minimum number of dimensions needed by the autoencoders to capture all the relevant physical scales of the data that is more apt for highly multiscale flows than previous criteria applied to lower-dimensional systems. We compare our architecture with two regularized variants as well as with linear methods, and find that manually fixing the dimension of the latent space produces the best results. We show how the estimated minimum dimension presents a sharp increase around $Ra\sim 10^7$, when the flow starts to transition to turbulence. Furthermore, we show how this dimension does not follow the same scaling as the physically relevant scales, such as the dissipation length scale and the thermal boundary layer.

We study the transition from the momentum- to buoyancy-dominated regime in temporal jets forced by gravity. From the conservation of the thermal content and of the volume flux, we develop a simple model which is able to describe accurately the transition between the two regimes in terms of a single parameter representing the entrainment coefficient. Our analytical results are validated against a set of numerical simulations of temporal planar forced plumes at different initial values of Reynolds and Froude numbers. We find that, although the the pure jet-scaling law is not clearly observed in simulations at finite Froude number, the model correctly describes the transition to the buoyancy-dominated regime which emerges at long times.

We present the first nonlinear results on the problem of non-rotating thermal convection in an internally heated full sphere. A nonlinear stability analysis by the energy method yields that, at least for no-slip boundary conditions, the critical Rayleigh numbers for linear stability and nonlinear stability coincide. We then explore different ranges of the parameter regime using direct numerical simulations. We first report on the system behaviour for a fixed Prandtl number of unity and both stress-free and no-slip boundary conditions up to very high forcing, reaching Rayleigh number $Ra=2\times 10^{12}$, approximately 250 million times the critical value ($Ra_c$) for the onset of convection under no-slip conditions. For both boundary conditions, we observe a scaling for the advective heat transfer measured by the Nusselt number $Nu$ close to $Nu \sim Ra^{1/4}$. This is consistent with a scaling prediction that we formulate analogously to the classical scaling in Rayleigh–Bénard convection. We then investigate the Prandtl number dependence at low to intermediate forcing for stress-free boundary conditions in the ranges $0.1 \leq Pr \leq 30$ and $Ra_c=3091\leq Ra \leq 3\times 10^5 \approx 100Ra_c$. We find five distinct dynamical regimes depending on the Prandtl number, describe each regime individually and issue heuristic interpretations of the system behaviour where possible.

We establish a theoretical framework for predicting friction and heat transfer coefficients in variable-property forced air convection. Drawing from concepts in high-speed wall turbulence, which also involves significant temperature, viscosity and density variations, we utilize the mean momentum balance and mean thermal balance equations to develop integral transformations that account for the impact of variable fluid properties. These transformations are then applied inversely to predict the friction and heat transfer coefficients, leveraging the universality of passive scalars transport theory. Our proposed approach is validated using a comprehensive dataset from direct numerical simulations (DNS), covering both heating and cooling conditions up to a friction Reynolds number $\textit {Re}_\tau \approx 3200$. The predicted friction and heat transfer coefficients closely match the DNS data with accuracy margin 1–2 %, representing a significant improvement over the current state of the art.

We provide scaling relations for the Nusselt number $Nu$ and the friction coefficient $C_{S}$ in sheared Rayleigh–Bénard convection, i.e. in Rayleigh–Bénard flow with Couette- or Poiseuille-type shear forcing, by extending the Grossmann & Lohse (J. Fluid Mech., vol. 407, 2000, pp. 27–56, Phys. Rev. Lett., vol. 86, 2001, pp. 3316–3319, Phys. Rev. E, vol. 66, 2002, 016305, Phys. Fluids, vol. 16, 2004, pp. 4462–4472) theory to sheared thermal convection. The control parameters for these systems are the Rayleigh number $Ra$, the Prandtl number $Pr$ and the Reynolds number $Re_S$ that characterises the strength of the imposed shear. By direct numerical simulations and theoretical considerations, we show that, in turbulent Rayleigh–Bénard convection, the friction coefficients associated with the applied shear and the shear generated by the large-scale convection rolls are both well described by Prandtl's (Ergeb. Aerodyn. Vers. Gött., vol. 4, 1932, pp. 18–29) logarithmic friction law, suggesting some kind of universality between purely shear-driven flows and thermal convection. These scaling relations hold well for $10^6 \leq Ra \leq 10^8$, $0.5 \leq Pr \leq 5.0$, and $0 \leq Re_S \leq 10^4$.

We employ direct numerical simulations to investigate the heat transfer and flow structures in turbulent Rayleigh–Bénard convection in both cylindrical cells and laterally periodic domains, spanning an unprecedentedly wide range of aspect ratios $0.075 \leqslant \varGamma \leqslant 32$. We focus on Prandtl number ${Pr}=1$ and Rayleigh numbers ${{Ra}}=2\times 10^7$ and ${{Ra}}=10^8$. In both cases, with increasing aspect ratio, the heat transfer first increases, then reaches a maximum (which is more pronounced for the cylindrical case due to confinement effects), and then slightly goes down again before it finally saturates at the large aspect ratio limit, which is achieved already at $\varGamma \approx 4$. Already for $\varGamma \gtrsim 0.75$, the heat transfers in both cylindrical and laterally periodic domains become identical. The large-$\varGamma$ limit for the volume-integrated Reynolds number and the boundary layer thicknesses are also reached at $\varGamma \approx 4$. However, while the integral flow properties converge at $\varGamma \approx 4$, the confinement of a cylindrical domain impacts the temperature and velocity variance distributions up to $\varGamma \approx 16$, as thermal superstructures cannot form close to the sidewall.

The seminal Bolgiano–Obukhov (BO) theory established the fundamental framework for turbulent mixing and energy transfer in stably stratified fluids. However, the presence of BO scalings remains debatable despite their being observed in stably stratified atmospheric layers and convective turbulence. In this study, we performed precise temperature measurements with 51 high-resolution loggers above the seafloor for 46 h on the continental shelf of the northern South China Sea. The temperature observation exhibits three layers with increasing distance from the seafloor: the bottom mixed layer (BML), the mixing zone and the internal wave zone. A BO-like scaling $\alpha =-1.34\pm 0.10$ is observed in the temperature spectrum when the BML is in a weakly stable stratified ($N\sim 0.0018$ rad s$^{-1}$) and strongly sheared ($Ri\sim 0.0027$) condition, whereas in the unstably stratified convective turbulence of the BML, the scaling $\alpha =-1.76\pm 0.10$ clearly deviated from the BO theory but approached the classical $-$5/3 scaling in isotropic turbulence. This suggests that the convective turbulence is not the promise of BO scaling. In the mixing zone, where internal waves alternately interact with the BML, the scaling follows the Kolmogorov scaling. In the internal wave zone, the scaling $\alpha =-2.12 \pm 0.15$ is observed in the turbulence range and possible mechanisms are provided.

Radial unstable stratification is a potential source of turbulence in the cold regions of accretion disks. To investigate this thermal effect, here we focus on two-dimensional Rayleigh–Bénard convection in an annulus subject to radially dependent gravitational acceleration $g \propto 1/r$. Next to the Rayleigh number $Ra$ and Prandtl number $Pr$, the radius ratio $\eta$, defined as the ratio of inner and outer cylinder radii, is a crucial parameter governing the flow dynamics. Using direct numerical simulations for $Pr=1$ and $Ra$ in the range from $10^7$ to $10^{10}$, we explore how variations in $\eta$ influence the asymmetry in the flow field, particularly in the boundary layers. Our results show that in the studied parameter range, the flow is dominated by convective rolls and that the thermal boundary-layer (TBL) thickness ratio between the inner and outer boundaries varies as $\eta ^{1/2}$. This scaling is attributed to the equality of velocity scales in the inner ($u_i$) and outer ($u_o$) regions. We further derive that the temperature drops in the inner and outer TBLs scale as $1/(1+\eta ^{1/2})$ and $\eta ^{1/2}/(1+\eta ^{1/2})$, respectively. The scalings and the temperature drops are in perfect agreement with the numerical data.

Geothermal gradients and heterogeneous permeability are commonly observed in natural geological formations for underground CO$_2$ sequestration. In this study, we conduct three-dimensional direct numerical simulations on the double-diffusive convection with both unstable temperature and concentration gradients in homogeneous and heterogeneous porous media. For homogeneous porous media, the root-mean-squared velocity increases linearly with density ratio defined as the buoyancy ratio by temperature and concentration differences. The flow structures show no remarkable changes when temperature Rayleigh number ${Ra}_T$ is less than its critical value, but alter from sheet-like to cellular structures as ${Ra}_T$ surpasses this threshold. The concentration wavenumber scales approximately as $k_{rS}\sim {Ra}_e^{0.47}$ with a defined effective Rayleigh number ${Ra}_e$. By using a scale analysis, the concentration flux exhibits a consistent linear relation with the total driving forces for all simulations. For heterogeneous porous media, where the Dykstra–Parsons coefficient $V_{DP}$ and correlation length $l_{r}$ determine the spatial distribution of the permeability field, the flow is strengthened in places with higher permeability. The velocity and concentration flux are less affected by $l_{r}$ than that by $V_{DP}$. For small correlation length, the flow structures coarsen and their characteristic width generally increases with increasing heterogeneity. For large correlation length, small structures emerge in the regions with large permeability, which can be attributed to the intensified local Rayleigh number triggering more vigorous convection there. The variations of concentration flux with $l_{r}$ and $V_{DP}$ can be explained by the portion of area covered by high concentration with large vertical velocity near the boundaries.

Turbulent convection in the interiors of the Sun and the Earth occurs at high Rayleigh numbers $Ra$, low Prandtl numbers $Pr$, and different levels of rotation rates. To understand the combined effects better, we study rotating turbulent convection for $Pr = 0.021$ (for which some laboratory data corresponding to liquid metals are available), and varying Rossby numbers $Ro$, using direct numerical simulations in a slender cylinder of aspect ratio 0.1; this confinement allows us to attain high enough Rayleigh numbers. We are motivated by the earlier finding in the absence of rotation that heat transport at high enough $Ra$ is similar between confined and extended domains. We make comparisons with higher aspect ratio data where possible. We study the effects of rotation on the global transport of heat and momentum as well as flow structures (a) for increasing rotation at a few fixed values of $Ra$, and (b) for increasing $Ra$ (up to $10^{10}$) at the fixed, low Ekman number $1.45 \times 10^{-6}$. We compare the results with those from unity $Pr$ simulations for the same range of $Ra$ and $Ro$, and with the non-rotating case over the same range of $Ra$ and low $Pr$. We find that the effects of rotation diminish with increasing $Ra$. These results and comparison studies suggest that for high enough $Ra$, rotation alters convective flows in a similar manner for small and large aspect ratios, so useful insights on the effects of high thermal forcing on convection can be obtained by considering slender domains.

We studied the transport and deposition behaviour of point particles in Rayleigh–Bénard convection cells subjected to Couette-type wall shear. Direct numerical simulations (DNSs) are performed for Rayleigh number ($Ra$) in the range $10^{7} \leq Ra \leq ~10^9$ with a fixed Prandtl number $Pr = 0.71$, while the wall-shear Reynolds number ($Re_w$) is in the range $0 \leq Re_w \leq ~12\,000$. With the increase of $Re_w$, the large-scale rolls expanded horizontally, evolving into zonal flow in two-dimensional simulations or streamwise-oriented rolls in three-dimensional simulations. We observed that, for particles with a small Stokes number ($St$), they either circulated within the large-scale rolls when buoyancy dominated or drifted near the walls when shear dominated. For medium $St$ particles, pronounced spatial inhomogeneity and preferential concentration were observed regardless of the prevailing flow state. For large $St$ particles, the turbulent flow structure had a minor influence on the particles’ motion; although clustering still occurred, wall shear had a negligible influence compared with that for medium $St$ particles. We then presented the settling curves to quantify the particle deposition ratio on the walls. Our DNS results aligned well with previous theoretical predictions, which state that small $St$ particles settle with an exponential deposition ratio and large $St$ particles settle with a linear deposition ratio. For medium $St$ particles, where complex particle–turbulence interaction emerges, we developed a new model describing the settling process with an initial linear stage followed by a nonlinear stage. Unknown parameters in our model can be determined either by fitting the settling curves or using empirical relations. Compared with DNS results, our model also accurately predicts the average residence time across a wide range of $St$ for various $Re_w$.

We discuss the applicability of quasilinear-type approximations for a turbulent system with a large range of spatial and temporal scales. We consider a paradigm fluid system of rotating convection with vertical and horizontal temperature gradients. In particular, the interaction of rotation with the horizontal temperature gradient drives a ‘thermal wind’ shear flow whose strength is controlled by the horizontal temperature gradient. Varying this parameter therefore systematically alters the ordering of the shearing time scale, the convective time scale and the correlation time scale. We demonstrate that quasilinear-type approximations work well when the shearing time scale or the correlation time scale is sufficiently short. In all cases, the generalised quasilinear approximation systematically outperforms the quasilinear approximation. We discuss the consequences for statistical theories of turbulence interacting with mean gradients.

Turbulent shear flows driven by a combination of a pressure gradient and buoyancy forcing are investigated using direct numerical simulations. Specifically, we consider the set-up of a differentially heated vertical channel subject to a Poiseuille-like horizontal pressure gradient. We explore the response of the system to its three control parameters: the Grashof number $Gr$, the Prandtl number $Pr$, and the Reynolds number $Re$ of the pressure-driven flow. From these input parameters, the relative strength of buoyancy driving to the pressure gradient can be quantified by the Richardson number $Ri=Gr/Re^2$. We compare the response of the mixed vertical convection configuration to that of mixed Rayleigh–Bénard convection, and find a nearly identical behaviour, including an increase in wall friction at higher $Gr$, and a drop in the heat flux relative to natural convection for $Ri=O(1)$. This closely matched response is despite vastly different flow structures in the systems. No large-scale organisation is visible in visualisations of mixed vertical convection – an observation that is confirmed quantitatively by spectral analysis. This analysis, combined with a statistical description of the wall heat flux, highlights how moderate shear suppresses the growth of small-scale plumes and reduces the likelihood of extreme events in the local wall heat flux. Vice versa, starting from a pure shear flow, the addition of thermal driving enhances the drag due to the emission of thermal plumes.

We study the effects of Prandtl number $\mathit {Pr}$ and Rayleigh number $\mathit {Ra}$ in two-dimensional Rayleigh–Bénard convection without boundaries, i.e. with periodic boundary conditions. For Prandtl numbers in the range $10^{-3} \leqslant \mathit {Pr} \leqslant 10^2$, the viscous dissipation scales as $\epsilon _\nu \propto \mathit {Pr}^{1/2}\mathit {Ra}^{-1/4}$, which is based on the observation that enstrophy $\langle {\omega ^2}\rangle \propto \mathit {Pr}^0 \mathit {Ra}^{1/4}$, and the Nusselt number tends to follow the ‘ultimate’ scaling $\mathit {Nu} \propto \mathit {Pr}^{1/2}\mathit {Ra}^{1/2}$ for all values of $\mathit {Pr}$ considered. The inverse cascade of kinetic energy forms the power-law spectrum $\hat {E}_u(k) \propto k^{-2.3}$, which is close to $k^{-11/5}$ proposed by the Bolgiano–Obukhov (BO) scaling. The potential energy flux is not constant, in contrast to one of the main assumptions underlying the BO phenomenology. So, the direct cascade of potential energy forms the power-law spectrum $\hat {E}_\theta (k) \propto k^{-1.2}$, which deviates from the expected $k^{-7/5}$. Finally, at $\mathit {Pr} \to 0$ and $\infty$, we find that the dynamics is dominated by vertically oriented elevator modes that grow without bound, even at high Rayleigh numbers and with large-scale dissipation present.

We consider the two-dimensional Rayleigh–Bénard convection problem between Navier-slip fixed-temperature boundary conditions, and present a new upper bound for the Nusselt number (${\textit {Nu}}$). The result, based on a localization principle for the Nusselt number and an interpolation bound, exploits the regularity of the flow. On one hand our method yields a shorter proof of the celebrated result of Whitehead & Doering (Phys. Rev. Lett., vol. 106, 2011, 244501) in the case of free-slip boundary conditions. On the other hand, its combination with a new, refined estimate for the pressure gives a substantial improvement of the interpolation bounds in Drivas et al. (Phil. Trans. R. Soc. A, vol. 380, issue 2225, 2022, 20210025) for slippery boundaries. A rich description of the scaling behaviour arises from our result: depending on the magnitude of the Prandtl number (${\textit {Pr}}$) and slip length, our upper bounds indicate five possible scaling laws (where ${\textit {Ra}}$ is the Rayleigh number): ${\textit {Nu}}\sim (L_s^{-1}\,{\textit {Ra}})^{{1}/{3}}$, ${\textit {Nu}}\sim (L_s^{-2/5}\,{\textit {Ra}})^{{5}/{13}}$, ${\textit {Nu}}\sim {\textit {Ra}}^{{5}/{12}}$, ${\textit {Nu}}\sim {\textit {Pr}}^{-1/6}(L_s^{-4/3}\,{\textit {Ra}})^{{1}/{2}}$ and ${\textit {Nu}}\sim {\textit {Pr}}^{-1/6}(L_s^{-1/3}\,{\textit {Ra}})^{{1}/{2}}$.

We experimentally and numerically characterize rapidly rotating radiatively driven thermal convection, beyond the sole heat transport measurements reported by Bouillaut et al. (Proc. Natl Acad. Sci., vol. 118, 2021, e2105015118). Based on a suite of direct numerical simulations (DNS) and additional processing of the experimental data collected by Bouillaut et al. (Proc. Natl Acad. Sci., vol. 118, 2021, e2105015118), we report the simultaneous validation of the scaling predictions of the ‘geostrophic turbulence’ regime – the diffusivity-free or ‘ultimate’ regime of rapidly rotating convection – for the heat transport and the temperature fluctuations. Following such cross-validation between DNS and laboratory experiments, we further process the numerical data to validate the ‘geostrophic turbulence’ scaling predictions for the flow velocity and horizontal scale. Radiatively driven convection thus appears as a versatile set-up for the laboratory observation of the diffusivity-free regimes of various convective flows of geophysical and/or astrophysical interest.

We report on an experimental study of turbulent Rayleigh–Bénard convection with asymmetric top and bottom plates. The plates are covered with pyramid-shaped roughness elements whose aspect ratios are $\lambda =1$ or $\lambda =4$. In the low-Rayleigh-number regime ($Ra<1.9\times 10^9$), the heat transport efficiencies in the asymmetric cells, characterized by the Nusselt number, are smaller than those measured in a symmetric $\lambda = 1$ cell and are greater than those for a symmetric $\lambda = 4$ cell, whereas in the high-Rayleigh-number regime ($Ra>1.9\times 10^9$), the Nusselt numbers of the asymmetric cells are, in turn, greater than those for the symmetric cell with $\lambda = 1$ and smaller than those for the symmetric cell with $\lambda = 4$. In addition, the heat transports of individual plates are studied based on the temperature drops across both halves of the cell. In the low-$Ra$ regime, the $\lambda =1$ plate shows higher heat transfer than the $\lambda =4$ plate, while for the high-$Ra$ regime, the $\lambda =4$ plate shows a higher heat transport ability. In both regimes, the individual Nusselt number of the plate with lower heat transfer is insensitive to the topology of the other plate. Besides, it is found that the symmetry of the centre temperature distribution is robust to the symmetry breaking of boundary topographies. For the $Ra$ range explored, a weak temperature inversion is observed in the bulk of asymmetric rough cells. Finally, we remark that the temperature fluctuation at the cell centre and the Reynolds number associated with the large-scale circulation show universal power laws in terms of the flux Rayleigh number as $\sigma _{T_{c}}\sim Ra_F^{0.68}$ and $Re_{LSC}\sim Ra_F^{0.36}$, respectively.

We report a numerical investigation of a previously noticed but less explored flow state transition in two-dimensional turbulent Rayleigh–Bénard convection. The simulations are performed in a square domain over a Rayleigh number range of $10^7 \leq Ra \leq 2 \times 10^{11}$ and a Prandtl number range of $0.25 \leq Pr \leq 20$. The transition is characterized by the emergence of multiple satellite eddies with increasing $Ra$, which orbit around and interact with the main vortex roll in the system. Consequently, the main roll is squeezed to a smaller size compared with the domain and wanders around in the bulk region irregularly and extensively. This is in sharp contrast to the flow state before the transition, which is featured by a domain-sized circulatory roll with its vortex centre ‘condensed’ near the domain's centre. Detailed velocity field analysis reveals that there exists an abrupt increase in the energy fluctuations of the Fourier modes during the transition. Based on this phase-transition-like signal, the critical condition for the transition is found to follow a scaling relation as $Ra_t \sim Pr^{1.41}$ where $Ra_t$ is the critical Rayleigh number for the transition. This scaling relation is quantitatively explained by a phenomenological model grounded on the bistability behaviour (i.e. spontaneous and stochastic switching between the two flow states) observed at the edge of the transition. The model can also account for the effects of aspect ratio on the transition reported in the literature (van der Poel et al., Phys. Fluids, vol. 24, 2012).