scholarly journals Rayleigh–Bénard stability and the validity of quasi-Boussinesq or quasi-anelastic liquid approximations

2017 ◽  
Vol 817 ◽  
pp. 264-305 ◽  
Author(s):  
Thierry Alboussière ◽  
Yanick Ricard
Keyword(s):  
Stability Analysis ◽  
Equation Of State ◽  
Rayleigh Number ◽  
Equations Of State ◽  
Critical Rayleigh Number ◽  
Ideal Gas ◽  
Dimensionless Temperature ◽  
Boussinesq Model ◽  
Onset Of Convection ◽  
Rayleigh Benard

The linear stability threshold of the Rayleigh–Bénard configuration is analysed with compressible effects taken into account. It is assumed that the fluid under investigation obeys a Newtonian rheology and Fourier’s law of thermal transport with constant, uniform (dynamic) viscosity and thermal conductivity in a uniform gravity field. Top and bottom boundaries are maintained at different constant temperatures and we consider here mechanical boundary conditions of zero tangential stress and impermeable walls. Under these conditions, and with the Boussinesq approximation, Rayleigh (Phil. Mag., vol. 32 (192), 1916, pp. 529–546) first obtained analytically the critical value $27\unicode[STIX]{x03C0}^{4}/4$ for a dimensionless parameter, now known as the Rayleigh number, at the onset of convection. This paper describes the changes of the critical Rayleigh number due to the compressibility of the fluid, measured by the dimensionless dissipation parameter ${\mathcal{D}}$ and due to a finite temperature difference between the hot and cold boundaries, measured by a dimensionless temperature gradient $a$. Different equations of state are examined: ideal gas equation, Murnaghan’s model (often used to describe the interiors of solid but convective planets) and a generic equation of state with adjustable parameters, which can represent any possible equation of state. In the perspective to assess approximations often made in convective models, we also consider two variations of this stability analysis. In a so-called quasi-Boussinesq model, we consider that density perturbations are solely due to temperature perturbations. In a so-called quasi-anelastic liquid approximation model, we consider that entropy perturbations are solely due to temperature perturbations. In addition to the numerical Chebyshev-based stability analysis, an analytical approximation is obtained when temperature fluctuations are written as a combination of only two modes, one being the original symmetrical (between top and bottom) mode introduced by Rayleigh, the other one being antisymmetrical. The analytical solution allows us to show that the antisymmetrical part of the critical eigenmode increases linearly with the parameters $a$ and ${\mathcal{D}}$, while the superadiabatic critical Rayleigh number departs quadratically in $a$ and ${\mathcal{D}}$ from $27\unicode[STIX]{x03C0}^{4}/4$. For any arbitrary equation of state, the coefficients of the quadratic departure are determined analytically from the coefficients of the expansion of density up to degree three in terms of pressure and temperature.

Non-linear Rayleigh-Benard Magnetoconvection in Temperature-sensitive Newtonian Liquids with Variable Heat Source

2021 ◽  
Vol 88 (1-2) ◽  
pp. 08
Author(s):  
A. S. Aruna ◽  
V. Ramachandramurthy ◽  
N. Kavitha
Keyword(s):  
Rayleigh Number ◽  
Heat Source ◽  
Variable Viscosity ◽  
Critical Rayleigh Number ◽  
Lorenz Model ◽  
Onset Of Convection ◽  
Variable Heat ◽  
Non Linear ◽  
Source Sink ◽  
Rayleigh Benard

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.

An analysis of convection in a mushy layer with a deformable permeable interface

2008 ◽  
Vol 596 ◽  
pp. 333-352 ◽  
Author(s):  
S. M. ROPER ◽  
S. H. DAVIS ◽  
P. W. VOORHEES
Keyword(s):  
Stability Analysis ◽  
Rayleigh Number ◽  
Free Boundary ◽  
Critical Rayleigh Number ◽  
Liquid Interface ◽  
Mushy Layer ◽  
Onset Of Convection ◽  
Weakly Nonlinear ◽  
Nonlinear Growth ◽  
Weakly Nonlinear Theory

We study the dynamics of a mushy layer in directional solidification for the case of a thin near-eutectic mush with a deformable and permeable mush–liquid interface. We examine the onset of convection using linear stability analysis, and the weakly nonlinear growth of liquid inclusions that signal the onset of chimneys. This analysis is compared to past analyses in which the mush–liquid interface is replaced by a rigid impermeable lid. We find qualitative agreement between the two models, but the rigid-lid approximation gives substantially different quantitative behaviour.In linear theory, the rigid-lid approximation leads to an over-estimate of the critical Rayleigh number and wavenumber of the instability. The condition for the onset of oscillatory instability is also changed by a factor of about 5 in composition number C. In the weakly nonlinear theory, the location of the onset of liquid inclusions is near the undisturbed front for the free-boundary analysis, whereas it lies at the centre of the mushy layer when the rigid-lid approximation is used. For hexagonal patterns, the boundary between regions of parameter space in which up and down hexagons are stable, shifts as a result of coupling between the liquid and mush regions.

On the Effect of Mechanical and Thermal Anisotropy on the Stability of Gravity Driven Convection in Rotating Porous Media

2005 ◽  
Author(s):  
Saneshan Govender ◽  
Peter Vadasz
Keyword(s):  
Rayleigh Number ◽  
Critical Rayleigh Number ◽  
Equation Model ◽  
Linear Stability Theory ◽  
Mechanical Anisotropy ◽  
Onset Of Convection ◽  
Thermal Anisotropy ◽  
Gravity Driven ◽  
The Stability ◽  
Rayleigh Benard

We investigate Rayleigh-Benard convection in a porous layer subjected to gravitational and Coriolis body forces, when the fluid and solid phases are not in local thermodynamic equilibrium. The Darcy model (extended to include Coriolis effects and anisotropic permeability) is used to describe the flow whilst the two-equation model is used for the energy equation (for the solid and fluid phases separately). The linear stability theory is used to evaluate the critical Rayleigh number for the onset of convection and the effect of both thermal and mechanical anisotropy on the critical Rayleigh number is discussed.

Heat transport by parallel-roll convection in a rectangular container

1987 ◽  
Vol 185 ◽  
pp. 205-234 ◽  
Author(s):  
R. W. Walden ◽  
Paul Kolodner ◽  
A. Passner ◽  
C. M. Surko
Keyword(s):  
Rayleigh Number ◽  
Heat Transport ◽  
Critical Rayleigh Number ◽  
Equation Model ◽  
Onset Of Convection ◽  
Infinite Layer ◽  
Lateral Extent ◽  
Prandtl Numbers ◽  
Rayleigh Numbers ◽  
Good Agreement

Heat-transport measurements are reported for thermal convection in a rectangular box of aspect’ ratio 10 x 5. Results are presented for Rayleigh numbers up to 35Rc, Prandtl numbers between 2 and 20, and wavenumbers between 0.6 and 1.0kc, where Rc and kc are the critical Rayleigh number and wavenumber for the onset of convection in a layer of infinite lateral extent. The measurements are in good agreement with a phenomenological model which combines the calculations of Nusselt number, as a function of Rayleigh number and roll wavenumber for two-dimensional convection in an infinite layer, with a nonlinear amplitude-equation model developed to account for sidewell attenuation. The appearance of bimodal convection increases the heat transport above that expected for simple parallel-roll convection.

Stability and Convection in Impulsively Heated Porous Layers

2008 ◽  
Vol 130 (11) ◽  
Author(s):  
M. J. Kohl ◽  
M. Kristoffersen ◽  
F. A. Kulacki
Keyword(s):  
Porous Medium ◽  
Rayleigh Number ◽  
Heated Surface ◽  
Critical Rayleigh Number ◽  
Axisymmetric Flow ◽  
Temperature Fluctuations ◽  
Lower Surface ◽  
Onset Of Convection ◽  
Upward Moving ◽  
The Stability

Experiments are reported on initial instability, turbulence, and overall heat transfer in a porous medium heated from below. The porous medium comprises either water or a water-glycerin solution and randomly stacked glass spheres in an insulated cylinder of height:diameter ratio of 1.9. Heating is with a constant flux lower surface and a constant temperature upper surface, and the stability criterion is determined for a step heat input. The critical Rayleigh number for the onset of convection is obtained in terms of a length scale normalized to the thermal penetration depth as Rac=83/(1.08η−0.08η2) for 0.02<η<0.18. Steady convection in terms of the Nusselt and Rayleigh numbers is Nu=0.047Ra0.91Pr0.11(μ/μ0)0.72 for 100<Ra<5000. Time-averaged temperatures suggest the existence of a unicellular axisymmetric flow dominated by upflow over the central region of the heated surface. When turbulence is present, the magnitude and frequency of temperature fluctuations increase weakly with increasing Rayleigh number. Analysis of temperature fluctuations in the fluid provides an estimate of the speed of the upward moving thermals, which decreases with distance from the heated surface.

scholarly journals Linear biglobal analysis of Rayleigh–Bénard instabilities in binary fluids with and without throughflow

2012 ◽  
Vol 713 ◽  
pp. 216-242 ◽  
Author(s):  
Jun Hu ◽  
Daniel Henry ◽  
Xie-Yuan Yin ◽  
Hamda BenHadid
Keyword(s):  
Reynolds Number ◽  
Rayleigh Number ◽  
Aspect Ratio ◽  
Critical Rayleigh Number ◽  
Two Dimensional ◽  
Binary Fluids ◽  
Transverse Modes ◽  
Aspect Ratios ◽  
Rayleigh Benard ◽  
Mode A

AbstractThree-dimensional Rayleigh–Bénard instabilities in binary fluids with Soret effect are studied by linear biglobal stability analysis. The fluid is confined transversally in a duct and a longitudinal throughflow may exist or not. A negative separation factor $\psi = \ensuremath{-} 0. 01$, giving rise to oscillatory transitions, has been considered. The numerical dispersion relation associated with this stability problem is obtained with a two-dimensional Chebyshev collocation method. Symmetry considerations are used in the analysis of the results, which allow the classification of the perturbation modes as ${S}_{l} $ modes (those which keep the left–right symmetry) or ${R}_{x} $ modes (those which keep the symmetry of rotation of $\lrm{\pi} $ about the longitudinal mid-axis). Without throughflow, four dominant pairs of travelling transverse modes with finite wavenumbers $k$ have been found. Each pair corresponds to two symmetry degenerate left and right travelling modes which have the same critical Rayleigh number ${\mathit{Ra}}_{c} $. With the increase of the duct aspect ratio $A$, the critical Rayleigh numbers for these four pairs of modes decrease and closely approach the critical value ${\mathit{Ra}}_{c} = 1743. 894$ obtained in a two-dimensional situation, one of the mode (a ${S}_{l} $ mode called mode A) always remaining the dominant mode. Oscillatory longitudinal instabilities ($k\approx 0$) corresponding to either ${S}_{l} $ or ${R}_{x} $ modes have also been found. Their critical curves, globally decreasing, present oscillatory variations when the duct aspect ratio $A$ is increased, associated with an increasing number of longitudinal rolls. When a throughflow is applied, the symmetry degeneracy of the pairs of travelling transverse modes is broken, giving distinct upstream and downstream modes. For small and moderate aspect ratios $A$, the overall critical Rayleigh number in the small Reynolds number range studied is only determined by the upstream transverse mode A. In contrast, for larger aspect ratios as $A= 7$, different modes are successively dominant as the Reynolds number is increased, involving both upstream and downstream transverse modes A and even the longitudinal mode.

Numerical simulation of two-dimensional compressible convection

1975 ◽  
Vol 70 (4) ◽  
pp. 689-703 ◽  
Author(s):  
Eric Graham
Keyword(s):  
Rayleigh Number ◽  
Numerical Solutions ◽  
Fluid Velocity ◽  
Critical Rayleigh Number ◽  
Finite Amplitude ◽  
Two Dimensions ◽  
Onset Of Convection ◽  
Local Sound ◽  
Constant Viscosity ◽  
Convection Problem

A procedure for obtaining numerical solutions to the equations describing thermal convection in a compressible fluid is outlined. The method is applied to the case of a perfect gas with constant viscosity and thermal conductivity. The fluid is considered to be confined in a rectangular region by fixed slippery boundaries and motions are restricted to two dimensions. The upper and lower boundaries are maintained at fixed temperatures and the side boundaries are thermally insulating. The resulting convection problem can be characterized by six dimension-less parameters. The onset of convection has been studied both by obtaining solutions to the nonlinear equations in the neighbourhood of the critical Rayleigh number Rc and by solving the linear stability problem. Solutions have been obtained for values of the Rayleigh number up to 100Rc and for pressure variations of a factor of 300 within the fluid. In some cases the fluid velocity is comparable to the local sound speed. The Nusselt number increases with decreasing Prandtl number for moderate values of the depth parameter. Steady finite amplitude solutions have been found in all the cases considered. As the horizontal dimension A of the rectangle is increased, the length of time needed to reach a steady state also increases. For large values of A the solution consists of a number of rolls. Even for small values of A, no solutions have been found where one roll is vertically above another.

scholarly journals Rayleigh–Bénard convection in a creeping solid with melting and freezing at either or both its horizontal boundaries

2018 ◽  
Vol 846 ◽  
pp. 5-36 ◽  
Author(s):  
Stéphane Labrosse ◽  
Adrien Morison ◽  
Renaud Deguen ◽  
Thierry Alboussière
Keyword(s):  
Boundary Condition ◽  
Rayleigh Number ◽  
Phase Change ◽  
Time Scale ◽  
Critical Rayleigh Number ◽  
Classical Case ◽  
Bénard Convection ◽  
Benard Convection ◽  
Rayleigh Benard Convection ◽  
Rayleigh Benard

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.

Three-dimensional instability of axisymmetric buoyant convection in cylinders heated from below

1996 ◽  
Vol 326 ◽  
pp. 399-415 ◽  
Author(s):  
M. Wanschura ◽  
H. C. Kuhlmann ◽  
H. J. Rath
Keyword(s):  
Rayleigh Number ◽  
Linear Stability ◽  
Critical Rayleigh Number ◽  
Three Dimensional ◽  
Base Flow ◽  
Onset Of Convection ◽  
Buoyant Convection ◽  
Thermal Boundary Conditions ◽  
Aspect Ratios ◽  
Prandtl Numbers

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.

Thermal Convective Instabilities in a Porous Medium

1984 ◽  
Vol 106 (1) ◽  
pp. 137-142 ◽  
Author(s):  
M. Kaviany
Keyword(s):  
Porous Medium ◽  
Steady State ◽  
Temperature Distribution ◽  
Rayleigh Number ◽  
Saturated Porous Medium ◽  
Critical Rayleigh Number ◽  
Linear Stability Theory ◽  
Onset Of Convection ◽  
Long Time ◽  
Nonlinear Temperature

The onset of convection due to a nonlinear and time-dependent temperature stratification in a saturated porous medium with upper and lower free surfaces is considered. The initial parabolic temperature distribution is due to uniform internal heating. The medium is then cooled by decreasing the upper surface temperature linearly with time. Linear stability theory is applied to the more formally developed governing equations. In order to obtain an asymptotic solution for transient problems involving very long time scales, the critical Rayleigh number for steady-state, nonlinear temperature distribution is also obtained. The effects of porosity, permeability, and Prandtl number on the time of the onset of convection are examined. The steady-state results show that the critical Rayleigh number depends only on the ratio of porosity to permeability and when this ratio exceeds a value of one thousand, the critical Rayleigh number is directly proportional to this ratio.

Sign in / Sign up

Export Citation Format

Share Document