Direct Numerical Simulation of a Low Prandtl Number Rayleigh–Bénard Convection in a Square Box

2019 ◽  
Vol 11 (6) ◽  
Author(s):  
Ojas Satbhai ◽  
Subhransu Roy ◽  
Sudipto Ghosh
Keyword(s):  
Nusselt Number ◽  
Steady State ◽  
Rayleigh Number ◽  
Prandtl Number ◽  
Dimensionless Temperature ◽  
Bénard Convection ◽  
Benard Convection ◽  
Wide Range ◽  
Rayleigh Benard Convection ◽  
Rayleigh Benard

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.

Exploring the severely confined regime in Rayleigh–Bénard convection

2016 ◽  
Vol 805 ◽  
Author(s):  
Kai Leong Chong ◽  
Ke-Qing Xia
Keyword(s):  
Rayleigh Number ◽  
Flow Topology ◽  
Thermal Plumes ◽  
Bénard Convection ◽  
Benard Convection ◽  
Wide Range ◽  
Classical Boundary ◽  
Global And Local ◽  
Rayleigh Benard Convection ◽  
Rayleigh Benard

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

The effect of distant side-walls on the evolution and stability of finite-amplitude Rayleigh-Bénard convection

1981 ◽  
Vol 378 (1775) ◽  
pp. 539-566 ◽  
Keyword(s):  
Rayleigh Number ◽  
Prandtl Number ◽  
Finite Amplitude ◽  
Two Dimensional ◽  
Bénard Convection ◽  
Benard Convection ◽  
Side Walls ◽  
Rayleigh Benard Convection ◽  
Rayleigh Numbers ◽  
Rayleigh Benard

A recent study by Cross et al . (1980) has described a class of finite-amplitude phase-winding solutions of the problem of two-dimensional Rayleigh-Bénard convection in a shallow fluid layer of aspect ratio 2 L (≫ 1) confined laterally by rigid side-walls. These solutions arise at Rayleigh numbers R = R 0 + O ( L -1 ) where R 0 is the critical Rayleigh number for the corresponding infinite layer. Nonlinear solutions of constant phase exist for Rayleigh numbers R = R 0 + O ( L -2 ) but of these only the two that bifurcate at the lowest value of R are stable to two-dimensional linearized disturbances in this range (Daniels 1978). In the present paper one set of the class of phase-winding solutions is found to be stable to two-dimensional disturbances. For certain values of the Prandtl number of the fluid and for stress-free horizontal boundaries the results predict that to preserve stability there must be a continual readjustment of the roll pattern as the Rayleigh number is raised, with a corresponding increase in wavelength proportional to R - R 0 . These solutions also exhibit hysteresis as the Rayleigh number is raised and lowered. For other values of the Prandtl number the number of rolls remains unchanged as the Rayleigh number is raised, and the wavelength remains close to its critical value. It is proposed that the complete evolution of the flow pattern from a static state must take place on a number of different time scales of which t = O(( R - R 0 ) -1 ) and t = O(( R - R 0 ) -2 ) are the most significant. When t = O(( R - R 0 ) -1 ) the amplitude of convection rises from zero to its steady-state value, but the final lateral positioning of the rolls is only completed on the much longer time scale t = O(( R - R 0 ) -2 ).

Turbulent Rayleigh–Bénard convection in a cylindrical container with aspect ratio Γ = 0.50 and Prandtl number Pr = 4.38

2011 ◽  
Vol 676 ◽  
pp. 5-40 ◽  
Author(s):  
STEPHAN WEISS ◽  
GUENTER AHLERS
Keyword(s):  
Nusselt Number ◽  
Aspect Ratio ◽  
Prandtl Number ◽  
State Transitions ◽  
Small Scale ◽  
Bénard Convection ◽  
Benard Convection ◽  
Run Time ◽  
Rayleigh Benard Convection ◽  
Rayleigh Benard

Measurements of the Nusselt number and properties of the large-scale circulation (LSC) are presented for turbulent Rayleigh–Bénard convection in water-filled cylindrical containers (Prandtl number Pr = 4.38) with aspect ratio Γ = 0.50. They cover the range 2 × 108 ≲ Ra ≲ 1 × 1011 of the Rayleigh number Ra. We confirm the occurrence of a double-roll state (DRS) of the LSC and focus on the statistics of the transitions between the DRS and a single-roll state (SRS). The fraction of the run time when the SRS existed varied continuously from about 0.12 near Ra = 2 × 108 to about 0.8 near Ra = 1011, while the fraction of the run time when the DRS could be detected changed from about 0.4 to about 0.06 over the same range of Ra. We determined separately the Nusselt number of the SRS and the DRS, and found the former to be larger than the latter by about 1.6% (0.9%) at Ra = 1010 (1011). We report a contribution to the dynamics of the SRS from a torsional oscillation similar to that observed for cylindrical samples with Γ = 1.00. Results for a number of statistical properties of the SRS are reported, and some are compared with the cases Γ = 0.50, Pr = 0.67 and Γ = 1.00, Pr = 4.38. We found that genuine cessations of the SRS were extremely rare and occurred only about 0.3 times per day, which is less frequent than for Γ = 1.00; however, the SRS was disrupted frequently by roll-state transitions and other less well-defined events. We show that the time derivative of the LSC plane orientation is a stochastic variable which, at constant LSC amplitude, is Gaussian distributed. Within the context of the LSC model of Brown & Ahlers (Phys. Fluids, vol. 20, 2008b, art. 075101), this demonstrates that the stochastic force due to the small-scale fluctuations that is driving the LSC dynamics has a Gaussian distribution.

Has the ultimate state of turbulent thermal convection been observed?

2015 ◽  
Vol 785 ◽  
pp. 270-282 ◽  
Author(s):  
L. Skrbek ◽  
P. Urban
Keyword(s):  
Nusselt Number ◽  
Rayleigh Number ◽  
Ultimate State ◽  
Logarithmic Factor ◽  
Aspect Ratios ◽  
Bénard Convection ◽  
Benard Convection ◽  
Rayleigh Benard Convection ◽  
Rayleigh Numbers ◽  
Rayleigh Benard

An important question in turbulent Rayleigh–Bénard convection is the scaling of the Nusselt number with the Rayleigh number in the so-called ultimate state, corresponding to asymptotically high Rayleigh numbers. A related but separate question is whether the measurements support the so-called Kraichnan law, according to which the Nusselt number varies as the square root of the Rayleigh number (modulo a logarithmic factor). Although there have been claims that the Kraichnan regime has been observed in laboratory experiments with low aspect ratios, the totality of existing experimental results presents a conflicting picture in the high-Rayleigh-number regime. We analyse the experimental data to show that the claims on the ultimate state leave open an important consideration relating to non-Oberbeck–Boussinesq effects. Thus, the nature of scaling in the ultimate state of Rayleigh–Bénard convection remains open.

Rayleigh–Bénard convection in the presence of spatial temperature modulations

2011 ◽  
Vol 673 ◽  
pp. 318-348 ◽  
Author(s):  
G. FREUND ◽  
W. PESCH ◽  
W. ZIMMERMANN
Keyword(s):  
Rayleigh Number ◽  
Numerical Solutions ◽  
Critical Rayleigh Number ◽  
Lower Boundary ◽  
Spatial Modulation ◽  
Bénard Convection ◽  
Benard Convection ◽  
Wide Range ◽  
Rayleigh Benard Convection ◽  
Rayleigh Benard

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.

Strongly turbulent Rayleigh–Bénard convection in mercury: comparison with results at moderate Prandtl number

1997 ◽  
Vol 335 ◽  
pp. 111-140 ◽  
Author(s):  
S. CIONI ◽  
S. CILIBERTO ◽  
J. SOMMERIA
Keyword(s):  
Nusselt Number ◽  
Prandtl Number ◽  
Scaling Laws ◽  
Convective Cell ◽  
Zone Model ◽  
Turbulent Regime ◽  
Bénard Convection ◽  
Benard Convection ◽  
Rayleigh Benard Convection ◽  
Rayleigh Benard

An experimental study of Rayleigh–Bénard convection in the strongly turbulent regime is presented. We report results obtained at low Prandtl number (in mercury, Pr = 0.025), covering a range of Rayleigh numbers 5 × 106 < Ra < 5 × 109, and compare them with results at Pr∼1. The convective chamber consists of a cylindrical cell of aspect ratio 1.Heat flux measurements indicate a regime with Nusselt number increasing as Ra0.26, close to the 2/7 power observed at Pr∼1, but with a smaller prefactor, which contradicts recent theoretical predictions. A transition to a new turbulent regime is suggested for Ra ≃ 2 × 109, with significant increase of the Nusselt number. The formation of a large convective cell in the bulk is revealed by its thermal signature on the bottom and top plates. One frequency of the temperature oscillation is related to the velocity of this convective cell. We then obtain the typical temperature and velocity in the bulk versus the Rayleigh number, and compare them with similar results known for Pr∼1.We review two recent theoretical models, namely the mixing zone model of Castaing et al. (1989), and a model of the turbulent boundary layer by Shraiman & Siggia (1990). We discuss how these models fail at low Prandtl number, and propose modifications for this case. Specific scaling laws for fluids at low Prandtl number are then obtained, providing an interpretation of our experimental results in mercury, as well as extrapolations for other liquid metals.

Bounds on Rayleigh–Bénard convection with an imposed heat flux

2002 ◽  
Vol 473 ◽  
pp. 191-199 ◽  
Author(s):  
JESSE OTERO ◽  
RALF W. WITTENBERG ◽  
RODNEY A. WORTHING ◽  
CHARLES R. DOERING
Keyword(s):  
Heat Flux ◽  
Nusselt Number ◽  
Rayleigh Number ◽  
Heat Transport ◽  
Absolute Constant ◽  
Temperature Drop ◽  
Bénard Convection ◽  
Benard Convection ◽  
Rayleigh Benard Convection ◽  
Rayleigh Benard

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.

scholarly journals Bounds for Rayleigh–Bénard convection between free-slip boundaries with an imposed heat flux

2018 ◽  
Vol 837 ◽  
Author(s):  
Giovanni Fantuzzi
Keyword(s):  
Heat Transfer ◽  
Heat Flux ◽  
Rayleigh Number ◽  
Prandtl Number ◽  
Temperature Drop ◽  
Three Dimensional ◽  
Bénard Convection ◽  
Benard Convection ◽  
Rayleigh Benard Convection ◽  
Rayleigh Benard

We prove the first rigorous bound on the heat transfer for three-dimensional Rayleigh–Bénard convection of finite-Prandtl-number fluids between free-slip boundaries with an imposed heat flux. Using the auxiliary functional method with a quadratic functional, which is equivalent to the background method, we prove that the Nusselt number $\mathit{Nu}$ is bounded by $\mathit{Nu}\leqslant 0.5999\mathit{R}^{1/3}$ uniformly in the Prandtl number, where $\mathit{R}$ is the Rayleigh number based on the imposed heat flux. In terms of the Rayleigh number based on the mean vertical temperature drop, $\mathit{Ra}$, we obtain $\mathit{Nu}\leqslant 0.4646\mathit{Ra}^{1/2}$. The scaling with Rayleigh number is the same as that of bounds obtained with no-slip isothermal, free-slip isothermal and no-slip fixed-flux boundaries, and numerical optimisation of the bound suggests that it cannot be improved within our bounding framework. Contrary to the two-dimensional case, therefore, the $\mathit{Ra}$-dependence of rigorous upper bounds on the heat transfer obtained with the background method for three-dimensional Rayleigh–Bénard convection is insensitive to both the thermal and the velocity boundary conditions.

Sign in / Sign up

Export Citation Format

Share Document