Noise source identification has been a long-standing challenge for decades. Although it is known that sound sources are closely related to flow structures, the underlying physical mechanisms remain controversial. This study develops a sound source identification method based on longitudinal and transverse process decomposition (LTD). Large-eddy simulations were performed on the flow around a cylinder at a Reynolds number of 3900. Using the new LTD method, sound sources in the cylinder flow were identified, and the mechanisms linking flow structures with noise generation were discussed in detail. Identifying the physical sound sources from two levels, low-order theory and high-order theory, the physical mechanism of wall sound sources was also analysed. Results indicate that the sound sources in the flow field mainly come from the leading edge, shear layer and wake region of the cylinder. The high-order theory reveals that sound sources are correlated with the spatio-temporal evolution of enstrophy, vortex stretching and surface deformation processes, this reflecting the coupling between transversal and longitudinal flow fields. The boundary thermodynamic flux and boundary dilatation flux distribution of the cylinder were analysed. Results indicate that the wall sound sources mainly come from the separation point and have a disorderly distribution on the leeward side of the cylinder, which is the main region where longitudinal variables enter the fluid from the wall surface, and the wall sound source is related to the boundary enstrophy flux.

We perform direct numerical simulations of centrifugal convection with an oscillating rotational velocity of small amplitude to study the effects of oscillatory boundary motion. The oscillation period is the main control parameter, with its range corresponding to a Womersley number in the range $1\lt Wo\lt 300$. Oscillating boundaries generate a circumferential shear flow, which significantly inhibits heat transfer, with maximum suppression $87\,\%$ observed in the present parameter space. Through analysis of the background flow, we find that as the oscillation period increases, the increasing penetration depth of the oscillation and weakening local shear strength result in non-monotonic changes in heat transfer. Under high-frequency oscillation, the characteristic length scale of the viscous layer induced by the oscillation is smaller than the convective length scale, and shear manifests primarily as a continuous suppression of the boundary layer. In contrast, under low-frequency oscillation, the shear flow covers the entire region but with weak strength. The suppression effect of such shear flow exhibits periodicity, leading to alternating phases of convection inhibition and convection generation. The present findings explore the physical mechanisms behind the suppression of convective heat transfer by oscillation, and offer a new strategy for controlling convection systems, with potential implications for both fundamental research and industrial applications.

This study presents a novel investigation into the vortex dynamics of flow around a near-wall rectangular cylinder based on direct numerical simulation at $Re=1000$, marking the first in-depth exploration of these phenomena. By varying aspect ratios ($L/D = 5$, $10$, $15$) and gap ratios ($G/D = 0.1$, $0.3$, $0.9$), the study reveals the vortex dynamics influenced by the near-wall effect, considering the incoming laminar boundary layer flow. Both $L/D$ and $G/D$ significantly influence vortex dynamics, leading to behaviours not observed in previous bluff body flows. As $G/D$ increases, the streamwise scale of the upper leading edge (ULE) recirculation grows, delaying flow reattachment. At smaller $G/D$, lower leading edge (LLE) recirculation is suppressed, with upper Kelvin–Helmholtz vortices merging to form the ULE vortex, followed by instability, differing from conventional flow dynamics. Larger $G/D$ promotes the formation of an LLE shear layer. An intriguing finding at $L/D = 5$ and $G/D = 0.1$ is the backward flow of fluid from the downstream region to the upper side of the cylinder. At $G/D = 0.3$, double-trailing-edge vortices emerge for larger $L/D$, with two distinct flow behaviours associated with two interactions between gap flow and wall recirculation. These interactions lead to different multiple flow separations. For $G/D = 0.9$, the secondary vortex (SV) from the plate wall induces the formation of a tertiary vortex from the lower side of the cylinder. Double-SVs are observed at $L/D = 5$. Frequency locking is observed in most cases, but is suppressed at $L/D = 10$ and $G/D = 0.9$, where competing shedding modes lead to two distinct evolutions of the SV.

We study the melting process of a solid under microgravity, driven solely by lateral vibrations that are perpendicular to the applied temperature gradient due to the absence of gravity-induced convection. Using direct numerical simulations with the phase-field method, we examine two-dimensional vibration-induced melting in a square cavity over four orders of magnitude of vibrational Rayleigh numbers, $10^5\le Ra_{{vib}}\le 10^9$. Our results show that as melting progresses, the flow structure transitions from a periodic-circulation regime with diffusion-dominated heat transfer to a columnar regime with vibroconvection. The mean height of the liquid–solid interface follows a power-law dependency with time, $\bar {\xi } \sim \tilde t^{1/(2-2\alpha )}$, where $\alpha = 0$ in the periodic-circulation regime and $\alpha = 1/2$ in the columnar regime. We further observe that within the columnar regime, the morphological evolution of the liquid–solid interface is influenced by the interaction of columnar thermal plumes in the central regions and the peripheral flow near the sidewalls. Specifically, we offer a comprehensive analysis of the plume merging behaviour, which is governed by the aspect ratio ($\bar {\xi }$) of the liquid layer and the intensity of vibration, quantified by the effective vibrational Rayleigh number $Ra_{vib}^{eff}$. We identify the relationship between the number of columnar plumes $K_m$ and $Ra_{vib}^{eff}$, finding that $K_m \sim \bar {\xi }^{-1} (Ra_{vib}^{eff})^{\gamma }$ with the fitting scaling exponent $\gamma = 0.150 \pm 0.025$. We subsequently quantify the characteristics of the interface roughness amplitude evolution in microgravity vibroconvection. Our results indicate that the roughness amplitude exhibits a power-law dependence on the mean height of the liquid layer. Drawing from the Stefan boundary condition, we theoretically deduce this dependence under the assumption of a non-uniform heat flux distribution at the interface, where the theory is corroborated by our numerical simulations.

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.

This Forum reflects the stimulating conversations we had as participants in a roundtable at the 2023 Boston meeting of the Association for Asian Studies. The discussion centred on the position of women (as both contributors and subjects) in competing visions of modernity in Southeast Asia. Our overall goals were to explore how women contributed to these visions and the extent to which their hopes were realised. In the process, numerous questions arose: What were women's experiences in these modernities? How were their roles influenced by class, ethnicity, religion, and age? And, finally, how do analyses of women-centred involvement inform broader understandings of social, cultural, economic, and political transformations in the region? The contributors reflect on their own trajectory in the study of women's history, the sources they have been able to leverage, the intersection between ‘modernity’ and nationalism, the related issue of women's position in postcolonial nation-states, and the question of how ideas of sexuality (in the family, society, politics) were changing. Throughout the Forum, we argue that the experience of women, and their (often but not necessarily agentic) contribution, is necessary to fully grasp the region's multiple entanglements with ‘modernity/ies’ in a variety of fields.

We report the unified constitutive law of vibroconvective turbulence in microgravity, i.e. $Nu \sim a^{-1} Re_{os}^\beta$ where the Nusselt number $Nu$ measures the global heat transport, $a$ is the dimensionless vibration amplitude, $Re_{os}$ is the oscillational Reynolds number and $\beta$ is the universal exponent. We find that the dynamics of boundary layers plays an essential role in vibroconvective heat transport and the $Nu$-scaling exponent $\beta$ is determined by the competition between the thermal boundary layer (TBL) and vibration-induced oscillating boundary layer (OBL). Then a physical model is proposed to explain the change of scaling exponent from $\beta =2$ in the TBL-dominant regime to $\beta = 4/3$ in the OBL-dominant regime. Our finding elucidates the emergence of universal constitutive laws in vibroconvective turbulence, and opens up a new avenue for generating a controllable effective heat transport under microgravity or even microfluidic environment in which the gravity effect is nearly absent.

Soft drink consumption has become a highly controversial public health issue. Given the pattern of consumption in China, sugar-sweetened beverage is the main type of soft drink consumed. Due to containing high levels of fructose, a soft drink may have a deleterious effect on handgrip strength (HGS) due to oxidative stress, inflammation and insulin resistance. However, few studies show an association between soft drink consumption and HGS in adults. We aimed to investigate the association between soft drink consumption and longitudinal changes in HGS among a Chinese adult population. A longitudinal population-based cohort study (5-year follow-up, median: 3·66 years) was conducted in Tianjin, China. A total of 11 125 participants (56·7 % men) were enrolled. HGS was measured using a handheld digital dynamometer. Soft drink consumption (mainly sugar-containing carbonated beverages) was measured at baseline using a validated FFQ. ANCOVA was used to evaluate the association between soft drink consumption and annual change in HGS or weight-adjusted HGS. After adjusting for multiple confounding factors, the least square means (95 % CI) of annual change in HGS across soft drink consumption frequencies were −0·70 (–2·49, 1·09) for rarely drinks, −0·82 (–2·62, 0·97) for < 1 cup/week and −0·86 (–2·66, 0·93) for ≥ 1 cup/week (Pfor trend < 0·05). Likewise, a similar association was observed between soft drink consumption and annual change in weight-adjusted HGS. The results indicate that higher soft drink consumption was associated with faster HGS decline in Chinese adults.

The numerical investigation focuses on the flow patterns around a rectangular cylinder with three aspect ratios ($L/D=5$, $10$, $15$) at a Reynolds number of $1000$. The study delves into the dynamics of vortices, their associated frequencies, the evolution of the boundary layer and the decay of the wake. Kelvin–Helmholtz (KH) vortices originate from the leading edge (LE) shear layer and transform into hairpin vortices. Specifically, at $L/D=5$, three KH vortices merge into a single LE vortex. However, at $L/D=10$ and $15$, two KH vortices combine to form a LE vortex, with the rapid formation of hairpin vortex packets. A fractional harmonic arises due to feedback from the split LE shear layer moving upstream, triggering interaction with the reverse flow. Trailing edge (TE) vortices shed, creating a Kármán-like street in the wake. The intensity of wake oscillation at $L/D=5$ surpasses that in the other two cases. Boundary layer transition occurs after the saturation of disturbance energy for $L/D=10$ and $15$, but not for $L/D=5$. The low-frequency disturbances are selected to generate streaks inside the boundary layer. The TE vortex shedding induces the formation of a favourable pressure gradient, accelerating the flow and fostering boundary layer relaminarization. The self-similarity of the velocity defect is observed in all three wakes, accompanied by the decay of disturbance energy. Importantly, the decrease in the shedding frequency of LE (TE) vortices significantly contributes to the overall decay of disturbance energy. This comprehensive exploration provides insights into complex flow phenomena and their underlying dynamics.

This study investigates the effect of vibration on the flow structure transitions in thermal vibrational convection (TVC) systems, which occur when a fluid layer with a temperature gradient is excited by vibration. Direct numerical simulation (DNS) of TVC in a two-dimensional enclosed square box is performed over a range of dimensionless vibration amplitudes $0.001 \le a \le 0.3$ and angular frequencies $10^{2} \le \omega \le 10^{7}$, with a fixed Prandtl number of 4.38. The flow visualisation shows the transition behaviour of flow structure upon the varying frequency, characterising three distinct regimes, which are the periodic-circulation regime, columnar regime and columnar-broken regime. Different statistical properties are distinguished from the temperature and velocity fluctuations at the boundary layer and mid-height. Upon transition into the columnar regime, columnar thermal coherent structures are formed, in contrast to the periodic oscillating circulation. These columns are contributed by the merging of thermal plumes near the boundary layer, and the resultant thermal updrafts remain at almost fixed lateral position, leading to a decrease in fluctuations. We further find that the critical point of this transition can be described nicely by the vibrational Rayleigh number ${{Ra}}_{vib}$. As the frequency continues to increase, entering the so-called columnar-broken regime, the columnar structures are broken, and eventually the flow state becomes a large-scale circulation (LSC), characterised by a sudden increase in fluctuations. Finally, a phase diagram is constructed to summarise the flow structure transition over a wide range of vibration amplitude and frequency parameters.

Intertemporal choices involve tradeoffs between outcomes that occur at different times. Most of the research has used pure gains tasks and the discount rates yielding from those tasks to explain and predict real-world behaviors and consequences. However, real decisions are often more complex and involve mixed outcomes (e.g., sooner-gain and later-loss or sooner-loss and later-gain). No study has used mixed gain-loss intertemporal tradeoff tasks to explain and predict real-world behaviors and consequences, and studies involving such tasks are also scarce. Considering that tasks involving a combination of gains and losses may yield different discount rates and that existing pure gains tasks do not explain or predict real-world outcomes well, this study conducted two experiments to compare the discount rates of mixed gain-loss intertemporal tradeoffs with those of pure gains or pure losses (Experiment 1) and to examine whether these tasks predicted different real-world behaviors and consequences (Experiment 2). Experiment 1 suggests that the discount rate ordering of the four tasks was, from highest to lowest, pure gains, sooner-loss and later-gain, pure losses, and sooner-gain and later-loss. Experiment 2 indicates that the evidence supporting the claim that the discount rates of the four tasks were related to different real-world behaviors and consequences was insufficient.

We investigate the energy transport and heat transfer efficiency in turbulent Rayleigh–Bénard (RB) convection laden with radiatively heated inertial particles. Direct numerical simulations combined with the Lagrangian point-particle mode were carried out in the range of density ratio $854.7\le \rho _p/\rho _0 \le 8547$ and radiation intensity $1\le \phi /\phi _{solar}\le 100$ for both two-dimensional (2-D) and three-dimensional (3-D) simulations. The Rayleigh number ranges from $2\times 10^6$ to $10^8$ for 2-D cases, and is $10^7$ for 3-D cases for $Pr=0.71$. It is found that particles with small density ratio that encounter strong radiation significantly alter the flow momentum transport and fluid heat transfer, so the fluid temperature of bulk is remarkably heated. We then derived the theoretical relation of the Nusselt number for interphase heat transfer in the heated particle-laden RB convection, which reveals that the heat transfer difference between the top and bottom plates stems from the interphase heat transfer. We further found that both the interphase heat transfer and the interphase thermal energy transport exhibit universal properties. They are both increased linearly with the reciprocal of the normalized density ratio. Additionally, both the interphase heat transfer and the interphase thermal energy transport increase linearly with the increase of radiation intensity. The growth rates exhibit specific scaling relations versus Rayleigh number and density ratio. Two different regimes distinguished by the critical density ratio, i.e. the exothermic particle regime and the endothermic particle regime, are observed. We further derived the power-law relation of the critical density ratios versus Rayleigh number and radiation intensity, i.e. $\rho _p/\rho _c \sim (\phi /\phi _{solar})^{1/2}\,Ra^{1/3}$, which is in remarkable agreement with the 3-D simulations.

The flow past a cylinder in proximity to a plane wall is investigated numerically for small gap ratios. Three vortex dynamic processes associated with different hairpin vortex generation mechanisms are identified for the first time, and the wake-induced turbulent transition is analysed. The vortex shedding is suppressed at $G/D = 0.1$, while the spanwise vortex is generated via a Kelvin–Helmholtz instability and evolves into hairpin vortices. For $G/D= 0.3$, the upper and lower rollers alternatively shedding from the cylinder, interact with the secondary vortex. The split secondary vortex merges with the upper roller and results in a new vortex downstream, which develops into hairpin vortices. When $G/D = 0.9$, the secondary vortex interacts with the lower roller and then evolves into hairpin vortices. A tertiary vortex induced by the secondary vortex is observed, rotating in the opposite direction to the secondary vortex the wake-induced transitions share the same route. The velocity fluctuations deviate from the optimal growth theory in the pre-transitional region. In the transitional region low-frequency disturbances penetrate the sheltering edge to generate streaks where the disturbance energy declines. In the turbulent region the logarithmic layer is formed, indicating that the turbulent equilibrium is established.

We report that vertical vibration with small amplitude and high frequency can tame convective heat transport in Rayleigh–Bénard convection in a turbulent regime. When vertical vibration is applied, a dynamically averaged ‘anti-gravity’ results that stabilizes the thermal boundary layer and inhibits the eruption of thermal plumes. This eventually leads to the attenuation of the intensity of large-scale mean flow and a significant suppression of turbulent heat transport. Accounting for both the thermally led buoyancy and the vibration-induced anti-gravitational effects, we propose an effective Rayleigh number that helps to extend the Grossmann–Lohse theory to thermal vibrational turbulence. The prediction of the reduction on both the Nusselt and Reynolds numbers obtained by the extended model is found to agree well with the numerical data. In addition, vibrational influences on the mean flow structure and the temporal evolution of Nusselt and Reynolds numbers are investigated. The non-uniform characteristic of vibration-induced ‘anti-gravity’ is discussed. The present findings provide a powerful basis for studying thermal vibrational turbulence and put forward a novel strategy for actively controlling thermal turbulence.

In this paper, we report that reversals of the large-scale circulation in two-dimensional Rayleigh–Bénard (RB) convection can be suppressed by imposing sinusoidally distributed heating to the bottom plate. We examine how the frequency of flow reversals depends on the dimensionless wavenumber $k$ of the spatial temperature modulation with various modulation amplitude $A$. For sufficiently large $k$, the flow reversal frequency is close to that in the standard RB convection under uniform heating. However, when $k$ decreases, the frequency of flow reversal gradually becomes lower and can even be largely reduced. Furthermore, we examine the growth of the corner roll and the global flow structure based on Fourier mode decomposition, and reveal that the size of the corner roll diminishes as the wavenumber decreases. The reason is that the regions occupied by the cold phase can absorb heat from the hot plumes and thus lower their temperature, which reduces the corner roll's kinetic energy input provided by the buoyancy force, and weakens the feeding process of the corner rolls. This results in the locking of the corner roll into a smaller region near the corner, making it harder for a reversal to occur. Using the concept of horizontal convection caused by non-uniform heating, we find a relevant parameter $k/A$ to describe briefly how the reversal frequency depends on wavenumber and modulation amplitude. The present work provides a new way to control the flow reversals in RB convection through modifying temperature boundary conditions.

Let $\mathbb {N}$ be the set of all nonnegative integers. For $S\subseteq \mathbb {N}$ and $n\in \mathbb {N}$ , let $R_S(n)$ denote the number of solutions of the equation $n=s_1+s_2$ , $s_1,s_2\in S$ and $s_1<s_2$ . Let A be the set of all nonnegative integers which contain an even number of digits $1$ in their binary representations and $B=\mathbb {N}\setminus A$ . Put $A_l=A\cap [0,2^l-1]$ and $B_l=B\cap [0,2^l-1]$ . We prove that if $C \cup D=[0, m]\setminus \{r\}$ with $0<r<m$ , $C \cap D=\emptyset $ and $0 \in C$ , then $R_{C}(n)=R_{D}(n)$ for any nonnegative integer n if and only if there exists an integer $l \geq 1$ such that $m=2^{l}$ , $r=2^{l-1}$ , $C=A_{l-1} \cup (2^{l-1}+1+B_{l-1})$ and $D=B_{l-1} \cup (2^{l-1}+1+A_{l-1})$ . Kiss and Sándor [‘Partitions of the set of nonnegative integers with the same representation functions’, Discrete Math. 340 (2017), 1154–1161] proved an analogous result when $C\cup D=[0,m]$ , $0\in C$ and $C\cap D=\{r\}$ .

External modulation on thermal convection has been studied extensively to achieve the control of flow structures and heat-transfer efficiency. In this paper, we carry out direct numerical simulations on Rayleigh–Bénard convection accounting for both the modulation of wall shear and roughness over the Rayleigh number range $1.0 \times 10^6 \le Ra \le 1.0 \times 10^8$, the wall shear Reynolds number range $0 \le Re_w \le 5000$, the aspect-ratio range $2 \le \varGamma \le 4{\rm \pi}$, and the dimensionless roughness height range $0 \le h \le 0.2$ at fixed Prandtl number $Pr = 1$. Under the combined actions of wall shear and roughness, with increasing $Re_w$, the heat flux is initially enhanced in the buoyancy-dominant regime, then has an abrupt transition near the critical shear Reynolds number $Re_{w,cr}$, and finally enters the purely diffusion regime dominated by shear. Based on the crossover of the kinetic energy production between the buoyancy-dominant and shear-dominant regimes, a physical model is proposed to predict the transitional scaling behaviour between $Re_{w,cr}$ and $Ra$, i.e. $Re_{w,cr} \sim Ra^{9/14}$, which agrees well with our numerical results. The reason for the observed heat-transport enhancement in the buoyancy-dominant regime is further explained by the fact that the moving rough plates introduce an external shear to strengthen the large-scale circulation (LSC) in the vertical direction and serve as a conveyor belt to increase the chances of the interaction between the LSC and secondary flows within cavities, which triggers more thermal plumes, efficiently transports the trapped hot (cold) fluids outside cavities.

We investigate the dynamic couplings between particles and fluid in turbulent Rayleigh–Bénard (RB) convection laden with isothermal inertial particles. Direct numerical simulations combined with the Lagrangian point-particle mode were carried out in the range of Rayleigh number $1\times 10^6 \le {Ra}\le 1 \times 10^8$ at Prandtl number ${Pr}=0.678$ for three Stokes numbers ${St_f}=1 \times 10^{-3}$, $8 \times 10^{-3}$ and $2.5 \times 10^{-2}$. It is found that the global heat transfer and the strength of turbulent momentum transfer are altered a small amount for the small Stokes number and large Stokes number as the coupling between the two phases is weak, whereas they are enhanced a large amount for the medium Stokes number due to strong coupling of the two phases. We then derived the exact relation of kinetic energy dissipation in the particle-laden RB convection to study the budget balance of induced and dissipated kinetic energy. The strength of the dynamic coupling can be clearly revealed from the percentage of particle-induced kinetic energy over the total induced kinetic energy. We further derived the power law relation of the averaged particles settling rate versus the Rayleigh number, i.e. $S_p/(d_p/H)^2{\sim} Ra^{1/2}$, which is in remarkable agreement with our simulation. We found that the settling and preferential concentration of particles are strongly correlated with the coupling mechanisms.

In this paper, we designed two different configurations with locally isothermal sidewalls, where the temperature is set to be the bulk temperature, to control the large-scale circulation in turbulent Rayleigh–Bénard convection, namely two-point control and four-point control. At fixed Rayleigh number $Ra=10^8$ and Prandtl number $Pr=2$, a series of direct numerical simulations are performed on both two-dimensional (2-D) and quasi-two-dimensional (quasi-2-D) cavities with both types of control, where the width of the control area is fixed at $\delta _c=0.05$ and the vertical distance from the cavity centre $h_c$ varies from 0 to 0.45 with an interval of 0.05. Our results show that the control effect depends on $h_c$, the control configurations as well as the flow dimensions. For 2-D cavities, both two-point control and four-point control suppress the flow reversal when $h_c \geq 0.05$, accompanied by the enhancement of vertical heat transfer and the strength of the large-scale circulation. For quasi-2-D cavities, the suppression of the flow reversals is obvious with two-point control and $h_c\geq 0.05$, while the effect is rather limited with four-point control. Further experiments with $Pr=5.7$ and $Ra$ up to $7.36\times10^8$ show that two-point control with $h_c=0.15$ can effectively suppress the flow reversal, while two-point control with $h_c=0$ can suppress the reversals at low $Ra=1.93\times 10^8$ and activate them at higher $Ra=7.36\times 10^8$, which agrees well with our numerical simulations.

We carry out direct numerical simulations of turbulent Rayleigh–Bénard convection in a square box with rough conducting plates over the Rayleigh number range $10^7\leqslant Ra\leqslant 10^9$ and the Prandtl number range $0.01\leqslant Pr\leqslant 100$. In Zhang et al. (J. Fluid Mech., vol. 836, 2018, R2), it was reported that while the measured Nusselt number $Nu$ is enhanced at large roughness height $h$, the global heat transport is reduced at small $h$. The division between the two regimes yields a critical roughness height $h_c$, and we now focus on the effects of the Prandtl number ($Pr$) on $h_c$. Based on the variations of $h_c$, we identify three regimes for $h_c(Pr)$. For low $Pr$, thermal boundary layers become thinner with increasing $Pr$. This makes the boundary layers easier to be disrupted by rough elements, leading to the decrease of $h_c$ with increasing $Pr$. For moderate $Pr$, the corner-flow rolls become much more pronounced and suppress the global heat transport via the competition between the corner-flow rolls and the large-scale circulation (LSC). As a consequence, $h_c$ increases with increasing $Pr$ due to the intensification of the corner–LSC competition. For high $Pr$, the convective flow transitions to the plume-controlled regime. As the rough elements trigger much stronger and more frequent plume emissions, $h_c$ again decreases with increasing $Pr$.