A numerical and experimental study of natural convective heat transfer from an inclined isothermal square cylinder with an exposed top surface

Natural convective heat transfer from an isothermal inclined cylinder with a square cross-section which have an exposed top surface and is, in general, inclined at an angle to the vertical has been numerically and experimentally studied. The cylinder is mounted on a flat adiabatic base plate, the cylinder being normal to the base plate. The numerical solution has been obtained by solving the dimensionless governing equations subject to the boundary conditions using the commercial cfd solver, FLUENT. The flow has been assumed to be symmetrical about the vertical center-plane through the cylinder. Results have only been obtained for Prandtl number of 0.7. Values of inclination angle between 0° and 180° and a wide range of Rayleigh number and the dimensionless cylinder width, W = w/h, have been considered. The effects of Dimensionless widths, Rayleigh numbers, and inclination angles on the mean Nusselt number for the entire cylinder and for the mean Nusselt numbers for the various surfaces that make up the cylinder have been examined. Empirical equations for the heat transfer rates from the entire cylinder have been derived.


List of symbols
The angle of inclination of the cylinder relative to the vertical (°)

Introduction
Natural convective heat transfer from an isothermal inclined cylinder with a square cross-section which have an exposed top surface and is, in general, inclined at an angle to the vertical has been numerically and experimentally studied. Some electrical and electronic component cooling problems can be approximately modeled as involving natural convective heat transfer from an isothermal inclined cylinder with a square cross-section mounted on a flat adiabatic base plate. The cylinder has an exposed top surface which is also isothermal and has the same temperature as that of the vertical heated side surfaces. This flow situation, which is considered here, is shown in Fig. 1. The heat transfer from the surface of the cylinder at various angles of inclination between vertically upwards and vertically downwards, as in Fig. 2, has been numerically and experimentally investigated. One of the main aims of the present work was to determine how the cross-sectional size-to-height ratio of the square cylinder, i.e., w/h, influences the mean heat transfer rate from the cylinder at various angles of inclination. Natural convection from vertical slender circular cylinders has been studied for many years. Reviews of early work were provided by Ede [1], Burmeister [2], and Jaluria [3]. Cebeci [4] gives numerical results for Prandtl numbers from 0.01 to 100, extending earlier work that had been for Pr = l and 0.72. This work was further extended by Lee et al. [5]. Integral method results were obtained by LeFevre and Ede [6]. The effect of Prandtl number was further studied by Crane [7]. There have been several more recent Fig. 1 Flow situation considered studies that have further extended this earlier work typical of these being that of Kimura et al. [8] and Oosthuizen [9,10]. Typical of the experimental studies of natural convective heat transfer from vertical circular cylinders are those of Jarall and Campo [11], Welling et al. [12] and Fukusawa and Iguchi [13]. Cylinders with exposed top inclined at an angle to the vertical have received attention by Oosthuizen [14,15].
Some studies of natural convective heat transfer from inclined cylinders have been undertaken. In very early study such as Koch [16] and Farber and Rennat [17]. Few studies however have been concerned with short inclined cylinders such as Oosthuizen [18], Eckert and Jackson [19], Al-Arabi and Salman [20], Al-Arabi and Khamis [21], Chand and Vir [22], Oosthuizen and Mansingh [23] and Oosthuizen and Paul [24].
Some studies of natural convective heat transfer from inclined square cylinders have been undertaken. An experimental study of natural convective heat transfer from an isothermal vertical short and slender square cylinder to air, where the top and bottom surfaces were insulated, were obtained using a lumped capacitance method by Popiel and Wojtkowiak [25]. Even with their widest square cylinder, the mean Nusselt numbers were higher than given by the correlation equation for the wide flat plate, this being especially true at low Rayleigh numbers. Experiments were carried out by Al-Arabi and Sarhan [26] to measure the average natural convection heat transfer rate from the outside surface of isothermal square cylinders of different length and side-lengths at different angle of inclination between 0°and 90°to the vertical in both laminar and turbulent regions. The heat transfer rates from the side surfaces of the square cylinder where the top and bottom ends are insulated were measured. The heat transfer rate was found to depend on side-length and the inclination angle. They found that due to the presence of edge effects, the heat transfer coefficient always exceeded that which would occur under the same conditions with a vertical plate. They also suggested general correlation equations. Natural convective heat transfer from an isothermal vertical square cylinder with exposed upper surface in air mounted on a flat adiabatic base was studied by Oosthuizen [27]. His solution has values of the parameters w/h between 0.1 and 1, Rayleigh number between 1 9 10 4 and 1 9 10 7 , and Prandtl number of 0.7. His results showed that at the lower values of Rayleigh number and w/h considered the mean Nusselt number increases with decreasing w/h, while at larger values of Rayleigh number and w/h the Nusselt number is essentially independent of w/h. Also the mean Nusselt number for the heated top surface was shown to be much lower than that for the vertical side surfaces. He provide a correlation equations for the top and the side surfaces.
Heat transfer from cylinders with a square or rectangular cross-section has mainly been studied from the point of view of their use in pin-fin heat sinks. Typical of these studies are those of Starner [37], Yazicioglu, and Yancu [38], and Sahray, et al. [39]. The present study is an extension of that described in Refs. [40][41][42].

Experimental apparatus and procedures
To study the natural convective heat transfer from a square cylinder which is inclined at an angle to the vertical, a free convection test chamber in which the model could be mounted was designed and constructed. Figure 3 shows a schematic front view of the experimental apparatus. The experimental test chamber had a size of 120 cm height 9 25 cm width 9 30 cm depth of Plexiglass. It was opened from the top and bottom mounted on a horizontal floor so it surrounds the test section. The volume of the test chamber was large enough to ensure that it did not interfere with the flow over and heat transfer rates from the model. The test chamber was constructed using transparent acrylic plates and constructed in such a way that it can be rotated around a fixed horizontal axis, in this way allowing the inclination angle of the test section to be changed. Furthermore the test chamber is placed in a bigger box 190 cm height 9 135 cm width 9 125 cm depth which was vented to the room. This arrangement ensures that no external disturbances in the room air or short term temperature changes in the room interfere with the experiments. The active component of the test apparatus was the solid aluminum square cylinder models used in the experiments, the model sizes are given in Table 1. The square cylinder model has a height of h and cross-sectional width w. The aspect ratio of the models used w/h range from 0.25 to 1. The bottom surface of the model is attached to a large base made of plexiglass 29 cm in length, 23.5 cm in width, and 1.2 cm thick. The end of the model in contact with this base being internally chamfered in order to reduce the contact area between the model and the sheet in order to reduce the conduction heat transfer from the model to the base. A series of small diameter holes were drilled longitudinally to various depths into the models. Seven thermocouples type-T 0.25 mm in diameter were inserted into these holes in various locations were used to measure temperatures. The thermocouples were monitored using a data acquisition system manufactured by TechmaTron Instrument (model USB-TC) which was self calibrating and that in turn was connected to a personal computer. The thermocouples with the data acquisition system were calibrated in a digital temperature controller water bath manufactured by VWR model 1166D, with the calibrated reference thermometer manufactured by Guildline Instruments Company model 9535/01,02,03 between 20 and 100°C and their output are less than ±0.5°C of the actual temperature. As noted before the model assemblies could be mounted inside the test chamber as shown in Fig. 3. This allowed the models to be set at any angle to the vertical.

Derivation of the experimental results
The heat transfer rates were determined using the transient method, i.e., by heating a model being tested and then measuring its temperature-time variation while it cooled. In an actual test, the model being tested was set at a required angle and then heated in an oven to a temperature of about 140°C. The model temperature variation with time was then measured while it cooled from 95-40°C. Tests and approximate calculations indicated that, because the Biot numbers existing during the tests were very small less than 1 9 10 -4 which meets the general requirement for applicability of the lumped system method of Bi \ 10 -1 . As a result of low Biot number, the temperature of the aluminum models remained effectively uniform at any given instant of time during the cooling process. The overall heat transfer coefficient could then be determined from the measured temperature-time variation using the usual procedure, i.e., the temperature variation of the model was therefore related to the time by: Hence, h t could be determined from the measured variation of ln(T i -T ? )/(T e -T ? ) with t, and known value of (A/mc).
The value of h t so determined is, of course, made up of the convective heat transfer to the surrounding air, the radiant heat transfer to the surroundings and the conduction from the model to the base. The radiant heat transfer could be allowed for by calculation using the known emissivity The conduction heat transfer to the base was calculated by fully covering the models with Styrofoam insulation and using the transient method as described before. The convective rate of the heat coefficient by convection is calculated as: The estimated values of convection, conduction and radiation heat transfer coefficients depend on the size of the experimental model, orientation and Rayleigh number, e.g. when the Rayleigh number is equal to 5. The dimensional analysis indicates that natural convection heat transfer from short square cylinders inclined at an angle to the vertical depends on Nusselt number, Nu, Rayleigh number, Ra, the ratio of the width to the height of the square cylinder, W, and the inclination angle to the vertical, u, i.e.; where Ra is defined as, and the average Nusselt number is defined as usual as: In expressing the experimental results in dimensionless form, all air properties in the Nusselt and Rayleigh numbers have been evaluated at the mean film temperature ðT w avg þ T F Þ=2 existing during the test.
The uncertainty in the present experimental values of Nusselt number arises due to uncertainties in the temperature measurements, non-uniform heating, temperature differences in the model during cooling, and uncertainties in the corrections applied for conduction heat transfer through the base. Therefore, an uncertainty analysis was performed by applying the estimation method proposed by Moffat [43,44]. The uncertainty in the average Nusselt number was found to be a function of temperature, the angle of inclination and the size of the experimental model used. Overall, the uncertainty in the average Nusselt number was estimated to be less than ±12 % for selected time interval, and the scatter in the experimental data is between ±3 % of the mean values.
Tests were performed with models mounted at various angles of inclination between vertically upwards and vertically downwards as shown Fig. 2. This work will focus on the effect of inclination angle and the aspect ratio of the square cylinder with exposed top surface on the mean heat transfer rate.

Numerical solution procedure
The solution domain used in obtaining the numerical solution is shown in Fig. 4. The flow has been assumed to be symmetric about the vertical center plane DEFOTS-IND shown in Fig. 4 and to be steady and laminar. It has also been assumed that the fluid properties are constant except for the density change with temperature which gives rise to the buoyancy forces, this being treated here by using the Boussinesq approach. The governing equations have been written in dimensionless form using the height, h, of the cylinder as the length scale and T w -T F as the temperature scale, T w being the temperature of surface of cylinder and T F being the fluid temperature far from the cylinder. The following reference velocity has been introduced: where Pr is the Prandtl number and Ra is the Rayleigh number based on h, i.e.: The following dimensionless variables have then been defined: where T is the temperature, T w is the temperature of the hot cylinder walls and T F is the fluid temperature far from the cylinder. As shown in Fig. 4, the z-coordinate is measured along the axis of the cylinder and the x-and y-coordinate are measured in the plane of the flat adiabatic base plate. In terms of these dimensionless variables and if u is the angle from the vertical at which the plate is inclined the governing equations are the following: Because the flow has been assumed to be symmetric about the vertical center-plane of the cylinder, the domain used in obtaining the solution is as shown in Fig. 4. The assumed boundary conditions on the solution in terms of the dimensionless variables are that for all the cylinder surfaces the dimensionless velocity components are all equal to zero and the dimensionless temperature is equal 1. On the adiabatic base the assumed boundary conditions are that all dimensionless velocities are equal to zero and that the dimensionless temperature gradient normal to all surfaces is equal to 0. On the plane of symmetry the dimensionless velocity component normal to the plane is zero and the gradients normal to this plane of the remaining velocity components and temperature are zero. On the outer planes of the solution domain the dimensionless temperature is equal to zero and the velocity components in the plane parallel to the surface are zero. On the upper surface HTSR the boundary conditions are the dimensionless temperature is equal to zero and the pressure is equal to zero. The outer boundary conditions of the solution domain were changed as the domain rotated from 0°to 180°.
The above dimensionless governing equations subject to the boundary conditions discussed above have been numerically solved using the commercial finite-volume solver, FLUENT. In obtaining the numerical results the pressure-based solver was used with the second-order upwind scheme for convective terms in the mass, momentum and energy equations. For pressure discretization, the Presto scheme has been employed while the SIMPLE-algorithm has been used for pressure-velocity coupling discretization. Extensive grid and convergence criterion independence testing was undertaken for all cases studied in the present work. A non-uniform grid distribution is used in the plane perpendicular and parallel to the main flow direction. This indicated that the heat transfer results presented here are to less than 1 % independent of the number of grid points and of the convergence-criterion used. Some typical results obtained here in the grid-independence testing are shown in Tables 2 and 3. These results are for the case of natural convective heat transfer from vertical and inclined square cylinders when dimensionless width, W, is equal to 1. In obtaining the results the number of nodes used was 1,191,641. Close to the square cylinder, the number of grid points or control volumes is increased to enhance the resolution and accuracy. The results in Tables 2 and 3 show that the mean Nusselt number for the surfaces that make up the cylinder for the selected number of nodes are to less than 1 % independent of the number of grid points. The effect of the distances of the outer surfaces of the solution domain (i.e., surfaces BRSN, AHTOF, ABRH, and HRST in Fig. 4) from the heated surfaces was also examined and the positions used in obtaining the results discussed here were chosen to ensure that the heat transfer results were independent of this positioning to within 1 %.
The mean Nusselt number for the heated cylinder is defined by: where " q 0 is the mean heat transfer rate per unit area from the entire surface of the cylinder. Similarly the mean Nusselt numbers for the left, right, front and top side surfaces of the cylinder are defined as follows: where q 0 left , q 0 right , q 0 front and q 0 top are the mean heat transfer rates per unit area from the left, right, front and top side surfaces of the cylinder respectively.

Results
The solution has the following parameters: 1. The Rayleigh number, Ra, based on the 'height' of the heated cylinder, h, and the overall temperature difference T w -T F 2. The dimensionless width of the cylinder, W = w/h 3. The Prandtl number, Pr 4. The inclination angle u of the cylinder Because of the applications that motivated this study, results have only been obtained for Pr = 0.7. A wide range of the other governing parameters have been considered.
In order to validate the results a comparison between the present experimental results and the correlation equations by Churchill and Chu [45] and Oosthuizen [27] are shown  in Fig. 5. Oosthuizen [27] correlation equation is valid for a vertical isothermal square cylinder with an exposed top surface where the dimensionless width of the cylinder is less than 0.6. Churchill and Chu [45] correlation equation are valid for a wide plate. Good agreement between the correlation equations and the experimental results for isothermal vertical square cylinder for W = 0.5 where the heat transfer from the top surface is small compared to the heat transfer from the side surfaces. Typical variations of numerical and experimental values results of the mean Nusselt number for the entire cylinder, Nu, with Rayleigh number, Ra, for various values of the dimensionless cylinder width, W, at different angles of inclination, u, between vertically upward and vertically downward are shown in Figs. 6, 7, 8, 9 and 10. These figures show that the mean Nusselt number increases as the Rayleigh number increases and the dimensionless cylinder width decreases for all values of inclination angles considered in this study. The increase in the Nusselt number with decreasing dimensionless cylinder width, W, arises from the fact that there is an induced inflow from the edges towards the center of the faces of the square cylinder and this causes the heat transfer rate to be higher near the vertical edges of the square cylinder than it is in the center region of the cylinder. The experimental results given in these Figures seen to be in good agreement with the numerical results and shows the same trends as the numerical results. However it will be seen from Figs. 9 and 10 that for angle of inclination of 135°and 180°the experimental results have lower values than the numerical results. This is because as the angle of inclination increases above 90°the base of the square cylinder which is made of plexiglass as mentioned before gets heated by the upward buoyant flow over the square cylinder and can no longer be treated as adiabatic as assumed in the numerical work. Typical variations of the mean Nusselt number for the entire cylinder, Nu, with dimensionless cylinder width, W, for various values of Ra and u are shown in Fig. 11. It will be seen from this figure that, the mean Nusselt number increases with decreasing W essentially at the lower values of W considered. At the lower values of Ra considered and when the dimensionless cylinder width, W, approximately greater than 0.6 the mean Nusselt number is essentially independent of the angle of inclination, u. However at larger values of Ra the mean Nusselt number depends on the angle of inclination u. Figures 12 and 13 shows the variation of the mean Nusselt number with angle of inclination for different values of W. Figure 12 shows a good agreement between the numerical and the experimental results for inclination angles below 90°and that the experimental results have lower values than the numerical results as the angle of inclination gets larger than 90°as discussed before. Considering the numerical results given in Fig. 13 which shows the variation of the mean Nusselt number for the square cylinder with u for various values of Rayleigh number and of the dimensionless cylinder width, W, it will clearly be seen that the mean Nusselt number at the lower values of Rayleigh numbers varies differently than at the higher values of the Rayleigh numbers. As already mentioned at the lower values of the Rayleigh number the mean Nusselt number is almost independent of the angle of inclination. However at higher values of Rayleigh number this is not the case. At the higher values of the Rayleigh number when dimensionless square cylinder width, W, is bigger than 0.25 the lowest mean Nusselt number occurs when the cylinder is in a horizontal position, i.e., when u is equal 90°while the highest mean Nusselt number occurs when the cylinder is at inclination angles, u, of 45°and 135°. When the dimensionless square cylinder width, W, is equal to 0.25 the lowest mean Nusselt number occurs when the cylinder is in a vertical position, i.e., when u is equal 0°and 180°while  the highest mean Nusselt number occurs when the cylinder is at inclination angles, u, of between 45°and 135°. In some situations the mean heat transfer rates from the individual surfaces, i.e., the left side, the right side, the front, and the top surfaces of the cylinder are of importance and these will be considered next. Typical variations of the mean Nusselt numbers for the front, right side, left side, and top surfaces with angles of inclination of various values of W and Ra are shown in Figs. 14, 15, 16 and 17, respectively.  Figure 14 shows that Nu front increases as the dimensionless cylinder width W decreases and that the variation of Nu front with u at lower values of Rayleigh number differs from the variation at higher values of Rayleigh number. At lower values of Rayleigh number Nu front has a maximum value when the inclination angle u is equal to 90°and a minimum value when u is equal to 180°whereas at higher values of Rayleigh number the maximum value of Nu front occurs when the angle of inclination u is equal to 45°while the lowest value again occurs when the angle of inclination u is equal to 180°. It will be noted that at the highest value of Ra considered in Fig. 14, Nu front increases as u increases from 0°to 45°. Nu front then decreases, as u increases from 45°to 90°, then increases as u increases from 90°to 135°, and then again decreases as u increases from 135°to 180°. When dimensionless square cylinder width, W, is equal to 0.25 the maximum mean Nusselt number occurs when the cylinder is in a horizontal position, i.e., when u is equal 90°. The results given in Fig. 15 shows that the variation of Nu right with u exhibits the same basic form of behavior as Nu front (Fig. 14), except that at the highest Rayleigh number considered when dimensionless square cylinder width,W, is equal to 0.25 the maximum mean Nusselt number occurs when u is equal 45°and 135°. Figure 16 shows the variation of Nu left with u for various values of Ra and W. It will be seen that at lower Rayleigh numbers Nu left increases as the dimensionless plate width W decreases at all angles of inclination. However at higher Rayleigh numbers a different form of variation of Nu left with W exists. At high Rayleigh numbers Nu left increases as the dimensionless plate width W decreases for values of u near 0°and 180°but decreases as the dimensionless plate width W decreases for values of u near 90°This is because as the angle of inclination of the cylinder increases 0°the left side of the cylinder changes from being an effective vertical flat plate to being horizontal plate for u near 90°and then again become an effective vertical flat plate for u near 180°.
The changes in the Nusselt number of the individual surfaces is further illustrated by the results given in Fig. 17 which shows that the Nu Top increases as the dimensionless cylinder width W decreases at all angles of inclination except at the highest Rayleigh number considered when at angles of inclination u near to 0°where Nu Top increases as the dimensionless plate width W increases. It should also be noted that at the lower values of Rayleigh number considered Nu Top increases continuously as the angle of inclination increases. However at higher values of Rayleigh number Nu Top increases as the angle of inclination increases but then passes through a maximum at an angle of inclination u of about 120°and then decreases with a further increase in u.
Typical variations of Nu Top , Nu front , Nu right , and Nu left , with an angle of inclination u and for two values of Ra and W equal 0.5 are shown in Figs. 18 and 19. Nu front , Nu right , and Nu left , have the same values when the angle of inclination is equal to 0°and to 180°. As expected, it will be seen that Nu Top is much less than Nu front , Nu right , and Nu left when the angle of inclination is equal to 0°. However as the angle of inclination increases this difference decreases and for angles greater than about 75°the Nu Top is greater than Nu front , Nu right , and Nu left . These changes arise from the fact that as the angle of inclination of each of the surfaces relative to the vertical changes, e.g., when the angle of inclination of the cylinder is 0°the 'top' surface faces upwards The changes in the three-dimensional flow patterns and the flow interaction between the surfaces that make up the square cylinder with angle of inclination on the top surface of the square cylinder are illustrated by the local Nusselt number distributions over the heated plate as shown in Fig. 20, where the dimensionless local Nusselt number are defined as: The variation of the mean Nusselt number for the top surface relative to the mean Nusselt numbers for the other surfaces has already been discussed. The relative importance of the heat transfer rate from the top surface compared to that from the vertical side surfaces of the cylinder will however depend both on the mean Nusselt numbers for the surfaces and the relative surface areas. Now since: where A total , A side , and A Top are the dimensionless surface areas of the entire cylinder, of all of the 'sides' of the cylinder, and of the 'top' surface of the cylinder respectively, it follows that: This equation indicates that over the range of values of W considered here, i.e., 0.25 to 1, A Top /A side is relatively small having a maximum value of 0.25 when W = 1, i.e., the area of the top surface relative to the area of the vertical side surfaces remains comparatively small. Now: This equation allows Q total /Q Top to be found using the calculated Nusselt number values. Some typical results obtained are shown in Fig. 21. From these results it will be seen that while Q total /Q Top is always greater than approximately 4 it is only at low values of u that the contribution of Q Top to Q total can be neglected. The assumption that the contribution of Q Top to Q total can be neglected has sometimes been used in deriving correlation equations for the heat transfer rate from vertical square cylinders using: and neglecting the second term on the right hand side. The present results show that this approach will not normally be applicable to the inclined cylinder case. The surfaces that make up the square cylinder are of the form of flat plates connected to each other. The correlation equation for the case of a wide flat plate with a uniform surface temperature has the form: The parameter B depending only on the Prandtl number of the fluid involved. When the width of the plate becomes relatively small and inclined at an angle to the vertical edge effects become important (see Kalendar and Oosthuizen [46]) and the Nusselt number depends on the boundary layer thickness to plate width ratio, i.e., since the boundary layer thickness will be dependent on the value of: The Nusselt number will of course also be dependent on the inclination angle. Therefore when the surfaces that make up the square cylinder are narrow and inclined at an angle to the vertical, the heat transfer rate from the square cylinder is given by an equation of the form: Using an equation of this form it has been found that the present results for the case of a heated square cylinder with an exposed top surface and inclined at an angle, u, of between u = 0°and u = 180°, and with different dimensionless widths can be approximately described by: The index on the left hand side of Eq. (28) has a Rayleigh number to the power of 0.28 instead of 0.25 since the interaction of flow over the plates that make up the square cylinder become more important as the plate width decreases. A comparison of the results given by this correlation Eq. (28) with the numerical and experimental results is shown in Fig. 22. It will be seen that the equation describes the computed numerical results to an accuracy of better than 92 %. Figure 23 shows the comparison of the experimental results with the correlation Eq. (28), the Fig. 21 Variation of the ratio of the total heat transfer rate from the surface of the cylinder to the heat transfer rate from the top surface with angle of inclinations for Ra values of 10 4 and 10 6 and for W = 1 Fig. 22 Comparison of correlation equation with the numerical and experimental for inclined square cylinder with exposed top surface Fig. 23 Comparison of correlation equation with the experimental results for inclined square cylinder with exposed top surface agreement being to an accuracy of better than 86 %, where the 12 % in Fig. 23 shows the 95 % prediction band of the experimental results.

Conclusions
Natural convective heat transfer from an isothermal inclined square cylinder which have an exposed top surface and is, in general, inclined at an angle to the vertical between vertically upwards and vertically downwards has been numerically and experimentally studied. The effects of Dimensionless widths, Rayleigh numbers, and inclination angles on the mean Nusselt number for the entire cylinder and for the mean Nusselt numbers for the various surfaces that make up the cylinder were conducted in this study. The results of the present study indicate that: 1. The experimental results are in a good agreement with the numerical results for square cylinder for all angle of inclination considered in this study.

The mean Nusselt number for the cylinder increases
with decreasing W under all conditions considered. 3. The relative magnitudes of the mean Nusselt numbers for the various faces that make up the cylinder vary considerably with the inclination angle and as a result it is only possible to neglect the heat transfer from the 'top' surface compared to that from the other surfaces when the dimensionless width is small and the angle of inclination is near 0°. 4. At lower values of Ra (approximately less than 10 4 ) for all dimensionless plate widths, W, considered the mean Nusselt number is independent of angle of inclination u. At larger values of Ra the dependence of the mean Nusselt number on the angle of inclination u becomes significant. 5. At the higher values of Rayleigh number the lowest mean Nusselt number occurs when the cylinder is in a horizontal position, i.e., when u is equal 90°while the highest mean Nusselt number occurs when the cylinder is at inclination angles, u, of 45°and 135°. 6. The results for the mean Nusselt number based on the mean heat transfer rate from vertical and inclined square cylinder can be adequately correlated by the Eq. (28).