Mechanism of mode selection in Rayleigh—Bénard convection with free—rigid boundaries

Higher harmonic resonances with wavenumber ratio of 1:2, 1:2:3 and so on are shown to take place in Rayleigh—Bénard convection under free—rigid boundary condition. Bifurcation diagrams for two-dimensional motion are obtained for the Prandtl number P = 7. The subharmonic instability is explained by a couple of amplitude equations obtained from weakly nonlinear stability theory. A straightforward extension of the coupled amplitude equations leads to a model which consists of n amplitude equations. The mechanism of mode selection is illustrated by numerical simulations of the model equations.     

Application of a non-linear local analysis method for the problem of mixed convection instability

We consider the problem of laminar mixed convection flow between parallel, vertical and uniformly heated plates where the governing dimensionless parameters are the Prandtl, Rayleigh and Reynolds numbers. Using the method based on the centre manifold theorem which was derived from the general theory of dynamical systems, we reduce a three-dimensional simplified model of ordinary differential amplitude equations emanating from the original Navier-Stokes system of the problem in the vicinity of a trivial stationary solution. We have found that when the forcing parameter, the Rayleigh number, increases beyond the critical value Ra_s, the stationary solution is a pitchfork bifurcation point of the system.     

Combined effect of thermal modulation and rotation on the onset of stationary convection in a porous layer

The stability of a fluid-saturated horizontal rotating porous layer subjected to time-periodic temperature modulation is investigated when the condition for the principle of exchange of stabilities is valid. The linear stability analysis is used to study the effect of infinitesimal disturbances. A regular perturbation method based on small amplitude of applied temperature field is used to compute the critical values of Darcy–Rayleigh number and wavenumber. The shift in critical Darcy–Rayleigh number is calculated as a function of frequency of modulation, Taylor number, and Darcy–Prandtl number. It is established that the convection can be advanced by the low frequency in-phase and lower-wall temperature modulation, where as delayed by the out-of-phase modulation. The effect of Taylor number and Darcy–Prandtl number on the stability of the system is also discussed. We found that by proper tuning of modulation frequency, Taylor number, and Darcy–Prandtl number it is possible to advance or delay the onset of convection.     

The Effect of Rotation on the Onset of Double Diffusive Convection in a Sparsely Packed Anisotropic Porous Layer

The effect of rotation on the onset of double diffusive convection in a sparsely packed anisotropic porous layer, which is heated and salted from below, is investigated analytically using the linear and nonlinear theories. The Brinkman model that includes the Coriolis term is employed for the momentum equation. The critical Rayleigh number, wavenumber for stationary and oscillatory modes and a dispersion relation are obtained analytically using linear theory. The effect of anisotropy parameters, Taylor number, Darcy number, solute Rayleigh number, Lewis number, Darcy–Prandtl number, and normalized porosity on the stationary, oscillatory and finite amplitude convection is shown graphically. It is found that contrary to its usual influence on the onset of convection in the absence of rotation, the mechanical anisotropy parameter show contrasting effect on the onset criterion at moderate and high rotation rates. The nonlinear theory based on the truncated representation of Fourier series method is used to find the heat and mass transfers. The effect of various parameters on heat and mass transfer is shown graphically. Some of the convection systems previously reported in the literature is shown to be special cases of the system presented in this study.     

Numerical study of steady magneto-convection around an adiabatic body inside a square enclosure in low Prandtl numbers

In this paper, the effects of Prandtl number on the steady magneto-convection around a centrally located adiabatic body inside a square enclosure are numerically investigated. Two-dimensional nonlinear governing equations are discretized using the control volume method and hybrid scheme. The equations are solved using SIMPLER algorithm. The results are displayed in the form of streamlines and isotherms when the Rayleigh number varies between 10³ and 10⁶, the Hartmann number changes between 0 and 100 and the Prandtl number ranges between 0.005 and 0.1. The ratio of the buoyancy force to the Lorentz force (Ra/Ha ²) is introduced as an index to compare the contribution of natural convection and magnetic field strength on heat transfer. The results obtained from numerical modeling show that the Prandtl number has not considerable effect on heat transfer at low Rayleigh numbers. The effect of magnetic field strength on convection is increased by increasing Prandtl number. The effect of Prandtl number on the average Nusselt number in the presence of a magnetic field is less than the case without a magnetic field.     

The onset of double diffusive convection in a viscoelastic fluid layer

The onset of double diffusive convection in a viscoelastic fluid layer is studied using a linear and a weak nonlinear stability analyses. The onset criterion for stationary, oscillatory and finite amplitude convection is derived analytically. There is a competition between the processes of thermal diffusion, solute diffusion and viscoelasticity that causes the convection to set in through oscillatory mode rather than stationary. The effect of Deborah number, retardation parameter, solutal Rayleigh number, Prandtl number, Lewis number on the stability of the system is investigated. It is shown that the critical frequency increases with Deborah number and solutal Rayleigh number while it decreases with retardation parameter and Lewis number. The nonlinear theory based on the truncated representation of Fourier series method is used to find the heat and mass transfers. The transient behaviour of the Nusselt number and Sherwood number is investigated by solving the finite amplitude equations using Runge-Kutta method. The effect of viscoelastic parameters on heat and mass transfer is brought out.     

The Prandtl number effect near the onset of Bénard convection in a porous medium

This note focuses on Kladias and Prasad's claim that the critical Rayleigh number for the onset of Bénard convection in an infinite horizontal porous layer increases as the Prandtl number decreases, and argues that the critical Rayleigh number (Ra_c) depends only on the Darcy number (Da), as linear stability analysis indicates. The two-dimensional steady-convection problem is then solved numerically to document the convection heat transfer effect of the Rayleigh number, Darcy number, Prandtl number, and porosity. The note concludes with an empirical correlation for the overall Nusselt number, which shows the effect of Prandtl number at above-critical Rayleigh numbers. The correlation is consistent with the corresponding correlation known for Bénard convection in a pure fluid.     

矩形倾斜腔体中的两种对流斑图和分区

为了研究矩形倾斜腔体中普朗特数Pr=0.72的流体对流斑图和分区,本文基于流体力学方程组进行了数值模拟。在相对瑞利数r=6.0的情况下,观察了倾角θ=10°和θ=60°时对流斑图随着时间的发展,发现系统存在单圈型对流和多圈型对流两种斑图。流线随着倾角的变化说明:随着倾角增加,对流圈数逐渐减少,对流波长逐渐增加,对流波数减小;然后,随着对流圈数显著地减少,系统由多圈型对流过渡到单圈型对流。根据模拟计算结果,给出了多圈型对流过渡到单圈型对流的临界倾角θc随着相对瑞利数r变化的关系曲线。对流在θ-r平面上分为两个区域:θ<θc时系统是单圈型对流,θ>θc时系统是多圈型对流。对流最大振幅A和努塞尔数Nu随着倾角θ的变化曲线被临界倾角θc分成两段,它们有不同的变化规律。因此,临界倾角也可以由对流最大振幅A或努塞尔数Nu的变化曲线来确定。     

Coriolis effect on convection for a low Prandtl number fluid

The problem of finite-amplitude thermal convection in a horizontal layer of a low Prandtl number heated from below and rotating about a vertical axis is studied. Linear stability and weak non-linear theories are used to investigate analytically the Coriolis effect on gravity-driven convection. The non-linear steady problem is solved by perturbation techniques, and the preferred mode of convection is determined by a stability analysis. Finite-amplitude results, obtained by using a weak amplitude, correspond to both stationary and oscillatory convections. These amplitude equations permit to identify from the post-transient conditions that the fluid is subject to Pitchfork bifurcation in the stationary convection and Hopf bifurcation in the oscillatory convection. Thereafter, in the small perturbations hypothesis, an amplitude solution is evaluated and drawn in time and space scales.     

Coriolis effect on thermal convection in a couple-stress fluid-saturated rotating rigid porous layer

Both linear and weakly nonlinear stability analyses are performed to study thermal convection in a rotating couple-stress fluid-saturated rigid porous layer. In the case of linear stability analysis, conditions for the occurrence of possible bifurcations are obtained. It is shown that Hopf bifurcation is possible due to Coriolis force, and it occurs at a lower value of the Rayleigh number at which the simple bifurcation occurs. In contrast to the nonrotating case, it is found that the couple-stress parameter plays a dual role in deciding the stability characteristics of the system, depending on the strength of rotation. Nonlinear stability analysis is carried out by constructing a set of coupled nonlinear ordinary differential equations using truncated representation of Fourier series. Sub-critical finite amplitude steady motions occur depending on the choice of physical parameters but at higher rotation rates oscillatory convection is found to be the preferred mode of instability. Besides, the stability of steady bifurcating equilibrium solution is discussed using modified perturbation theory. Heat transfer is calculated in terms of Nusselt number. Also, the transient behavior of the Nusselt number is investigated by solving the nonlinear differential equations numerically using the Runge–Kutta–Gill method. It is noted that increase in the value of Taylor number and the couple-stress parameter is to dampen the oscillations of Nusselt number and thereby to decrease the heat transfer.     

Convective heat transfer and temperature stratification in a sphere completely filled with a liquid,with a given heat flux

A number of articles have been devoted to the theoretical and experimental investigation of natural convection in spherical vessels completely filled with a liquid [1–6]. Analytical solutions are known, obtained by the expansion of the sought function in series in powers of the Rayleigh number (see, for example, [1]), valid for very small values of this number. A numerical solution of the nonlinear Boussinesq equations can be used to obtain solutions with larger Rayleigh numbers, but the existing data for spherical regions [2, 3] embrace a relatively narrow range of Rayleigh numbers. The experimental data with a given heat flux, published in [4–6], were obtained with relatively large Rayleigh numbers (Ra*=10⁹?10¹¹) and Prandtl numbers (P= 3?1500). Data on the characteristics of convection in spherical vessels are still not very numerous and, in a number of cases, contradictory. This relates, in particular, to the boundaries of unsteady-state conditions. The present article, continuing [7–9], expounds a method and gives the results of a calculation of convection in a sphere with a thinwalled shell, in a range of Rayleigh and Fourier numbers embracing the principal conditions of unsteady-state laminar convection with a given heat flux.     

Effect of Thermal/Gravity Modulation on the Onset of Convection in a Maxwell Fluid Saturated Porous Layer

The effect of thermal/gravity modulation on the onset of convection in a Maxwell fluid saturated porous layer is investigated by a linear stability analysis. Modified Darcy–Maxwell model is used to describe the fluid motion. The regular perturbation method based on the small amplitude of modulation is employed to compute the critical Rayleigh number and the corresponding wavenumber. The stability of the system characterized by a correction Rayleigh number is calculated as a function of the viscoelastic parameter, Darcy–Prandtl number, normalized porosity, and the frequency of modulation. It is found that the low frequency symmetric thermal modulation is destabilizing while moderate and high frequency symmetric modulation is always stabilizing. The asymmetric modulation and lower wall temperature modulations are, in general, stabilizing while the system becomes unstable for large values of Darcy–Prandtl number and for small frequencies. It is shown that in general the gravity modulation produces a stabilizing effect on the onset of convection for moderate and high frequency. The small frequency gravity modulation is found to have destabilizing effect on the stability of the system.     

Route to Chaos for Moderate Prandtl Number Convection in a Porous Layer Heated from Below

The route to chaos for moderate Prandtl number gravity driven convection in porous media is analysed by using Adomian's decomposition method which provides an accurate analytical solution in terms of infinite power series. The practical need to evaluate numerical values from the infinite power series, the consequent series truncation, and the practical procedure to accomplish this task, transform the otherwise analytical results into a computational solution achieved up to a desired but finite accuracy. The solution shows a transition to chaos via a period doubling sequence of bifurcations at a Rayleigh number value far beyond the critical value associated with the loss of stability of the convection steady solution. This result is extremely distinct from the sequence of events leading to chaos in low Prandtl number convection in porous media, where a sudden transition from steady convection to chaos associated with an homoclinic explosion occurs in the neighbourhood of the critical Rayleigh number (unless mentioned otherwise by 'the critical Rayleigh number' we mean the value associated with the loss of stability of the convection steady solution). In the present case of moderate Prandtl number convection the homoclinic explosion leads to a transition from steady convection to a period-2 periodic solution in the neighbourhood of the critical Rayleigh number. This occurs at a slightly sub-critical value of Rayleigh number via a transition associated with a period-1 limit cycle which seem to belong to the sub-critical Hopf bifurcation around the point where the convection steady solution looses its stability. The different regimes are analysed and periodic windows within the chaotic regime are identified. The significance of including a time derivative term in Darcy's equation when wave phenomena are being investigated becomes evident from the results.     

Convective Instability in Slip Flow in a Vertical Circular Porous Microchannel

The paper represents an analysis of convective instability in a vertical cylindrical porous microchannel performed using the Galerkin method. The dependence of the critical Rayleigh number on the Darcy, Knudsen, and Prandtl numbers, as well as on the ratio of the thermal conductivities of the fluid and the wall, was obtained. It was shown that a decrease in permeability of the porous medium (in other words, increase in its porosity) causes an increase in flow stability. This effect is substantially nonlinear. Under the condition Da?>?0.1, the effect of the porosity on the critical Rayleigh number practically vanishes. Strengthening of the slippage effects leads to an increase in the instability of the entire system. The slippage effect on the critical Rayleigh number is nonlinear. The level of nonlinearity depends on the Prandtl number. With an increase in the Prandtl number, the effect of slippage on the onset of convection weakens. With an increase in the ratio of the thermal conductivities of the fluid and the wall, the influence of the Prandtl number decreases. At high values of the Prandtl numbers (Pr?>?10), its influence practically vanishes.

     

Rayleigh–Benard convection subject to time dependent wall temperature/gravity in a fluid-saturated anisotropic porous medium

 The effect of time-periodic temperature/gravity modulation at the onset of convection in a Boussinesq fluid-saturated anisotropic porous medium is investigated by making a linear stability analysis. Brinkman flow model with effective viscosity larger than the viscosity of the fluid is considered to give a more general theoretical result. The perturbation method is applied for computing the critical Rayleigh and wave numbers for small amplitude temperature/gravity modulation. The shift in the critical Rayleigh number is calculated as a function of frequency of the modulation, viscosity ratio, anisotropy parameter and porous parameter. We have shown that it is possible to advance or delay the onset of convection by time-periodic modulation of the wall temperature and to advance convection by gravity modulation. It is also shown that the small anisotropy parameter has a strong influence on the stability of the system. The effect of viscosity ratio, anisotropy parameter, the porous parameter and the Prandtl number is discussed. Received on 28 July 2000 / Published online: 29 November 2001     

Onset of unsteady axi-symmetric laminar natural convection in a vertical cylindrical enclosure heated at the wall

In the present study laminar transition to oscillatory convection of fluids having different Prandtl numbers in a laterally heated vertical cylindrical enclosure for different aspect ratios (melt height to crucible radius) of 2–4 is investigated numerically for 0.01 ≤ Pr ≤ 10. Numerical solution to two-dimensional axisymmetric transient Navier Stokes equations and energy equation were solved by finite volume method using SIMPLE algorithm. Numerical results illustrate that there exists a critical Rayleigh number for each Prandtl number beyond which sustained laminar oscillatory flow sets in. The oscillatory regime was characterised by the oscillation of the average kinetic energy and average thermal energy of the melt. For a given aspect ratio, critical Rayleigh number increases with Pr upto 1 and then flattens. It was observed that for low Prandtl number fluids, Pr < 1.0, critical Rayleigh number is found to increase with increase in aspect ratio while for high Prandtl number fluids, Pr ≥ 1.0, it is found to decrease with increase in aspect ratio. The influence of aspect ratio on the transient behaviour of the melt volume below and above the critical Rayleigh number was studied.     

具有强SORET效应的充分发展的混合流体行进波对流

混合流体Rayleigh-Benard对流是研究对流稳定性,时空结构和非线性特性的典型模型之一。本文利用流体力学扰动方程组的数值模拟,讨论了偏离传导状态具有强SORET效应的混合流体行进波对流的温度场和浓度场的成长过程,分析了充分发展对流情况下的对流振幅,Nusselt数及混合参数与相对瑞利数的关系。并给出了行进波相速度对相对瑞利数的依赖关系。结果说明混合参数的曲线与行进波相速度的分布曲线是类似的。文末,给出了垂直速度,温度和浓度场的分布并讨论了相对瑞利数对场的分布及不同场之间的相位差的影响。     

多孔介质内黏弹性流体的热对流稳定性研究

基于修正的Darcy模型, 介绍了多孔介质内黏弹性流体热对流稳定性研究的现状和主要进展. 通过线性稳定性理论, 分析计算多孔介质几何形状(水平多孔介质层、多孔圆柱以及多孔方腔)、热边界条件(底部等温加热、底部等热流加热、底部对流换热以及顶部自由开口边界)、黏弹性流体的流动模型(Darcy-Jeffrey, Darcy-Brinkman-Oldroyd以及Darcy-Brinkman -Maxwell模型)、局部热非平衡效应以及旋转效应对黏弹性流体热对流失稳的临界Rayleigh数的影响. 利用弱非线性分析方法, 揭示失稳临界点附近热对流流动的分叉情况, 以及失稳临界点附近黏弹性流体换热Nusselt数的解析表达式. 采用数值模拟方法, 研究高Rayleigh数下黏弹性流体换热Nusselt数和流场的演化规律,分析各参数对黏弹性流体热对流失稳和对流换热速率的影响.主要结果:(1)流体的黏弹性能够促进振荡对流的发生;(2)旋转效应、流体与多孔介质间的传热能够抑制黏弹性流体的热对流失稳;(3)在临界Rayleigh数附近,静态对流分叉解是超临界稳定的, 而振荡对流分叉可能是超临界或者亚临界的,主要取决于流体的黏弹性参数、Prandtl数以及Darcy数;(4)随着Rayleigh数的增加,热对流的流场从单个涡胞逐渐演化为多个不规则单元涡胞, 最后发展为混沌状态.      

Onset of Darcy–Brinkman convection in a binary viscoelastic fluid saturated porous layer

The onset of Darcy–Brinkman double-diffusive convection in a binary viscoelastic fluid-saturated porous layer is studied using both linear and weakly nonlinear stability analyses. The Oldroyd-B model is employed to describe the rheological behavior of the fluid. An extended form of Darcy–Oldroyd law incorporating the Brinkman’s correction and time derivative is used to describe the fluid flow and the Oberbeck–Boussinesq approximation is invoked. The onset criterion for stationary and oscillatory convection is derived analytically. The effects of rheological parameters, Darcy number, normalized porosity, Lewis number, solute Rayleigh number, and Darcy–Prandtl number on the stability of the system is investigated. The results indicated that there is a competition among the processes of thermal, solute diffusions and viscoelasticity that causes the convection to set in through the oscillatory modes rather than the stationary. The Darcy–Prandtl number has a dual effect on the threshold of oscillatory convection. The nonlinear theory based on the method of truncated representation of Fourier series is used to find the transient heat and mass transfer. Some existing results are reproduced as the particular cases of present study.     

Fluid convection in a rotating porous layer under modulated temperature on the boundaries

The linear stability of thermal convection in a rotating horizontal layer of fluid-saturated porous medium, confined between two rigid boundaries, is studied for temperature modulation, using Brinkman’s model. In addition to a steady temperature difference between the walls of the porous layer, a time-dependent periodic perturbation is applied to the wall temperatures. Only infinitesimal disturbances are considered. The combined effect of rotation, permeability and modulation of walls’ temperature on the stability of flow through porous medium has been investigated using Galerkin method and Floquet theory. The critical Rayleigh number is calculated as function of amplitude and frequency of modulation, Taylor number, porous parameter and Prandtl number. It is found that both, rotation and permeability are having stabilizing influence on the onset of thermal instability. Further it is also found that it is possible to advance or delay the onset of convection by proper tuning of the frequency of modulation of the walls’ temperature.