Numerical Simulation of Improving The Efficienty of Photovoltaic
Numerical Simulation of Improving The Efficienty of Photovoltaic
C.-G. Popovici a, G. Filipescua , N.-C. Cherecheș a*, S.-V. Hudișteanu a, E.-F.Țurcanu a, Ana
Diana Ancas a , Marina Verdeș a, P. V. Roșua
a
Technical University "Ghe. Asachi" of Iasi, Faculty of Civil Engineering and Building Servicies, Department Building Sevicies, D.
Mangeron 67 str.,700050, Romania; E-mails: [email protected], [email protected],
*[email protected], [email protected], [email protected],
[email protected], [email protected], [email protected]
DOI: 10.2478/jaes-2024-0010
KEY WORDS: green energy, photovoltaic thermal panels, numerical simulation, Ansys Software.
ABSTRACT:
The need to optimize the operation of photovoltaic modules inevitably arises with the development of green energy production
technology. In order to achieve a good technological yield, durability and efficiency in production, continuous studies and
innovations are required. This study focuses on simulating the operation of water-cooled and uncooled PV modules in order to
understand the temperature-dependent PV operation. This cooling module consists, in the first phase, of a copper coil through which
water circulates, and in the second phase, of a coolant distributor/collector system. The module is attached to the lower surface of the
photovoltaic panel, respectively to the teller foil layer. For this simulation we used the Ansys software package (Discovery, Fluent
and Space Claim).
viable solution, different means of addressing this temperature Meshing phase, where the calculation nodes and points cannot
issue must be found, resulting in an overall increase in conversion exceed 500,000 nodes. Generally, this is not a problem; the
efficiency (Wang Y. et.al. 2006; Bhubaneswari P. et. al.2011). simulation can proceed smoothly with a level of calculation
detail accepted in the research field. The geometry was created
using the SpaceClaim program, which is integrated into the
2. CASE OF STUDY same Ansys package.
The present study aims to calculate the temperature distribution Fig. 1 and Fig. 2 represent a cross-sectional view of the
on the surface of a photovoltaic panel under fixed ambient photovoltaic panel created using the SpaceClaim software.
temperature and solar radiation conditions. The solar radiation
has been chosen to be 1000W/m2K. Two types of active
cooling methods were used: one with a serpentine cooling
system and the other with a distributor/collector cooling system.
Additionally, the simulation also includes the baseline scenario
without any cooling. The inlet temperature of the cooling agent
is set to be 20°C, and the thermal agent's circulation velocity.
The system was simulated at 0.1; 0.3; 0.5; 0.7; and 0.9 m/s for
each of the two cooling systems, as well as for the non-cooling
scenario. This resulted in 11 simulations that were compared to Figure 1. PV Cross-Section, 3D Top View
understand which type of active cooling is more suitable under
these conditions compared to the system without cooling. The
data extracted from these simulations were provided a clearer
understanding of this phenomenon. The PV surface has a size of
2m² and was subjected to uniform solar radiation over its entire
surface.
For this simulation and numerical analysis study, a standard Fig.3 and Fig.4 represent the view in X-Z axis view of the PV
photovoltaic panel with dimensions of 2000/1000mm was used, with Serpentine Cooling respectively the PV with
composed of the following layers: Distributor/Collector Cooling.
Cover and protective layer made of glass
EVA layer 1 (ethylene-vinyl acetate)
PV cells
EVA layer 2
Tedlar film (polyvinyl fluoride (PVF) film)
82
JOURNAL OF APPLIED ENGINEERING SCIENCES VOL. 14(27), ISSUE 1/2024
ISSN: 2247-3769 / e-ISSN: 2284-7197 ART.NO. 398 pp. 81-88
Figure 6. View of the Component Layers of PV Figure 9. Geometry and Dimensions of the Distributor/Collector
Cooling System
The serpentine cooling system consists of a sequence of 9
copper pipes with a diameter of Dn10mm, Fig. 8. The distance
between the pipes is 101mm, and they are consecutively 3. NUMERICAL SIMULATION
connected to each other through bends with a radius of 55mm.
The total length of the serpentine is 18m. The numerical simulation was performed using the Ansys
Workbench 2023-R1 program, student version. The geometry
was created with the same CAD program, specifically the three-
dimensional (3D) subdivision in SpaceClaim, followed by the
Meshing phase, settings, solution, and results carried out with
Ansys Fluent, a computational fluid dynamics (CFD) tool
capable of physical modeling and simulation.
In this project, the numerical simulation for the operation of the
same photovoltaic panel includes 3 cases:
83
JOURNAL OF APPLIED ENGINEERING SCIENCES VOL. 14(27), ISSUE 1/2024
ISSN: 2247-3769 / e-ISSN: 2284-7197 ART.NO. 398 pp. 81-88
Photovoltaic panel 2000/1000 mm without cooling system 0,5 20 34,8 32,7 42,7 37,1
VH2O TH2Oint TH2Oout TminPV TmaxPV TmedPV 0,7 20 30,6 30,5 37,5 33,5
(m/s) (°C) (°C) (°C) (°C) (°C) 0,9 20 28,33 29,1 34,6 31,5
NO NO NO 64,5 112,3 88,4 Table 3. Data from Numerical Simulation of PV with
Serpentine Cooling
Table 2. Data from Numerical Simulation of PV Without
Cooling In Fig. 11, the temperature contour on the surface of the
photovoltaic panel can be observed for the color spectrum set
In Fig. 10, the temperature contour on the surface of the within the range of 20...110°C. The water circulation velocity
photovoltaic panel can be observed for the color spectrum set ranges from 0.1 to 0.9 m/s.
within the range of 20...110°C. The maximum temperature of
the PV panel is at 112.3°C. The average temperature is 88.4°C.
The minimum temperature is 64.5°C.
a.
b.
Figure 10. a. Front; b.Bottom View of PV in the Numerical
Simulation Without Cooling
84
JOURNAL OF APPLIED ENGINEERING SCIENCES VOL. 14(27), ISSUE 1/2024
ISSN: 2247-3769 / e-ISSN: 2284-7197 ART.NO. 398 pp. 81-88
85
JOURNAL OF APPLIED ENGINEERING SCIENCES VOL. 14(27), ISSUE 1/2024
ISSN: 2247-3769 / e-ISSN: 2284-7197 ART.NO. 398 pp. 81-88
In the domain of PV operation with serpentine cooling, as Regarding the PV cooling system with distributor/collector, Fig.
shown in Fig. 14, the graph of temperatures obtained from the 16 displays the graph of temperatures obtained from the
numerical simulation is presented. Similarly, I consider the numerical simulation. The maximum temperatures obtained
obtained maximum temperatures as preliminary data to ensure serve as preliminary data to ensure the logical operation of PV
the logical operation of PV under the most adverse exterior under the most extreme external climate conditions. At the
climatic conditions. At the beginning of each serpentine cooling beginning of each cooling cycle with the distributor/collector,
cycle, for each of the five allowed circulation velocities for each of the five allowable circulation velocitys considered, a
considered, a maximum temperature is observed. In the velocity maximum temperature is observed. In the range of velocitys
domain from 0.1 m/s to 0.5 m/s, the maximum temperature from 0.1 m/s to 0.3 m/s, the maximum temperature obtained
obtained does not match the STC efficiency calculation. In this does not match the STC efficiency calculation, indicating that
case, the resulting temperatures are too high. In the rest of the the resulting temperatures are too high. However, in the rest of
approached velocity domain, given the forced circulation of the the velocity range considered, due to the forced circulation of
cooling agent and better distribution on the PV surface, the the cooling agent and a uniform and constant distribution over
operating temperatures are lower and even optimal in terms of the entire PV surface, the operating temperatures are lower.
efficiency.
Velocity of
cooling agent 0,1 0,7 0,9
0,3 0,5
(m/s)
TminPV(°C) 55,9 37,5 32,7 30,5 29,1
STC(%) 86,095 94,375 96,535 97,525 98,155
TmedPV(°C) 78,8 44,9 37,1 33,5 31,5
STC(%) 75,79 45,79 94,555 96,175 97,075
TmaxPV(°C) 105,3 105,3 42,7 37,5 34,6 Figure 16. Distribution of temperatures for the PV system with
STC(%) 63,865 63,865 92,035 94,375 95,68 distributor/collector cooling
In Table 7, according to the STC test, we extracted the
Table 6. Calculation of PV efficiency with Serpentine Cooling following calculation for the efficiency of the PV system with
according to STC distributor/collector cooling.
TmedPV (°C)
TmaxPV (°C)
TminPV (°C)
Velocity of
STC2(%)
STC3(%)
(m/s)
85,6 73,36
45,8
64,9
84,2
0,1
45,8
57,0
0,3
39,9
49,2
over time.
0,5
86
JOURNAL OF APPLIED ENGINEERING SCIENCES VOL. 14(27), ISSUE 1/2024
ISSN: 2247-3769 / e-ISSN: 2284-7197 ART.NO. 398 pp. 81-88
97,84 97,175
95,725 94,78
92,125 90,91
31,3
36,6
45,2
0,7
system manages to exchange heat more efficiently, resulting in
lower temperatures and higher efficiency for the entire PV
29,8
34,5
42,5
system.
0,9
87
JOURNAL OF APPLIED ENGINEERING SCIENCES VOL. 14(27), ISSUE 1/2024
ISSN: 2247-3769 / e-ISSN: 2284-7197 ART.NO. 398 pp. 81-88
As already evident, in this simulation, we numerically identified Rauf A, Liu X, Amin W, Ozturk I, Rehman O.U, Sarwar S,
this thermal potential. The benefits of improved heat exchange Energy and ecological sustainability: Challenges and
are reflected in the efficient operation, durability, and reliability panoramas in belt and road initiative countries, 2018,
of the photovoltaic system, as well as the available thermal Sustainability, 10, 2743. [CrossRef].
potential. Thermal energy storage is possible, and there are
numerous approaches in this regard. The advantage lies in Robles Algari, C., Rodríguez Álvarez, O, An overview of the
reducing the environmental impact in the field of electricity and renewable energy in the World, Latin America and Colombia,
thermal energy production through this hybrid electro-thermal 2018, Espacios, 39, 1–16.
system, optimizing electricity and thermal consumption,
reducing fossil fuel consumption, minimizing the environmental Tahar, M., Rosli, N., Simulation Study of Cooling System for
impact, reducing assembly and maintenance periods for Photovoltaic Panel Using Ansys M.S.M, 2022, Journal of
photovoltaic-thermal (hybrid) systems, and ultimately reducing Modern Manufacturing Systems and Technology (JMMST) vol.
costs. 6, issue 2, 26 – 33.
Yau YH, Lim KS. Energy analysis of green office buildings in the
tropics-Photovoltaic system. Energy Build 2016; 126:177–93