The onset of penetrative double-diffusive convection

The onset of double-diffusive convection in a horizontal fluid layer is studied. The density is assumed to depend quadratically on the temperature and linearly on the solute concentration. Under the Boussinesq approximation, the linear stability of the conduction state is investigated with respect to the oscillatory and steady convection modes. For steady onset, the critical thermal Rayleigh number is found to be a double-valued function of the solutal Rayleigh number as long as the relative maximum of the density profile exists within the fluid layer. Driving mechanisms of the steady convections are discussed.     

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.     

Double diffusive convection in a porous medium with modulated temperature on the boundaries

The effect of temperature modulation on the onset of double diffusive convection in a sparsely packed porous medium is studied by making linear stability analysis, and using Brinkman-Forchheimer extended Darcy model. The temperature field between the walls of the porous layer consists of a steady part and a time dependent periodic part that oscillates with time. Only infinitesimal disturbances are considered. The effect of permeability and thermal modulation on the onset of double diffusive convection has been studied using Galerkin method and Floquet theory. The critical Rayleigh number is calculated as a function of frequency and amplitude of modulation, Vadasz number, Darcy number, diffusivity ratio, and solute Rayleigh number. Stabilizing and destabilizing effects of modulation on the onset of double diffusive convection have been obtained. The effects of other parameters are also discussed on the stability of the system. Some results as the particular cases of the present study have also been obtained. Also the results corresponding to the Brinkman model and Darcy model have been compared.     

A numerical simulation of combined radiation and natural convection heat transfer in a square enclosure heated by a centric circular cylinder

A numerical simulation of combined natural convection and radiation in a square enclosure heated by a centric circular cylinder and filled with absorbing-emitting medium is presented. The ideal gas law and the discrete ordinates method are used to model the density changes due to temperature differences and the radiation heat transfer correspondingly. The influence of Rayleigh number, optical thickness and temperature difference on flow and temperature fields along with the natural convection, radiation and total Nusselt number at the source surfaces is studied. The results reveal that the radiation heat transfer as well as the optical thickness of the fluid has a distinct effect on the fluid flow phenomena, especially at high Rayleigh number. The heat transfer and so the Nusselt number decreases with increase in optical thickness, while increases greatly with increase in temperature difference. The variation in radiation heat transfer with optical thickness and temperature difference is much more obvious as comparison with convection heat transfer.     

Turbulent thermal convection in wide horizontal fluid layers

Turbulence in thermal convection is investigated for flows in which the production of turbulence energy is due solely to buoyancy, and the statistics of the flow are homogeneous in horizontal planes. New experimental results for high Rayleigh number unsteady turbulent convection in a horizontal layer heated from below and insulated from above are presented and compared to turbulent Rayleigh convection, convection in the planetary boundary layer, and laboratory penetrative convection. Mean temperature fields are correlated in terms of wall layer scales and convection scales. Joint statistics of turbulent temperature and horizontal velocity and vertical velocity through fourth order are presented for the core region of the convection layer.This paper was presented at the Ninth Symposium on Turbulence, University of Missouri-Rolla, October 1–3, 1984     

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

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

侧加热腔内的自然对流

开展侧加热腔内自然对流的研究具有重大的环境及工业应用背景. 总结侧加热腔内水平温差驱动的自然对流的最新研究进展, 并概述相应的流动性质、动力机制和传热特性以及对不同无量纲控制参数的依赖也有重要的科学价值. 已取得的研究结果显示突然侧加热的腔内自然对流的发展可包括初始阶段、过渡阶段和定常或准定常阶段. 不同发展阶段的流动依赖于瑞利数、普朗特数及腔体的高宽比, 且定常或准定常阶段的流态可以是定常层流流动、非定常周期性流动或者湍流流动. 此外, 回顾了对流流动失稳机制的研究成果以及湍流自然对流方面的新进展. 最后, 展望了侧加热腔内的自然对流研究的前景.      

Thermal effects of the micropolar fluid parameters on the laminar free convection in spherical annuli with mixed boundary conditions

The effects of micro-rotation and vortex viscosity in micropolar fluids have been investigated numerically to determine heat transfer by natural convection between concentric and vertically eccentric spheres with specified mixed boundary conditions. Calculations were carried out systematically for several different eccentricities and a range of modified Rayleigh numbers to determine the average Nusslet numbers which are affected by the micropolar parameters (F) of the flow and temperature fields. The skin friction stress on the walls has also been studied and discussed. The governing equations, in terms of vorticity, stream function, temperature and angular momentum are expressed in a spherical polar coordinate system. Results were obtained for steady heat-transfer in spherical annuli at a Prandtl number of 0.7, with the modified Rayleigh number ranging from 10³ to 5 × 10⁵, for a radius ratio of 2.0 and eccentricities varying from −0.625 to +0.625. Comparisons are attempted between the Newtonian fluid and micropolar fluid.     

非定常流函数涡量方程的一种数值解法的研究

对非定常流函数涡量方程的数值求解方法进行了改进,其中流函数一阶导数即速度项采用四阶精度的Hermitian公式,对流项由一般二阶精度的中心差分提高到四阶精度离散差分,包含温度方程在内的离散方程组采用ADI迭代方法求得定常解．以无内热体及有一内热体的封闭方腔内自然对流为例,进行了不同瑞利数（Ra）条件下的数值研究．结果表明,该方法推导简单,求解精度高且计算稳定,适用于封闭腔内高瑞利数复杂混合对流的数值模拟．     

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.     

Nonlinear behaviour of the instability of steady natural convection in a vertical air layer

In this paper we deal with the flow of natural convection between two vertical planes with horizontal temperature gradient. We show that periodic steady flow occur when Rayleigh number increases. Critical values are obtained numerically and a nonlinear analysis is presented according to Center manifold method.     

Natural Convection Induced by a Centrifugal Force Field in a Horizontal Annular Porous Layer Saturated with a Binary Fluid

The Darcy Model with the Boussinesq approximation is used to study natural convection in a horizontal annular porous layer filled with a binary fluid, under the influence of a centrifugal force field. Neumann boundary conditions for temperature and concentration are applied on the inner and outer boundary of the enclosure. The governing parameters for the problem are the Rayleigh number, Ra, the Lewis number, Le, the buoyancy ratio, j{\varphi } , the radius ratio of the cavity, R, the normalized porosity, e{\varepsilon } , and parameter a defining double-diffusive convection (a = 0) or Soret induced convection (a = 1). For convection in a thin annular layer (R → 1), analytical solutions for the stream function, temperature and concentration fields are obtained using a concentric flow approximation and an integral form of the energy equation. The critical Rayleigh number for the onset of supercritical convection is predicted explicitly by the present model. Also, results are obtained from the analytical model for finite amplitude convection for which the flow and heat and mass transfer are presented in terms of the governing parameters of the problem. Numerical solutions of the full governing equations are obtained for a wide range of the governing parameters. A good agreement is observed between the analytical model and the numerical simulations.     

Weak Turbulence and Chaos for Low Prandtl Number Gravity Driven Convection in Porous Media

Low Prandtl number convection in porous media is relevant to modern applications of transport phenomena in porous media such as the process of solidification of binary alloys. The transition from steady convection to chaos is analysed by using Adomian's decomposition method to obtain an analytical solution in terms of infinite power series. The practical need to evaluate the solution and obtain numerical values from the infinite power series, the consequent series truncation, and the practical procedure to accomplish this task, transform the analytical results into a computational solution evaluated up to a finite accuracy. The solution shows a transition from steady convection to chaos via a Hopf bifurcation producing a 'solitary limit cycle which may be associated with an homoclinic explosion. This occurs at a slightly subcritical value of Rayleigh number, the critical value being associated with the loss of linear stability of the steady convection solution. Periodic windows within the broad band of parameter regime where the chaotic solution persists are identified and analysed. It is evident that the further transition from chaos to a high Rayleigh number periodic convection occurs via a period halving sequence of bifurcations.     

Effect of modulation on onset of thermal convection of a second-grade fluid layer

This study examines the stability of a horizontally extended second-grade fluid layer heated from below, when a steady temperature difference between the walls is superimposed on sinusoidal temperature perturbations. A linear stability analysis proposed by Venezian (J. Fluid Mech. 35 (1969) 243) is employed to obtain the critical Rayleigh numbers for different types of temperature modulation. The free–free and isothermal boundary conditions are considered so as to allow analytic solutions. The stability characterized by the shift in critical Rayleigh number R_2c is calculated as a function of the modulation frequency ω, the Prandtl number Pr, and the viscoelastic parameter Q. It is found that the onset of convection can be delayed or advanced by these parameters.     

Effect of Thermal Modulation on the Onset of Convection in a Porous Medium Layer Saturated by a Nanofluid

The effect of time-periodic temperature modulation at the onset of convection in a Boussinesq porous medium saturated by a nanofluid is studied analytically. The model used for the nanofluid incorporates the effects of Brownian motion. Three types of boundary temperature modulations are considered namely, symmetric, asymmetric, and only the lower wall temperature is modulated while the upper wall is held at constant temperature. The perturbation method is applied for computing the critical Rayleigh and wave numbers for small amplitude temperature modulation. The shift in the critical Rayleigh number is calculated as a function of frequency of modulation, concentration Rayleigh number, porosity, Lewis number, and thermal capacity ratio. It has been shown that it is possible to advance or delay the onset of convection by time-periodic modulation of the wall temperature. The nanofluid is found to have more stabilizing effect when compared to regular fluid. Low frequency is destabilizing, while high frequency is always stabilizing for symmetric modulation. Asymmetric modulation and only lower wall temperature modulation is stabilizing for all frequencies when concentration Rayleigh number is greater than one.     

Numerical study of convection of a liquid heated from below

The equilibrium of a liquid heated from below is stable only for small values of the vertical temperature gradient. With increase of the temperature gradient a critical equilibrium situation occurs, as a result of which convection develops. If the liquid fills a closed cavity, then there is a discrete sequence of critical temperature gradients (Rayleigh numbers) for which the equilibrium loses stability with respect to small characteristic disturbances. This sequence of critical gradients and motions may be found from the solution of the linear problem of equilibrium stability relative to small disturbances. If the temperature gradient exceeds the lower critical value, then (for steady-state heating conditions) there is established in the liquid a steady convective motion of a definite amplitude which depends on the magnitude of the temperature gradient. Naturally, the amplitude of the steady convective motion cannot be determined from linear stability theory; to find this amplitude we must solve the problem of convection with heating from below in the nonlinear formulation. A nonlinear study of the steady motion of a liquid in a closed cavity with heating from below was made in [1]. In that study it was shown that for Rayleigh numbers R which are less than the lower critical value R_c steady-state motions of the liquid are not possible. With R>R_c a steady convection arises, whose amplitude near the threshold is small and proportional to (R–R_c)^1/2 (the so-called soft instability)-this is in complete agreement with the results of the phenom-enological theory of Landau [2, 3].Primarily, various versions of the method of expansion in powers of the amplitude [4–8] have been used, and, consequently, the results obtained in those studies are valid only for values of R which are close to R_c, i. e., near the convection threshold.It is apparent that the study of developed convective motion far from the threshold can be carried out only numerically, with the use of digital computers. In [9, 10] the numerical methods have been successfully used for the study of developed convection in an infinite plane horizontal liquid layer.The present paper undertakes the numerical study of plane convective motions of a liquid in a closed cavity of square section. The complete nonlinear system of convection equations is solved by the method of finite differences on a digital computer for various values of the Rayleigh number, the maximal value exceeding by a factor of 40 the minimal critical value R_c. The numerical solution permits following the development of the steady motion which arises with R>R_c in the course of increase of the Rayleigh number and permits study of the oscillatory motions which occur at some value of the parameter R. The heat transfer through the cavity is studied. The corresponding linear problem on equilibrium stability is solved approximately by the Galerkin method.     

Linear and non-linear double diffusive convection in a rotating porous layer using a thermal non-equilibrium model

Double diffusive convection in a fluid-saturated rotating porous layer heated from below and cooled from above is studied when the fluid and solid phases are not in local thermal equilibrium, using both linear and non-linear stability analyses. The Darcy model that includes the time derivative and Coriolis terms is employed as momentum equation. A two-field model that represents the fluid and solid phase temperature fields separately is used for energy equation. The onset criterion for stationary, oscillatory and finite amplitude convection is derived analytically. It is found that small inter-phase heat transfer coefficient has significant effect on the stability of the system. There is a competition between the processes of thermal and solute diffusions that causes the convection to set in through either oscillatory or finite amplitude mode rather than stationary. The effect of solute Rayleigh number, porosity modified conductivity ratio, Lewis number, diffusivity ratio, Vadasz number and Taylor number on the stability of the system is investigated. The non-linear theory based on the truncated representation of Fourier series method predicts the occurrence of subcritical instability in the form of finite amplitude motions. The effect of thermal non-equilibrium on heat and mass transfer is also brought out.     

Interaction of turbulent natural convection and surface thermal radiation in inclined square enclosures

Two-dimensional numerical studies of flow and temperature fields for turbulent natural convection and surface radiation in inclined differentially heated enclosures are performed. Investigations are carried out over a wide range of Rayleigh numbers from 10⁸ to 10¹², with the angle of inclination varying between 0° and 90°. Turbulence is modeled with a novel variant of the k–ε closure model. The predicted results are validated against experimental and numerical results reported in literature. The effect of the inclination of the enclosure on pure turbulent natural convection and the latter’s interaction with surface radiation are brought out. Profiles of turbulent kinetic energy and effective viscosity are studied to observe the net effect on the intensity of turbulence caused by the interaction of natural convection and surface radiation. The variations of local Nusselt number and average Nusselt number are presented for various inclination angles. Marked change in the convective Nusselt number is found with the orientation of enclosure. Also analyzed is the influence of change in emissivity on the flow and heat transfer. A correlation relevant to practical applications in the form of average Nusselt number, as a function of Rayleigh number, Ra, radiation convection parameter, N _RC and inclination angle of the enclosure, φ is proposed.     

Aspect ratio effect on natural convection in water near its density maximum temperature

The effect of the aspect ratio on natural convection in water subjected to density inversion has been investigated in this study. Numerical simulations of the two-dimensional, steady state, incompressible flow in a rectangular enclosure with a variety of aspect ratios, ranging from 0.125 to 100, have been accomplished using a finite element model. Computations cover Rayleigh numbers from 10³ to 10⁶. Results reveal that the aspect ratio, A, the Rayleigh number, Ra, and the density distribution parameter, R, are the key parameters to determine the heat transfer and fluid flow characteristics for density inversion fluids in an enclosure. A new correlation for predicting the maximum mean Nusselt number is proposed in the form of , with the constants a and b depending on density distribution number R. It is demonstrated that the aspect ratio has a strong impact on flow patterns and temperature distributions in rectangular enclosures. The stream function ratio Ψ_inv/|Ψ_reg| is introduced to describe quantitatively the interaction between inversional and regular convection. For R=0.33, the density inversion enhancement is observed in the regime near A=3.     

A comparative study on low‐order,amplitude‐equation and perturbation approaches in thermal convection

A comparison among three weakly nonlinear approaches for thermo‐gravitational instability in a Newtonian fluid layer heated from below is presented. First, the dynamical systems describing the time evolution of the problem from different weakly nonlinear approaches, namely, the Lorenz model, the amplitude equations and the perturbation expansion approaches are obtained. Next, the steady states and their stability, as well as the transient behaviour are obtained from each dynamical system. The similarity and difference among the three models are emphasized. The role of each of the nondimensional groups, the Rayleigh number and the Prandtl number is compared for the three models. The different approaches lead to similar behaviours when the Rayleigh number is just above its critical value and Prandtl number is high. However, only the dynamical system obtained from the amplitude equations is able to reflect the role of the Prandtl number. On the other hand, the amplitude equations and perturbation expansion techniques are not suitable for predicting the uniform oscillatory behaviour observed frequently in Rayleigh–Bénard convection. The novelty of the current work lies in studying the critical differences in the findings of the three popular approaches to investigate weakly nonlinear thermal convection for the first time. Copyright © 2011 John Wiley & Sons, Ltd.