An official website of the United States government
We experimentally explored the effect of single-sidewall cooling on Rayleigh–Bénard (RB) convection. Canonical RB was also studied to aid insight. The scenarios shared tank dimensions and bottom and top wall temperatures; the single sidewall cooling had the top wall temperature. Turbulence was explored at two canonical Rayleigh numbers, $$Ra=1.6\times 10^{10}$$ and $$Ra=2\times 10^9$$ under Prandtl number $Pr=5.4$ . Particle image velocimetry described vertical planes parallel and perpendicular to the sidewall cooling. The two $Ra$ scenarios reveal pronounced changes in the flow structure and large-scale circulation (LSC) due to the sidewall cooling. The density gradient induced by the sidewall cooling led to asymmetric descending and ascending flows and irregular LSC. Flow statistics departed from the canonical case, exhibiting lower buoyancy effects, represented by an effective Rayleigh number with effective height dependent on the distance from the lateral cooling. Velocity spectra show two scalings, $$\varPhi \propto f^{-5/3}$$ Kolmogorov (KO41) and $$\varPhi \propto f^{-11/5}$$ Bolgiano (BO59) in the larger $Ra$ ; the latter was not present in the smaller set-up. The BO59 scaling with sidewall cooling appears at higher frequencies than its canonical counterpart, suggesting weaker buoyancy effects. The LSC core motions allowed us to identify a characteristic time scale of the order of vortex turnover time associated with distinct vortex modes. The velocity spectra of the vortex core oscillation along its principal axis showed a scaling of $$\varPhi _c \propto f^{-5/3}$$ for the single sidewall cooling, which was dominant closer there. It did not occur in the canonical case, evidencing the modulation of LSC oscillation on the flow.
more »
« less
This content will become publicly available on August 10, 2026
Particle dynamics and dune formation in Rayleigh–Bénard convection: a particle-resolved simulation study
This paper presents numerical results for Rayleigh–Bénard convection with suspended particles at Rayleigh numbers $Ra=10^7$ and $10^8$, and unit Prandtl number. Accounting for their finite size makes it possible to investigate in detail the mechanism by which the particles, which are 10% heavier than the fluid, get resuspended after settling, thus maintaining a two-phase circulating flow. It is shown that an essential component of this mechanism is the formation of particle accumulations, or ‘dunes’, on the bottom of the Rayleigh–Bénard cell. Ascending plumes become localised on these dunes. Particles are dragged up the dune slopes, and when they reach the top, are entrained into the rising plumes. Direct resuspension of particles from the cell bottom, if it happens at all, is very rare. For $Ra=10^7$, aspect ratios (width/height) $$\Gamma =1,2,4$$ are considered. It is found that in these and in the other cases simulated, at steady state, a single dune evolves, the largest linear dimension of which is comparable to the cell size. A remarkable consequence is that even at the low volume fraction considered here, 3.27%, the particles are able to structure the flow and to determine the size and position of the largest ascending plumes. Their effect on the Nusselt number, however, remains small. This and other results are explained on the basis of the ratio of the cell-bottom viscous boundary-layer thickness to the particle diameter.
more »
« less
- Award ID(s):
- 2053204
- PAR ID:
- 10630890
- Publisher / Repository:
- Cambridge University Press
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 1016
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Steady two-dimensional Rayleigh–Bénard convection between stress-free isothermal boundaries is studied via numerical computations. We explore properties of steady convective rolls with aspect ratios $${\rm \pi} /5\leqslant \varGamma \leqslant 4{\rm \pi}$$ , where $$\varGamma$$ is the width-to-height ratio for a pair of counter-rotating rolls, over eight orders of magnitude in the Rayleigh number, $$10^3\leqslant Ra\leqslant 10^{11}$$ , and four orders of magnitude in the Prandtl number, $$10^{-2}\leqslant Pr\leqslant 10^2$$ . At large $Ra$ where steady rolls are dynamically unstable, the computed rolls display $$Ra \rightarrow \infty$$ asymptotic scaling. In this regime, the Nusselt number $Nu$ that measures heat transport scales as $$Ra^{1/3}$$ uniformly in $Pr$ . The prefactor of this scaling depends on $$\varGamma$$ and is largest at $$\varGamma \approx 1.9$$ . The Reynolds number $Re$ for large- $Ra$ rolls scales as $$Pr^{-1} Ra^{2/3}$$ with a prefactor that is largest at $$\varGamma \approx 4.5$$ . All of these large- $Ra$ features agree quantitatively with the semi-analytical asymptotic solutions constructed by Chini & Cox ( Phys. Fluids , vol. 21, 2009, 083603). Convergence of $Nu$ and $Re$ to their asymptotic scalings occurs more slowly when $Pr$ is larger and when $$\varGamma$$ is smaller.more » « less
-
Linear and nonlinear stability analyses are performed to determine critical Rayleigh numbers ( $${Ra}_{cr}$$ ) for a Rayleigh–Bénard convection configuration with an imposed bottom boundary heat flux that varies harmonically in time with zero mean. The $${Ra}_{cr}$$ value depends on the non-dimensional frequency $$\omega$$ of the boundary heat-flux modulation. Floquet theory is used to find $${Ra}_{cr}$$ for linear stability, and the energy method is used to find $${Ra}_{cr}$$ for two different types of nonlinear stability: strong and asymptotic. The most unstable linear mode alternates between synchronous and subharmonic frequencies at low $$\omega$$ , with only the latter at large $$\omega$$ . For a given frequency, the linear stability $${Ra}_{cr}$$ is generally higher than the nonlinear stability $${Ra}_{cr}$$ , as expected. For large $$\omega$$ , $${Ra}_{cr} \omega ^{-2}$$ approaches an $O(10)$ constant for linear stability but zero for nonlinear stability. Hence the domain for subcritical instability becomes increasingly large with increasing $$\omega$$ . The same conclusion is reached for decreasing Prandtl number. Changing temperature and/or velocity boundary conditions at the modulated or non-modulated plate leads to the same conclusions. These stability results are confirmed by selected direct numerical simulations of the initial value problem.more » « less
-
null (Ed.)We use well resolved numerical simulations with the lattice Boltzmann method to study Rayleigh–Bénard convection in cells with a fractal boundary in two dimensions for $Pr = 1$ and $$Ra \in \left [10^7, 10^{10}\right ]$$ , where Pr and Ra are the Prandtl and Rayleigh numbers. The fractal boundaries are functions characterized by power spectral densities $S(k)$ that decay with wavenumber, $$k$$ , as $$S(k) \sim k^{p}$$ ( $p < 0$ ). The degree of roughness is quantified by the exponent $$p$$ with $p < -3$ for smooth (differentiable) surfaces and $$-3 \le p < -1$$ for rough surfaces with Hausdorff dimension $$D_f=\frac {1}{2}(p+5)$$ . By computing the exponent $$\beta$$ using power law fits of $$Nu \sim Ra^{\beta }$$ , where $Nu$ is the Nusselt number, we find that the heat transport scaling increases with roughness through the top two decades of $$Ra \in \left [10^8, 10^{10}\right ]$$ . For $$p$$ $= -3.0$ , $-2.0$ and $-1.5$ we find $$\beta = 0.288 \pm 0.005, 0.329 \pm 0.006$$ and $$0.352 \pm 0.011$$ , respectively. We also find that the Reynolds number, $Re$ , scales as $$Re \sim Ra^{\xi }$$ , where $$\xi \approx 0.57$$ over $$Ra \in \left [10^7, 10^{10}\right ]$$ , for all $$p$$ used in the study. For a given value of $$p$$ , the averaged $Nu$ and $Re$ are insensitive to the specific realization of the roughness.more » « less
-
This study presents direct numerical simulations of turbulent Rayleigh–Bénard convection in non-colloidal suspensions, with special focus on the heat transfer modifications in the flow. Adopting a Rayleigh number of $10^8$ and Prandtl number of 7, parametric investigations of the particle volume fraction $$0\leq \varPhi \leq 40\,\%$$ and particle diameter $$1/20\leq d^*_p\leq 1/10$$ with respect to the cavity height, are carried out. The particles are neutrally buoyant, rigid spheres with physical properties that match the fluid phase. Up to $$\varPhi =25\,\%$$ , the Nusselt number increases weakly but steadily, mainly due to the increased thermal agitation that overcomes the decreased kinetic energy of the flow. Beyond $$\varPhi =30\,\%$$ , the Nusselt number exhibits a substantial drop, down to approximately 1/3 of the single-phase value. This decrease is attributed to the dense particle layering in the near-wall region, confirmed by the time-averaged local volume fraction. The dense particle layer reduces the convection in the near-wall region and negates the formation of any coherent structures within one particle diameter from the wall. Significant differences between $$\varPhi \leq 30\,\%$$ and 40 % are observed in all statistical quantities, including heat transfer and turbulent kinetic energy budgets, and two-point correlations. Special attention is also given to the role of particle rotation, which is shown to contribute to maintaining high heat transfer rates in moderate volume fractions. Furthermore, decreasing the particle size promotes the particle layering next to the wall, inducing a similar heat transfer reduction as in the highest particle volume fraction case.more » « less
