Finite Element Simulation on Natural Convection Flow in a Triangular Enclosure Due to Uniform and Nonuniform Bottom Heating

2008 ◽  
Vol 130 (3) ◽  
Author(s):  
S. Roy ◽  
Tanmay Basak ◽  
Ch. Thirumalesha ◽  
Ch. Murali Krishna

A penalty finite element analysis with biquadratic elements has been carried out to investigate natural convection flows within an isosceles triangular enclosure with an aspect ratio of 0.5. Two cases of thermal boundary conditions are considered with uniform and nonuniform heating of bottom wall. The numerical solution of the problem is illustrated for Rayleigh numbers (Ra), 103⩽Ra⩽105 and Prandtl numbers (Pr), 0.026⩽Pr⩽1000. In general, the intensity of circulation is found to be larger for nonuniform heating at a specific Pr and Ra. Multiple circulation cells are found to occur at the central and corner regimes of the bottom wall for a small Prandtl number regime (Pr=0.026−0.07). As a result, the oscillatory distribution of the local Nusselt number or heat transfer rate is seen. In contrast, the intensity of primary circulation is found to be stronger, and secondary circulation is completely absent for a high Prandtl number regime (Pr=0.7–1000). Based on overall heat transfer rates, it is found that the average Nusselt number for the bottom wall is 2 times that of the inclined wall, which is well, matched in two cases, verifying the thermal equilibrium of the system. The correlations are proposed for the average Nusselt number in terms of the Rayleigh number for a convection dominant region with higher Prandtl numbers (Pr=0.7 and 10).

1970 ◽  
Vol 39 (1) ◽  
pp. 30-39 ◽  
Author(s):  
Tamanna Sultana ◽  
Sumon Saha ◽  
Goutam Saha ◽  
Md Quamrul Islam

A numerical study of natural convection in a tilted square cavity with heated horizontal base and cold upper wall is presented. The present study is based in such a configuration where the top triangular wall of two different shapes is maintained at a constant low temperature. A constant heat flux source whose length is 20% of the total length of the cavity is discretely embedded at the left corner of the bottom wall. The remaining part of the bottom wall and the two sidewalls are considered to be adiabatic. The study includes computations for inclination angles of the cavity from 0° to 45°, where the Grashof number, Gr varies from 103 to 106. The Penalty finite element method has been used to see the effects of inclination angles and Grashof number on heat transfer process in the cavity. Results are presented in the form of streamline and isotherm plots as well as the variation of the average Nusselt number. Observation shows the significant effect of different triangular top surface on the heat transfer characteristics at the higher Grashof number and inclination angle. Keywords: Natural convection, Penalty finite element, Nusselt number, Isoflux heating. doi:10.3329/jme.v39i1.1831 Journal of Mechanical Engineering, vol. ME39, No. 1, June 2008 30-39  


2012 ◽  
Vol 16 (1) ◽  
pp. 79-91 ◽  
Author(s):  
Varun Sharma ◽  
Kumar Dhiman

In this work, effects of Prandtl number on the heat transfer characteristics of an unconfined rotating circular cylinder are investigated for varying rotation rate (? = 0 - 5) in the Reynolds number range 1 - 35 and Prandtl numbers range 0.7 - 100 in the steady flow regime. The numerical calculations are carried out by using a finite volume method based commercial CFD solver FLUENT. The isotherm patterns are presented for varying values of Prandtl number and rotation rate in the steady regime. The variation of the local and the average Nusselt numbers with Reynolds number, Prandtl number and rotation rate are presented for the above range of conditions. The average Nusselt number is found to decrease with increasing value of the rotation rate for the fixed value of the Reynolds and Prandtl numbers. With increasing value of the Prandtl number, the average Nusselt number increases for the fixed value of the rotation rate and the Reynolds number; however, the larger values of the Prandtl numbers show a large reduction in the value of the average Nusselt number with increasing rotation rate.


2020 ◽  
pp. 190-190
Author(s):  
Houssem Laidoudi ◽  
Mustapha Helmaoui ◽  
Mohamed Bouzit ◽  
Abdellah Ghenaim

In this paper, we performed a numerical simulation of natural convection of Newtonian fluids between two cylinders of different cross-sectional form. The inner cylinder is supposed to be hot and the outer cylinder is assumed to be cold. The diameter of inner cylinder to the diameter of outer cylinder defines the radii ratio (RR= 2.5). The governing equations describing the physical behavior of fluid flow and heat transfer are solved using finite volume method. The effects of Prandtl number (Pr = 0.71 to 100), Rayleigh number (Ra = 103 to 105) and inclination angle of inner cylinder (? = 0? to 80?) on streamlines, isotherms and dimensionless velocity are presented and discussed. Also, the mean average Nusselt number of inner cylinder is plotted versus the governing parameters. All present simulations are considered in two-dimensions for steady laminar flow regime. The obtained results showed that the flow between cylinders is more stable for the inclination angle ? = 0?. Increase in Rayleigh number increases the heat transfer rate for all values of inclination angle. Furthermore, the effect of Prandtl number on the mean average Nusselt number becomes negligible when Pr is over the value 7.01. For example at Pr = 0.71 and Ra =105, increase in inclination from 0? to 40?decreases the average Nusselt number by 5.4%. A new correlation is also provided to describe the average Nusselt number as function of Pr and Ra at ? = 0?.


2021 ◽  
Vol 25 (01) ◽  
pp. 131-142
Author(s):  
Hassan H. Hatem ◽  
◽  
Luma F. Ali ◽  

Steady conjugate natural convection heat transfer in a two-dimensional enclosure filled with fluid saturated porous medium is studied numerically. The two vertical boundaries of the enclosure are kept isothermally at same temperature, the horizontal upper wall is adiabatic, and the horizontal lower wall is partially heated. The Darcy extended Brinkman Forcheimer model is used as the momentum equation and Ansys Fluent software is utilized to solve the governing equations. Rayleigh number (1.38 ≤ Ra ≤ 2.32), Darcy number (3.9 * 10-8), the ratio of conjugate wall thickness to its height (0.025 ≤ W ≤ 0.1) , heater length to the bottom wall ratio (1/4 ≤  ≤ 3/4) and inclination angle (0°, 30° and 60°) are the main considered parameters. The presented results show the effect of these parameters on the heat transfer and fluid flow characteristics. These results include streamlines, isotherm patterns, and local and average Nusselt number for different values of the governing parameters. It is found that either increasing the Rayleigh number and the ratio of conjugate wall thickness to its height (d/H) or decreasing the ratio of heat source width to bottom wall (l/L ), the average Nusselt number is increased. Also, it was observed that the average Nusselt number does not change substantially with inclination angle.


1970 ◽  
Vol 39 (1) ◽  
pp. 1-7 ◽  
Author(s):  
Sumon Saha ◽  
Noman Hasan ◽  
Chowdhury Md Feroz

A numerical study has been carried out for laminar natural convection heat transfer within a two-dimensional modified square enclosure having a triangular roof. The vertical sidewalls are differentially heated considering a constant flux heat source strip is flush mounted with the left wall. The opposite wall is considered isothermal having a temperature of the surrounding fluid. The rest of the walls are adiabatic. Air is considered as the fluid inside the enclosure. The solution has been carried out on the basis of finite element analysis by a non-linear parametric solver to examine the heat transfer and fluid flow characteristics. Different heights of the triangular roof have been considered for the present analysis. Fluid flow fields and isotherm patterns and the average Nusselt number are presented for the Rayleigh numbers ranging from 103 to 106 in order to show the effects of these governing parameters. The average Nusselt number computed for the case of isoflux heating is also compared with the case of isothermal heating as available in the literature. The outcome of the present investigation shows that the convective phenomenon is greatly influenced by the inclined roof height. Keywords: Natural convection, triangular roof, Rayleigh number, isoflux heating. Doi:10.3329/jme.v39i1.1826 Journal of Mechanical Engineering, vol. ME39, No. 1, June 2008 1-7


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.


2019 ◽  
Vol 29 (11) ◽  
pp. 4130-4141 ◽  
Author(s):  
Abdulmajeed Mohamad ◽  
Mikhail A. Sheremet ◽  
Jan Taler ◽  
Paweł Ocłoń

Purpose Natural convection in differentially heated enclosures has been extensively investigated due to its importance in many industrial applications and has been used as a benchmark solution for testing numerical schemes. However, most of the published works considered uniform heating and cooling of the vertical boundaries. This paper aims to examine non-uniform heating and cooling of the mentioned boundaries. The mentioned case is very common in many electronic cooling devices, thermal storage systems, energy managements in buildings, material processing, etc. Design/methodology/approach Four cases are considered, the left-hand wall’s temperature linearly decreases along the wall, while the right-hand wall’s temperature is kept at a constant, cold temperature. In the second case, the left-hand wall’s temperature linearly increases along the wall, while the right-hand wall’s temperature is kept a constant, cold temperature. The third case, the left-hand wall’s temperature linearly decreases along the wall, while the right-hand wall’s temperature linearly increases along the wall. In the fourth case, the left-hand and the right-hand walls’ temperatures decrease along the wall, symmetry condition. Hence, four scenarios of natural convection in enclosures were covered. Findings It has been found that the average Nusselt number of the mentioned cases is less than the average Nusselt number of the uniformly heated and cooled enclosure, which reflects the physics of the problem. The work quantifies the deficiency in the rate of the heat transfer. Interestingly one of the mentioned cases showed two counter-rotating horizontal circulations. Such a flow structure can be considered for passively, highly controlled mechanism for species mixing processes application. Originality/value Previous works assumed that the vertical boundary is subjected to a constant temperature or to a sinusoidal varying temperature. The subject of the work is to examine the effect of non-uniformly heating and/or cooling vertical boundaries on the rate of heat transfer and flow structure for natural convection in a square enclosure. The temperature either linearly increases or decreases along the vertical coordinate at the boundary. Four scenarios are explored.


Author(s):  
T. W. Lin ◽  
M. C. Wu ◽  
C. H. Peng ◽  
P. L. Chen ◽  
Y. H. Hung

In the present study, an experimental setup with stringent measurement methods for performing the natural convection from a horizontal heated chip mounted with a silicon heat spreader coated with diamond film has been successfully established. The parametric studies on the local and average effective heat transfer characteristics for natural convection from a horizontal smooth silicon wafer, rough silicon wafer or silicon wafer coated with diamond film spreader have been explored. The influencing parameters and conditions include Grashof number and spreader material with different surface treatment conditions. From the results, an axisymmetric bowl-shaped Nu profile is achieved, and the highest heat transfer performance occurs at the location near the rim of the heated surfaces for various heat spreaders. The local Nusselt number for a specified convective heat flux decreases along the distance from the disk rim toward the center. The local or average Nusselt number increases with increasing Grashof number for various heat spreaders. As compared with the average Nusselt number for smooth water surface (Ra=5.69nm), the heat transfer enhancements for rough silicon surface (Ra=516.61nm) and rough diamond surface (Ra=319.51nm) are 10.42% and 7.69%, respectively. Furthermore, new correlations for local and average Nusselt numbers for various heat spreaders are presented, respectively. As compared with the smooth silicon surface, the external thermal resistance for rough silicon surface and rough diamond surface are reduced to 91.18% and 90.73%, respectively; and the maximum thermal resistances for rough silicon wafer and silicon wafer coated with diamond film are reduced to 90.43% and 92.61%, respectively.


Author(s):  
Dean Ferley ◽  
Scott J. Ormiston

Numerical analysis of steady, two-dimensional, laminar forced convection in corrugated-plate channels is performed using a commercial CFD code: ANSYS CFX. The flow domain consists of six modules in each of three wall corrugations: sinusoidal-wavy-shaped (SWS), rounded-ellipse-shaped (RES), and rounded-vee-shaped (RVS). One ratio of minimum-to-maximum plate spacings and one module length-to-height ratio is considered. Fluid flow and heat transfer are repeating in the modules and the results are examined in a typical module in the fully-developed region for Reynolds numbers in the range of 25 to 300 for Prandtl numbers of 0.7 (air), 2.29 (water), and 34.6 (ethylene glycol). The RES corrugation produced the highest peak value of local Nusselt number as well as the highest friction factor. The SWS corrugation produced the highest average Nusselt number, except at a Prandtl number of 34.6 at higher Reynolds number where the RES corrugation had the highest value. The RVS corrugation had the lowest friction factor for the geometric configuration considered. The highest heat transfer rate per unit pumping power was found at the highest Prandtl number for the RES corrugation.


1981 ◽  
Vol 103 (2) ◽  
pp. 226-231 ◽  
Author(s):  
G. S. Shiralkar ◽  
C. L. Tien

Heat transfer by natural convection in a horizontal cavity with adiabatic horizontal walls and isothermal side walls is investigated numerically for high aspect ratios (width/height). Comparison is made with existing analytical and experimental results. Agreement is generally good at moderate and high Prandtl numbers to which most previous works have been restricted. Improvements of the existing correlation have been proposed in regions of discrepancy. Extension to the low Prandtl number case, including the range of liquid metals, has been made on the basis of an analytical model for high Rayleigh numbers as well as by numerical solution of the full equations. The agreement between the two is found to be very good. A correlation for the heat transfer is proposed for each of the two different cases of high and low Prandtl number.


Sign in / Sign up

Export Citation Format

Share Document