Abstract

An increase in the operating temperature of photovoltaic (PV) panels caused by high levels of solar irradiation can affect the efficiency and lifespan of PV panels. This study uses numerical and experimental analyses to investigate the reduction in the operating temperature of PV panels with an air-cooled heat sink. The proposed heat sink was designed as an aluminum plate with perforated fins that is attached to the back of the PV panel. A comprehensive computational fluid dynamics (CFD) simulation was conducted using the software ANSYS Fluent to ensure that the heat sink model worked properly. The influence of heat sinks on the heat transfer between a PV panel and the circulating ambient air was investigated. The results showed a substantial decrease in the operating temperature of the PV panel and an increase in its electrical performance. The CFD analysis in the heat sink model with an air flow velocity of 1.5 m/s and temperature of 35°C under a heat flux of 1000 W/m2 showed a decrease in the PV panel’s average temperature from 85.3°C to 72.8°C. As a consequence of decreasing its temperature, the heat sink increased the open-circuit photovoltage () and maximum power point () of the PV panel by 10% and 18.67%, respectively. Therefore, the use of aluminum heat sinks could provide a potential solution to prevent PV panels from overheating and may indirectly lead to a reduction in CO2 emissions due to the increased electricity production from the PV system.

1. Introduction

The use of renewable energy resources is of interest to researchers and governments around the world due to increasing energy consumption and climate change issues caused by the exploitation of conventional energy sources [1]. Solar energy is the most abundant renewable energy resource on Earth and could therefore be the solution for the growing demand of global energy consumption [13]. In general, the utilization of solar energy is divided into two methods: photothermal and photovoltaic (PV). Photothermal systems utilize the heat energy from solar radiation for various purposes such as crop drying, solar stoves, and solar water heaters [4]. The PV system converts photons from solar radiation directly into electrical energy using solar cell technology [5]. Solar cell technology is widely applied in both small-scale uses, such as street lighting and providing residential electricity, and large-scale uses such as in national power plants [57].

Although solar energy harvesting using solar cell technology is highly promising, there are challenges to maximizing the efficiency of solar cells [810]. PV modules based on silicon material can convert 8–20% of solar radiation into electrical energy [1113] but the rest of the solar radiation is reflected back to the surrounding environment and converted into heat energy that can increase the temperature of the PV device and consequently reduce the total power output [1416]. Radziemska [12] adjusted the working temperature of solar cells and showed an inverse linear relationship between the temperature and the power output of solar cells. Cotfas and Cotfas [8] found that increasing the temperature of solar cells would degrade the voltage because the electron excitation of the thermal energy was higher than the electrical energy of the semiconductor material. As a result, the efficiency of solar cells with crystalline silicon structures decreases by ±0.4–0.5% for every 1°C increase in temperature [17]. Hence, decreasing the operating temperature of solar cells will be important for achieving a higher efficiency and longer lifespan of these cells [18].

Researchers have developed cooling systems to decrease the operating temperature of solar cells; these are categorized as active and passive cooling systems. Active cooling systems operate mechanical or electrical devices such as fans or pumps, which require external power input, while passive cooling does not require additional power to operate [1921]. Passive cooling may require additional parts, such as a heat pipe or heat sink, that rely on convective heat transfer to drive coolant flow to the solar cells [22]. A heat sink with thermal conductive material attached to the bottom of a solar cell will increase the area of heat transfer from solar cell to its surrounding environment [17, 22]. Thus, because they are relatively simple and inexpensive to manufacture, heat sinks have a high potential as devices to cool PV panels and should be developed further.

Popovici et al. [23] conducted a simulation using the ANSYS-Fluent software to study the characteristics of heat transfer in a heat sink under turbulent flow conditions. The results showed that the cooling rate of a solar cell was proportional to its fin height and inversely proportional to the configuration of the inclination angle. The temperature of the solar cell was reduced to 10°C and its electrical power capability was increased to 7.55%. Modeling and simulations of heat sinks conducted by Zhu and Sun [24] showed that a simple heat sink can affect the longevity of semiconductor material. Research by Luo et al. [25] used an experimental method to demonstrate several heat sinks with flared fin configurations and obtained an overall reduction in thermal resistance of up to 10%. Research by Grubišić-Čabo et al. [26, 27] examined the aluminum fins mounted on the backside of a PV panel (Si-poly, 50 Wp). The results showed that, under low wind conditions, the electrical efficiency of solar cells increases by 0.3% and 0.2% under solar irradiation of 850 W/m2 and 500 W/m2, respectively. Based on the various studies currently being conducted on heat sink applications for PV solar cells, there is a great potential to combine simulation and experimental methods to produce an integrated study. This applied research will aid the development of PV cooling systems by providing a complete theoretical and analytical overview of the methods to decrease the temperature of solar cells.

This study investigated the application of heat sinks to a PV module performance with a simple combination of the computational fluid dynamics (CFD) approach and experimental testing. The CFD study was conducted using the ANSYS Fluent software, and the three-dimensional (3D) modeling was carried out to identify the working temperature of the PV module. The results of the simulation were presented as the temperature contour of the module and the velocity profile of the air flowing through the heat sink. The air flow characteristics shown by the velocity profile could provide valuable information regarding the relation between turbulent air flow and the cooling performance of the panel by using an air-cooled heat sink. Additionally, an experimental analysis was performed to study the effects of heat sinks on the electrical characteristics of PV modules. This study also discusses comprehensive physical phenomena in the silicon material of solar cells caused by temperature changes. This heat sink cooling system is expected to decrease operating temperature effectively and increase the efficiency of PV modules.

2. Materials and Methods

2.1. Numerical Setup

Figure 1 shows the design and geometry of a PV module with an aluminum heat sink. The aluminum heat sink was mounted on the back of a vertical solar panel; the fins of the panel were perforated to improve air circulation around them and allow the absorption of more heat from the PV panel. In the modeling program, the PV panel was assumed to be a unique composite layer [2830]. Table 1 shows the properties of each layer in the solar module. The properties of these layers are assumed as independent properties with temperature and pressure change.

The base and fins of the heat sink are both 2 mm thick. The fluid domain for air has a width of 0.1 m, with the resulting hydraulic diameter () of 0.166 m. The air flow velocity () at the inlet was 1.5 m/s, with a temperature of 35°C according to the average temperature and velocity of wind in Indonesia [3133]. The Reynolds (Re) number and turbulence intensity () were calculated using Equations (1) and (2) below: where and are density (kg/m3) and viscosity (kg/m s) of fluid, respectively [23]. The Reynolds number for the imposed air velocity was approximately 13100 and the estimated turbulence intensity was 4.8%. The simulation was carried out under steady-state conditions using the k-ε re-normalization group (RNG) turbulence model. The semi-implicit method for pressure-linked equation (SIMPLE) pressure–velocity coupling and the second-order upwind scheme were used to solve the equations, with convergence criteria of 10–6 for energy and 10–3 for pressure, velocity, and continuity equations.

The PV panel was irradiated with 1000 W/m2 of solar energy in standard test conditions; it converted this into electrical energy through the mechanism of PV effects [34, 35]. In general, the electrical efficiency of the crystalline silicon solar cells ranges from 11 to 20%, while tempered glass reflects as much as 3 to 10% of the solar radiation to the surroundings of the cell. This study assumed that 25% of solar energy is converted into electrical energy and reflected to the environment; therefore, the remaining solar radiation was treated as input heat flux to the solar module [23, 35].

2.2. Experimental Procedure

This study used 50 Wp polycrystalline solar cells with the dimensions of  mm. The experimental set up is shown in Figure 2, with the aluminum heat sink mounted on the back of a PV panel with fins and a base that are both 2 mm thick. The fins were connected to the baseplate with a rectangular hole to accommodate the PV panel junction box. Before attaching the heat sink to the bottom of the solar panel, the thermal grease HT-GY260 ( W/m K and  C-in2/W) was applied on the contact areas to minimize thermal contact resistance.

The electrical characteristics of the PV panels were measured using a variable resistance method. This method provides an electrical resistance load that varies from zero resistance (to measure the maximum current of a solar cell) to infinite resistance (to measure the maximum voltage of a solar cell). The variable resistance was achieved by arranging several resistors with electrical resistances of 2.5, 3.5, 4.7, 5.4, 5.8, 6.0, 6.4, 6.6, 6.9, 7.4, 8.5, 13.8, 19.6, 42.5, 111, and 330 Ω adjusted to the size of the PV panel. A LI-COR LI-250A pyranometer with an accuracy of ±0.6% was used to measure incoming solar irradiation. The temperature of the PV panel was measured with a K-type thermocouple with an accuracy of ±1.2°C. The thermocouples were placed on top of the PV panel to measure its average temperature. The wind speed passing through the underside of the PV panel was measured using an anemometer. The position and distance between the 35 W fan blower and the PV panel was adjusted to obtain a uniform wind speed of approximately 1.5 m/s. Figure 3 shows a photograph of the experimental procedure and heat sink configuration on the back of the PV panel.

3. Results and Discussion

3.1. Computational Fluid Dynamics (CFD) Analysis

In the CFD study, the meshing process divided domains into many cells. The cell size varied from big to small such that the meshing process yielded different numbers of elements. This was done to generate independent simulation results not influenced by the number of the cells in the mesh. Meshing was performed with 1 and 2 mm cell size for heat sinks and 4, 6, and 8 mm cell size for the fluid domain. The mesh independence study result had an error of 0.13%.

Figure 4 shows the temperature contour of the PV panel with and without aluminum heat sinks with a heat flux of 1000 W/m2. The average temperature for the PV panel without the aluminum heat sinks was 85.3°C and the average temperature for the PV panel with the aluminum heat sinks was 72.8°C. The results showed a significant decrease in temperature down to 12.5°C, and an increase in the heat transfer performance from the panel to the air when aluminum heat sinks are used.

Figure 5 shows the velocity vector of the circulating air which flows through the heat sink. Fins in the heat sink caused turbulence of air flow, so the heat transfer from the PV panel to the environment increased. The high turbulence level can reduce the thermal boundary layer thickness and have an impact on the high temperature gradient in the back of the PV panel. As a result, the average temperature of PV panels with heat sinks was lower than in PV panels without heat sinks. Research by Omeroglu [36] revealed a similar result where the heat sink configuration affected the air velocity and heat transfer rate from a PV panel to the surrounding environment. Popovici et al. [23] demonstrated that the orientation and direction of the fin in the heat sink also affected the air temperature distribution. The 0.03 m height of the fin and 1.5 m/s air inlet velocity decreased the base temperature of the PV panel from 56°C to 42.35°C ( up to 13.65°C), which is similar to this result.

3.2. Experimental Results

Figure 6 shows the average temperature of a PV panel as a function of solar irradiation. As previously noted, the average temperature of a PV panel without a heat sink was higher than that of a PV panel with a heat sink. We also observed that increasing the intensity of solar radiation would consequently increase the temperature of a solar cell. The higher the intensity of the radiation, the greater the amount of photon energy that hits the solar cell. In turn, a higher amount of photon energy will increase electron excitation in solar panels, which will then result in increasing the temperature of the solar cell [12, 37]. However, the installation of a heat sink increases the heat transfer area and heat transfer rate of a PV panel, thereby reducing the temperature of the panel.

The performance of solar cells can be characterized by a photocurrent-voltage curve (I-V curve) as shown in Figure 7. The measurement results of the I-V curve produce several important parameters, including open-circuit photovoltage (), short-circuit photocurrent (), fill factor (FF), and efficiency (). This study used a variable resistor system to measure the solar cell characteristics.

is determined at the open-circuit condition or when electrical resistance is very high, with no flow of current and the voltage at maximum, whereas is determined during the short-circuit condition or when electrical resistance approaches zero and generates a maximum current. The product of current and voltage produces electrical power (), and the maximum product of current and voltage is called the maximum power point (). Another important parameter of solar cell performance is the fill factor (FF), which is the ratio of to the product of and . The efficiency of solar cells () is defined as the ratio of to the solar irradiation energy received by the cells. The electrical characteristics of a PV panel are shown in Table 2.

The heat sink increased the by 10% and reduced the by 1.18%. As the increase in was higher than the decrease in , the electrical power of solar cells also increased. The from solar cells increased by 18.67%, with the improvement in the electrical performance being primarily caused by the decrease in the operating temperature. The electrical efficiency increase gained with the heat sink could reach 2.6% under 1000 W/m2 solar irradiation. Research by Grubišić-Čabo [26, 27] showed a similar trend of increasing efficiency with aluminum fins mounted on the backside of a PV panel (Si-poly, 50 Wp).

Semiconductor properties in solar cells are similar to temperature-sensitive electronic devices. Increasing temperature can decrease the band gap value of the semiconductor, which in turn increases the energy of the electrons in the semiconductor material. As a result, the electron bonds in the semiconductor will be broken easily only by the entry of energy in low quantities. The parameter that is most affected by fluctuations in temperature values is [38, 39]. decreases as operating temperature increases due to the influence of the reverse saturation current () parameter, as shown in Equations (3), (4), and (5). where is the intrinsic carrier concentration (cm–3), is the Boltzman constant (1.381 × 10–23 J/K or 8.617 × 10–5 eV/K), is the temperature (K), is the electric charge (1.602 × 10–19 Coulomb), is the current of the solar cell, is the area of the solar cell, is the diffusivity of the minority carrier for silicon as a function of doping material, is the length of diffusivity of the minority carrier, and is the doping material. All of these parameters could change with an increase in temperature, but the most significant change is in the intrinsic carrier concentration as shown in Equations (6) and (7): where is the band gap energy, is Planck’s constant (–34 J.s), is the effective mass of the electron, and is the effective mass of the hole. The constant is not a function of temperature; therefore, the parameter depends on the operating temperature of the semiconductor. Semiconductors with low band gap values will produce a high intrinsic carrier concentration; thus, they will also produce high values at high temperatures. The intrinsic carrier concentration value of silicon material can be determined as a function of temperature as shown in Equation (8) [40]:

Based on Equations (3) and (4), a high value produces a large reverse saturation current () value, thereby reducing the value of the solar cell. Additionally, Equation (5) indicates that the increase in the values of and will result in an increase in the current. However, the increase in the current is negligible compared to the decrease in voltage from the solar cell [41].

In addition to the changing current and voltage values, changes in the operating temperature of the solar cell also affect the fill factor (FF). The rise in temperature increases the internal resistance of solar cells, and these obstacles contribute to the charge flow at the p-n junction, resulting in a loss in the form of charge recombination in the cell. This decreases FF when the operating temperature of solar cells increases. Therefore, because the output power of a solar cell is the product of voltage and current, an increase in temperature will result in a decrease in power [11, 12].

4. Conclusions

We conducted a CFD simulation for PV panels with and without aluminum heat sinks installed. The results showed a reduction of up to 10°C in the average temperature of the PV panels with a heat sink. A physical experiment was also conducted with a PV module that had a heat sink installed, and various values of solar irradiation were applied to PV module to observe their influence on the temperature distribution of the PV panel. The results showed that the installed heat sink could reduce the panel temperature due to the increased heat transfer area and heat transfer performance. Because of its decreasing temperature, the heat sink increased the and of the PV panel by 10% and 18.67%, respectively. The simple design of this heat sink model provides a potential solution to prevent PV panels from overheating and may indirectly lead to a reduction in CO2 emissions due to the increased electricity production from the PV system. However, the analysis could be developed further from the economic perspective, to produce a cost-based and performance-based optimum design for heat sinks.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors state that there are no conflicts of interest to declare.

Acknowledgments

This research was funded by the Direktorat Penelitian dan Pengabdian kepada Masyarakat (DP2M), Direktorat Jenderal Pendidikan Tinggi (DIKTI) (research grant number 719/UN27.21/PN/2019). The authors are grateful to Ian Guardian for his assistance during the research.