Conjugate turbulent natural convection with surface radiation in air filled rectangular enclosures

Conjugate turbulent natural convection and surface radiation in rectangular enclosures heated from below and cooled from other walls, typically encountered in Liquid Metal Fast Breeder Reactor (LMFBR) subsystems, have been investigated by a finite volume method for various aspect ratios. The formulation comprises the standard two equation k–ε turbulence model with physical boundary conditions (no wall functions), along with the Boussinesq approximation, for the flow and heat transfer. As far as radiation is concerned, the radiosity – irradiation formulation for a transparent fluid of Prandtl number 0.7 has been employed. The conjugate coupling on the walls has been handled by using a fin type formulation. The Rayleigh number based on the width of the enclosure is varied from 10⁸ to 10¹² and the aspect ratio is varied from 0.5 to 2.0. Detailed results including stream lines, temperature profiles, and convective, radiative and overall Nusselt numbers are presented. A correlation for the mean convection Nusselt number in terms of Rayleigh number and aspect ratio is proposed for design purposes. The influence of the wall emissivity and the external heat transfer coefficient on the heat transfer from the enclosure has also been investigated.     

Conjugate Wall Conduction-Fluid Natural Convection in a Three-Dimensional Inclined Enclosure

Laminar conjugate conduction-natural convection heat transfer in a 3-D inclined cubic enclosure comprised of finite thickness conductive walls and central cavity is numerically investigated. The dimensionless governing equations describing the convective flow and wall heat conduction are solved by the high accuracy multidomain pseudospectral method. Computations are performed for different Rayleigh numbers (10³ ≤ Ra* ≤ 10⁶), thermal conductivity ratios (1 ≤ k ≤ 100), dimensionless wall thickness (0 ≤ s ≤ 0.25), and enclosure inclinations (?30° ≤ α ₁ ≤ 30°, 0° ≤ α ₂ ≤ 45°). The effects of the above controlling parameters on the heat transfer performances of the enclosure system are investigated in detail, with emphases on the variations of wall conduction and fluid convection heat transfer, and the interactive heat transfer conditions between solid walls and fluid in the central cavity. Numerical results reveal that the existence of enclosure walls reduces the temperature gradient across the cavity and alters the temperature distribution within the solid walls; thus, the fluid convection is complexly determined by the combined effects of k and s, and is greatly affected by enclosure inclinations at high Rayleigh numbers. Moreover, the temperature distributions and solid-fluid interactive heat transfer conditions are provided for further interpretation and demonstration of the effects of the solid walls.     

Effect of Darcy,fluid rayleigh and heat generation parameters on natural convection in a porous square enclosure: A Brinkman-extended Darcy model

A Pressure-velocity solution for natural convection for fluid saturated heat generating porous medium in a square enclosure is analysed by finite element method. The numerical solutions obtained for wide range of fluid Rayleigh number, Ra_f, Darcy number, Da, and heat generating number, Q_d. The justification for taking these non-dimensional parameters independently is to establish the effect of individual parameters on flow patterns. It has been observed that peak temperature occurs at the top central part and weaker velocity prevails near the vertical walls of the enclosure due to the heat generation parameter alone. On comparison, the modified Rayleigh number used by the earlier investigators[4,6], can not explain explicitly the effect of heat generation parameter on natural convection within an enclosure having differentially heated vertical walls. At higher Darcy number, the peak temperature and peak velocity are comparatively more, resulting in better enhancement of heat transfer rate.     

Free convection in porous media filled right-angle triangular enclosures

Steady-state free convection heat transfer in a right-angle triangular enclosure, whose vertical wall insulated and inclined and bottom walls are differentially heated, is performed in this study. The governing equations are obtained using Darcy model. In this study, the governing equations were solved by finite difference method and solution of algebraic equations was made via Successive Under Relaxation method. The effect of aspect ratios ranging from 0.25 to 1.0 and Rayleigh numbers 50 ≤ Ra ≤ 1000 is investigated as governing parameters on heat transfer and flow field. It is observed that heat transfer is increased with the decreasing of aspect ratio and multiple cells are formed at high Rayleigh numbers.     

Natural convection heat transfer by heated partitions within enclosure

In this study, natural convection heat transfer and fluid flow of two heated partitions. within an enclosure have been analysed numerically. The right side wall and the bottom wall of the enclosure were insulated perfectly while the left side wall and top wall were maintained at the same uniform temperature. The partitions were placed on the bottom of the enclosure and their temperatures were kept higher than the non-isolated walls. The effects of position and heights of the partitions on heat transfer and flow field have been investigated. Computations for Rayleigh number in the range of 10⁴ and 10⁶ have been conducted. Using the control volume approach, finite difference equations are obtained with non-staggered grid arrangement, a computer program based on the SIMPLEM algorithm was developed. The finite difference equations were solved iteratively with a line-by-line Thomas algorithm.     

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.     

Influences of Boundary Condition on Laminar Natural Convection of Bingham Fluids in Square Cross-Sectioned Cylindrical Annular Enclosures with Differentially Heated Vertical Walls

Numerical simulations have been carried out to analyze steady-state laminar natural convection of yield stress fluids obeying Bingham model in square cross-sectioned cylindrical annular enclosures with differentially heated vertical walls for both constant wall temperature and constant wall heat flux boundary conditions for active walls. The simulations have been performed under the assumption of axisymmetry for a nominal Rayleigh number range of 10³ to 10⁶ and nominal Prandtl number range of 10 to 10³ for different ratio of internal cylinder radius to cylinder height range of 0.125 to 16. The mean Nusselt number on the inner periphery for the constant wall heat flux configuration has been found to be smaller than that in the case of constant wall temperature configuration for a given set of values of nominal Rayleigh and Prandtl numbers for both Newtonian and Bingham fluid cases. The mean Nusselt number normalized by the corresponding value obtained for pure conductive transport increases with increasing internal radius before approaching the corresponding mean Nusselt number for square enclosures regardless of the boundary conditions. Detailed physical explanations have been provided for the effects of the aforementioned parameters on the mean Nusselt number on the inner periphery. Finally, the new Nusselt number correlations have been proposed for laminar natural convection of both Newtonian and Bingham fluids in square cross-sectioned cylindrical annular enclosures for both constant wall temperature and constant wall heat flux boundary conditions.     

NATURAL CONVECTION IN A DIFFERENTIALLY HEATED SQUARE CAVITY WITH INTERNAL HEAT GENERATION

A high-resolution, finite-difference numerical study is reported on natural convection in a square cavity. The vertical sidewatts of the cavity are differentially heated, and a uniform internal heat generation is also present. Two principal parameters are considered, the internal Rayleigh number Ra_I, which represents the strength of the internal heat generation, and the external Rayleigh number Rag, which denotes the effect due to the differential heating of the side walls. The internal Rayleigh number varies in the range 10¹⁰ Ra_I ≤ 10⁷, while the external Rayleigh number is set at Ra_E = 5 x 10⁷ for most computations. As the relative strength of the internal heat generation increases, the flows near the tap portion of the heated sidewall are directed downward. When the effect of the internal heat generation is dominant, the thermal energy leaves the system for the surroundings over the top portion of the heated wall. Only in the bottom pari of the heated wall is heat transfer directed into the system. These numerical solutions are in qualitative agreement with the available experimental measurements.     

Entropy generation due to conjugate natural convection in enclosures bounded by vertical solid walls with different thicknesses

Entropy generation due to conjugate natural convection heat transfer and fluid flow has been studied inside an enclosure with bounded by two solid massive walls from vertical sides at different thicknesses. Enclosure is differentially heated from vertical walls and horizontal walls are adiabatic. Governing equations which are written in streamfunction-vorticity form solved by finite difference technique for the governing parameters as Rayleigh number, 10³ ≤ Ra ≤ 10⁶, length ratio of solid walls as ₁ (for left vertical wall) and ₂ (for right vertical wall) and thermal conductivity ratio of solid to fluid (k), 1 ≤ k ≤ 10. Entropy generation contours due to fluid friction and heat transfer irreversibility, isotherms, streamlines, Nusselt numbers and velocity profiles were obtained. It is found that entropy generation increases with increasing of thermal conductivity ratio and thicknesses of the walls. Entropy generation due to heat transfer is more significant than that of fluid flow irreversibility for all values of thickness of the solid vertical walls.     

Experimental study of natural convection heat transfer in a semicircular enclosure

An experimental study of natural convection heat transfer in a differentially heated semicircular enclosure was carried out. The flat surface was heated and the radial surface was cooled isothermally. The effects of angle of enclosure inclination on the heat transfer across semicircular regions of several radii were measured for Rayleigh numbers Ra_R ranging from 6.72 × 10⁶ to 2.33 × 10⁸, using water as the working fluid. The angle of inclination varied from −90 degrees to 90 degrees with radii R of 50, 40, and 30 mm. The flow patterns were sketched from the results of a visualization experiment using aluminum powder. The temperature measurements in the enclosure were carried out using liquid crystals and thermocouples. The results indicate that different flow patterns were encountered as the angle of inclination varied, and the heat transfer rate was largely dependent on the flow pattern. In particular, enhanced heat transfer rates can be obtained when plume-like flow occurs along both hot and cold walls in the case of an upward-facing hot wall. Heat transfer for the inclined enclosure can be predicted using the equation for a vertical enclosure presented in this paper. © 1998 Scripta Technica, Inc. Heat Trans Jpn Res, 26(2): 131–142, 1997     

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.     

Lattice Boltzmann simulation of natural convection and heat transfer from multiple heated blocks

This study is aimed to investigate the natural convection heat transfer from discrete heat sources (similar to heated microchips) using Bhatnagar‐Gross‐Krook lattice Boltzmann method via graphics process unit computing. The simulation is carried out separately for three and six heated blocks model for different Rayleigh numbers and fixed Prandtl number, $Pr=0.71$ (air). The uniformly heated blocks are placed at the bottom wall inside a rectangular enclosure. The enclosure is maintained by the cold temperature at its left and right walls. The top and bottom surface is maintained by adiabatic conditions apart from the regions where blocks are attached to the bottom wall. The numerical code is validated with the benchmark heat transfer problem of side‐heated square cavity as well as with an experimental study for one discrete heat source. The rate of heat transfer is presented in terms of the local Nusselt and average Nusselt number for each block. It is found that the heat transfer rate becomes maximized in the leftmost and rightmost blocks due to the adjacent cold walls. It is found that the number of blocks and their positions play a substantial role in determining their collective performance on the heat transfer rate.     

NUMERICAL EVALUATION OF WEAKLY TURBULENT FLOW PATTERNS OF NATURAL CONVECTION IN A SQUARE ENCLOSURE WITH DIFFERENTIALLY HEATED SIDE WALLS

This article presents the results of numerical evaluation of weakly turbulent natural convection of air in a rectangular enclosure with differentially heated side walls and adiabatic horizontal walls. The turbulence in the natural convection was described by k–ε equations, which were solved by Strang splitting, while average thermal and fluid flow fields were described by statistically averaged equations, which were solved by the projection method PmIII. The combined application of projection method and the Strang splitting characterizes the numerical method in this study. Numerical results for Rayleigh number 1.58 × 10⁹ have revealed reasonable agreement with the existing experimental ones, with some discrepancy attributable to the adiabatic boundary conditions on the horizontal walls. The results are also in good agreement with some published numerical results, particularly at higher Rayleigh numbers. However, comparison with the latest experimental data reveals that the turbulent heat flux model is not quite capable of giving satisfactory temperature distribution.     

Natural convective flow in an inclined lid-driven enclosure with a heated thin plate in the middle

Natural convection heat transfer from a heated thin plate located in the middle of a lid-driven inclined square enclosure has been analyzed numerically. Left and right of the cavity are adiabatic, the two horizontal walls have constant temperature lower than the plate’s temperature. The study is formulated in terms of the vorticity-stream function procedure and numerical solution was performed using a fully higher-order compact (FHOC) finite difference scheme on the 9-point 2D stencil. Air was chosen as a working fluid (Pr = 0.71). Two cases are considered depending on the position of heated thin plate (Case I, horizontal position; Case II, vertical position). Governing parameters, which are effective on flow field and temperature distribution, are Rayleigh number values (Ra) ranging from 10³ to 10⁵ and inclination angles γ (0° ? γ < 360°). The fluid flow, heat transfer and heat transport characteristics were illustrated by streamlines, isotherms and Nusselt number (Nu). It is found that fluid flow and temperature fields strongly depend on Rayleigh numbers and inclination angles. Further, for the vertical located position of thin plate heat transfer becomes more enhanced with lower γ at various Rayleigh numbers.     

Lattice Boltzmann Method Based Simulation of Natural Convection in Partially Open Enclosure

A thermal lattice Boltzmann method‐based analysis was performed to numerically investigate the heat transfer by natural convection from an enclosure with a large vertical side opening. The height of the opening was less than the enclosure height and the vertical wall opposite to the opening was maintained at constant temperature. A parametric study was carried out for different values of Rayleigh number (Ra) ranging from 10³ to 10⁵ with air as the working fluid for three opening sizes and three opening locations. The Prandtl number was fixed at 0.71 and the enclosure aspect ratio was also fixed at 2 in all calculations. With Boussinesq approximation, the temperature distribution and stream functions in the enclosure were predicted. The profile of the normal velocity component at the opening location was determined. The opening size affects the stratification and recirculation pattern within the enclosure. The average Nusselt number at the heated wall was determined for all cases. © 2013 Wiley Periodicals, Inc. Heat Trans Asian Res; Published online in Wiley Online Library ( wileyonlinelibrary.com/journal/htj ). DOI 10.1002/htj.21110     

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.     

Thermal transport analysis for natural convection in a porous corrugated rhombic enclosure

A numerical study of fluid flow and heat transfer, applying natural convection is carried out in a porous corrugated rhombic enclosure. A uniform heating source is applied from the bottom boundary wall while the inclined side walls are maintained to a constant cold temperature and the top corrugated wall is retained at insulated condition inside the enclosure. The heat transfer and flow features are presented for a wide spectrum of Rayleigh numbers (Ra), 10⁴ ≤ Ra ≤ 10⁶, and Darcy numbers (Da), 10^?3 ≤ Da ≤ 10^?2. The number of undulations (n) for the top and bottom walls have been varied from 1 to 13 keeping the amplitude of undulation fixed. It is revealed that the characteristics of heat transfer are conceivably modulated by changing the parameter of the undulation number on the enclosure walls, specifically at the bottom and top. The influencing control of n in altering the heat transfer rate is felt maximum on the left wall and minimum for the right wall, and there is a strong interplay between Ra and Da together with n on dictating the heat transfer characteristics. The critical value, where heat transfer rate is observed as maximum is at n = 11 and thereafter the values decrease.     

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.     

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.     

Lattice Boltzmann simulation of natural convection heat transfer in an elliptical-triangular annulus

A numerical study for steady-state, laminar natural convection in a horizontal annulus between a heated triangular inner cylinder and cold elliptical outer cylinder was investigated using lattice Boltzmann method. Both inner and outer surfaces are maintained at the constant temperature and air is the working fluid. Study is carried out for Rayleigh numbers ranging from 1.0 × 10³ to 5.0 × 10⁵. The effects of different aspect ratios and elliptical cylinder orientation were studied at different Rayleigh numbers. The local and average Nusselt numbers and percentage of increment heat transfer rate were presented. The average Nusselt number was correlated. The results show that by decreasing the value of aspect ratio and/or increasing the Rayleigh number, the Nusselt number increases. Also the heat transfer rate increases when the ellipse positioned vertically.