Nonlinear Rayleigh–Bénard Convection With Variable Heat Source

2013 ◽  
Vol 135 (12) ◽  
Author(s):  
P. G. Siddheshwar ◽  
P. Stephen Titus

Linear and nonlinear Rayleigh–Bénard convections with variable heat source (sink) are studied analytically using the Fourier series. The strength of the heat source is characterized by an internal Rayleigh number, RI, whose effect is to decrease the critical external Rayleigh number. Linear theory involving an autonomous system (linearized Lorenz model) further reveals that the critical point at pre-onset can only be a saddle point. In the postonset nonlinear study, analysis of the generalized Lorenz model leads us to two other critical points that take over from the critical point of the pre-onset regime. Classical analysis of the Lorenz model points to the possibility of chaos. The effect of RI is shown to delay or advance the appearance of chaos depending on whether RI is negative or positive. This aspect is also reflected in its effect on the Nusselt number. The Lyapunov exponents provide useful information on the closing in and opening out of the trajectories of the solution of the Lorenz model in the cases of heat sink and heat source, respectively. The Ginzburg-Landau models for the problem are obtained via the 3-mode and 5-mode Lorenz models of the paper.

2021 ◽  
Vol 88 (1-2) ◽  
pp. 08
Author(s):  
A. S. Aruna ◽  
V. Ramachandramurthy ◽  
N. Kavitha

The present paper aims at weak non-linear stability analysis followed by linear analysis of nite-amplitude Rayleigh-Benard magneto convection problem in an electrically conducting Newtonian liquid with heat source/sink. It is shown that the internal Rayleigh number, ther- morheological parameter, and the Chandrasekhar number in uence the onset of convection. The generalized Lorenz model derived for the prob- lem is essentially the classical Lorenz model but with some coecient depending on the variable heat source (sink), viscosity, and the applied magnetic eld. The result of the parameters' in uence on the critical Rayleigh number explains their in uence on the Nusselt number. It is found that an increasing strength of the magnetic eld is to stabilize the system and diminishes heat transport whereas the heat source and variable viscosity in-tandem to work system unstable and enhances heat transfer.


2020 ◽  
Vol 16 (5) ◽  
pp. 1175-1188
Author(s):  
Keerthi R ◽  
B. Mahanthesh ◽  
Smita Saklesh Nagouda

PurposeThe study of instability due to the effects of Maxwell–Cattaneo law and internal heat source/sink on Casson dielectric fluid horizontal layer is an open question. Therefore, in this paper, the impact of internal heat generation/absorption on Rayleigh–Bénard convection in a non-Newtonian dielectric fluid with Maxwell–Cattaneo heat flux is investigated. The horizontal layer of the fluid is cooled from the upper boundary, while an isothermal boundary condition is utilized at the lower boundary.Design/methodology/approachThe Casson fluid model is utilized to characterize the non-Newtonian fluid behavior. The horizontal layer of the fluid is cooled from the upper boundary, while an isothermal boundary condition is utilized at the lower boundary. The governing equations are non-dimensionalized using appropriate dimensionless variables and the subsequent equations are solved for the critical Rayleigh number using the normal mode technique (NMT).FindingsResults are presented for two different cases namely dielectric Newtonian fluid (DNF) and dielectric non-Newtonian Casson fluid (DNCF). The effects of Cattaneo number, Casson fluid parameter, heat source/sink parameter on critical Rayleigh number and wavenumber are analyzed in detail. It is found that the value Rayleigh number for non-Newtonian fluid is higher than that of Newtonian fluid; also the heat source aspect decreases the magnitude of the Rayleigh number.Originality/valueThe effect of Maxwell–Cattaneo heat flux and internal heat source/sink on Rayleigh-Bénard convection in Casson dielectric fluid is investigated for the first time.


1997 ◽  
Vol 55 (4) ◽  
pp. R3824-R3827 ◽  
Author(s):  
Omar Osenda ◽  
Carlos B. Briozzo ◽  
Manuel O. Cáceres

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.


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.


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.


Sign in / Sign up

Export Citation Format

Share Document