1. Introduction
The current global energy context is characterized by the energy consumption of humanity, which has reached impressive values in recent years. Thus, in the current year (2021) the well-known Overshoot Day was “celebrated” on 29 July. This meant that humanity had already used up all of the restorable resources of 2021. From 30 July we started consuming more resources than the planet could regenerate in a year [
1]. Therefore, by the end of the year 2021, the energy consumption is projected to reach about 1.73× Earth’s resource regeneration capacity. It is a global emergency and every small step in reducing this tendency to over consume energy will become decisive in this fight.
Solar energy is a clean and sustainable source of energy. It is involved directly or indirectly in the vast majority of other renewable energy sources, and it is also related to the conventional fuels—as most of them are the result of photosynthesis. The solar energy incident on the Earth could cover the annual energy demand for the entire world in only 1.5 h [
2,
3].
Regarding the conversion of solar radiation into electricity, this can be achieved by using photovoltaic (PV) cells or it can be generated through solar thermal technologies [
1].
Photovoltaic conversion has important advantages: the energy source is renewable and free, and the processes of electricity production are environmentally friendly (in terms of CO2 and NOx emissions, waste, noise etc.). The generation of electricity is achieved without the use of moving parts, so the cost of maintenance is reduced. Furthermore, the energy produced is consumed locally, which results in reduced power losses.
By implementing efficient PV systems, the safety of the users of the electricity supply is increased [
4,
5]. These systems also reduce the vulnerability of European space and could eliminate future energy crises or the lack of an energy supply. Considering the long service life and the low CO
2 emissions of PV systems, they represent a sustainable alternative to electricity production. In addition, the energy market in Europe, which lately has been affected by internal and external factors, could become more stable and competitive by increasing the power it produces and the number of users of photovoltaic systems [
6].
The main environmental parameters, which are influencing the efficiency and the lifetime of PV panels are solar radiation, wind velocity and direction, dust, ambient temperature and their operating temperature. Most of these meteorological parameters are random and uncontrollable. From the point of view of the installation methods, the tilt angle, orientation and positioning are very difficult to manage in the case of typical PV systems [
7,
8]. However, the management of the operating temperature of PV cells has been demonstrated to be a good solution for enhancing their efficiency [
9].
Almost 80% of the solar energy incident on the PV panel’s surface is converted into heat, and it is well known that the performance of PV panels is temperature-dependent [
9,
10]. Skoplaki, E. et al. [
10] analyzed different methods and relations for calculating the dependence between conversion efficiency and photovoltaic cell temperature; a linear dependence between efficiency and temperature has been proved. The decrease in the efficiency of the produced power related to the increase in the temperature has values between −0.46–−0.50%/°C and −0.47–−0.50%/°C, respectively [
11,
12,
13].
Taking into account that the power of a photovoltaic cell is influenced by the temperature and radiation levels, standard test condition (STC) parameters were defined:
tPV = 25 °C,
G = 1000 W/m
2, AM1.5 [
14]. In STCs the PV panels produce the watt peak power [Wp], but a real concern is that in regular operation, at 1000 W/m
2, the photovoltaic panels can reach extreme temperatures of 70–80 °C [
11,
15,
16], leading to a significant decrease in the efficiency.
Complementary cooling solutions are recommended in order to achieve high productivity [
17,
18,
19,
20,
21]. The literature presents a large number of papers on different solutions for improving the conversion efficiency by decreasing the temperature of the photovoltaic panels [
13,
15,
20,
22,
23,
24,
25,
26,
27].
Concepts frequently used in analyzing PV systems are the NOCT (nominal operating cell temperature) conditions. In contrast to STCs, the NOCT conditions are closer to reality, taking into account the following operating conditions: tilt of PV panel = 45°;
G = 800 W/m
2;
tair = 20 °C;
vair = 1 m/s towards the backside of PV panel [
15]. According to these conditions, every producer must supply, via the data sheet, the average operating temperature of the PV panel—
tNOCT.
The aim of this study was to overcome technological and conceptual barriers that were currently attributed implicitly to photovoltaic systems by improving their efficiency during months with maximum energy production but lower efficiencies due to high operating temperatures. For example, in hot summer days it turned out that the most important cause of the production potential not being reached, when the intensity of solar radiation had the highest values of the year, was the uncontrolled increase in the operating temperatures of the photovoltaic panels.
The improvement in the efficiency of PV panels when air-cooled heat sinks were used, was reasonable for the residential sector, urban areas and power plants [
28,
29].
The novelty of the present study consisted of analyzing the impact of longitudinal and transversal heat sinks with perforated and non-perforated fins on the performance of PV panels. The impact was analyzed using more realistic wind conditions and tilt of the panel, compared to existing studies in the literature.
2. Related Research
Experimental tests and numerical analyses on monocrystalline or polycrystalline photovoltaic panels show that the decrease in efficiency, as the operating temperature increases, and the necessity of cooling solutions are current concerns in the literature [
11,
12,
13,
16,
30,
31].
Different technologies to improve the efficiency of photovoltaic panels by passive cooling, using air, are studied in the literature [
8,
13,
32,
33,
34,
35,
36,
37,
38,
39]. Detailed studies on air cooling of photovoltaic panels are carried out in the following works [
11,
40]. Air cooling with heat sinks and perforated fins is studied both numerically [
11,
32] and experimentally [
34]. The characteristic values for multiple photovoltaic panels cooled with aluminum and copper heat sinks were compared through experimental tests under different conditions [
34]. The results of the simulation and from the experiment showed that the heat sink with perforated fins can significantly reduce the temperature of the panel and increase its performance; there was an increase in the heat extracted from the photovoltaic cell and a decrease in its temperature, with values between 7 °C and 10 °C [
34].
Also, numerical simulations on a particular case of backside ventilation [
11,
16,
41] revealed that the temperature was reduced by 10 °C below the value obtained in the base case. This aspect was favorable for conversion efficiency, determining a generation of approximately 90% of the maximum power produced in STCs [
14], while the average power production of PV panels was about 70–80% of
Pmp in STCs [
42].
3. Materials and Methods
3.1. Problem Description
The goal of this work was to develop a numerical model regarding the passive cooling of the PV panels, through perforated and non-perforated heat sinks, and determine the enhancement of operating temperatures compared to the base case (without cooling). A photovoltaic panel with an area of approximately
S = 1.6 m
2 was studied in a fixed position and tilted at 45° from horizontal, as well as in variable solar radiation, wind velocity and ambient temperature. This setup, together with the wind direction towards the backside of the PV panel, was established according to NOCT conditions [
14] in order to validate the model.
A common monocrystalline Si-based PV panel (AE Solar AE320HM6-60), with an advantageous position of the junction box, was analyzed with the following technical data (
Table 1) [
43].
A step-by-step analysis of heat sink usage for cooling PV panels was conducted. After the initial validation of the base case—PV panel without heat sink—a total of six photovoltaic panels with attached heat sinks were analyzed. Two of them were non-perforated, while the other four had perforated fins.
The heat sink consisted of a metal plate with copper attached to the photovoltaic panel and perforated or non-perforated metal fins,
Figure 1. The role of the designed holes on the fins was assumed to be both quantitative (by reducing the amount of material necessary for device construction) and qualitative (by amplifying the turbulence and heat transfer in the immediate vicinity of the PV panel).
Due to preliminary analysis and the literature review [
11,
32,
34,
41], the 30 mm and 60 mm diameter circular perforations distributed on a single row turned out to be feasible. It emerged that other shapes, dimensions or multiple rows of perforations were heading towards lower performance of the heat sink (the operating temperatures were at least 2–3 °C higher). It must be specified that the literature comprises a diversity of heat sinks with different approaches of fins [
42,
44]. The present study is focused on using heat sinks with continuous fins.
Two types of heat sinks were analyzed:
Type 1—heat sinks with horizontal (transversal) fins of 100 mm height (with 3 subtypes: non-perforated, Φ = 30 mm perforations, Φ = 60 mm perforations);
Type 2—heat sinks with vertical (longitudinal) fins of 100 mm height (with 3 subtypes: non-perforated, Φ = 30 mm perforations, Φ = 60 mm perforations).
3.2. Numerical Simulation
The numerical simulation was performed using the ANSYS-Fluent 2021-R1 software. The geometry and discretization of the computing domain was performed within the platform using the Fluent-Meshing and SpaceClaim programs.
For simulation the control volume method, SIMPLE pressure–velocity coupling and second-order accuracy for equation solving were used, taking into account that the flow was not aligned with the mesh [
45]. This method used a relationship between velocity and pressure corrections to enforce mass conservation and to obtain the pressure field. A fully implicit numerical scheme was employed, in which upwind differences were used for the convective terms and central differences for the diffusion terms [
11]. The calculation was iterative, with a convergence criteria of 10
−6 for the energy equation and 10
−3 for the pressure, velocities,
k,
ε and continuity equations.
The simulation of solar radiation was realized using the Solar Ray Tracing model of the Fluent program. The absorption of solar radiation by the front of the PV panel was imposed as a global coefficient α = 0.7, taking into account all the layers of the PV panel [
11,
46].
The heat transfer balance was set in order to include the convection and radiation with the surroundings [
46]. The numerical simulation was realized using the
k–ε realizable model [
45]. For all studied configurations, the distance between fins was set to 150 mm. The fins were perpendicular to the heat sink plate—
Figure 1. For the case of perforated fins, they had circular holes of 30 mm or 60 mm diameter, which were placed at a distance of 100 mm to one to another [
34]. Their role was to improve the air circulation near the heat sink and to extract more heat from the PV panel.
The accuracy and the solution time were both highly depend on the characteristic of the mesh. Computational fluid dynamics (CFDs) solvers depended on the highly orthogonal mesh that could be difficult to achieve on different types of geometries. For the present study, a polyhedral mesh was used for the surfaces and a poly-hex core for one of the volumes. The choice of these types of cells was made due to an increase in speed of simulation from 20 to 50% while maintaining the high level of accuracy. In addition, this approach in the simulation opened up the possibility of improving the solver speed and accuracy compared with the tetrahedral cells.
The mesh was realized with different refinements for the heat sink and air volume. The mesh quality achieved for the numerical analysis was as follows:
- -
Aspect ratio was in the range 5–60;
- -
Skewness was in the range 0.4–0.7;
- -
Orthogonal quality was in the range 0.1–0.3.
A mesh independence study was carried out, taking into account the large number of elements that resulted due to modeling the air volume that flowed over the PV panel and heat sinks.
Figure 2 presents the mesh independence study, which was achieved for approximately 3 million cells in the base case and over 5.5 million cells for cases with heat sinks and perforated fins.
The numerical analysis was conducted on a workstation that had 48 cores, 256 GB of RAM and a graphical computing video card Nvidia Quadro P6000.
Typical approaches from the literature considered the position of the PV panel as either horizontal [
31] or vertical [
11,
34], with a channel simulating the wind behind the PV panel. There were only a few studies on cooling tilted PV panels [
16,
40].
During this study the modelling was realized by taking into account the air circulation in the vicinity of a tilted PV panel cooled with heat sinks. In this way, by using an optimized resource consumer mesh, the airflow and the heat transfer were more realistic and easier to visualize and understand. The material properties used in the simulations are shown in
Table 2, and the sketch of the model is presented in
Figure 3.
During the study, the non-perforated heat sinks were compared with perforated heat sinks in order to determine the optimum configuration for implementation in NOCT conditions of wind.
Simulations were realized in steady state conditions with the following input parameters:
- -
ambient temperature, tair = 20, 25, 30 and 35 °C;
- -
velocity of the air towards the back of the photovoltaic panel, vair = 1–5 m/s;
- -
the normal component of solar radiation varied between 200–1000 W/m2, assuming the typical conditions of operation for the PV panel.
It should be mentioned that heat sinks were studied on some of the worst scenarios of temperature operation for the base case. Therefore, some cases with low levels of solar radiation, low ambient temperatures or high wind velocities were excluded from the results presented during the paper.
The main outputs of the simulations were:
- -
average temperature of the PV panel;
- -
convective heat transfer coefficients, achieved for both the front and the back of the PV panel or heat sink.
Studied Cases
A total of six configurations of heat sinks were analyzed, plus the base case, these are shown in
Table 3.
The position of the PV panel and air circulation was set according to NOCT conditions: the PV panel was tilted at 45° and wind direction toward the backside of the PV panel (
Figure 4 and
Figure 5).
Heat sinks consist of a metal plate and perforated or non-perforated longitudinal or transversal fins arranged at equal distances of 150 mm from each other (
Table 4). Both the base plate and the fins were made of copper. The base plate was in direct contact with the back of the photovoltaic panel. The fins had circular perforations with a diameter of 30 or 60 mm (
Table 4).
During the study, the fins were arranged longitudinally (vertically—parallel to the length of the PV panel) or transversally (horizontally—parallel to the width of the PV panel). The distance between fins was s = 150 mm, while their height was hfins = 100 mm.
The approximate length of the fins was
Lfins = 1520 mm (for longitudinal ones) and
Lfins = 860 mm (for transversal ones). The number of fins was 11 for the horizontal case and 5 for vertical case (
Table 4).
The perforations had a circular shape and a diameter of Φ = 30 mm or Φ = 60 mm. The distance measured on the fins between the centers of the perforations was d = 100 mm. The perforations were staggered by 50 mm between two successive fins.
The back, side, top and bottom view of the heat sinks attached to the PV panel are presented in
Table 4. In addition, the main components of the geometry conceptions of the photovoltaic panel, mounting system and heat sinks with perforated and non-perforated fins were easy to follow in 3D perspective views.
5. Results and Discussion
Most simulations in the literature were realized by imposing a constant convective heat transfer coefficient (
hc), which determined a quasi-uniform temperature distribution [
50]. In the present work, the simulation of the airflow near the PV panel determined a more realistic temperature distribution on its surface compared with the typical approaches [
42], while the convective heat transfer coefficient was varying on the PV panel’s surface.
The results of the simulations consisted of the photovoltaic panel temperature and convective heat transfer coefficient variation, depending on the type of the heat sink. The spectra of the operating temperatures of the PV panels and velocity path lines on the cross-sectional plane were obtained (
Table 5). The following images present the qualitative visualization of the results for cases simulated in NOCT configuration:
tair = 20 °C,
vair = 1 m/s,
G = 800 W/m
2 (
Table 5).
The visual results clarify the temperature distribution achieved for each case analyzed, together with the air velocity field on the cross-sectional plane. The effect of the heat sinks was determining a decrease in maximum operating temperature of approximately 10 °C. In addition, from a qualitative point of view, for this particular wind direction (towards the backside of the PV panel), the air turbulence was increased compared to the base case, determining an improved heat transfer and a better cooling of the PV panel. In addition, each type of heat sink determined a slightly different distribution of wind velocities and recirculation areas. It can be remarked that the cooling effect was decreased in the zones with air recirculation in front of the PV panel compared to the rest of the area.
The results for
vair = 1 m/s, various inlet temperatures (
tair = 20–35 °C) and solar radiation (
G = 600–1000 W/m
2) were merged in the following plots:
Figure 8,
Figure 9 and
Figure 10.
An interesting aspect was recorded on convective heat transfer coefficients. Therefore, for horizontal heat sinks the
hc on the front of the PV panel were much closer to the
hc on the back, compared to vertical ones, where an important difference was recorded. The ambient temperature variation did not seem to change the best configuration (Case 2—horizontal 30 mm). This configuration determined a drop of temperature as presented in
Table 6.
The results showed that for the studied wind conditions, the perforations did not improve the temperature of the PV panel when the fins were vertical. In addition, in the case of vertical fins, the larger perforations (Φ = 60 mm) were better, while for the horizontal fins the smaller ones (Φ = 30 mm) were better.
In terms of efficiency and power production, according to the assumption that the efficiency dropped linearly with the rise in temperature of −0.37%/°C [
43], the following power production could be achieved—
Table 7.
The power production was calculated as the percentage of maximum power possible, by assuming that at 25 °C the ratio achieved was 100% (320 Wp in STCs). From here on, we considered that the rise of temperature determined a drop in power production for each degree over
tSTC = 25 °C, according to the PV panel’s data sheet [
43].
The most important effect of air cooling heat sinks was achieved in low wind conditions and high levels of solar radiation. It could be observed,
Table 7, that for
vair = 1 m/s,
G = 1000 W/m
2 and ambient temperature
tair = 35 °C, the percentage of maximum power production achieved was 83.33% for the base case, while the best cooling scenario (Case 3—horizontal 30 mm) could reach 88.74%. This meant a rise in power production of 6.49% compared to the base case.
For the position of the PV panel and wind direction studied in the paper, the best solution of heat sink for various conditions was the following: Case 3—Horizontal 30 mm. For this type of heat sink, the enhancement of power production was also reached for higher wind velocities (2–5 m/s), but with lower impact: from 3.60–4.97%, compared to the base case.
When lower ambient temperatures were involved, the effect of the heat sinks was decreased, this was due to inferior operating temperatures of the PV panel that were reached for base case. Therefore, the improvement of the power production reached values in the range 3.31–4.99%.