Theoretical and experimental analysis of the performances of a heat sink with vertical orientation in natural convection

In the various areas in which electrical components are used, the problem of heat dissipation generated due to the absorption of electrical energy assumes great interest and is worthy of an in-depth study. In steady state conditions, the thermal power generated can equal the electrical power absorbed and leads to an alteration in the physical properties of electrical components compromising their performance and correct functioning. One of the most frequently adopted solutions consists in the application of a heat sink on the surface to be cooled. Experimental tests were conducted using an infrared thermal camera, an internal climate control unit for the recording of the thermo hygrometric conditions of the environment and a finite element software (ProENGINEER) to simulate the thermal behaviour of the heat sink in order to analyse the modalities of thermal exchange of the heat sink. The results obtained were subsequently compared with the heat sink properties provided by the manufacturer. The main objective of the work is that of providing a methodology that blends the use of thermographic and simulation techniques with finite elements, in order to render the development of a theoretical–experimental correlation possible for any physical condition and geometrical configuration taken into consideration. This methodology is confirmed in the field of technological development of electrical components, where at each stage of the planning process exists a marked intertwining of computing, electronics, mechanics and heat transmission.

Abstract In the various areas in which electrical components are used, the problem of heat dissipation generated due to the absorption of electrical energy assumes great interest and is worthy of an in-depth study. In steady state conditions, the thermal power generated can equal the electrical power absorbed and leads to an alteration in the physical properties of electrical components compromising their performance and correct functioning. One of the most frequently adopted solutions consists in the application of a heat sink on the surface to be cooled. Experimental tests were conducted using an infrared thermal camera, an internal climate control unit for the recording of the thermo hygrometric conditions of the environment and a finite element software (ProENGINEER) to simulate the thermal behaviour of the heat sink in order to analyse the modalities of thermal exchange of the heat sink. The results obtained were subsequently compared with the heat sink properties provided by the manufacturer. The main objective of the work is that of providing a methodology that blends the use of thermographic and simulation techniques with finite elements, in order to render the development of a theoretical-experimental correlation possible for any physical condition and geometrical configuration taken into consideration. This methodology is confirmed in the field of technological development of electrical components, where at each stage of the planning process exists a marked intertwining of computing, electronics, mechanics and heat transmission. Greek symbols DT exper Experimental temperature difference (°C) e Percentage error

Introduction
The passing of current in semiconductors causes, due to the Joule effect, a generation of heat which determines an increase in the temperature of the element. In a steady state condition, the dissipated thermal power Q d is equal to: The equivalent thermal resistance between the junction and the external environment R th is the parameter which characterises the efficiency of a heat sink. Such resistance is supplied by the manufacturer and is defined as the increase in temperature undergone due to the application of an active electrical power of 1 W.
The component to be cooled is placed in contact with the heat sink, which presents high thermal conductivity and a high external surface due to its particular geometry (finning), with the aim of favouring convective and radiative heat exchange with the external environment. It is therefore possible to argue that the heat sink does nothing but reduce the thermal resistance between the component and the external environment.
In order to optimise heat conduction through the component and the heat sink, it is opportune to smooth the contact surface to obtain the greatest contact area between the two, thus reducing roughness to a minimum. However, since it is not possible to guarantee a perfect result in polishing and therefore in the elimination of any possible air interspace (increase in thermal resistance), use is made of highly conductive greases, such as silicone grease with silver, zinc or copper oxides.
To evaluate the radiative heat exchange of a heat sink, it is necessary to take into consideration different factors which characterise its functioning, such as the surface finish or the surface colour. Rough or shiny surface reduces heat sink efficiency and increase thermal resistance, while a dark colour consents heat sink thermal behaviour that is close to that of a 'black body', thus increasing thermal exchange efficiency.
The parametrical study of the heat sink was extensively addressed by Kraus and Bar-Cohen [1]. Elenbaas [2] was the first to create a detailed study on heat sink fins. Instead, Starner and McManus [3], Welling and Wooldridge [4] carried out several experimental studies, while Van de Pol and Tierney [5] developed a correlation with Welling and Woolridge's experimental results. Finally, Bilitsky [6] has carried out a complete parametrical analysis of heat sinks in relation to the height and thickness of the fins.
A simple and commonly used fin configuration is that obtained by using different parallel fins on a baseplate that exchanges mainly by conduction (Fig. 1). In this case, the viscous resistances in the cavity between the two adjacent fins can assume, according to the distance between the two surfaces, notable importance and reduce the value of the convective heat exchange coefficient. The optimal spacing assumes importance in the case of large surfaces in order to obtain the best equilibrium value between the opposing effects provided by the number of fins (greater exchange surface) and by the viscosity. The optimal choice is that which maximises the dissipated power for a given base area.
In Ref. [1], Kraus and Bar-Cohen provide relations for the obtainment of the optimal distance S opt for isotherm fins, whose thickness D is much smaller than the distance S between the fins (Fig. 2).
where Ra is the dimensionless Rayleigh number and H is the fin length. The average convective heat exchange coefficient for a vertical plate is equal to: where K is the thermal conductance coefficient of the baseplate. The optimal number of fins n is defined by the relation: where B is the width of the baseplate.

Description of the heat sink and instruments used
The heat sink used for the test (Fig. 3) is in black anodised aluminium in order to favour radiative heat exchange and has two series of continuous fins and a central cavity.
It is characterised by a thermal resistance (linearised) towards the environment, declared by the manufacturer to be equal to 1.1 (K W -1 ), from a baseplate of 102 mm 9 120 mm and by a total of 10 fins placed symmetrically to the central cavity. In the tests carried out, the thermal load of the electrical component was simulated with a resistance positioned centrally on the surface opposite the fins (Fig. 4), whose theoretical power can vary from 25 to 60 W.
The thermal camera used for the thermographic investigations was supplied by NEC, model TH7102MV ( Fig. 5a), with a resolution of 0.08°C at 30°C and a measurement range from -20 to 100°C and with a 48 mm objective. Emissivity can be set on the thermal camera and can vary from 0.1 to 1.0 with an increase of 0.01. The software used for the elaboration of the thermographic images ( Fig. 5b) is the MikroSpec v1.2 furnished by NEC. Such images contain all the information regarding the various measurement points (temperature, position and emissivity) and allow visualisation of the temperature distribution on the heat sink.
The power supply used to deliver electrical current to the resistance has three outputs for which the tension and current (DC) can be adjusted, and it is equipped with a digital display for the visualisation of the tension and  current values (Fig. 5c). The maximum tension, which can be supplied, is 32 V and a current equal to 2 A DC. There are 2 U with output 0-32 V, 2 A and 1 U 0-5.5 V, 5 A. The power output is therefore equal to 64 W for the two canals 0-32 V and equal to 27.5 W for the output 0-5.5 V. It is possible to obtain greater power output operating in a serial or parallel mode. In the case of the parallel mode of the two units with a potential difference equal to 32 V, the resulting current is 4 A for a total power of 128 W. For the measurement of the environmental hygrothermal conditions, a portable climatic data logger, LSI Babuc/A ( Fig. 5d), equipped with opportune probes and sensors, such as a psychrometer and an anemometer to measure air velocity ( Fig. 5e), was used.

Methodology used
Once the electrical resistance and heat sink simulation models were created with ProENGINEER ( Fig. 6), laboratory testing was started using the equipment as described in the previous paragraph. Resistance was set on the rear and central part of the heat sink and was connected to the power supply. The heat sink was positioned vertically and was placed frontally on a black surface with an approximate emissivity of 0.95 and at a distance to avoid local thermal exchange phenomenon, isolating as much as possible from unwanted thermal exchanges, such as conductive and radiative exchanges with other bodies.
The emissivity value for the heat sink used has been certified by the manufacturer. However, during the tests performed in the laboratory, the emissivity value has been set on the thermal camera and, using a thermocouple, was measured the temperature of various points of the heat sink and compared with the values resulting from the camera by varying the value of emissivity so long as the two values was coincident.
The experimental tests were carried out over following days since the attainment of steady state conditions has required a long time, since even though the environmental conditions were stable, some small temperature oscillations occurred. The meeting of stationary conditions is at the basis of the methodology used in that all the electrical energy absorbed must be dissipated in thermal energy. For each test, by means of the climatic data logger, the dry bulb and relative humidity temperature values were recorded. Furthermore, through the use of the hot-wire anemometer, fluid speed in the various cavities was measured, and it was found that the flow on all the plate surfaces and for all the range of powers applied resulted as being laminar.
By means of thermal camera monitoring, it was possible to observe the transitory phenomenon trend, and in the moment in which the stationary condition was met, a thermographic image of the heat sink temperature map was recorded. By analysing the image, the temperature of the different heat sink surfaces was obtained which is necessary in order to evaluate the flow and calculate all the fluid- The film temperature is defined as: where t p is the fin surface temperature and t a is the temperature of the fluid in the cavities. Once the film temperature is determined, it was possible to estimate the convective heat exchange coefficient, utilising the well-known McAdams correlation: By using the MicroSpec software, an analysis of the temperature maps for each individual test was carried out in order to find the correlation that provided the convective heat exchange coefficient which best approximates the obtained results. Due to the particular geometry and the natural convection conditions, the use of a single convective heat exchange value for the total area would not allow the obtainment of the real distribution of the temperatures on the heat sink by means of a simulation with ProENGINEER.
The convective coefficient value calculated with the theoretical method on the surface taken into consideration was the starting point for each individual test. Subsequently, proceeding with attempts and diverse simulations, the condition in which the convective coefficient values calculated provided the temperature distribution which best approximated those recorded with the thermal camera, was met. Once such coefficients were set on the different surfaces of the model created with ProENGINEER software and the surface thermal load produced by the resistance on the rear of the heat sink was applied, a finite elements analysis was carried out thus obtaining the experimental trend of the temperatures.
The method used is an iterative method and its control criterion is that to compare the real map of temperatures obtained through thermographic measurement with the map of temperatures obtained through the use of variable convective coefficients up to since in all points discretised using the FEM method, the temperature difference does not exceed 0.1°C, which corresponds to the minimum temperature difference that the thermal camera is able to estimate.
The algorithm that allows to obtain this result checks that the absolute error, e ¼ T calc;i À T exper;i , is \0.1°C. To simplify the problem, however, the heat sink has been divided into four zones that, from the thermographic measurement, showed different temperatures and, for these areas, were calculated the four values of the convective heat exchange coefficient. The algorithm implemented on Matlab allows to assign to the four areas, in the next step, the values of the convective heat exchange coefficients for which the whole temperature distribution presented a minimum RMS value.
The dimensions of the heat sink, those taken into consideration in the calculations of the convective heat exchange coefficient with the correlation, are reported in Table 1. The heat sink was divided into four critical zones (Fig. 7), in correspondence of which different convective coefficient values were used. The radiative effect was overlooked due to the low operating temperatures reached by the heat sink.
The upper and lower surfaces of the heat sink, the rounded summit of the fins, the rear surface and the internal fin surfaces which face the central cavity are indicated with S esp (Fig. 7a). S int indicates the internal surface of the fins, which includes all the interstices between the heat sink fins being subject to convective exchange of different entities and air speed different from the rest of the heat sink (Fig. 7b). The central cavity surface, S cav (Fig. 7c), is the zone in which the highest temperatures are recorded, while the external surface of the fins at the extremities is indicated with S estr (Fig. 7d).
This type of numerical analysis was carried out by Culham et al. [7] and [8], by Narasimhan [9], Higuera and Ryazantsev [10], Da Silva et al. [11] and Floryan and Novak [12], even with different hypotheses than those considered in this study. Moreover, in the literature, there are studies on natural convention in heat sinks varying the thickness of the fins and the width of the cavity (Kim et al. [13]) with unrealistic hypotheses, in that isothermal and insulated heating surfaces were used in order to obtain thermal conditions similar to those discussed in the literature and for which it is possible to obtain an analytical solution. In the literature, there are studies which characterise the performance of heat sinks in forced convention (Hussam and Axcell [14]), always based on the hypothesis of classical operating.
In scientific literature, the inverse problem of determining the convective heat exchange coefficient from temperature distribution has been widely investigated as well as reported in the references [15][16][17][18][19][20].
In references [15] and [16], the thermographic measurements were carried out with reference to a thin plate of copper. In Ref. [17] are compared various calculation methods that allow the measurement of the temperatures in transitory regime and further traced back to the convective heat exchange coefficient. In Ref. [19] are applied, with reference to the composite materials, calculation methods similar to those described above in the presence of  [20], the method of Tikhonov is used, in addition to the method of Nelder-Mead, for a two-dimensional case.

Experimental tests
Five experimental tests were conducted with powers to be dissipated of between 9.2283 and 30.861 W, while the dry bulb temperature of the environment air varied between 18.2 and 20.84°C. The time required by the heat sink to reach stationary conditions, from an initial starting temperature equal to the environment temperature, was on average 20 min, after which it was possible to acquire a thermographic image of the heat sink temperature map.
For the first test, a tension of 9.6 V was applied to the resistance, which simulates the electrical component, and a continual current equal to 2.12 A was supplied, obtaining a total power equal to 20.352 W. The heat sink took around 20 min to reach the stationary regime, cooled down until the dry bulb environment temperature equal to 18.2°C. The parameters used during the first test for the calculation of the initial heat exchange convective coefficients for the four surfaces are reported in Table 2.
The maximum temperature reached by the heat sink in the centre of the cavity is equal to t cav = 48.6°C (point 1), while the exposed surface reaches a temperature of t esp = 45.4°C (point 2). In the interstice between the two most distant fins from the centre, the temperature is equal to t int = 44.2°C (point 3), and on the lateral surface of the fin closest to the cavity, a temperature equal to 46.4°C (point 4) is observed. The average temperature of the heat sink resulted as being equal to t mean = 45°C (Fig. 8).
Utilising such temperatures for the calculation of the convective heat exchange coefficients in the ProENGI-NEER simulation, the map reported in Fig. 9 is obtained. It is possible to note that the temperature distribution obtained is very different from the experimental values recorded, with differences in the region of 10-13°C; therefore, the correlation used is not appropriate.
The convective coefficients obtained after several attempts and for each surface examined are the following: h esp = 8.6 (Wm -2 K), h int = 8.4 (Wm -2 K), h cav = 8.4 (Wm -2 K) and h estr = 8.6 (Wm -2 K). The final result obtained with ProENGINEER comparing it with the starting thermographic image is reported in Fig. 10.
During the course of the following days, four further tests were conducted using the same procedure as the first test. The most relevant data relating to the five tests, temperature of the four surfaces and the average heat sink temperature, are reported in Table 3.
The comparison between the experimental values (h exper ) of the convective coefficient and those calculated with the correlation (h calc ) for the four surfaces and for all the experimental tests is reported in Table 4.
The final column (h mean ) refers to the average values of the convective coefficient for the four surfaces and therefore for the entire heat sink.
Comparisons, for the different heat sink surfaces, between the theoretical convective coefficient calculated with the correlation and the experimental ones obtained by   Comparison of the manufacturer's data and the correlation proposed for average values The possibility of using a mean heat convective exchange coefficient on all the surfaces of the heat sink, in such a way as to approximate the different point values calculated, was evaluated at this stage. From the comparison between the temperature maps obtained using the point values of the heat exchange coefficient and those obtainable using its average value, with similar power and environmental temperature of the test in consideration, it is possible to observe how the variation in temperature which occurs locally compared to the effective values, only varies for some tenths of centigrade degrees. The average value of the convective heat exchange coefficient, for the five tests, is reported in Table 4.
The trend of the average values of the convective heat exchange coefficient in relation to the average temperature, both for experimental and theoretical values, is reported in Fig. 15. From the experimental trend values, it is possible to obtain the trend lines for the average convective coefficient in relation to the power dissipated: The power that is effectively dissipated by the electrical component can be calculated in relation to the experimental DT exper , with the correlation: Subsequently, the trend of the power dissipated in relation to the difference between the average temperature of the heat sink t mean , and the environment temperature t a was analysed: The value of the typical resistance towards the environment supplied by the manufacturer, valid in stationary conditions, is equal to R manufact = 1.1 KW -1 .  Therefore, the dissipated thermal power, indicated with Q manufact (Q manufact = DT exper /R manufact ), and the experimentally calculated dissipated thermal power, indicated with Q exper (Table 5), were calculated for the different tests. The thermal resistance values calculated experimentally (R exper = DT exper /Q exper ) are reported in Table 5.
The trends of the two powers are reported in Fig. 16. The power obtained with the resistance supplied by the manufacturer provides values for the dissipated power which are much higher than the total thermal power supplied to the system. Such conditions during the project phase would mean that an electrical component set on the heat sink could dissipate much more power.
It is important to specify that the temperature range investigated in this paper is such as to consider negligible the effects of radiative heat exchange and the use of a mean convective heat exchange coefficient; the case of the heat sink under consideration can be considered as a good solution due to the fact that for the calculation of the dissipated power, the local DT is not used but the mean DT, obtained from the mean temperature of the heat sink and the air temperature. Moreover, given the low variability of the convective heat exchange coefficient, even the areas at a higher temperature having small heat exchange areas dissipate power as the surfaces are with larger areas but with lower temperature gradient.    The uncertainty of the experimental measurement of the surface temperature T (x, y) of the measuring point is linked to the two variables x and y [21].
Indicating the absolute uncertainty of the temperature at a point P (x ? dx, y) with u x , and the absolute uncertainty of the temperature at a point P (x, y ? dy) with u y , was analysed the combined uncertainty on the temperature u T with the expression: The percentage uncertainty is given by the equation: The maximum percentage error of the camera is equal to l %, both along the x-direction and the y-direction, because the emissivity does not present spatial variations. Applying the Eq. (12), a value of the uncertainty on the measurement of the temperature T (x, y) equal to 1.41 % is obtained.
The error was defined with the relation: and the mean error as: Instead, the root mean square was defined with the following expression: With reference to the experimental and theoretical data relating to the dissipated powers reported in Table 5 and in the graphic in Fig. 16, a mean error of 23.4 % and RMS of 25.45 % are obtained.

Conclusions
The objective of the present work is that of providing the guidelines, methodology and the instruments for the observation and characterisation of the convective phenomena of heat dissipation of electrical components and how the complexity of thermal exchange follows trends, which are not linear. However, this is not just a case study but has as its goal the finding of a possible link between thermographic measuring techniques and the analysis of finite elements, which cannot solely be based on the use of correlations, but in specific cases, it almost always requires an experimental confirmation. In this case, the study of new surfaces and complex geometrical forms or the properties of new materials for the removal of heat becomes possible.  In the field of electronics, but also in other fields in which heat to be dissipated is involved, the main focus of research must not be the removal of heat produced in a manner that is collateral to the useful process, but rather its reduction with the aim of limiting thermal risks. In fact, the creation of components capable of working without an excessive production of heat would certainly eliminate costs linked to heat dissipation, but it would also be great progress from a viewpoint of environmental impact which, as is well known, has become the most urgent theme of the modern age in all aspects of daily life.
A further consideration regards the use of classical correlations from the literature, which when applied locally do not describe in an exact way the phenomenon of thermal exchange, indicating the flow conditions, even if it is evident from the trend of convective coefficient that after the third test, there is a discontinuity step in the calculation of the convective coefficient.
Moreover, since the values suggested by the manufacturer's characteristic curve are always greater than those recorded experimentally, it is inferred that in the planning phase, it is necessary to oversize the heat sink so as to avoid encountering overheating problems.
Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.