Numerical Investigation of Unsteady Natural Convection in an Inclined Square Enclosure With Heat-Generating Porous Medium

The unsteady laminar natural convection in an inclined square enclosure with heat-generating porous medium whose heat varies by a cosine function is investigated by a thermal equilibrium model and the Brinkman–Darcy–Forchheimer model numerically, with the four cooled walls of closure as isothermal. The numerical code based on the finite-volume method has been validated by reference data before it was adopted. Influence of dimensionless frequency and inclination angle on heat transfer characteristics in a square enclosure, such as flow distribution, isotherm, averaged Nusselt number on each wall, and time-averaged Nusselt number, are discussed, with specified value for Rayleigh number = 10⁸, Darcy number = 10^?4, Prandtl number = 7, porosity = 0.4, and specific heat ratio = 1. It is found that when the internal heat source varies by cosine, the Nusselt numbers of the four walls oscillate with the same frequency as the internal heat source; however, phase difference occurs. Moreover, frequency has little impact on time-averaged Nusselt number of the four walls, which is different from the phenomenon discovered in natural convection with suitable periodic varying wall temperature boundary condition. Moreover, inclination angle plays an important role in the heat transfer characteristics of the walls studied.     

Numerical Studies on Natural Convection in a Trapezoidal Enclosure With Discrete Heating

Abstract

A steady state laminar natural convection flow in a trapezoidal enclosure with discretely heated bottom wall, adiabatic top wall, and constant temperature cold inclined walls is performed. The finite volume based commercial code “ANSYS-FLUENT” is used to investigate the influence of discrete heating on natural convection flows in a trapezoidal cavity. The numerical solution of the problem covers various Rayleigh numbers ranging from 10³ to 10⁶, non-dimensional heating length ranging from 0.2 to 0.8 and Prandtl number is 0.7. The performance of the present numerical approach is represented in the form of streamfunction, temperature profile and Nusselt number. Heat transfer increases with increase of Rayleigh numbers at the corners of the cavity for same heating length from center of the bottom wall. However, the heat transfer rate is less and almost constant for the Rayleigh numbers considered. It is found that the average Nusselt number monotonically increases with increase of Rayleigh number and length of heat source. The variation of local and average Nusselt numbers is more significant for larger length of heating than smaller one. The heat transfer correlations useful for practical design problems have been predicted.     

Investigation of buoyancy-driven flow and heat transfer in a trapezoidal cavity filled with water–Cu nanofluid

A numerical study is conducted to investigate the transport mechanism of free convection in a trapezoidal enclosure filled with water–Cu nanofluid. The horizontal walls of the enclosure are insulated while the inclined walls are kept at constant but different temperatures. The numerical approach is based on the finite element technique with Galerkin's weighted residual simulation. Solutions are obtained for a wide range of the aspect ratio (AR) and Prandtl number (Pr) with Rayleigh number (Ra = 10⁵) and solid volume fraction (? = 0.05). The streamlines, isotherm plots and the variation of the average Nusselt number at the left hot wall are presented and discussed. It is found that both AR and Pr affect the fluid flow and heat transfer in the enclosure. A correlation is also developed graphically for the average Nusselt number as a function of the Prandtl number as well as the cavity aspect ratio.     

Conjugate natural convection in the roof cavity of heavy construction building during summer

This paper presents the results of a study of conjugate natural convection inside a building attic in the shape of a rectangular enclosure bounded by realistic walls made from composite construction materials under summer day boundary conditions. The effects of cavity aspect ratio, Rayleigh number (Ra), and orientation of the external surfaces on the flow and heat transfer characteristics were the main focus of the investigation. The problem was formulated in terms of the vorticity-stream function procedure, and the governing equations for steady, laminar, two-dimensional conjugate natural convection heat transfer were solved by employing the Alternating Direction Implicit (ADI) control volume method along with under-relaxation factors for temperature, vorticity, and stream functions. For Ra ranging from 10³ to 10¹⁰, steady state results of the streamline and temperature contours in addition to local and mean Nusselt number at all surfaces of the cavity were obtained. The results show that the values of Ra and the aspect ratio have significant effect on the temperature and stream function contours within the enclosure. Another important finding of the study is that heat flux into the room increases with the increase of both the aspect ratio and Rayleigh number.     

A Differential Quadrature Solution of Natural Convection in an Enclosure with a Partial Partition

Natural convection in a partially divided enclosure has been examined numerically using a differential quadrature method. Governing equations for the flow and heat transfer have been constructed by the vorticity-stream function formulation and computational results have been obtained for the Rayleigh numbers, 10⁴, 10⁵, 10⁶, the locations of the partition, 0.2, 0.4, 0.6, 0.8, and the partition ratios, 0.25, 0.50, 0.75, 1.0. The results show that circulation strength, and therefore heat transfer, decreases considerably with increasing partition height especially for the higher values of the Rayleigh number. As the distance of the partition from the hot wall increases, first a decrease and then an increase is observed in the average Nusselt number for the low Rayleigh numbers and first an increase and then a decrease turns out for the high Rayleigh numbers. As for the average Nusselt number on the partition, first it shows a decreasing trend and then an increasing trend as the partition is distanced from the hot wall towards the cold wall. The increase in the partition height brings about significant increase in the average Nusselt number on the partition.     

Three-Dimensional Numerical Study of Natural Convective Heat Transfer of Liquid in a Cubic Enclosure

ABSTRACT

Steady-state laminar natural convection in a cubic enclosure with a cold vertical wall and two hot square heaters with constant temperature on the opposite wall is studied numerically. The enclosure is filled with various liquids. Three-dimensional Navier–Stokes Equations are solved by employing the SIMPLE algorithm. Computations are performed for a range of Rayleigh number from 10³ to 10⁷ while enclosure aspect ratio varies from 0.05 to 1.6. The effects of Rayleigh number, enclosure aspect ratio, and Prandtl number on heat transfer characteristics are studied in detail. The results show that the flow field is very complex and heat transfer from the two heaters is not the same. The effects of Prandtl number are negligible in the range from 5 to 140 with other parameters kept constant. This allows the use of liquids such as water for studying other dielectric liquids, provided the flow geometry and other nondimensional parameters are similar. The overall Nusselt number increases markedly with Rayleigh number. It is also affected by enclosure aspect ratio. It attains the maximum value when aspect ratio is in the range of 0.1–0.2 and decreases as enclosure aspect ratio varies from 0.2 to 1.6. Also, various settings of cooling face and arrangement of heaters are investigated, and the results show that they have considerable effects on heat transfer of both heaters.     

Effect of aspect ratio on natural convection in an inclined rectangular enclosure with sinusoidal boundary condition

The aim of the present numerical study is to understand the natural convection flow and heat transfer in an inclined rectangular enclosure with sinusoidal temperature profile on the left wall. The top and bottom walls of the enclosure are kept to be adiabatic. The finite difference method is used to solve the governing equations with a range of inclination angles, aspect ratios and Rayleigh numbers. The results are presented in the form of streamlines, isotherms and Nusselt numbers. The heat transfer increases first then decreases with increasing the inclination of the enclosure for all aspect ratio and Rayleigh number. Increasing the aspect ratio shows a decreasing trend of the heat transfer for all Rayleigh numbers considered. A correlation equation is also introduced for the heat transfer analysis in this study.     

Free Convection in Confined Power-Law Fluids From Two Side-by-Side Cylinders in a Square Enclosure

Laminar free convection heat transfer in power-law fluids from two side-by-side cylinders (one hot and one cold) confined in a square duct has been studied numerically in the two-dimensional flow regime. For a fixed value of the ratio of cylinder radius to size of the enclosure, the effect of geometrical placement of the cylinders is studied on the resulting velocity and temperature fields in the laminar free convection regime by considering six asymmetric locations of the two cylinders. In particular, extensive results reported herein span the range of conditions of Grashof number, 10 to 10⁵; Prandtl number, 0.7 to 100, thereby yielding the range of the Rayleigh number as 7 to 10⁷; power-law index, 0.3 to 1.8; and the relative positions (dimensionless) of the cylinders with respect to the centerline, –0.25 to 0.25. The heat transfer characteristics are analyzed in terms of the local Nusselt number along the surfaces of the two cylinders and the enclosure walls. Overall, the average Nusselt number shows a positive dependence on both the Grashof number and the Prandtl number irrespective of the values of power-law index and relative positioning of the cylinders. Also, all else being equal, shear-thinning fluid behavior promotes heat transfer with reference to that in Newtonian fluids. When the two cylinders are situated close to the bottom wall, the rate of heat transfer is augmented with reference to that for the symmetric positioning of the cylinders along the horizontal mid-plane of the enclosure. Conversely, heat transfer deteriorates as the cylinders are located above the centerline of the enclosure. The present numerical results have been consolidated via the use of a modified Rayleigh number, thereby enabling the estimation of the average Nusselt number in a new application.     

Lattice Boltzmann Simulation of Natural Convection in a Differentially Heated Square Enclosure Containing Heat Generating Low‐Prandtl Number Fluid

Lattice Boltzmann simulations were conducted for the free convective flow of a low‐Prandtl number (Pr = 0.0321) fluid with internal heat generation in a square enclosure having adiabatic top and bottom walls and isothermal side walls. The problem of free convection with volumetric heat source has represented itself in connection with advanced engineering applications, such as water‐cooled lithium–lead breeder blankets for nuclear fusion reactors and liquid metal sources of spallation neutrons for subcritical fission systems. A single relaxation time (SRT) thermal lattice Boltzmann method (LBM) was employed. While applying SRT, a D2Q9 model was used to simulate the flow field and temperature field. Results have been obtained for various Rayleigh numbers characterizing internal and external heating from 10³ to 10⁶. Flow and temperature fields in terms of stream function and isotherms in the enclosure were predicted for these cases. The temperature of the fluid in the enclosure was found higher than the heated wall temperature at high values of internal Rayleigh numbers. The internal heat generation affected the rate of heat transfer significantly as two convection loops are observed in the enclosure. The average Nusselt number at the heated and cold wall was determined for all the cases.     

WALL HEAT CONDUCTION EFFECT ON NATURAL CONVECTION IN AN ENCLOSURE FILLED WITH A NON-DARCIAN POROUS MEDIUM

ABSTRACT

A numerical analysis has been made of the conjugate natural convection in a rectangular enclosure filled with a fluid-saturated porous medium and surrounded with four solid walls. The conductance of the walls is assumed to be much greater than that of the cavity filled with a porous medium. The main objective was to investigate the influences of the ratio of thermal conductivity of the wall to that of the fluid-porous matrix composite, the Darcy-modified Rayleigh number, the Prandtl number, and the aspect ratio. The streamlines and isotherms are presented; also, the local and average Nusselt numbers are presented along the interface between walls and cavity. A non-Darcian model was employed and the numerical method was SIMPLE-C. The numerical results indicate that the wall heat conduction effects decrease the heat transfer rate. When the wall heat conduction is considered, the greater the conductance of the solid walls surrounding the cavity, the greater is the rate of heat transfer.     

Study of buoyancy-induced flows subjected to partially heated sources on the left and bottom walls in a square enclosure

In this paper, numerical simulations of laminar, steady, two-dimensional natural convection flows in a square enclosure with discrete heat sources on the left and bottom walls are presented using a finite-volume method. Two different orientated wall boundary conditions are designed to investigate the natural convection features. The computational results are expressed in the form of streamlines and isothermal lines for Rayleigh numbers ranging from 10² to 10⁷ in the cavity. In the course of study, a combination of third-order and exponential interpolating profile based on the convective boundedness criterion is proposed and tested against the partially heated cavity flow up to the highest Rayleigh number 10⁷. The effects of thermal strength and heating length on the hydrodynamic and thermal fields inside the enclosure are also presented. Numerical results indicate that the average Nusselt number increases as Rayleigh number increases for both cases. Moreover, it is seen that the effect of the heat transfer rate due to the heating strength on the left wall is different from the one on the bottom. For the heater size effect, it is observed that by increasing the length of heat source segment, the heat transfer rate is gradually increased for both cases.     

Free convective heat transfer in an enclosure with an internal louvered blind

A numerical study has been conducted of free convection in a tall vertical enclosure with an internal louvered metal blind. The study considers the effects of Rayleigh number, enclosure aspect ratio, and blind geometry on the convective heat transfer. The numerical model has been validated against experimental measurements and the results have been presented in terms of an empirical correlation for the average Nusselt number. The correlation is applicable to an enclosure with an internal metal blind. It has been shown that the Nusselt number correlation can be combined with a simple one-dimensional model to closely predict the enclosure U-value.     

Combined Effect of Mixed Convection and Surface Radiation Heat Transfer for Thermally Developing Flow in Vertical Channels

Combined effect of laminar flow mixed convection and surface radiation heat transfer for thermally developing airflow in a vertical channel heated from a side has been experimentally examined with different thermal and geometric parameters. The channel boundary is made of two isothermal walls and two adiabatic walls, the isothermal parallel wall is heated uniformly and the opposite cold wall temperature is maintained equal to the inlet conditions. The heated wall temperature ranged from 55 to 100°C, Reynolds number ranged from 800 to 2900 and the heat flux was varied from 250 to 870 W/m². To cover the wide range of Reynolds numbers, two aspect ratios of square and rectangular section were used. Surface radiation from the internal walls is considered through two emissivities i.e. 0.05 and 0.85, to represent weak and strong radiation effects, respectively. From the experiments, surface temperature and Nusselt number distributions of convection and radiation heat transfer are obtained for different heat flux values. Flow structure inside the channel is visualized to observe the flow pattern. The results show the combined effect of laminar flow mixed convection and surface radiation on the total heat transfer rate within the channel. The accumulating buoyancy force and airflow moves together vertically in the upward direction to give significant heat transfer enhancement in the vertical orientation of the channel.     

Natural convective heat transfer in square enclosures heated from below

The paper deals with the results of an experimental and numerical study of free convective heat transfer in a square enclosure characterized by a discrete heater located on the lower wall and cooled from the lateral walls.The study analysed how the heat transfer develops inside the cavity at the increasing of the heat source length.The experimental data are obtained by measuring the temperature distribution in the air layer by real-time and double-exposure holographic interferometry while the commercial finite volumes code Fluent 6.0 is used for the numerical study. Convection has been studied for Rayleigh number from 10³ to 10⁶. Different convective forms are obtained depending on Ra and on the heat source length.The local Nusselt number is evaluated on the heat source surface and it shows a symmetrical form raising near the heat source borders. Graphs of the local Nusselt number on the heat source and of the average Nusselt number at several Ra are finally presented.     

Free convection in a shallow wavy enclosure

This paper presents the simulation of free convection heat transfer and fluid flow in a horizontal and shallow wavy enclosure. Its bottom wall is varied with a sinusoidal function while the top wall is flat. The wavy wall has constant hot temperature and ceiling is cold as verticals are insulated. The simulations are performed using CFDRC software. Results are presented in terms of streamlines, isotherms and Nusselt number for different aspect ratios, non-dimensional wave length and Rayleigh number. The obtained results showed that heat transfer is increased with the decreasing non-dimensional wave length and it is increased with the increasing of aspect ratio and Rayleigh number.     

Conjugate turbulent natural convection in the roof enclosure of a heavy construction building during winter

This paper presents the results of a study of conjugate turbulent natural convection inside a building attic in the shape of a rectangular enclosure bounded by realistic walls made from composite construction materials under winter day boundary conditions. The effects of cavity aspect ratio, Rayleigh number (Ra), depth of the external concrete beam, and external wall construction materials on the flow and heat transfer characteristics were the main focus of the investigation. The Shear stress transport k–ω turbulence model is implemented to calculate air-flow velocities and temperatures in a steady, turbulent, two-dimensional conjugate natural convection heat transfer inside an attic. The governing equations were solved by employing the line-by-line tri-diagonal matrix algorithm (TDMA) control volume method. For Ra ranging from 10⁷ to 10¹⁰, steady-state results of the streamline and temperature contours in addition to local and mean Nusselt numbers at all surfaces of the cavity were obtained. The results show that the values of Ra, attic aspect ratio and the composite wall materials have significant effect on the temperature and stream function contours within the enclosure, and the heat flux out of the room through the enclosure.     

Natural convection of nanofluids in enclosures with low aspect ratios

This paper analyzes heat transfer and fluid flow of natural convection in inclined cavity filled with CuO-water nanofluid heated from one side and cooled from the ceiling. The transport equations for the flow are solved numerically by the finite volume element method using the SIMPLER algorithm Based on numerical predictions. The effects of Rayleigh number and aspect ratio on flow pattern and energy transport are investigated for Rayleigh numbers ranging from 10⁴ to 10⁷ volume fraction of solid varied to 0%–4% and for five different aspect ratios of 0.08, 0.1, 0.125, 0.25 and 0.5. It is found that the effect of Rayleigh number on heat transfer is less significant when the enclosure is shallow (AR = 0.5) and the influence of aspect ratio is stronger when the enclosure is tall and the Rayleigh number is high.     

Effects of thermal boundary conditions on natural convection flows within a square cavity

A numerical study to investigate the steady laminar natural convection flow in a square cavity with uniformly and non-uniformly heated bottom wall, and adiabatic top wall maintaining constant temperature of cold vertical walls has been performed. A penalty finite element method with bi-quadratic rectangular elements has been used to solve the governing mass, momentum and energy equations. The numerical procedure adopted in the present study yields consistent performance over a wide range of parameters (Rayleigh number Ra, 10³ ⩽ Ra ⩽ 10⁵ and Prandtl number Pr, 0.7 ⩽ Pr ⩽ 10) with respect to continuous and discontinuous Dirichlet boundary conditions. Non-uniform heating of the bottom wall produces greater heat transfer rates at the center of the bottom wall than the uniform heating case for all Rayleigh numbers; however, average Nusselt numbers show overall lower heat transfer rates for the non-uniform heating case. Critical Rayleigh numbers for conduction dominant heat transfer cases have been obtained and for convection dominated regimes, power law correlations between average Nusselt number and Rayleigh numbers are presented.     

Lattice Boltzmann simulation of natural convection in tall enclosures using water/SiO2 nanofluid

Natural convection in enclosures using water/SiO₂ nanofluid is simulated with Lattice Boltzmann method (LBM). This investigation compared with other numerical methods and found to be in excellent agreement. This study has been carried out for the pertinent parameters in the following ranges: the Rayleigh number of base fluid, Ra = 10³-10⁵, the volumetric fraction of nanoparticles between 0 and 4% and aspect ratio (A) of the enclosure between 0.5 and 2. The thermal conductivity of nanofluids is obtained on basis of experimental data. The comparisons show that the average Nusselt number increases with volume fraction for the whole range of Rayleigh numbers and aspect ratios. Also the effect of nanoparticles on heat transfer augments as the enclosure aspect ratio increases.     

Three-dimensional natural convection in an enclosure with a sphere at different vertical locations

Numerical calculations are carried out for the three-dimensional natural convection induced by a temperature difference between a cold outer cubic enclosure and a hot inner sphere. The immersed-boundary method (IBM) to model a sphere based on the finite volume method is used to study a three-dimensional natural convection for different Rayleigh numbers varying in the range of 10³–10⁶. This study investigates the effect of the inner sphere location on the heat transfer and fluid flow. The flow and thermal fields eventually reach the steady state for all Rayleigh numbers regardless of the sphere location. For Rayleigh numbers of 10⁵ and 10⁶, the variation of local Nusselt number of the sphere along the circumferential direction is large, showing the strong three dimensionality of the natural convection in the enclosure unlike to the cases of lower Rayleigh numbers of 10³ and 10⁴. For the highest Rayleigh number, the local peaks of the Nusselt number on the top wall of the enclosure shows the sinusoidal distribution along the circumferential direction. The flow and thermal fields, and the local and surface-averaged Nusselt numbers on the sphere and the enclosure are highlighted in detail.