Available potential energy in Rayleigh–Bénard convection

2013 ◽  
Vol 729 ◽  
Author(s):  
Graham O. Hughes ◽  
Bishakhdatta Gayen ◽  
Ross W. Griffiths

AbstractThe mechanical energy budget for thermally equilibrated Rayleigh–Bénard convection is developed theoretically, with explicit consideration of the role of available potential energy, this being the form in which all the mechanical energy for the flow is supplied. The analysis allows derivation for the first time of a closed analytical expression relating the rate of mixing in symmetric fully developed convection to the rate at which available potential energy is supplied by the thermal forcing. Only about half this supplied energy is dissipated viscously. The remainder is consumed by mixing acting to homogenize the density field. This finding is expected to apply over a wide range of Rayleigh and Prandtl numbers for which the Nusselt number is significantly greater than unity. Thus convection at large Rayleigh number involves energetically efficient mixing of density variations. In contrast to conventional approaches to Rayleigh–Bénard convection, the dissipation of temperature or density variance is shown not to be of direct relevance to the mechanical energy budget. Thus, explicit recognition of available potential energy as the source of mechanical energy for convection, and of both mixing and viscous dissipation as the sinks of this energy, could be of further use in understanding the physics.

2018 ◽  
Vol 841 ◽  
pp. 825-850 ◽  
Author(s):  
Chong Shen Ng ◽  
Andrew Ooi ◽  
Detlef Lohse ◽  
Daniel Chung

Previous numerical studies on homogeneous Rayleigh–Bénard convection, which is Rayleigh–Bénard convection (RBC) without walls, and therefore without boundary layers, have revealed a scaling regime that is consistent with theoretical predictions of bulk-dominated thermal convection. In this so-called asymptotic regime, previous studies have predicted that the Nusselt number ($\mathit{Nu}$) and the Reynolds number ($\mathit{Re}$) vary with the Rayleigh number ($\mathit{Ra}$) according to $\mathit{Nu}\sim \mathit{Ra}^{1/2}$ and $\mathit{Re}\sim \mathit{Ra}^{1/2}$ at small Prandtl numbers ($\mathit{Pr}$). In this study, we consider a flow that is similar to RBC but with the direction of temperature gradient perpendicular to gravity instead of parallel to it; we refer to this configuration as vertical natural convection (VC). Since the direction of the temperature gradient is different in VC, there is no exact relation for the average kinetic dissipation rate, which makes it necessary to explore alternative definitions for $\mathit{Nu}$, $\mathit{Re}$ and $\mathit{Ra}$ and to find physical arguments for closure, rather than making use of the exact relation between $\mathit{Nu}$ and the dissipation rates as in RBC. Once we remove the walls from VC to obtain the homogeneous set-up, we find that the aforementioned $1/2$-power-law scaling is present, similar to the case of homogeneous RBC. When focusing on the bulk, we find that the Nusselt and Reynolds numbers in the bulk of VC too exhibit the $1/2$-power-law scaling. These results suggest that the $1/2$-power-law scaling may even be found at lower Rayleigh numbers if the appropriate quantities in the turbulent bulk flow are employed for the definitions of $\mathit{Ra}$, $\mathit{Re}$ and $\mathit{Nu}$. From a stability perspective, at low- to moderate-$\mathit{Ra}$, we find that the time evolution of the Nusselt number for homogenous vertical natural convection is unsteady, which is consistent with the nature of the elevator modes reported in previous studies on homogeneous RBC.


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 }$.


2020 ◽  
Vol 887 ◽  
Author(s):  
Gerrit Green ◽  
Dimitar G. Vlaykov ◽  
Juan Pedro Mellado ◽  
Michael Wilczek


2018 ◽  
Vol 22 (Suppl. 2) ◽  
pp. 535-545
Author(s):  
Weihua Cai ◽  
Zhifeng Zheng ◽  
Changye Huang ◽  
Yue Wang ◽  
Xin Zheng ◽  
...  

In order to clarify the controversies for the role of nanoparticles on heat transfer in natural convection, lattice Boltzmann method is used to investigate Rayleigh-Benard convection heat transfer in differentially-heated enclosures filled with Al2O3-water nanofluids. The results for streamline and isotherm contours, vertical velocity, and temperature profiles as well as the local and average Nusselt number are discussed for a wide range of Rayleigh numbers and nanoparticle volume fractions (0 ? ? ? 5%). The results show that with the increase of Rayleigh number and nanoparticles loading, Nuave increases. It is suggested that the addition of nanoparticles can enhance the heat transfer in Rayleigh-Benard convection.


2021 ◽  
Vol 930 ◽  
Author(s):  
Hugo N. Ulloa ◽  
Juvenal A. Letelier

Thermally driven flows in fractures play a key role in enhancing the heat transfer and fluid mixing across the Earth's lithosphere. Yet the energy pathways in such confined environments have not been characterised. Building on Letelier et al. (J. Fluid Mech., vol. 864, 2019, pp. 746–767), we introduce novel expressions for energy transfer rates – energetics – of geometrically constrained Rayleigh–Bénard convection in Hele-Shaw cells (HS-RBC) based on two different conceptual frameworks. First, we derived the energetics following the well-established framework introduced by Winters et al. (J. Fluid Mech., vol. 289, 1995, pp. 115–128), in which the gravitational potential energy, $E_{\textit {p}}$ , is decomposed into its available, $E_{\textit {ap}}$ , and background, $E_{\textit {bp}}$ , components. Secondly, we derived the energetics considering a new decomposition for $E_{\textit {p}}$ , named dynamic, $E_{\textit {dp}}$ , and reference, $E_{\textit {rp}}$ , potential energies; $E_{\textit {dp}}$ is defined as the departure of the system's potential energy from the reference state $E_{\textit {rp}}$ , determined by the ‘energy’ of the scalar fluctuations. For HS-RBC, both frameworks lead to the same energy transfer rates at a steady state, satisfying the relationship $\langle E_{\textit {ap}} \rangle _{\tau } = \langle E_{\textit {dp}} \rangle _{\tau } + 1/6$ . Consistent with the work by Hughes et al. (J. Fluid Mech., vol. 729, 2013) on three-dimensional Rayleigh–Bénard convection, we report analytical expressions for the energetics and efficiencies of HS-RBC in terms of the Rayleigh number and the global Nusselt number. Additionally, we performed numerical experiments to illustrate the application of the energetics for the analysis of HS-RBC. Finally, we discuss the impact of the thermal forcing and the geometrical control exerted by Hele-Shaw cells on the development of boundary layers, protoplumes and the self-organisation of large-scale flows.


2015 ◽  
Vol 764 ◽  
pp. 349-361 ◽  
Author(s):  
Chong Shen Ng ◽  
Andrew Ooi ◽  
Detlef Lohse ◽  
Daniel Chung

AbstractResults from direct numerical simulations of vertical natural convection at Rayleigh numbers $1.0\times 10^{5}$–$1.0\times 10^{9}$ and Prandtl number $0.709$ support a generalised applicability of the Grossmann–Lohse (GL) theory, which was originally developed for horizontal natural (Rayleigh–Bénard) convection. In accordance with the GL theory, it is shown that the boundary-layer thicknesses of the velocity and temperature fields in vertical natural convection obey laminar-like Prandtl–Blasius–Pohlhausen scaling. Specifically, the normalised mean boundary-layer thicknesses scale with the $-1/2$-power of a wind-based Reynolds number, where the ‘wind’ of the GL theory is interpreted as the maximum mean velocity. Away from the walls, the dissipation of the turbulent fluctuations, which can be interpreted as the ‘bulk’ or ‘background’ dissipation of the GL theory, is found to obey the Kolmogorov–Obukhov–Corrsin scaling for fully developed turbulence. In contrast to Rayleigh–Bénard convection, the direction of gravity in vertical natural convection is parallel to the mean flow. The orientation of this flow presents an added challenge because there no longer exists an exact relation that links the normalised global dissipations to the Nusselt, Rayleigh and Prandtl numbers. Nevertheless, we show that the unclosed term, namely the global-averaged buoyancy flux that produces the kinetic energy, also exhibits both laminar and turbulent scaling behaviours, consistent with the GL theory. The present results suggest that, similar to Rayleigh–Bénard convection, a pure power-law relationship between the Nusselt, Rayleigh and Prandtl numbers is not the best description for vertical natural convection and existing empirical relationships should be recalibrated to better reflect the underlying physics.


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.


Sign in / Sign up

Export Citation Format

Share Document