scholarly journals Effects of modulation on Rayleigh-Benard convection. Part I

2004 ◽  
Vol 2004 (19) ◽  
pp. 991-1001 ◽  
Author(s):  
B. S. Bhadauria ◽  
Lokenath Debnath

The linear stability of a horizontal layer of fluid heated from below and above is considered. In addition to a steady temperature difference between the walls of the fluid layer, a time-dependent periodic perturbation is applied to the wall temperatures. Only infinitesimal disturbances are considered. Numerical results for the critical Rayleigh number are obtained at various Prandtl numbers and for various values of the frequency. Some comparisons have been made with the known results.

2003 ◽  
Vol 58 (2-3) ◽  
pp. 176-182
Author(s):  
B. S. Bhadauria

The linear stability of a horizontal fluid layer, confined between two rigid walls, heated from below and cooled from above is considered. The temperature gradient between the walls consists of a steady part and a periodic part that oscillates with time. Only infinitesimal disturbances are considered. Numerical results for the critical Rayleigh number are obtained for various Prandtl numbers and for various values of the frequency. Some comparisons with known results have also been made.


2004 ◽  
Vol 59 (4-5) ◽  
pp. 266-274
Author(s):  
B. S. Bhadauria

The linear thermal instability of a horizontal fluid layer with time-periodic temperature distribution is studied with the help of the Floquet theory. The time-dependent part of the temperature has been expressed in Fourier series. Disturbances are assumed to be infinitesimal. Only even solutions are considered. Numerical results for the critical Rayleigh number are obtained at various Prandtl numbers and for various values of the frequency. It is found that the disturbances are either synchronous with the primary temperature field or have half its frequency. - 2000 Mathematics Subject Classification: 76E06, 76R10.


2021 ◽  
Vol 2057 (1) ◽  
pp. 012012
Author(s):  
A I Fedyushkin

Abstract The paper presents the results of a numerical study of convective heat transfer in a long horizontal layer heated from below with and without the vibration effect of the lower wall. The simulation was carried out on the basis of solving the Navier-Stokes 2D equations for an incompressible fluid in the Boussinesq approximation. It is shown that the influence of vibrations of the lower heated wall on the wave number of the convective flow roll structure, on the time and on the critical Rayleigh number of convection. The influence of controlled harmonic vibrations of wall on the structure of convective flow in the Rayleigh-Benard problem has been investigated. It is shown that the wave number of the periodic convective structure, the critical Rayleigh number, and the time of occurrence of Rayleigh-Benard convection under the vertical vibration effect on the horizontal layer from the lower wall are reduced.


1984 ◽  
Vol 143 ◽  
pp. 125-152 ◽  
Author(s):  
P. G. Daniels

This paper considers the temporal evolution of two-dimensional Rayleigh–Bénard convection in a shallow fluid layer of aspect ratio 2L ([Gt ] 1) confined laterally by rigid sidewalls. Recent studies by Cross et al. (1980, 1983) have shown that for Rayleigh numbers in the range R = R0 + O(L−1) (where R0 is the critical Rayleigh number for the corresponding infinite layer) there exists a class of finite-amplitude steady-state ‘phase-winding’ solutions which correspond physically to the possibility of an adjustment in the number of rolls in the container as the local value of the Rayleigh number is varied. It has been shown (Daniels 1981) that in the temporal evolution of the system the final lateral positioning of the rolls occurs on the long timescale t = O(L2) when the phase function which determines the number of rolls in the system satisfies a one-dimensional diffusion equation but with novel boundary conditions that represent the effect of the sidewalls. In the present paper this system is solved numerically in order to determine the precise way in which the roll pattern adjusts after a change in the Rayleigh number of the system. There is an interesting balance between, on the one hand, a tendency for the number of rolls to change by the least number possible and, on the other, a tendency for the even or odd nature of the initial configuration to be preserved during the transition. In some cases this second property renders the natural evolution susceptible to arbitrarily small external disturbances, which then dictate the form of the final roll pattern.The complete transition involves an analysis of the motion on three timescales, a conductive scale t = O(1), a convective growth scale t = O(L) and a convective diffusion scale t = O(L2).


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.


1996 ◽  
Vol 326 ◽  
pp. 399-415 ◽  
Author(s):  
M. Wanschura ◽  
H. C. Kuhlmann ◽  
H. J. Rath

The stability of steady axisymmetric convection in cylinders heated from below and insulated laterally is investigated numerically using a mixed finite-difference/Chebyshev collocation method to solve the base flow and the linear stability equations. Linear stability boundaries are given for radius to height ratios γ from 0.9 to 1.56 and for Prandtl numbers Pr = 0.02 and Pr = 1. Depending on γ and Pr, the azimuthal wavenumber of the critical mode may be m = 1, 2, 3, or 4. The dependence of the critical Rayleigh number on the aspect ratio and the instability mechanisms are explained by analysing the energy transfer to the critical modes for selected cases. In addition to these results the onset of buoyant convection in liquid bridges with stress-free conditions on the cylindrical surface is considered. For insulating thermal boundary conditions, the onset of convection is never axisymmetric and the critical azimuthal wavenumber increases monotonically with γ. The critical Rayleigh number is less then 1708 for most aspect ratios.


1970 ◽  
Vol 42 (1) ◽  
pp. 161-175 ◽  
Author(s):  
D. R. Caldwell

Observations of temperature drop as a function of heat flow in Rayleigh–Bénard convection with curved density profiles show: (1) reversal of slope in the heating curve, (2) oscillations with time, (3) history dependence, and (4) an increase in critical Rayleigh number as the curvature of the density profile is increased. Some of the results are quite similar to the predictions of Busse.


2013 ◽  
Vol 736 ◽  
pp. 464-494 ◽  
Author(s):  
P. Pearce ◽  
J. Daou

AbstractWe investigate the Rayleigh–Bénard convection problem within the context of a diffusion flame formed in a horizontal channel where the fuel and oxidizer concentrations are prescribed at the porous walls. This problem seems to have received no attention in the literature. When formulated in the low-Mach-number approximation the model depends on two main non-dimensional parameters, the Rayleigh number and the Damköhler number, which govern gravitational strength and reaction speed respectively. In the steady state the system admits a planar diffusion flame solution; the aim is to find the critical Rayleigh number at which this solution becomes unstable to infinitesimal perturbations. In the Boussinesq approximation, a linear stability analysis reduces the system to a matrix equation with a solution comparable to that of the well-studied non-reactive case of Rayleigh–Bénard convection with a hot lower boundary. The planar Burke–Schumann diffusion flame, which has been previously considered unconditionally stable in studies disregarding gravity, is shown to become unstable when the Rayleigh number exceeds a critical value. A numerical treatment is performed to test the effects of compressibility and finite chemistry on the stability of the system. For weak values of the thermal expansion coefficient $\alpha $, the numerical results show strong agreement with those of the linear stability analysis. It is found that as $\alpha $ increases to a more realistic value the system becomes considerably more stable, and also exhibits hysteresis at the onset of instability. Finally, a reduction in the Damköhler number is found to decrease the stability of the system.


2003 ◽  
pp. 29-40 ◽  
Author(s):  
Jürgen Zierep

We discuss the solution of the small perturbation equations for a horizontal fluid layer heated from below with an applied magnetic field either in vertical or in horizontal direction. The magnetic field stabilizes, due to the Lorentz force, more or less Rayleigh-B?nard convective cellular motion. The solution of the eigenvalue problem shows that the critical Rayleigh number increases with increasing Hartmann number while the corresponding wave length decreases. Interesting analogies to solar granulation and black spots phenomena are obvious. The influence of a horizontal field is stronger than that of a vertical field. It is easy to understand this by discussing the influence of the Lorentz force on the Rayleigh-B?nard convection. This result corrects earlier calculations in the literature.


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