Experimental Study For Optimum Fin Spacing of Rectangular Fin Arrangements Under The Influences of Free Convection
Experimental Study For Optimum Fin Spacing of Rectangular Fin Arrangements Under The Influences of Free Convection
                                        Abstract
    Heat transfer enhancement systems are widely used in various industries, such
    as thermal power plants, heating and air conditioning systems, and electronic
    equipment. In this research, the thermal performance of a vertical heat sink with
    rectangular fins geometry was investigated experimentally. The present study
    mainly aims to find the optimum fin spacing that provides the highest heat
    transfer rate and reduces the heat sink base temperature. An Aluminum heat sink
    of six 300 mm length fins was considered in the test. The fin thickness and height
    were kept constant at 4 mm and 45 mm, respectively. In the analysis, a range of
    fins spacing (i.e., 10, 11, 12, 13, and 15 mm) was used in the experiments, and a
    variety of Rayleigh number (i.e., from 7.6*106 to 1.48*108) was considered. The
    difference of heat sink and surrounding temperatures were obtained at various
    heat input and all the different parameters above. These temperature differences
    are then used for the estimation of the heat transfer rate of the fin array. The
    results showed that the optimum fin spacing was S = 12 mm, which gives the
    highest heat transfer rate. Besides, a reduction of 18 % of the heat sink weight is
    achieved due to the use of the optimum fin spacing compared to the weight of
    using 15 mm as a fin spacing. The reduction of heat sink weight and size can no
    doubt lead to a decrease in manufacturing cost, which consequently lead to the
    spread of heat sink applications. Furthermore, the decline of the heat sink base
    temperature will improve the working life span of the electronic devices that use
    the heat sink. The results of the present study also show the development of a
    new empirical equation predicting the Nusselt number in terms of Rayleigh
    number and the fin spacing to fin height ratio. The developed equation shows the
    high accuracy of prediction (e.g. about 90%).
    Keywords: Fins, Fin spacing, Heat sink, Heat transfer rate, Natural convection.
                                           4075
4076     A. J. J. Al-Jassani
1. Introduction
Natural convection phenomena in enclosures are necessary for an excellent
performance of high-power density electrics. Buoyancy drove flows have many
applications in widely preferred phenomena, such as air-cooled car engines, cooling
of generators, motors, refrigerators, transformers, and cooling of computer
processors. Natural convection is a method of heat dissipation that is relatively
inexpensive, discreet, and most dependable. The heat generated by electronic devices
can be controlled using fins. In the design of an efficient cooling system, generally,
55% of failure mechanisms in electronic devices are related to thermal effects Pascoe
[1]. However, many parameters, such as fins height and fins spacing, affect the rate
of convection heat dissipation. These parameters should be considered when the fins
best design or the enhancement of heat dissipation rate is targeted.
    The subject of investigation of the thermal performance of finned systems still
attracts the interest of many researchers. A considerable number of researches have
been published. Part of these publications focuses on the development of
mathematical formulas that considered the fin arrangement specification. Some
other parts of the papers focus on the investigation on fins system arrangements
such as fins number, shapes, dimension, and fin spacing. For examples of the
former part of studies above (i.e., studies that develop mathematical correlation
formula), Elenbaas [2] studied the effect of small gap width between parallel plates
on heat dissipation under natural convection condition experimentally.
    The results developed a proportional relation between the Nusselt number and
the Rayleigh number. Yazicioğlu and Yüncü [3] investigated the thermal
performance of the vertical heat sink of rectangular fins under different parameters.
The fin base–to-surrounding temperature difference ranges from 30 to 150 K. The
fin length is taken as 250 mm and 340mm. The space between the fins and the
height of the fin was from 4.5 to 85.5 mm and 5 to 25 mm, respectively. The authors
reported that the optimum value of fin spacing varies with the height of the fin,
which was from 6.1 to 11.9 mm. Besides, the rate of heat transfer enhanced
depends on the average base temperature, fin length, and fin spacing.
    The influence of fin height, fin spacing, and fin orientation on the natural heat
transfer coefficient were studied experimentally by Walunj, et al. [4]. Bar-Cohen
and Rohsenow [5] conducted an analytical investigation for the natural convective
heat transfer from two parallel plates. Their work includes the development of a
correlation formula describing Nusselt number in terms of Rayleigh number for
isothermal and isoflux plates. The optimum fin spacing was also reported in the
same study.
    Bodoia and Osterle [6] developed a numerical technique for the investigation
of flow in channels and heat transfer between symmetrically heated, isothermal
plate to predict the channel length required to achieve fully developed flow as a
function of the channel width and wall temperature. The heat transfer in inclined
interrupted fin channels was investigated in Fujii [7]. A mathematical correlation
was developed for fitting the experimental results. As was mentioned previously,
some publications investigated the influence of some fin geometrical properties
such as fins number, fins shapes, fins dimension, and fin spacing. Experimental
studies on vertical rectangular fin arrays ware carried out by Güvenç and Yüncü
[8] and Leung and Probert [9].
heat dissipation of the heat sink. Consequently, this optimum spacing will no doubt
reduce the weight and manufacturing cost of the heat sink. In addition, an empirical
equation was developed for describing the relation between Nusselt number,
Rayleigh number, and dimensionless parameter of fins spacing and fins height.
3. Methodology
The methodology contains constructing a test rig of the finned heat sink for
conducting several experiments to obtain an average base and ambient temperature
at different heat inputs. The heat input can be determined by knowing the supplied
electrical voltage and current provided to the electrical heater. Then, the convection
heat transfer rate can be determined by subtracting the radiation and loss heat rates
from the generated heat transfer (e.g., input heat). Different convection heat
transfer coefficients can then be obtained at different input voltage and current,
different heat sink fins spacing. Eventually, an optimum fin spacing can be found
as that fin spacing providing heights heat transfer coefficients. Besides, an
empirical equation predicting the relation of the Nusselt number in terms of
Rayleigh number and the fin spacing to fin height ratio. The reduction of heat sink
weight and size can no doubt lead to a decrease in manufacturing cost, which
consequently lead to the spread of heat sink applications.
4. Experimental Procedure
The experiments were conducted at various heat flux. This variation is made through
controlling the heater voltage by a voltage regulator. The provided voltage to the
heater is varied from 50 to 150 V by 25 V step. The surrounding temperature was
measured via a thermometer. The fin base plate temperature was acquired at eight
positions to calculate the average temperature along with the base plate. The
temperature was measured using a temperature scanner thermometer (8 channel
thermocouple) manufactured by Altop Industries Ltd, Gujarat, India (Si. 1005164.
Model ADT 5003) with a resolution of 0.1. The average of these eight positions
represents the fin base temperature. Figure 2 shows the experimental rig with the
measuring instrument actions. The predefined heater inputs were modified with the
assistance of dimmer stat. The temperature of an aggregated fin heat sink array of
various locations and surrounding temperatures was registered during time intervals
of 30 min. It requires close to four hours to reach a steady condition. Consequently,
when the temperature reading does not vary by 0.3°C, in this case, the system reaches
a steady-state condition. The experimental tested was carry out in the Fluid Dynamics
Laboratory of the mechanical engineering department/Wasit University.
5. Analytical Analysis
The heat input of the heating base was dissipated based on two approaches,
primarily through the heat sink and, slightly by the insulation fins and heat sink
structure. A voltage V is provided through the voltage regulator to the heater, which
is fixed on the heat sink base. The temperature from the eight positions (predefined
in the previous section) as well as the ambient temperature is acquired. The
provided voltage and the received temperature are used for the calculation of
Qconv, as in the following equations Jubear [31].
𝑄𝑄𝑔𝑔𝑔𝑔𝑔𝑔. = 𝑄𝑄𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐. + 𝑄𝑄𝑟𝑟𝑟𝑟𝑟𝑟. + 𝑄𝑄𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙                                 (1)
where
𝑄𝑄𝑔𝑔𝑔𝑔𝑔𝑔. = 𝑉𝑉 × 𝐼𝐼                                                               (2)
     And the heat transfer by radiation was:
                                         4       4
𝑄𝑄𝑟𝑟𝑟𝑟𝑟𝑟. = σ × 𝜀𝜀 × 𝑆𝑆𝑠𝑠𝑠𝑠𝑠𝑠 × 𝐴𝐴𝑡𝑡 (𝑇𝑇av. − 𝑇𝑇𝑎𝑎𝑎𝑎𝑎𝑎 )                         (3)
    The loss is estimated according to the temperature difference between the heat
sink base and ambient, as follow Feng, et al. [32]:
𝑄𝑄𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙 = −9 ∗ 10−5 ∗ ∆𝑇𝑇 2 + 0.0202 ∗ ∆𝑇𝑇                                     (4)
   In the present study, the heat loss for all tested was less than 4% of the overall
power input. Consequently, for these ranges of heat input, the heat loss is ignored.
Therefore, the convection heat transfer was calculated as follows:
     𝑄𝑄𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐. = 𝑄𝑄𝑔𝑔𝑔𝑔𝑔𝑔. − 𝑄𝑄𝑟𝑟𝑟𝑟𝑟𝑟. − 𝑄𝑄𝑙𝑙𝑙𝑙𝑙𝑙𝑙𝑙                            (5)
   Hence, the heat transfer coefficient in free convection was calculated from the
equation which, is known as Newton’s cooling law, as follows:
                  𝑄𝑄𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐.
     ℎ=                                                                          (6)
            𝐴𝐴𝑡𝑡 (𝑇𝑇𝑎𝑎𝑎𝑎. − 𝑇𝑇𝑎𝑎𝑎𝑎𝑎𝑎)
   The Rayleigh Number was defined as the ratio between buoyancy and viscosity
and written as follows:
               𝐿𝐿3 𝑔𝑔 𝛽𝛽 (𝑇𝑇𝑎𝑎𝑎𝑎. −𝑇𝑇𝑎𝑎𝑎𝑎𝑎𝑎)
     𝑅𝑅𝑅𝑅 =                                    𝑃𝑃𝑃𝑃                              (9)
                             𝜈𝜈2
Fig. 7. The effect of heat transfer               Fig. 8. The effect of heat transfer
coefficient on heat transfer rate.                 coefficient on Nusselt number.
Figs. 9 to 11. The effect of heat input on the fin base-to- surrounding temperature
difference can be seen more clearly in Fig. 9. In this figure, the fin base-to-
surrounding temperature difference is plotted as a function of fin spacing. It can be
deduced that the fin base-to-surrounding temperature difference from a fin array
decreases with the space of fins, and reach to a minimum level when the fin spacing
was 12 mm, and then increasing. The corresponding fin spacing value of the
minimum base-to-ambient temperature difference is named the optimum fin
spacing. It was noted that the optimum fin spacing of this study was S = 12 mm.
This result reveals that the optimum value of fin spacing was hypersensitive to the
variations in heat generated and temperature difference parameters.
    The values of the natural convection heat transfer coefficient obtained for
different heat inputs are shown in Fig. 10. It can be seen that when the space
between two fins (i.e., gap) increases from 10 to 12 mm, the free convection heat
transfer coefficient ultimately increases. Besides, a further increase in the fin
spacing leads to a decrease in this coefficient. If the gaps between the fins are
evenly matched or small, the free convection heat transfer coefficient is minimized
as mixing of the thermal boundary layer occurs. Figure 10 clearly shows that the
heat transfer coefficient decreases as the gap between fins decreases. However, if
the fins are closely spaced, there is also a more dissipating surface area (more fins
for a given volume). The additional surface area can counteract the reduced heat
transfer coefficient.
    More precisely, a variation of the convective heat transfer coefficient with
different heat input is studied in this figure. Initially, the heat transfer coefficient
enhanced when the space among the fins increases. It reaches a maximum value at
fin spacing (S = 12 mm), and with further additions of fin spacing, it begins to
decrease. The space between the fins is called optimum value when the convection
heat transfer coefficient has maximized. Inspection of Fig. 10, reveals that the
optimum fin spacing value was (S = 12 mm) at the different heat input. One of the
essential parameters in designing a heat sink is the fin spacing, S. The enhancement
   of the heat transfer coefficient due to optimum fin spacing was 27%. The same
   behaviour is shown in Fig. 11.
 Fig. 10. The effect of fin spacing on heat    Fig. 11. The effect of fin spacing on Nusselt
transfer coefficient of various heat input.           number of various heat input.
 Fig. 12. The relation between heat                  Fig. 13. The relation between
 transfer rate and Nusselt number.               Nusselt number and Rayleigh number.
      It is essential to mention that 18% of heat sink weight is achieved due to the
   use of the optimum fin spacing compared to use the 15 mm. The reduction
   percentage can be calculated either by 1) weighing the heat sink at using fin spacing
of 15mm and optimum fin spacing (e.g.12 mm) or 2) calculating the heat sink
volume and multiplying by the material density to fins the mass.
    The accuracy of the empirical equation can be statistically assessed using the
following measures Zubaidi, et al. [36]:
                       𝑦𝑦𝑜𝑜 −𝑦𝑦p 2
%𝑅𝑅𝑅𝑅 = ∑𝑁𝑁
         𝑚𝑚=1 �                 � ∗ 100                                          (11)
                         𝑦𝑦𝑜𝑜
                                            2
                         ∑𝑁𝑁
                          𝑚𝑚=1�𝑦𝑦𝑜𝑜 −𝑦𝑦p �
𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁𝑁 = 1 − �                                                               (12)
                                     𝑁𝑁
            ∑𝑁𝑁
             𝑚𝑚=1(𝑦𝑦𝑜𝑜 −𝑦𝑦
                        ����)       �����
                          𝑜𝑜 �𝑦𝑦𝑝𝑝 −𝑦𝑦𝑝𝑝
𝑅𝑅 =                                            2
                                                                                 (13)
       �∑𝑁𝑁
         𝑚𝑚=1(𝑦𝑦𝑜𝑜 −𝑦𝑦
                          2 𝑁𝑁
                      𝑜𝑜 ∑𝑚𝑚=1�𝑦𝑦𝑝𝑝 −𝑦𝑦
                    ����)            �����
                                       𝑝𝑝
where:
RE     Relative error.
RMSE Root mean square error
NRMSE normalized root means square error
  The statistical measures of the predicted and experimental results show that
%RE = 10, NRMSE= 0.43 and R= 0.84
    For the present work, the uncertainty of the experimental results is obtained
using mathematical formulas in Moffat [37]. The uncertainty parameter for the
measured base temperature was ∓1.2℃. The correlation coefficient of the
innovative and predicted Nusselt values are computed and found 0.82.
7. Conclusions
In this study, the effect of different fin spacing, heat flux, and Rayleigh number
(Ra) on the performance of vertical rectangular fins under free convection heat
transfer was investigated experimentally. The heat transfer coefficient and base
temperature are obtained analytically. From the different results, it was found that
the optimum fin spacing is S=12 mm. This optimum S was selected in terms of
providing maximum heat transfer coefficient, maximum Nusselt number, and heat
transfer rate. Additionally, the fin base-surrounding temperature increases by 25%
at this fin spacing. The selection of optimum S also leads to a18% of fin weight
reduction. In the present study, a precise empirical correlation equation was also
developed for describing the Nusselt in terms of Rayleigh number and the fin
spacing to fin height (S/H). The maximum error of the equation prediction was 10
%. The studies on the selection of fin spacing are still of great of importance as it
provides practical guidelines for manufacturing heat sinks with lighter weights,
less expensive cost and smaller size.
Nomenclatures
     Greek Symbols
     σ        Steven-Boltzmann, W/m2 K4
     β        The volume coefficient of expansion, K-1
     ε        Emissivity
     ν        Kinematic viscosity, m2/s
References
1.     Pascoe, N. (2011) Reliability Technology: Principles and Practice of Failure
       Prevention in Electronic Systems: John Wiley & Sons.
2.     Elenbaas, W. (1942). Heat dissipation of parallel plates by free convection,
       Physica, 9(1), 1-28.
3.     Yazicioğlu, B.; and Yüncü, H. (2007). Optimum fin spacing of rectangular
       fins on a vertical base in free convection heat transfer, Heat and Mass
       Transfer, 44(1), 11-21.
4.     Walunj, A.; Daund, V.; and Palande, D. (2013). Parametric analysis of plate-
       fin heat sink over heat transfer, International Journal of Research, 1(4).
5.     Bar-Cohen, A.; and Rohsenow, W. (1984). Thermally optimum spacing of
       vertical, natural convection cooled, parallel plates.
6.     Bodoia, J.; and Osterle, J. (1962). The development of free convection
       between heated vertical plates.
7.     Fujii, M. (2007). Enhancement of natural convection heat transfer from a
       vertical heated plate using inclined fins, Heat Transfer—Asian Research: Co‐
       sponsored by the Society of Chemical Engineers of Japan and the Heat
       Transfer Division of ASME, 36(6), 334-344.
23. Leung, C.; and Probert, S. (1987). Natural-convective heat exchanger with
    vertical rectangular fins and base: design criteria, Proceedings of the
    Institution of Mechanical Engineers, Part C: Journal of Mechanical
    Engineering Science, 201(5), 365-372.
24. Leung, C.; and Probert, S. (1989). Thermal effectiveness of short-protrusion
    rectangular, heat-exchanger fins, Applied energy, 34(1), 1-8.
25. Van de Pol, D.; and Tierney, J. (1974). Free convection heat transfer from
    vertical fin-arrays, IEEE Transactions on Parts, Hybrids, and Packaging,
    10(4), 267-271.
26. Chaddock, J. (1970). Free convection heat transfer from vertical rectangular
    fin arrays, ASHRAE journal, 12(8), 53-&.
27. Ahmed, S.E.; Hussein, A.K.; Mohammed, H.; Adegun, I.; Zhang, X.; Kolsi,
    L. (2014). Viscous dissipation and radiation effects on MHD natural
    convection in a square enclosure filled with a porous medium, Nuclear
    Engineering and Design, 266(34-42.
28. Ahmed, S.E.; Hussein, A.K.; El-Aziz, M.A.; and Sivasankaran, S. (2016).
    Conjugate natural convection in an inclined square porous enclosure with
    finite wall thickness and partially heated from its left sidewall, Heat Transfer
    Research, 47(4).
29. Hoseinzadeh, S.; Heyns, P.; Chamkha, A.; and Shirkhani, A. (2019). Thermal
    analysis of porous fins enclosure with the comparison of analytical and
    numerical methods, Journal of Thermal Analysis and Calorimetry, 1-9.
30. Hoseinzadeh, S.; Moafi, A.; Shirkhani, A.; and Chamkha, A.J. (2019).
    Numerical validation heat transfer of rectangular cross-section porous fins,
    Journal of Thermophysics and Heat Transfer, 1-7.
31. Jubear, A.J. (2017). An experimental investigation of heat transfer
    enhancement for vertical interrupted fin in free convection, International
    Journal of Recent Scientific Research, 8(6), 17297-17284.
32. Feng, S.; Li, F.; Zhang, F.; and Lu, T.J. (2018). Natural convection in metal
    foam heat sinks with open slots, Experimental Thermal and Fluid Science,
    91(354-362.
33. Tari, I.; and Mehrtash, M. (2013). Natural convection heat transfer from
    inclined plate-fin heat sinks, International Journal of Heat and Mass Transfer,
    56(1), 574-593.
34. Tari, I.; and Mehrtash, M. (2013). Natural convection heat transfer from
    horizontal and slightly inclined plate-fin heat sinks, Applied Thermal
    Engineering, 61(2), 728-736.
35. Mehrtash, M.; and Tari, I. (2013). A correlation for natural convection heat
    transfer from inclined plate-finned heat sinks, Applied Thermal Engineering,
    51(1-2), 1067-1075.
36. Zubaidi, S.L.; Dooley, J.; Alkhaddar, R.M.; Abdellatif, M.; Al-Bugharbee, H.;
    and Ortega-Martorell, S. (2018). A Novel approach for predicting monthly
    water demand by combining singular spectrum analysis with neural networks,
    Journal of hydrology, 561(136-145.
37. Moffat, R.J. (1988). Describing the uncertainties in experimental results,
    Experimental thermal and fluid science, 1(1), 3-17.