scholarly journals A RANS K-? SIMULATION OF 2D TURBULENT NATURAL CONVECTION IN AN ENCLOSURE WITH HEATING SOURCES

2019 ◽  
Vol 20 (1) ◽  
pp. 229-244
Author(s):  
Mehdi Ahmadi ◽  
Seyed Ali Agha Mirjalily ◽  
Seyed Amir Abbas Oloomi

ABSTRACT: This study is conducted to investigate turbulent natural convection flow in an enclosure with thermal sources using the low-Reynolds number (LRN) k-? model. This enclosure has a cold source with temperature Tc and a hot source with temperature Th as thermal sources, other walls of the enclosure are adiabatic. The aim of this study is to predict the effect of change in Rayleigh number, repositioning of cold and hot sources, and thermal sources aspect ratio on the flow field, temperature, and rate of heat transfer. To achieve this aim, the equations of continuity, momentum, energy, turbulent kinetic energy, and kinetic energy dissipation are employed in the case of 2D turbulence with constant thermo-physical properties except the density in the buoyancy term (Boussinesq approximation). To numerically solve these equations, the finite volume method and SIMPLE algorithm are used. According to the modeling results, the most optimal temperature distribution in the enclosure is seen when the hot source is below the cold source. With decreasing distance between hot and cold sources, heat transfer rate increases. The maximal heat transfer rate is derived via study of the heating sources aspect ratio. In constant positions of cold and hot sources on a wall, the heat transfer rate increases with increasing Rayleigh number (Ra=109-1011). ABSTAK: Kajian ini dijalankan bagi mengkaji perubahan semula jadi aliran perolakan dalam tempat tertutup dengan sumber haba menggunakan model k-? nombor Reynolds-rendah (LRN). Bekas tertutup ini mempunyai dua sumber haba iaitu sumber sejuk dengan suhu Tc dan sumber panas dengan suhu Th, manakala dinding lain bekas ini adalah adiabatik. Tujuan kajian ini adalah bagi mengesan perubahan nombor Rayleigh, mengubah sumber sejuk dan panas dan nisbah sumber haba kepada kawasan aliran, suhu dan halaju perubahan haba. Bagi mencapai tujuan tersebut, persamaan sambungan, momentum, tenaga, tenaga kinetik perolakan, dan pengurangan tenaga kinetik telah dilaksanakan dalam kes perolakan 2D dengan sifat fizikal-haba berterusan (malar) kecuali isipadu terma keapungan (anggaran Boussinesq). Bagi menyelesaikan persamaan ini secara berangka, kaedah isipadu terhad dan algorithma MUDAH telah digunakan. Berdasarkan keputusan model, suhu distribusi optimal dalam bekas tertutup dilihat apabila sumber panas adalah kurang daripada sumber sejuk. Dengan pengurangan jarak antara sumber panas dan sejuk, kadar pertukaran haba meningkat. Kadar pertukaran haba maksima telah diperoleh melalui kajian nisbah  aspek sumber pemanasan. Kadar pertukaran haba bertambah dengan bertambahnya nombor Rayleigh  (Ra=109-1011), pada posisi tetap sumber sejuk dan panas pada dinding bekas.

2020 ◽  
Vol 330 ◽  
pp. 01004
Author(s):  
Abdennacer Belazizia ◽  
Smail Benissaad ◽  
Said Abboudi

Steady, laminar, natural convection flow in a square enclosure with partially active vertical wall is considered. The enclosure is filled with air and subjected to horizontal temperature gradient. Finite volume method is used to solve the dimensionless governing equations. The physical problem depends on three parameters: Rayleigh number (Ra =103-106), Prandtl number (Pr=0.71), and the aspect ratio of the enclosure (A=1). The active location takes two positions in the left wall: top (T) and middle (M). The main focus of the study is on examining the effect of Rayleigh number on fluid flow and heat transfer rate. The results including the streamlines, isotherm patterns, flow velocity and the average Nusselt number for different values of Ra. The obtained results show that the increase of Ra leads to enhance heat transfer rate. The fluid particles move with greater velocity for higher thermal Rayleigh number. Also by moving the active location from the top to the middle on the left vertical wall, convection and heat transfer rate are more important in case (M). Furthermore for high Rayleigh number (Ra=106), Convection mechanism in (T) case is principally in the top of the enclosure, whereas in the remaining case it covers the entire enclosure.


1990 ◽  
Vol 112 (3) ◽  
pp. 648-652 ◽  
Author(s):  
K. S. Chen ◽  
A. C. Ku ◽  
C. H. Chou

Experimental results are presented for steady natural convection in a two-dimensional, partially divided, rectangular enclosure, in which two of the vertical walls were maintained at different uniform temperatures and the top and bottom walls were insulated. The partition plate was adiabatic, and the experiment was carried out both with and without an opening in the partition. Rayleigh numbers ranging from 106 to 108 and opening ratios of 0, 1/8, and 1/4 were investigated for an enclosure aspect ratio (length/height) of 2 and Prandtl number of 7 (for water). Local velocity and temperature measurements were made with a laser-Doppler velocimeter and thermocouple probes. Flow visualization using colored dye was also performed. Results show that there was a recirculation zone in the upper and left quadrant of the enclosure when there was no opening in the partition plate. With an opening in the partition, the recirculation zone was absent and the heat transfer rate increased. An unopened partial obstruction would reduce the heat transfer rate by an amount of 12 to 30 percent depending on the Rayleigh number. However, the opening seems to have little effect on the velocity and temperature profiles of the left-moving fluid on the bottom wall. A correlation of the Nusselt number is derived, which shows that the heat transfer rate increases as the Rayleigh number or opening ratio increases.


Author(s):  
Ridha Jmai ◽  
Brahim Ben Beya ◽  
Taieb Lili

Natural convection in a rectangular cavity with aspect ratio (Ax), partially heated and filled with a nanofluid (Cu-Water) has been studied numerically. Two heat sources with length (B) are placed on the opposite vertical walls; the remainder of the walls is maintained adiabatic while the horizontal walls are brought to a cold temperature. The equations governing the flow are solved using a finite volume home code using a multigrid technique. Among the parameters governing the flow, a detailed study on the effects of the aspect ratio (Ax) and the length of the source (B) on flow and heat transfer rate is given. The results are shown in terms of streamlines and isotherms. It was found that the transfer of heat significantly increases with the aspect ratio (Ax) and the length of the source (B). A correlation expressing the Nusselt number as a function of (Ax) and d is established.


2020 ◽  
Vol 17 (2) ◽  
pp. 89-99
Author(s):  
Houssem Laidoudi

The governing equations of continuity, momentum and energy are numerically solved to study the laminar natural convection heat transfer of Newtonian fluid confined within two concentric cylinders. The inner cylinder is elliptical cross-section with different aspect ratio E = 0.1 to 0.5 and it is considered to be hot, whereas the outer cylinder is circular and it is supposed to be cold.    The annular spacing between the cylinders is defined based on radii ratio (RR = 2.5). Also, the inner cylinder is inclined with an inclination angle (θ = 0 to 90). The main purpose of this study is to determine the effects of inclination angle (θ = 0° to 90°), aspect ratio of inner cylinder (E = 0.1 to 0.5), Prandtl number (Pr = 0.71 and 7.01) and Rayleigh number (Ra = 103 to 105) on fluid flow and heat transfer rate. The flow patterns and temperature distributions are potted in terms of streamlines and isotherms respectively. The obtained results showed that increase in inclination angle enhances the heat transfer rate of inner cylinder for all values of aspect ratio. Also, for the inclination angle          (θ = 90°), the decrease in aspect ratio (E) improves the heat transfer rate of inner cylinder.


1999 ◽  
Vol 122 (1) ◽  
pp. 176-179
Author(s):  
Q. W. Wang ◽  
W. Q. Tao ◽  
Z. Lin ◽  
T. T. Chow

Natural convection in a cubic inclined enclosure with three isolated plates was investigated experimentally. The influences of the plates’ spacing, the inclination angles of the enclosure, and the Rayleigh number on the heat transfer of the plate group were obtained. It was found that under a fixed Rayleigh number, there is a plate spacing at which the heat transfer rate of the three plates is approximately the same for the horizontal plate group. Under the range of Rayleigh numbers considered, the heat transfer rate of the plate group is less than that of the natural convection of the plate group in infinite space. The heat transfer rate of the plates increases with the plate inclined angle tilting from θ=90 deg to θ=0 deg, with the most steep increase occurring in the range of θ=90 deg to 45 deg. A global correlation of the heat transfer results for all the inclination between 0 deg and 60 deg can be obtained as Nul,m=0.5360Ralcos θ0.25 with a spread of ±8.9 percent. [S0022-1481(00)00501-6]


Energies ◽  
2019 ◽  
Vol 12 (14) ◽  
pp. 2807 ◽  
Author(s):  
Hamed Bagheri ◽  
Mohammadali Behrang ◽  
Ehsanolah Assareh ◽  
Mohsen Izadi ◽  
Mikhail A. Sheremet

In the present investigation, the free convection energy transport was studied in a C-shaped tilted chamber with the inclination angle α that was filled with the MWCNT (MultiWall Carbon Nanotubes)-Fe3O4-H2O hybrid nanofluid and it is affected by the magnetic field and thermal flux. The control equations were numerically resolved by the finite element method (FEM). Then, using the artificial neural network (ANN) combined with the particles swarm optimization algorithm (PSO), the Nusselt number was predicted, followed by investigating the effect of parameters including the Rayleigh number (Ra), the Hartmann number (Ha), the nanoparticles concentration (φ), the inclination angle of the chamber (α), and the aspect ratio (AR) on the heat transfer rate. The results showed the high accuracy of the ANN optimized by the PSO algorithm in the prediction of the Nusselt number such that the mean squared error in the ANN model is 0.35, while in the ANN model, it was optimized using the PSO algorithm (ANN-PSO) is 0.22, suggesting the higher accuracy of the latter. It was also found that, among the studied parameters with an effect on the heat transfer rate, the Rayleigh number and aspect ratio have the greatest impact on the thermal transmission intensification. The obtained data also showed that a growth of the Hartmann number illustrates a reduction of the Nusselt number for high Rayleigh numbers and the heat transfer rate is almost constant for low Rayleigh number values.


Author(s):  
Basma Souayeh ◽  
Nader Ben-Cheikh ◽  
Brahim Ben-Beya

Purpose The purpose of this paper is to examine numerically the three natural convection of air induced by temperature difference between a cold outer cubic enclosure and a hot inner cylinder. Simulations have been carried out for Rayleigh numbers ranging from 103 to 107 and titled angle of the enclosure from 0° to 90°. The developed mathematical model is governed by the coupled equations of continuity, momentum and energy, and is solved by finite volume method. The effects of cylinder inclination and Rayleigh number on fluid flow and heat transfer are presented. The distribution of isocontours of temperature and isosurfaces of velocity eventually reaches a steady state in the range of Rayleigh numbers between 103 and 107 for titled inclination of 90°; however, for the remaining inclinations, Rayleigh number must be in the range 103-106 to avoid unsteady state, which is manifested by the division of the area containing the maximum local heat transfer rate into three parts for a Rayleigh number equal to 107 and an inclination of 90°. We mention that instability study is not included in the present paper, which is solely devoted to three-dimensional calculations. Results also indicate that optimal average heat transfer rate is obtained for both high Rayleigh number of 106 and high inclination of 90° for the two cases of the inner cylinder and cubical enclosure. Design/methodology/approach The manuscript deals with prediction of the three-dimensional natural convection phenomena in a cubical cavity induced by an isothermal cylinder at the center with different inclinations by simulating the flow using highly numerical methods such as finite volume method. Findings It is found that the local Nusselt number through active walls for titled inclination set at 90°, the symmetry of the flow is conserved and the area containing the maximum heat transfer is divided into three smaller areas situated near the upper portion of the wall, taking the maximum value. That may be due to the preparation of local occurrence of instabilities and bifurcation phenomena that appear for Ra > 107, which is not included in the present paper to save journal space. It was found also that an optimal heat transfer appears when the cylinder orientation becomes vertical (a = 90°). For this inclination, buoyancy forces act upward, corresponding to an aiding situation. In addition, heat transfer rate is increasing with Rayleigh numbers, so correlations of average Nusselt through the cubical cavity and the cylinder are established as function of two parameters (Ra, a). Comparisons of the numerical results with those obtained from all correlations show good agreements. Originality/value To the author’s knowledge, studies have thus far adressed three-dimensional cuboids enclosures induced by an inner shape which the location is changed. However, no study has examined three-dimensional natural convection between the inner isothermal cylinder and outer cooled cubical enclosure when the outer enclosure is tilted.


2021 ◽  
Vol 13 (8) ◽  
pp. 168781402110391
Author(s):  
Ben Abdelmlek Khaoula ◽  
Ben Nejma Fayçal

This paper deals with a numerical study of mixed convection heat transfer in horizontal eccentric annulus. The inner cylinder is supposed hot and rotating, however the outer one is kept cold and motionless. The numerical problem was solved using COMSOL Multiphysics® which is based on finite element method. The resolution of the partial differential equations was conducted through an implicit scheme with the use of the damped Newton’s method. The present numerical analysis concerns the effect of eccentricity, rotation speed and Rayleigh number on the flow patterns, heat transfer rate, and energy efficiency of the process. It was found that the heat transfer rate increases with the increase of Rayleigh number. In addition, the heat transfer rate drops with the increase of rotation speed. Finally, we have demonstrated that maximum energy efficiency is achieved not only with higher Rayleigh number but also it is maximum with small eccentricity.


2009 ◽  
Vol 131 (7) ◽  
Author(s):  
Ashok Kumar ◽  
P. Bera

A comprehensive numerical investigation on the natural convection in a hydrodynamically anisotropic porous enclosure is presented. The flow is due to nonuniformly heated bottom wall and maintenance of constant temperature at cold vertical walls along with adiabatic top wall. Brinkman-extended non-Darcy model, including material derivative, is considered. The principal direction of the permeability tensor has been taken oblique to the gravity vector. The spectral element method has been adopted to solve numerically the governing conservative equations of mass, momentum, and energy by using a stream-function vorticity formulation. Special attention is given to understand the effect of anisotropic parameters on the heat transfer rate as well as flow configurations. The numerical experiments show that in the case of isotropic porous enclosure, the maximum rates of bottom as well as side heat transfers (Nub and Nus) take place at the aspect ratio, A, of the enclosure equal to 1, which is, in general, not true in the case of anisotropic porous enclosures. The flow in the enclosure is governed by two different types of convective cells: rotating (i) clockwise and (ii) anticlockwise. Based on the value of media permeability as well as orientation angle, in the anisotropic case, one of the cells will dominate the other. In contrast to isotropic porous media, enhancement of flow convection in the anisotropic porous enclosure does not mean increasing the side heat transfer rate always. Furthermore, the results show that anisotropy causes significant changes in the bottom as well as side average Nusselt numbers. In particular, the present analysis shows that permeability orientation angle has a significant effect on the flow dynamics and temperature profile and consequently on the heat transfer rates.


Sign in / Sign up

Export Citation Format

Share Document