Laminar natural convection in a partially divided rectangular cavity at high Rayleigh number

Finite element predictions of two-dimensional laminar natural convection in a partially divided rectangular cavity at high Rayleigh number are presented. The walls are differentially heated, the horizontal surfaces are insulated and the cavity contains a partial vertical divider which is centrally located and whose height is varied. Detailed results are presented for an aluminium half-divider in water for Rayleigh number up to 10¹¹ and compared directly with recent experiments in a cavity of aspect ratio 1/2. The predicted flow and heat transfer are in good agreement with the measurements and confirm the existence of a high Rayleigh number regime with characteristic behaviour that differs significantly from that found at lower Rayleigh number. In addition, the effects of the divider height, the divider conductivity, the fluid Prandtl number and the cavity aspect ratio are studied. The results show that a direct simulation of the complex flow and heat transfer that occurs in partially divided cavities is possible for realistic physical conditions.     

Natural convection flow of a viscous incompressible fluid in a rectangular porous cavity heated from below with cold sidewalls

We consider the flow, which is induced by differential heating on the boundaries of a porous cavity heated from below. In particular we allow the sidewalls to have the same cold temperature as the upper surface, and thus the problem is a variant of the Darcy-Bénard convection problem, but one where there is flow at all non-zero Grashof numbers. Attention is focused on how the flow and heat transfer is affected by variations in the cavity aspect ratio, the Grashof number and the Darcy number. The flow becomes weaker as the Darcy number decreases from the pure fluid limit towards the Darcy-flow limit. In addition the number of cells which form in the cavity varies primarily with the aspect ratio and is always even due to the symmetry imposed by the cold sidewalls.     

Time-dependent natural convection in a square cavity: Application of a new finite volume method

A new finite volume (FV) approach with adaptive upwind convection is used to predict the two-dimensional unsteady flow in a square cavity. The fluid is air and natural convection is induced by differentially heated vertical walls. The formulation is made in terms of the vorticity and the integral velocity (induction) law. Biquadratic interpolation formulae are used to approximate the temperature and vorticity fields over the finite volumes, to which the conservation laws are applied in integral form. Image vorticity is used to enforce the zero-penetration condition at the cavity walls. Unsteady predictions are carried sufficiently forward in time to reach a steady state. Results are presented for a Prandtl number (Pr) of 0-71 and Rayleigh numbers equal to 10³, 10⁴ and 10⁵. Both 11 × 11 and 21 × 21 meshes are used. The steady state predictions are compared with published results obtained using a finite difference (FD) scheme for the same values of Pr and Ra and the same meshes, as well as a numerical bench-mark solution. For the most part the FV predictions are closer to the bench-mark solution than are the FD predictions.     

Transient natural convection in a rectangular enclosure with one heated side wall

This paper describes a numerical and theoretical study of the transient natural convection heating of a two-dimensional rectangular enclosure filled with fluid. The heating is applied suddenly along one of the side walls, while the remaining three walls are maintained insulated. It is shown that the process has two distinct phases, an early period dominated by conduction and a late period dominated by convection. The scaling laws for the heat transfer rate and the effectiveness (energy storage fraction) are determined based on scale analysis. These theoretical results are confirmed by numerical experiments conducted in the domain Ra = 10³−10⁶, Pr = 7, A = 1, where Ra is the Rayleigh number based on height and initial temperature difference, Pr is the Prandtl number, and A is the height/length ratio of the enclosure. Correlations for heat transfer rate and effectiveness are constructed by comparing the theoretical scaling laws with the numerical results.     

Transient natural convection of low-Prandtl-number fluids in a differentially heated cavity

The transient convective motion in a two-dimensional square cavity driven by a temperature gradient is analysed. The cavity is filled with a low-Prandtl-number fluid and the vertical walls are maintained at constant but different temperatures, while the horizontal boundaries are adiabatic. A control volume approach with a staggered grid is employed to formulate the finite difference equations. Numerically accurate solutions are obtained for Prandtl numbers of 0·001, 0·005 and 0·01 and for Grashof numbers up to 1 × 10⁷. It was found that the flow field exhibits periodic oscillation at the critical Grashof numbers, which are dependent on the Prandtl number. As the Prandtl number is decreased, the critical Grashof number and the frequency of oscillation decrease. Prior to the oscillatory flow, steady state solutions with an oscillatory transient period were predicted. In addition to the main circulation, four weak circulations were predicted at the corners of the cavity.     

Transient natural convection in a cavity with heat input and a constant temperature wall on opposite sides

The transient natural convection of a fluid with Prandtl number of order 200 in a two-dimensional square cavity has been numerically studied. One of the vertical walls of the cavity is kept at a constant (ambient) temperature and a constant heat flux is applied on the opposite wall. The other walls are adiabatic. Initially, a boundary layer is formed near the heated wall; subsequently, a large vortical structure is generated, together with an upper intrusion layer. As time progresses, the average temperature in the cavity increases, and a descending boundary layer is formed near the constant temperature wall. During the transition to the steady-state regime, a thermal stratification pattern is formed. The results are compared with the scale analysis presented by Patterson and Imberger (1980).     

Three-dimensional simulations of natural convection in a sidewall-heated cube

A high-resolution, finite difference numerical study is reported on three-dimensional natural convection of air in a differentially heated cubical enclosure over an extensive range of Rayleigh number from 10³ to 10¹⁰. The maximum number of grid points is 122 × 62 × 62. Solutions to the primitive variable formulation of the incompressible Navier-Stokes and energy equations are acquired by a control-volume-based procedure together with a higher-order upwind-differencing technique. The field characteristics at large-time limits are examined in detail by state-of-the-art numerical visualizations of the three-dimensional results. The emergence of the well-defined boundary layers and the interior core at high Rayleigh numbers is captured by elaborate numerical visualizations. Both the similarities and discrepancies between the three- and two-dimensional computations are pointed out. These emphasize the need for three-dimensional calculations to accurately determine the flow characteristics and heat transfer properties in realistic, high-Rayleigh-number situations.     

Experiments on turbulent natural convection in an enclosed tall cavity

Experiments have been undertaken to investigate the natural convection of air in a tall differentially heated rectangular cavity (2.18 m high by 0.076 m wide by 0.52 m in depth). They were performed with temperature differentials between the vertical plates of 19.6°C and 39.9°C, giving Rayleigh numbers based on the width of 0.86×10⁶ and 1.43×10⁶. Under these conditions the flow in the core of the cavity is fully turbulent and property variations with temperature are comparatively small. A previously used experimental rig has been modified, by fitting partially conducting top and bottom walls and outer guard channels, to provide boundary conditions which avoid the inadequately defined sharp changes in temperature gradient and other problems associated with insufficient insulation on nominally adiabatic walls. Mean and turbulent temperature and velocity variations within the cavity have been measured, together with heat fluxes and turbulent shear stresses. The temperature and flow fields were found to be closely two-dimensional, except close to the front and back walls, and anti-symmetric across the diagonal of the cavity. The partially conducting roof and floor provide locally unstable thermal stratification in the wall jet flows there, which enhances the turbulence as the flow moves towards the temperature controlled plates. The results provide a greatly improved benchmark for the testing of turbulence models in this low turbulence Reynolds number flow.     

Natural convection in a cavity with time-dependent flux boundary

Numerical simulations have been carried out to investigate the unsteady natural convection flow in a cavity subjected to a sidewall heat flux varying sinusoidally with time. With all walls non-slip and the upper and lower boundaries and the other sidewall adiabatic, the heating and cooling produces an alternating direction natural convection boundary layer that discharges hot fluid to the top and cold fluid to the bottom of the cavity, generating a time-varying thermal stratification in the cavity interior. Scaling analysis has been conducted for different flow regimes based on the forcing frequency, with the characteristic time scales being the forcing period and the boundary layer development time. The scaling relations are then verified using the simulations, with the results showing overall good agreement with the derived scaling relations.     

Combined forced and free flow in a vertical rectangular duct with prescribed wall heat flux

In this paper, combined forced and free convection is studied in a vertical rectangular duct with a prescribed uniform wall heat flux (H2 boundary condition). A different heat flux value for each plane wall is considered; the condition of a uniform wall heat flux throughout the duct results as a special case. The local momentum and energy balance equations are written in a dimensionless form and solved numerically, by means of a Galerkin finite element method. The numerical solution gives the dimensionless velocity and temperature distributions, together with the values of the Fanning friction factor, of the Nusselt number, of the momentum flux correction factor and of the kinetic energy correction factor. These dimensionless parameters are reported as functions of the aspect ratio and of the ratio between the Grashof number, Gr, and the Reynolds number, Re. The threshold values of Gr/Re for the onset of flow reversal are evaluated.     

Heat transfer by conduction,natural convection and radiation across a rectangular cellular structure

This paper describes results on the effects of wall conduction and radiation heat exchange among surfaces on laminar natural convection heat transfer in a two-dimensional rectangular cavity modelling a cellular structure. Parametric heat transfer calculations have been performed, and numerical results are presented in graphical and tabular form. Local and average Nusselt numbers along the cavity walls are reported for a range of parameters of physical interest. The findings suggest that the local or the average Nusselt number is one of many parameters that control conjugate heat transfer problems. The results indicate that natural convection heat transfer in the cavity is reduced by heat conduction in the walls and radiation exchange among surfaces. The results obtaibed for the total heat transfer rate through the system using the two-dimensional model are compared with those based on a one-dimensional model.     

Laminar natural convection in a fully partitioned enclosure containing fluid with nonlinear thermophysical properties

The effects of a heat conducting partition on the laminar natural convection heat transfer and fluid flow were obtained by comparing the numerical and experimental results for a cubic enclosure without and with a partition. The two opposite vertical walls of the enclosure were isothermal at different temperatures. The working fluid was glycerol. The complete vertical partition, made of Plexiglass, was positioned in the middle of the enclosure. The visualizations of the velocity and temperature fields were obtained by using respectively, Plexiglass and liquid crystal particles as tracers. A middle plane perpendicular to the partition was numerically modeled. The steady two-dimensional model accounted for the variable thermophysical properties of the fluid. The finite volume method based on the finite difference approach was applied. The convective terms were approximated using a deferred correction central difference scheme. The velocity and temperature fields and the distribution of the local and average Nusselt numbers were found as a function of the Rayleigh (38 000 <Ra <369 000) and Prandtl (2700 < Pr < 7000) numbers.     

Natural convection flow in a square cavity revisited: Laminar and turbulent models with wall functions

Numerical simulations have been undertaken for the benchmark problem of natural convection flow in a square cavity. The control volume method is used to solve the conservation equations for laminar and turbulent flows for a series of Rayleigh numbers (Ra) reaching values up to 10¹⁰. The k-? model has been used for turbulence modelling with and without logarithmic wall functions. Uniform and non-uniform (stretched) grids have been employed with increasing density to guarantee accurate solutions, especially near the walls for high Ra-values. ADI and SIP solvers are implemented to accelerate convergence. Excellent agreement is obtained with previous numerical solutions, while some discrepancies with others for high Ra-values may be due to a possibly different implementation of the wall functions. Comparisons with experimental data for heat transfer (Nusselt number) clearly demonstrates the limitations of the standard k-? model with logarithmic wall functions, which gives significant overpredictions.     

Lattice Boltzmann simulation of natural convection in open end cavity with inclined hot wall

Natural convection in an open end cavity with a hot inclined wall is simulated based on the lattice Boltzmann method (LBM). The physics of flow and energy transfer in open end cavities are addressed when the hot wall is inclined. The combination of the two topics (open cavity and inclined walls) is the main novelty of the present study. The effects of the angle of the hot inclined wall on the flow field and heat transfer are thoroughly investigated. The Prandtl number is fixed to 0.71 (air). The Rayleigh number and the angle of the hot inclined wall are varied in the range of 10⁴ to 10⁶ and 60? to 85?, respectively. The results are presented for two different aspect ratios, i.e., A = 1 and 2. The results obtained with the LBM are also compared with those of the finite volume method (FVM). The predicted results of the LBM conform to those of the FVM. The results show that by increasing the angle of the hot inclined wall and the aspect ratio of the cavity, the average Nusselt number decreases. The trend of the local Nusselt number on the inclined wall is also discussed.     

Double diffusive convection of anomalous density fluids in a porous cavity

A numerical study has been performed to analyze the combined effect of temperature and species gradients on the buoyancy-driven natural convection flow of cold water near its density extremum contained in a porous cavity. The governing equations are descretized using the finite volume method. The results of the investigation are presented in the form of steady-state streamlines, velocity vectors, isotherms, and isoconcentrationlines. The results are discussed for different porosities, Darcy numbers, and Grashof numbers. The heat and mass transfer rates calculated are found to behave nonlinearly with hot wall temperature. The heat and mass transfer are increased with increasing Darcy number and porosity. It is found that the convective heat and mass transfer rate are greatly affected by the presence of density maximum.     

Non-linear approximations for natural convection in a horizontal annulus

Two different problems are proposed as approximations of the usual system modelling natural convection under the Oberbeck-Boussinesq assumptions. The error is evaluated by means of the norm of its gradient in the Hilbert space. The average Nusselt number is also estimated.     

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

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

Melting process with solid-liquid density change and natural convection in a rectangular cavity

This paper presents a numerical method that simulates the melting process in the presence of solid-liquid density change and natural convection in the melt. The physical model concerned is two-dimensional melting of a phase-change material, initially at its fusion temperature, charged in a rectangular cavity with isothermally heated side walls and an adiabatic bottom wall. The presence of the density change brings no change into the basic form of governing equation, so it is considered through the reformulation of boundary conditions. Difficulties associated with the complex time-dependent melt region, whose shape is also a part of the solutions, are overcome by employing the boundary-fitted coordinate system. Comparison with other works validates the present numerical model and reveals the effects of density change qualitatively. Also, it is confirmed that the present method is preferable to others with natural convection only. Computed results for interesting cases are shown in forms of transient position of the interface, temperature distribution, flow pattern, heat transfer coefficient, and melting fraction as a function of time. Closer examination on melting patterns allows a correlation to be made between the melting fraction and a new independent variable Ste·Fo·Ra^1/4.     

Mass and heat transfer by natural convection in a vertical cavity

This paper reports a fundamental study of laminar natural convection in a rectangular enclosure with heat and mass transfer from the side, when the bouyancy effect is due to density variations caused by either temperature or concentration variations. In the first part of the study scale analysis is used to determine the scales of the flow, temperature and concentration fields in boundary layer flow for all values of Prandtl and Lewis numbers. In particular, scale analysis shows that in the extreme case where the flow is driven by bouyancy due to temperature variations, the ratio of mass transfer rate divided by heat transfer rate scales as $\text{Le}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$ only if (Pr > 1, Le < 1) or (Pr < 1, Sc < 1), and as $\text{Le}{}^{\mathrm{<mtext>1</mtext><mtext>3</mtext>}}$ if (Pr > 1, Le > 1) or (Pr < 1, Sc > 1). In the second part of the study, the boundary layer scales derived in the first part are used to determine the heat and mass transport characteristics of a vertical slot filled with fluid. Criteria for the existence of distinct thermal and concentration boundary layers in the slot are determined. Numerical solutions for the flow and concentration fields in a slot without distinct thermal boundary layers are reported. These solutions support further the method of scale analysis employed in the first part of the study