ISSN: 2319-9873

Reach Us
+44 7456 035580

^{1}Assistant Professor, Department of Mechanical Engineering, Coimbatore Institute of Technology, Coimbatore – 641 014, Tamilnadu, India

^{2}Post Graduate Student, Heat Power Engineering, Coimbatore Institute of Technology, Coimbatore – 641 014, Tamilnadu, India

- *Corresponding Author:
- Balamurugan S

Assistant Professor

Department of Mechanical Engineering

Coimbatore Institute of Technology

Coimbatore – 641 014, Tamilnadu, India

**Tel:**+91-(0)-9865062725

**E-mail:**[email protected]

**Received:** 27 June 2015 **Accepted:** 13 September 2015 **Published:** 25 September 2015

**Visit for more related articles at** Research & Reviews: Journal of Engineering and Technology

In this paper, natural convection around a square and triangle bar kept in a square enclosure has been studied. Aim of this research has analysis, the heat transfer based on natural convection by using air as a medium. The problems taken for study were a varying square enclosure with different aspect ratio of sources of both square and triangle bar. Numerical simulations and boundary conditions are carried out for analysis these systems by using Fluent Software. The square and triangle bar are considered as a hot source and the square enclosure of vertical walls are cold surfaces and horizontal walls are insulated. The size of the square enclosure are taken 20 mm, 40 mm and 80 mm corresponding to each enclosure were taken as the aspect ratio of both square and triangle bar of 0.2, 0.3 and 0.4 respectively. The heat transfer results are obtained from simulation with the help of temperature distribution, Nusselt Number and flow stream function. From this analysis the heat transfer rate was compared each other and concludes that the higher heat transfer rate in square source than triangle source varying enclosures and different aspect ratio.

Natural convection, Numerical simulation, Comparison.

The curiosity of natural convection in enclosures filled with air having cold vertical walls and adiabatic horizontal walls has been theme of research over the past years. The reason of considering this geometry is, it has application in various fields such as building and thermal insulation systems [1-7], solar engineering applications [8,9], geophysical fluid mechanics, etc...

The study of numerically the natural convection around tilted square cylinders in the range of (0° ≤ θ ≤ 45°) inside an enclosure having horizontal adiabatic wall and cold vertical wall is studied. The investigated of the two dimensional natural convective flow and heat transfer around a heated cavity kept in a square enclosure in the range of 103≤ Ra ≤106 [10-12]. Oosthuizen [13] has studied numerically the natural convective air flow in an enclosure with a horizontal lower wall, vertical side-walls and a straight inclined top wall. Sun [14] have investigated laminar natural convection heat transfer from a horizontal triangular cylinder to its concentric cylindrical enclosure. Fan [15] studied the effect of Prandtl number on the heat transfer in a horizontal cylindrical enclosure with a coaxial triangular cylinder inside it. Hussain and Hussein [16] investigated numerically the natural convection in a uniformly heated circular cylinder at different locations inside a square enclosure.

In the present work, the effect of size of a square and triangle bar in the square enclosure on the flow and heat transfer rate is compared for both the square and triangle sources according to the aspect ratio.

The work is done systematically in a sequential manner to fulfill the research objective. Based on the objective of this study, literature survey was done. Available modeling techniques have been used for this purpose. The applied modeling techniques were of two types, mathematical and 2D modeling. The obtained results were keenly analyzed. Depending on the analysis output, results were elaborated and the final conclusion is provided (**Figure 1**).

**Figures** **8** and **9** shows the geometry in which natural convective heat transfer is studied in the present work. Natural convection is the phenomenon where the movement of the fluid is characterised by the density changes. Here the natural convection is modelled with a heat source of square and triangle bar is placed inside an enclosure. Air is present between each of the heat source and the enclosure. Given the conditions of the different enclosure walls, the heat transfer is modelled. The aim is to study the best method that ensures the best heat transfer rate. The procedure is to be carried out with a different shape of the heat source and the heat transfer rate is studied. The problem is carried out in two-dimensional form.

Initially, the heat tranfer was studied for heated square bar domain and heated triangle bar domain by applying hot and cold sources (**Figures** **2** and **3**). Then the respective surface meshes of fluid around these bars were checked (**Figures** **4** and **5**).

The heat transfer was studied for the square enclosure by using different walls types, boundaries and conditions, which are represented in the following **Table 1**.

A details simulation study was conducted for the square bar containing various lengths (L) and aspect ratios (A). These different length parameters and associated aspect ratios that were used for the square bar is shown individually as follows:

The simulation with length L= 20 mm and Aspect ratio (A) of 0.2 was studied and the results is provided in **Figures** **6** and **7**, where **Figure** **6** shows the contours of stream function and **Figure** **7** shows velocity of fluid flow around square bar following the condition of A 0.2 and L 20 mm.

The simulation study with length L= 20 mm and Aspect ratio (A) of 0.3 and the observed results is shown in **Figures** **8** and **9**, where **Figure** **8** depicts the contours of stream function and **Figure** **9** represents velocity of fluid flow around square bar maintaining a value of A 0.2 and L 20 mm.

In continuation, similar simulation study was conducted with length L= 20 mm and Aspect ratio (A) = 0.4. The results are shown in **Figures** **10** and **11**, where **Figure** **10** shows the contours of stream function and **Figure** **11** shows velocity of fluid flow around square bar of A 0.4 and L 20 mm.

Another parameter considered for the simulation was length L= 40 mm and Aspect ratio (A) = 0.2 and the results are shown in **Figures** **12** and **13**, where **Figure** **12** shows the contours of stream function and **Figure** **13** shows velocity of fluid flow around square bar of A 0.2 and L 40 mm.

The simulation with length L= 40 mm and Aspect ratio (A) = 0.3 was conducted whereas the secured results is shown in **Figures** **14** and **15**, where **Figure** **14** shows the contours of stream function and **Figure** **15** shows velocity of fluid flow around square bar.

The simulation with length L= 40 mm and Aspect ratio (A) = 0.4 was conducted whereas the secured results were shown in **Figures** **16** and **17**, where **Figure** **16** shows the contours of stream function and **Figure** **17** shows velocity of fluid flow around square bar of A 0.4 and L 40 mm.

The simulation with length L= 80 mm and Aspect ratio (A) = 0.2 were conducted whereas the secured results were shown in **Figures** **18** and **19**, where **Figure** **18** shows the contours of stream function and **Figure** **19** shows velocity of fluid flow around square bar of A 0.2 and L 80 mm.

The simulation with length L= 80 mm and Aspect ratio (A) = 0.3 were conducted and the secured results shown in **Figures** **20** and **21**, where **Figure** **20** shows the contours of stream function and **Figure** **21** shows velocity of fluid flow around square bar of A 0.3 and L 80 mm.

The simulation with length L= 80 mm and Aspect ratio (A) = 0.4 were conduct6ed and the secured results were shown in **Figures** **22** and **23**, where **Figure** **22** shows the contours of stream function and **Figure** **23** shows velocity of fluid flow around square bar of A 0.4 and L 80 mm.

The simulations for a heated triangle bar was studied with various lengths and aspect ratios. The results that were obtained, were shown individually as follows:

The simulation of heated triangle bar of length L = 20 mm and Aspect ratio (A) = 0.2 were conducted and the secured results were shown in **Figures** **23** and **24**. **Figure** **23** shows the contours of stream function and **Figure** **24** shows the velocity of fluid around the triangle bar.

The simulation of heated triangle bar of length L = 20 mm and Aspect ratio (A) = 0.3 were conducted and the secured results were shown in **Figures** **25**-**27**. **Figure** **26** shows the contours of stream function and **Figure** **27** shows the velocity of fluid around the triangle bar.

The simulation of heated triangle bar of length L = 20 mm and Aspect ratio (A) = 0.4 were studied and the secured results were shown in **Figures** **28** and **29**. **Figure** **28** shows the contours of stream function and **Figure** **29** shows the velocity of fluid around the triangle bar.

The simulation of heated triangle bar of length L = 40 mm and Aspect ratio (A) = 0.2 were conducted and the secured results were shown in **Figures** **30** and **31**. **Figure** **30** shows the contours of stream function and **Figure** **31** shows the velocity of fluid around the triangle bar.

The simulation of heated triangle bar of length L = 20 mm and Aspect ratio (A) = 0.2 were conducted and the secured results were shown in **Figures** **32** and **33**. **Figure** **32** shows the contours of stream function and **Figure** **33** shows the velocity of fluid around the triangle bar.

The simulation of heated triangle bar of length L = 40 mm and Aspect ratio (A) = 0.4 were conducted and the secured results were shown in **Figures** **34** and **35**. **Figure** **34** shows the contours of stream function and **Figure** **35** shows the velocity of fluid around the triangle bar.

The simulation of heated triangle bar of length L = 80 mm and Aspect ratio (A) = 0.2 were conducted and the secured results were shown in **Figures** **36** and **37**. **Figure** **36** shows the contours of stream function and **Figure** **37** shows the velocity of fluid around the triangle bar.

The simulation of heated triangle bar of length L = 80 mm and Aspect ratio (A) = 0.3 were conducted and the secured results were shown in **Figures** **38** and **39**. **Figure** **38** shows the contours of stream function and **Figure** **39** shows the velocity of fluid around the triangle bar.

The simulation of heated triangle bar of length L = 80 mm and Aspect ratio (A) = 0.4 were conducted and the secured results were shown in **Figures** **40** and **41**. **Figure** **40** shows the contours of stream function and **Figure** **41** shows the velocity of fluid around the triangle bar.

From the above work, the results that were obtained is shown in **Tables** **2** and **3**. The obtained results for a square bar of varying lengths and aspect ratios are represeted in an understandable manner. From the above results obtained, variations in Nusselt no. and enclosure left wall were estimated.

The Nusselt number (Nu) is the ratio of convective to conductive heat transfer across the boundary. In this context, convection includes both advection and diffusion. It is calculated for a heat transfer across a boundary within a fluid. The conductive component is measured under the same conditions as the heat convection but with a (hypothetically) stagnant (or motionless) fluid.

The variations of Nusselt number and enclosure side wall of square and triangle bars of various lengths were estimated as shown in following **Figures** **42**-**56**. The results obtained for a triangle bar with varying lengths and aspect ratios can be shown in **Table 4**. For determining the variation of Nusselt number with enclosure of wall was determined for both square and triangular bars by using various lengths of 20, 40 and 80 mm.

**1. For the square bar enclosure of length L=20 mm (Figures 43-45):**

**2. For the square bar enclosure of length L=40 mm (Figures 46-49):**

**3. For the square bar enclosure of length L=80 mm (Figures 48-50) (Table 4):**

**1. For the triangle bar enclosure of length L=20mm (Figures 51 and 52):**

**2. For the triangle bar enclosure of length L=40mm (Figure 53 and 54):**

**3. For the triangle bar enclosure of length L=80mm (Figure 55 and 56):**

Heat transfer and fluid flow due to natural convection in air around heated square cylinders of different sizes inside an enclosure having adiabatic horizontal and diathermic vertical walls of size 20 mm, 40 mm and 80 mm with respect to aspect ratio of 0.2, 0.3 and 0.4 are analyzed, and results are exhibited in the dimension form of Stream function and Velocity vector diagram.

The fluid motion and circulation rate increase with increase in enclosure size. The fluid motion is almost uniform for lower for smaller enclosure size and the fluid motion is prominent near the walls and the fluid is almost stagnant in the core region for higher larger enclosure size. The heat transfer rate is comparatively higher at the upper portions of the vertical enclosure walls and from the base of the triangle bar. For smaller enclosures, the circulation rate decreases with increase in bar size, but for larger enclosure, the circulation rate increases with bar size.

It’s proved that Nusselt number dependent of temperature and heated source bar can absorb more energy by increasing the size and can transfer more heat to the diathermic walls for various purposes.

- Akinsete VA and Coleman TA. Heat transfer by steady laminar free convection in triangular enclosures, Int. J. Heat Mass Transfer 1982; 25: 991–998.
- Holtzman GA, et al. Laminar natural convection in isosceles triangular enclosures heated from below and symmetrically cooled from above, J. Heat Transfer – Trans. ASME, 2000; 485–491.
- Haese PM andTeubner MD. Heat exchange in an attic space, Int. J. Heat Mass Transfer 2002;45:4925–4936.
- SahaSC. Unsteady natural convection in a triangular enclosure under isothermal heating, Energy Build 2011; 43:695–703.
- Saha SC. Scaling of free convection heat transfer in a triangular cavity for Pr> 1, Energy Build 2011;43: 2908–2917.
- Saha SC, et al. Natural convection in attics subject to instantaneous and ramp cooling boundary conditions, Energy Build 2010; 42:1192–1204.
- Saha SC, et al. Natural convection in attic-shaped spaces subject to sudden and ramp heating boundary conditions, Heat Mass Transfer 2010;46:621–638.
- Joudi KA, et al. Computational model for a prism shaped storage solar collector with a right triangular cross section, Energy Convers. Manage 2004; 45:337–342.
- Kaushik SC, et al. Transient analysis of a triangular built-in-storage solar water-heater under winter conditions, Heat Recovery Syst. 1994; 14:391–409.
- Roychowdhury DG, et al. Numerical simulation of natural convective heat transfer and fluid flow around a heated cylinder inside an enclosure, Heat Mass Transfer 2002;38: 565–576.
- Dalal A, et al. Natural convection around a heated square cylinder placed in different angles inside an enclosure, Heat and mass transfer conference 2008; 45.
- De K and DalalA. A numerical study of natural convection around a square horizontal heated cylinder placed in an enclosure, International Journal of Heat and Mass Transfer 2006;49:4608-4623.
- Oosthuizen PH.Free convective flow in an enclosure with a cooled inclined upper surface, Computational mechanics 1994;14:420-430.
- Xu X, et al. Numerical investigation of laminar natural convective heat transfer from a horizontal triangular cylinder to its concentric cylindrical enclosure, International Journal of Heat and Mass Transfer 2009;52:3176–3186.
- Yu Z, et al. Prandtl number dependence of laminar natural convection heat transfer in a horizontal cylindrical enclosure with an inner coaxial triangular cylinder, International Journal of Heat and Mass Transfer 2010; 53: 1333–1340.
- Hussain S and Hussein K. Numerical investigation of natural convection phenomena in a uniformly heated circular cylinder immersed in square enclosure filled with air at different vertical locations, International Communications in Heat and Mass Transfer 2010;37:1115–1126.