Comparison of the quasi-steady-state heat transport in phase-change and classical Rayleigh-Bénard convection for a wide range of Stefan number and Rayleigh number

2019 ◽  
Vol 31 (9) ◽  
pp. 096605 ◽  
Author(s):  
Ojas Satbhai ◽  
Subhransu Roy ◽  
Sudipto Ghosh ◽  
Suman Chakraborty ◽  
Rajaram Lakkaraju
Author(s):  
Ojas Satbhai ◽  
Subhransu Roy ◽  
Sudipto Ghosh

Direct numerical simulations for low Prandtl number fluid (Pr = 0.0216) are used to study the steady-state Rayleigh–Bénard convection (RB) in a two-dimensional unit aspect ratio box. The steady-state RB convection is characterized by analyzing the time-averaged temperature-field, and flow field for a wide range of Rayleigh number (2.1 × 105 ⩽ Ra ⩽ 2.1 × 108). It is seen that the time-averaged and space-averaged Nusselt number (Nuh¯) at the hot-wall monotonically increases with the increase in Rayleigh number (Ra) and the results show a power law scaling Nuh¯∝Ra0.2593. The current Nusselt number results are compared with the results available in the literature. The complex flow is analyzed by studying the frequency power spectra of the steady-state signal of the vertical velocity at the midpoint of the box for different Ra and probability density function of dimensionless temperature at various locations along the midline of the box.


2018 ◽  
Vol 846 ◽  
pp. 5-36 ◽  
Author(s):  
Stéphane Labrosse ◽  
Adrien Morison ◽  
Renaud Deguen ◽  
Thierry Alboussière

Solid-state convection can take place in the rocky or icy mantles of planetary objects, and these mantles can be surrounded above or below or both by molten layers of similar composition. A flow towards the interface can proceed through it by changing phase. This behaviour is modelled by a boundary condition taking into account the competition between viscous stress in the solid, which builds topography of the interface with a time scale $\unicode[STIX]{x1D70F}_{\unicode[STIX]{x1D702}}$, and convective transfer of the latent heat in the liquid from places of the boundary where freezing occurs to places of melting, which acts to erase topography, with a time scale $\unicode[STIX]{x1D70F}_{\unicode[STIX]{x1D719}}$. The ratio $\unicode[STIX]{x1D6F7}=\unicode[STIX]{x1D70F}_{\unicode[STIX]{x1D719}}/\unicode[STIX]{x1D70F}_{\unicode[STIX]{x1D702}}$ controls whether the boundary condition is the classical non-penetrative one ($\unicode[STIX]{x1D6F7}\rightarrow \infty$) or allows for a finite flow through the boundary (small $\unicode[STIX]{x1D6F7}$). We study Rayleigh–Bénard convection in a plane layer subject to this boundary condition at either or both its boundaries using linear and weakly nonlinear analyses. When both boundaries are phase-change interfaces with equal values of $\unicode[STIX]{x1D6F7}$, a non-deforming translation mode is possible with a critical Rayleigh number equal to $24\unicode[STIX]{x1D6F7}$. At small values of $\unicode[STIX]{x1D6F7}$, this mode competes with a weakly deforming mode having a slightly lower critical Rayleigh number and a very long wavelength, $\unicode[STIX]{x1D706}_{c}\sim 8\sqrt{2}\unicode[STIX]{x03C0}/3\sqrt{\unicode[STIX]{x1D6F7}}$. Both modes lead to very efficient heat transfer, as expressed by the relationship between the Nusselt and Rayleigh numbers. When only one boundary is subject to a phase-change condition, the critical Rayleigh number is $\mathit{Ra}_{c}=153$ and the critical wavelength is $\unicode[STIX]{x1D706}_{c}=5$. The Nusselt number increases approximately two times faster with the Rayleigh number than in the classical case with non-penetrative conditions, and the average temperature diverges from $1/2$ when the Rayleigh number is increased, towards larger values when the bottom boundary is a phase-change interface.


2016 ◽  
Vol 805 ◽  
Author(s):  
Kai Leong Chong ◽  
Ke-Qing Xia

We study the effect of severe geometrical confinement in Rayleigh–Bénard convection with a wide range of width-to-height aspect ratio $\unicode[STIX]{x1D6E4}$, $1/128\leqslant \unicode[STIX]{x1D6E4}\leqslant 1$, and Rayleigh number $Ra$, $3\times 10^{4}\leqslant Ra\leqslant 1\times 10^{11}$, at a fixed Prandtl number of $Pr=4.38$ by means of direct numerical simulations in Cartesian geometry with no-slip walls. For convection under geometrical confinement (decreasing $\unicode[STIX]{x1D6E4}$ from 1), three regimes can be recognized (Chong et al., Phys. Rev. Lett., vol. 115, 2015, 264503) based on the global and local properties in terms of heat transport, plume morphology and flow structures. These are Regime I: classical boundary-layer-controlled regime; Regime II: plume-controlled regime; and Regime III: severely confined regime. The study reveals that the transition into Regime III leads to totally different heat and momentum transport scalings and flow topology from the classical regime. The convective heat transfer scaling, in terms of the Nusselt number $Nu$, exhibits the scaling $Nu-1\sim Ra^{0.61}$ over three decades of $Ra$ at $\unicode[STIX]{x1D6E4}=1/128$, which contrasts sharply with the classical scaling $Nu-1\sim Ra^{0.31}$ found at $\unicode[STIX]{x1D6E4}=1$. The flow in Regime III is found to be dominated by finger-like, long-lived plume columns, again in sharp contrast with the mushroom-like, fragmented thermal plumes typically observed in the classical regime. Moreover, we identify a Rayleigh number for regime transition, $Ra^{\ast }=(29.37/\unicode[STIX]{x1D6E4})^{3.23}$, such that the scaling transition in $Nu$ and $Re$ can be clearly demonstrated when plotted against $Ra/Ra^{\ast }$.


2012 ◽  
Vol 710 ◽  
pp. 260-276 ◽  
Author(s):  
Quan Zhou ◽  
Bo-Fang Liu ◽  
Chun-Mei Li ◽  
Bao-Chang Zhong

AbstractWe report high-precision measurements of the Nusselt number $Nu$ as a function of the Rayleigh number $Ra$ in water-filled rectangular Rayleigh–Bénard convection cells. The horizontal length $L$ and width $W$ of the cells are 50.0 and 15.0 cm, respectively, and the heights $H= 49. 9$, 25.0, 12.5, 6.9, 3.5, and 2.4 cm, corresponding to the aspect ratios $({\Gamma }_{x} \equiv L/ H, {\Gamma }_{y} \equiv W/ H)= (1, 0. 3)$, $(2, 0. 6)$, $(4, 1. 2)$, $(7. 3, 2. 2)$, $(14. 3, 4. 3)$, and $(20. 8, 6. 3)$. The measurements were carried out over the Rayleigh number range $6\ensuremath{\times} 1{0}^{5} \lesssim Ra\lesssim 1{0}^{11} $ and the Prandtl number range $5. 2\lesssim Pr\lesssim 7$. Our results show that for rectangular geometry turbulent heat transport is independent of the cells’ aspect ratios and hence is insensitive to the nature and structures of the large-scale mean flows of the system. This is slightly different from the observations in cylindrical cells where $Nu$ is found to be in general a decreasing function of $\Gamma $, at least for $\Gamma = 1$ and larger. Such a difference is probably a manifestation of the finite plate conductivity effect. Corrections for the influence of the finite conductivity of the top and bottom plates are made to obtain the estimates of $N{u}_{\infty } $ for plates with perfect conductivity. The local scaling exponents ${\ensuremath{\beta} }_{l} $ of $N{u}_{\infty } \ensuremath{\sim} R{a}^{{\ensuremath{\beta} }_{l} } $ are calculated and found to increase from 0.243 at $Ra\simeq 9\ensuremath{\times} 1{0}^{5} $ to 0.327 at $Ra\simeq 4\ensuremath{\times} 1{0}^{10} $.


2011 ◽  
Vol 673 ◽  
pp. 318-348 ◽  
Author(s):  
G. FREUND ◽  
W. PESCH ◽  
W. ZIMMERMANN

Motivated by recent experiments, we study a rich variation of the familiar Rayleigh–Bénard convection (RBC), where the temperature at the lower boundary varies sinusoidally about a mean value. As usual the Rayleigh number R measures the average temperature gradient, while the additional spatial modulation is characterized by a (small) amplitude δm and a wavevector qm. Our analysis relies on precise numerical solutions of suitably adapted Oberbeck–Boussinesq equations (OBE). In the absence of forcing (δm = 0), convection rolls with wavenumber qc bifurcate only for R above the critical Rayleigh number Rc. In contrast, for δm≠0, convection is unavoidable for any finite R; in the most simple case in the form of ‘forced rolls’ with wavevector qm. According to our first comprehensive stability diagram of these forced rolls in the qm – R plane, they develop instabilities against resonant oblique modes at R ≲ Rc in a wide range of qm/qc. Only for qm in the vicinity of qc, the forced rolls remain stable up to fairly large R > Rc. Direct numerical simulations of the OBE support and extend the findings of the stability analysis. Moreover, we are in line with the experimental results and also with some earlier theoretical results on this problem, based on asymptotic expansions in the limit δm → 0 and R → Rc. It is satisfying that in many cases the numerical results can be directly interpreted in terms of suitably constructed amplitude and generalized Swift–Hohenberg equations.


2012 ◽  
Vol 707 ◽  
pp. 241-259 ◽  
Author(s):  
Jared P. Whitehead ◽  
Charles R. Doering

AbstractRigorous bounds on heat transport are derived for thermal convection between stress-free horizontal plates. For three-dimensional Rayleigh–Bénard convection at infinite Prandtl number ($\mathit{Pr}$), the Nusselt number ($\mathit{Nu}$) is bounded according to $\mathit{Nu}\leq 0. 28764{\mathit{Ra}}^{5/ 12} $ where $\mathit{Ra}$ is the standard Rayleigh number. For convection driven by a uniform steady internal heat source between isothermal boundaries, the spatially and temporally averaged (non-dimensional) temperature is bounded from below by $\langle T\rangle \geq 0. 6910{\mathit{R}}^{\ensuremath{-} 5/ 17} $ in three dimensions at infinite $\mathit{Pr}$ and by $\langle T\rangle \geq 0. 8473{\mathit{R}}^{\ensuremath{-} 5/ 17} $ in two dimensions at arbitrary $\mathit{Pr}$, where $\mathit{R}$ is the heat Rayleigh number proportional to the injected flux.


2017 ◽  
Vol 836 ◽  
Author(s):  
Yi-Zhao Zhang ◽  
Chao Sun ◽  
Yun Bao ◽  
Quan Zhou

Rough surfaces have been widely used as an efficient way to enhance the heat-transfer efficiency in turbulent thermal convection. In this paper, however, we show that roughness does not always mean a heat-transfer enhancement, but in some cases it can also reduce the overall heat transport through the system. To reveal this, we carry out numerical investigations of turbulent Rayleigh–Bénard convection over rough conducting plates. Our study includes two-dimensional (2D) simulations over the Rayleigh number range $10^{7}\leqslant Ra\leqslant 10^{11}$ and three-dimensional (3D) simulations at $Ra=10^{8}$. The Prandtl number is fixed to $Pr=0.7$ for both the 2D and the 3D cases. At a fixed Rayleigh number $Ra$, reduction of the Nusselt number $Nu$ is observed for small roughness height $h$, whereas heat-transport enhancement occurs for large $h$. The crossover between the two regimes yields a critical roughness height $h_{c}$, which is found to decrease with increasing $Ra$ as $h_{c}\sim Ra^{-0.6}$. Through dimensional analysis, we provide a physical explanation for this dependence. The physical reason for the $Nu$ reduction is that the hot/cold fluid is trapped and accumulated inside the cavity regions between the rough elements, leading to a much thicker thermal boundary layer and thus impeding the overall heat flux through the system.


2005 ◽  
Vol 19 (28n29) ◽  
pp. 1707-1710
Author(s):  
CHUANSHAN DAI

Two new methods for correlating experimental data of Rayleigh-Bénard convection in a fluid layer dispersed with phase-change-material particles have been proposed. Instead of the conventional arithmetic mean temperature method, two new methods for evaluating the intensity of convection in a horizontal fluid layer with strongly temperature-dependent specific heat were given. One method is by introducing a modified Stefan number in a Nu-Ra correlation or using the integrally averaged specific heat across the layer for calculating the Rayleigh number. The other method is by using the supercriticality (Ra/Rac) as an intensity parameter, in which the critical Rayleigh number Rac can be calculated by using the linear stability theory. Two models of non-dimensionalized apparent specific heat functions with respect to temperature were introduced. Results show that the two proposed methods are more rational comparing with the convectional 1/2-height rule method.


2002 ◽  
Vol 473 ◽  
pp. 191-199 ◽  
Author(s):  
JESSE OTERO ◽  
RALF W. WITTENBERG ◽  
RODNEY A. WORTHING ◽  
CHARLES R. DOERING

We formulate a bounding principle for the heat transport in Rayleigh–Bénard convection with fixed heat flux through the boundaries. The heat transport, as measured by a conventional Nusselt number, is inversely proportional to the temperature drop across the layer and is bounded above according to Nu [les ] cRˆ1/3, where c < 0.42 is an absolute constant and Rˆ = αγβh4/(νκ) is the ‘effective’ Rayleigh number, the non-dimensional forcing scale set by the imposed heat flux κβ. The relation among the parameter Rˆ, the Nusselt number, and the conventional Rayleigh number defined in terms of the temperature drop across the layer, is NuRa = Rˆ, yielding the bound Nu [les ] c3/2Ra1/2.


Sign in / Sign up

Export Citation Format

Share Document