scholarly journals ROBUST HETEROCLINIC CYCLES IN TWO-DIMENSIONAL RAYLEIGH–BÉNARD CONVECTION WITHOUT BOUSSINESQ SYMMETRY

2002 ◽  
Vol 12 (11) ◽  
pp. 2501-2522 ◽  
Author(s):  
ISABEL MERCADER ◽  
JOANA PRAT ◽  
EDGAR KNOBLOCH

The onset of convection in systems that are heated via current dissipation in the lower boundary or that lose heat from the top boundary via Newton's law of cooling is formulated as a bifurcation problem. The Rayleigh number as usually defined is shown to be inappropriate as a bifurcation parameter since the temperature difference across the layer depends on the amplitude of convection and hence changes as convection evolves at fixed external parameter values. A modified Rayleigh number is introduced that does remain constant even when the system is evolving, and solutions obtained with the standard formulation are compared with those obtained via the new one. Near the 1 : 2 spatial resonance in low Prandtl number fluids these effects open up intervals of Rayleigh number with no stable solutions in the form of steady convection or steadily traveling waves. Direct numerical simulations in two dimensions show that in such intervals the dynamics typically take the form of a nearly heteroclinic modulated traveling wave. This wave may be quasiperiodic or chaotic.

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.


2019 ◽  
Vol 881 ◽  
pp. 1073-1096 ◽  
Author(s):  
Andreas D. Demou ◽  
Dimokratis G. E. Grigoriadis

Rayleigh–Bénard convection in water is studied by means of direct numerical simulations, taking into account the variation of properties. The simulations considered a three-dimensional (3-D) cavity with a square cross-section and its two-dimensional (2-D) equivalent, covering a Rayleigh number range of $10^{6}\leqslant Ra\leqslant 10^{9}$ and using temperature differences up to 60 K. The main objectives of this study are (i) to investigate and report differences obtained by 2-D and 3-D simulations and (ii) to provide a first appreciation of the non-Oberbeck–Boussinesq (NOB) effects on the near-wall time-averaged and root-mean-squared (r.m.s.) temperature fields. The Nusselt number and the thermal boundary layer thickness exhibit the most pronounced differences when calculated in two dimensions and three dimensions, even though the $Ra$ scaling exponents are similar. These differences are closely related to the modification of the large-scale circulation pattern and become less pronounced when the NOB values are normalised with the respective Oberbeck–Boussinesq (OB) values. It is also demonstrated that NOB effects modify the near-wall temperature statistics, promoting the breaking of the top–bottom symmetry which characterises the OB approximation. The most prominent NOB effect in the near-wall region is the modification of the maximum r.m.s. values of temperature, which are found to increase at the top and decrease at the bottom of the cavity.


1988 ◽  
Vol 190 ◽  
pp. 451-469 ◽  
Author(s):  
D. R. Jenkins

The relationship between observations of cellular Rayleigh-Bénard convection using shadowgraphs and theoretical expressions for convection planforms is considered. We determine the shadowgraphs that ought to be observed if the convection is as given by theoretical expressions for roll, square or hexagonal planforms and compare them with actual experiments. Expressions for the planforms derived from linear theory, valid for low supercritical Rayleigh number, produce unambiguous shadowgraphs consisting of cells bounded by bright lines, which correspond to surfaces through which no fluid flows and on which the vertical component of velocity is directed downwards. Dark spots at the centre of cells, indicating regions of hot, rising fluid, are not accounted for by linear theory, but can be produced by adding higher-order terms, predominantly due to the temperature dependence of a material property of the fluid, such as its viscosity.


1981 ◽  
Vol 102 ◽  
pp. 85-100 ◽  
Author(s):  
D. E. Fitzjarrald

Convection flows have been systematically observed in a layer of fluid between two isothermal horizontal boundaries. The working fluid was a nematic liquid crystal, which exhibits a liquid–liquid phase change at which latent heat is released and the density changed. In addition to ordinary Rayleigh–Bénard convection when either phase is present alone, there exist two distinct types of convective motions initiated by the unstable density difference. When a thin layer of heavy fluid is present near the top boundary, hexagons with downgoing centres exist with no imposed thermal gradient. When a thin layer of light fluid is brought on near the lower boundary, the hexagons have upshooting centres. In both cases, the motions are kept going once they are initiated by the instability due to release of latent heat. Relation of the results to applicable theories is discussed.


1994 ◽  
Vol 272 ◽  
pp. 67-90 ◽  
Author(s):  
Michael D. Graham ◽  
Paul H. Steen

The classical boundary-layer scaling laws proposed by Howard for Rayleigh–Bénard convection at high Rayleigh number extend to the analogous case of convection in saturated porous media. We computationally study two-dimensional porous-media convection near the onset of this scaling behaviour. The main result of the paper is the observation and study of instabilities that lead to deviations from the scaling relations.At Rayleigh numbers below the scaling regime, boundary-layer fluctuations born at a Hopf bifurcation strengthen and eventually develop into thermal plumes. The appearance of plumes corresponds to the onset of the boundary-layer scaling behaviour of the oscillation frequency and mean Nusselt number, in agreement with the classical theory. As the Rayleigh number increases further, the flow undergoes instabilities that lead to ‘bubbles’ in parameter space of quasi-periodic flow, and eventually to weakly chaotic flow. The instabilities disturb the plume formation process, effectively leading to a phase modulation of the process and to deviations from the scaling laws. We argue that these instabilities correspond to parametric resonances between the timescale for plume formation and the characteristic convection timescale of the flow.


2016 ◽  
Vol 804 ◽  
pp. 370-386 ◽  
Author(s):  
David Goluskin ◽  
Charles R. Doering

We consider Rayleigh–Bénard convection in a layer of fluid between rough no-slip boundaries where the top and bottom boundary heights are functions of the horizontal coordinates with square-integrable gradients. We use the background method to derive an upper bound on the mean heat flux across the layer for all admissible boundary geometries. This flux, normalized by the temperature difference between the boundaries, can grow with the Rayleigh number ($Ra$) no faster than $O(Ra^{1/2})$ as $Ra\rightarrow \infty$. Our analysis yields a family of similar bounds, depending on how various estimates are tuned, but every version depends explicitly on the boundary geometry. In one version the coefficient of the $O(Ra^{1/2})$ leading term is $0.242+2.925\Vert \unicode[STIX]{x1D735}h\Vert ^{2}$, where $\Vert \unicode[STIX]{x1D735}h\Vert ^{2}$ is the mean squared magnitude of the boundary height gradients. Application to a particular geometry is illustrated for sinusoidal boundaries.


Entropy ◽  
2020 ◽  
Vol 22 (3) ◽  
pp. 283
Author(s):  
Laiyun Zheng ◽  
Bingxin Zhao ◽  
Jianqing Yang ◽  
Zhenfu Tian ◽  
Ming Ye

This paper studied the Rayleigh–Bénard convection in binary fluid mixtures with a strong Soret effect (separation ratio ψ = − 0.6 ) in a rectangular container heated uniformly from below. We used a high-accuracy compact finite difference method to solve the hydrodynamic equations used to describe the Rayleigh–Bénard convection. A stable traveling-wave convective state with periodic source defects (PSD-TW) is obtained and its properties are discussed in detail. Our numerical results show that the novel PSD-TW state is maintained by the Eckhaus instability and the difference between the creation and annihilation frequencies of convective rolls at the left and right boundaries of the container. In the range of Rayleigh number in which the PSD-TW state is stable, the period of defect occurrence increases first and then decreases with increasing Rayleigh number. At the upper bound of this range, the system transitions from PSD-TW state to another type of traveling-wave state with aperiodic and more dislocated defects. Moreover, we consider the problem with the Prandtl number P r ranging from 0.1 to 20 and the Lewis number L e from 0.001 to 1, and discuss the stabilities of the PSD-TW states and present the results as phase diagrams.


2019 ◽  
Vol 141 (6) ◽  
Author(s):  
P. G. Siddheshwar ◽  
K. M. Lakshmi

Unicellular Rayleigh–Bénard convection of water–copper nanoliquid confined in a high-porosity enclosure is studied analytically. The modified-Buongiorno–Brinkman two-phase model is used for nanoliquid description to include the effects of Brownian motion, thermophoresis, porous medium friction, and thermophysical properties. Free–free and rigid–rigid boundaries are considered for investigation of onset of convection and heat transport. Boundary effects on onset of convection are shown to be classical in nature. Stability boundaries in the R1*–R2 plane are drawn to specify the regions in which various instabilities appear. Specifically, subcritical instabilities' region of appearance is highlighted. Square, shallow, and tall porous enclosures are considered for study, and it is found that the maximum heat transport occurs in the case of a tall enclosure and minimum in the case of a shallow enclosure. The analysis also reveals that the addition of a dilute concentration of nanoparticles in a liquid-saturated porous enclosure advances onset and thereby enhances the heat transport irrespective of the type of boundaries. The presence of porous medium serves the purpose of heat storage in the system because of its low thermal conductivity.


Sign in / Sign up

Export Citation Format

Share Document