Numerical investigation of Rayleigh-Benard Natural convection and entropy generation in a cubic cavity with discrete heat sources

Author(s):  
Chemseddine Maatki ◽  
Kaouther Ghachem ◽  
Mohammed A. Almeshaal ◽  
Nidhal Ben Khedher ◽  
Lioua Kolsi

Abstract The Three-dimensional natural convection with isothermal discrete heat sources in a cubical cavity has been carefully studied using the 3D vector potential-vorticity formulation. Based on the finite volume method, the governing equations are solved with a home-made computational code (written in Fortran). Assuming that all cavity vertical walls are adiabatic, the upper wall of the cavity is kept at a cold temperature. However, in the bottom face, heat sources are placed under different configurations. The size of the discrete sources, their positions, and their numbers are varied for different Rayleigh numbers. The Prandtl number is fixed at 0.71. Three-dimensional distribution of the temperature iso-surfaces, the heat transfer rate, and entropy generation are evaluated. It is found that heat transfer and entropy generation are strongly affected by the arrangement of the discrete heated sources. In conclusion, the heat transfer rate is maximized, and the entropy generation is minimized for the inline arrangement of more than two heaters compared to the diagonal one.

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 ahead-of-print (ahead-of-print) ◽  
Author(s):  
Mahyar Ashouri ◽  
Mohammad Mehdi Zarei ◽  
Ali Moosavi

Purpose The purpose of this paper is to investigate the effects of geometrical parameters, eccentricity and perforated fins on natural convection heat transfer in a finned horizontal annulus using three-dimensional lattice Boltzmann flux solver. Design/methodology/approach Three-dimensional lattice Boltzmann flux solver is used in the present study for simulating conjugate heat transfer within an annulus. D3Q15 and D3Q7 models are used to solve the fluid flow and temperature field, respectively. The finite volume method is used to discretize mass, momentum and energy equations. The Chapman–Enskog expansion analysis is used to establish the connection between the lattice Boltzmann equation local solution and macroscopic fluxes. To improve the accuracy of the lattice Boltzmann method for curved boundaries, lattice Boltzmann equation local solution at each cell interface is considered to be independent of each other. Findings It is found that the maximum heat transfer rate occurs at low fin spacing especially by increasing the fin height and decreasing the internal-cylindrical distance. The effect of inner cylinder eccentricity is not much considerable (up to 5.2% enhancement) while the impact of fin eccentricity is more remarkable. Negative fin eccentricity further enhances the heat transfer rate compared to a positive fin eccentricity and the maximum heat transfer enhancement of 91.7% is obtained. The influence of using perforated fins is more considerable at low fin spacing although some heat transfer enhancements are observed at higher fin spacing. Originality/value The originality of this paper is to study three-dimensional natural convection in a finned-horizontal annulus using three-dimensional lattice Boltzmann flux solver, as well as to apply symmetry and periodic boundary conditions and to analyze the effect of eccentric annular fins (for the first time for air) and perforated annular fins (for the first time so far) on the heat transfer rate.


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.


Processes ◽  
2021 ◽  
Vol 9 (8) ◽  
pp. 1339
Author(s):  
Yacine Khetib ◽  
Ahmad Aziz Alahmadi ◽  
Ali Alzaed ◽  
Hamidreza Azimy ◽  
Mohsen Sharifpur ◽  
...  

In this paper, the free convective heat transfer of nanofluids in a square cavity is simulated using a numerical method. The angle of the cavity could be changed in the horizontal axis from 0 to 90 degrees. The cavity is exposed under a constant magnetic field. Two opposite walls of the cavity are cold and warm, and the rest of the walls are insulated. On the hot wall, there are two fins with the same wall temperature. The equations were discretized by the finite volume method (FVM) and then solved using the SIMPLE algorithm. Three different fin configurations (straight, inclined and curved) were studied in terms of heat transfer rate and generation of entropy. According to the simulation results, the heat transfer rate was improved by tilting the fins toward the top or bottom of the cavity. At Ra = 105 and Ha = 20, the maximum heat transfer rate was achieved at a cavity inclination of 90° and 45°, respectively, for straight and curved fins. In the horizontal cavity, heat transfer rate could be improved up to 6.4% by tilting the fins and up to 4.9% by warping them. Increasing the Hartmann number from 0 to 40 reduced the Nusselt number and entropy generation by 37.9% and 33.8%, respectively.


Author(s):  
M.A. Mansour ◽  
Sameh Elsayed Ahmed ◽  
Ali J. Chamkha

Purpose This paper aims to investigate the entropy generation due to magnetohydrodynamic natural convection flow and heat transfer in a porous enclosure filled with Cu-water nanofluid in the presence of viscous dissipation effect. The left and right walls of the cavity are thermally insulated. There are heated and cold parts, and these are placed on the bottom and top wall, respectively, whereas the remaining parts are thermally insulated. Design/methodology/approach The finite volume method is used to solve the dimensionless partial differential equations governing the problem. A comparison with previously published woks is presented and is found to be in an excellent agreement. Findings The minimization of entropy generation and local heat transfer according to different values of the governing parameters are presented in details. It is found that the presence of magnetic field has negative effects on the local entropy generation because of heat transfer and the local total entropy generation. Also, the increase in the heated part length leads to a decrease in the local Nusselt number. Originality/value This problem is original, as it has not been considered previously.


2001 ◽  
Author(s):  
Patrick H. Oosthuizen ◽  
Matt Garrett

Abstract Natural convective heat transfer from a wide isothermal plate which has a “wavy” surface, i.e., has a surface which periodically rises and falls, has been numerically studied. The surface waves run parallel to the direction of flow over the surface and have a relatively small amplitude. Two types of wavy surface have been considered here — saw-tooth and sinusoidal. Surfaces of the type considered are approximate models of situations that occur in certain window covering applications, for example, and are also sometimes used to try to enhance the heat transfer rate from the surface. The flow has been assumed to be laminar. Because the surface waves are parallel to the direction of flow, the flow over the surface will be three-dimensional. Fluid properties have been assumed constant except for the density change with temperature that gives rise to the buoyancy forces, this being treated by means of the Boussinesq type approximation. The governing equations have been written in dimensionless form, the height of the surface being used as the characteristic length scale and the temperature difference between the surface temperature and the temperature of the fluid far from the plate being used as the characteristic temperature. The dimensionless equations have been solved using a finite-element method. Although the flow is three-dimensional because the surface waves are all assumed to have the same shape, the flow over each surface thus being the same, and it was only necessary to solve for the flow over one of the surface waves. The solution has the following parameters: the Grashof number based on the height, the Prandtl number, the dimensionless amplitude of the surface waviness, the dimensionless pitch of the surface waviness, and the form of the surface waviness (saw-tooth or sinusoidal). Results have been obtained for a Prandtl number of 0.7 for Grashof numbers up to 106. The effects of Grashof number, dimensionless amplitude and dimensionless pitch on the mean heat transfer rate have been studied. It is convenient to introduce two mean heat transfer rates, one based on the total surface area and the other based on the projected frontal area of the surface. A comparison of the values of these quantities gives a measure of the effectiveness of the surface waviness in increasing the mean heat transfer rate. The results show that while surface waviness increases the heat transfer rate based on the frontal area, the modifications of the flow produced by the surface waves are such that the increase in heat transfer rate is less than the increase in surface area.


2016 ◽  
Vol 831 ◽  
pp. 83-91
Author(s):  
Lahoucine Belarche ◽  
Btissam Abourida

The three-dimensional numerical study of natural convection in a cubical enclosure, discretely heated, was carried out in this study. Two heating square sections, similar to the integrated electronic components, are placed on the vertical wall of the enclosure. The imposed heating fluxes vary sinusoidally with time, in phase and in opposition of phase. The temperature of the opposite vertical wall is maintained at a cold uniform temperature and the other walls are adiabatic. The governing equations are solved using Control volume method by SIMPLEC algorithm. The sections dimension ε = D / H and the Rayleigh number Ra were fixed respectively at 0,35 and 106. The average heat transfer and the maximum temperature on the active portions will be examined for a given set of the governing parameters, namely the amplitude of the variable temperatures a and their period τp. The obtained results show significant changes in terms of heat transfer, by proper choice of the heating mode and the governing parameters.


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.


2013 ◽  
Vol 465-466 ◽  
pp. 500-504 ◽  
Author(s):  
Shahrin Hisham Amirnordin ◽  
Hissein Didane Djamal ◽  
Mohd Norani Mansor ◽  
Amir Khalid ◽  
Md Seri Suzairin ◽  
...  

This paper presents the effect of the changes in fin geometry on pressure drop and heat transfer characteristics of louvered fin heat exchanger numerically. Three dimensional simulation using ANSYS Fluent have been conducted for six different configurations at Reynolds number ranging from 200 to 1000 based on louver pitch. The performance of this system has been evaluated by calculating pressure drop and heat transfer coefficient. The result shows that, the fin pitch and the louver pitch have a very considerable effect on pressure drop as well as heat transfer rate. It is observed that increasing the fin pitch will relatively result in an increase in heat transfer rate but at the same time, the pressure drop will decrease. On the other hand, low pressure drop and low heat transfer rate will be obtained when the louver pitch is increased. Final result shows a good agreement between experimental and numerical results of the louvered fin which is about 12%. This indicates the capability of louvered fin in enhancing the performance of heat exchangers.


Sign in / Sign up

Export Citation Format

Share Document