Photovoltaic-thermal system combined with wavy tubes, twisted tape inserts and a novel coolant fluid: energy and exergy analysis

To create a highly efficient photovoltaic-thermal (PV-T) system and maximise the energy and exergy efficiency, this study aims to propose an innovative configuration of a PV-T system comprising wavy tubes with twisted-tape inserts. Following the validation of a numerical model, a parametric study has been conducted to assess the geometrical effects of twisted tape and wavy tubes, as well as the coolant fluid type and velocity, on the overall performance of a PV-T system, located in Shiraz, Iran. It is found that employing twisted tape improves the energy and exergy efficiency by approx. 6.3%. The best configuration yields 12.4% and 16.8% increase in energy and exergy efficiency compared to conventional PV systems. This is achieved at 15% volumetric concentration of microencapsulated phase change material slurry. The monthly variation of global horizontal irradiance in Shiraz highly affected the energy and exergy efficiency of PV-T, with July and October exhibiting the most efficient months, corresponding to 90% and 11.3%, respectively. Highlights Energy and exergy of photovoltaic-thermal system with twisted tape and wavy tubes. Reduction in the temperature of the PV panel by 8.3°C compared to the base case. 12.4% and 16.8% increase in energy and exergy efficiency in the optimum configuration. Higher energy and exergy efficiency from microencapsulated phase change material.


Introduction
A key challenge associated with using photovoltaic (PV) modules to produce electricity is an increased panel temperature, leading to reduced electrical efficiency (Essam et al., 2022).To overcome this challenge, heat transfer fluid (HTF) is used underneath the PV panel to cool it down.The heated fluid can then be utilised for thermal applications (Huang et al., 2020).These photovoltaicthermal (PV-T) systems are cutting-edge renewablepowered technology that combines the power of PV cells and thermal collectors in one integrated system (Li et al., 2022).PV-T systems have gained increasing attention in recent years as an efficient and cost-effective way to meet energy demands while reducing greenhouse gas emissions.They can be utilised in a wide range of applications such as space heating, water heating, building, and industrial process heat.
The geometrical configuration of the PV-T system is crucial for achieving higher efficiency (Kim & Tae, 2022).Researchers have evaluated various configurations to improve the performance of PV-T systems.Shahsavar et al. (2020) experimentally investigated a PV-T system equipped with cooling tubes with different shapes.They found that electrical efficiency was increased by 12% by adding a coolant tube.Moreover, the use of nanoparticles enhanced the exergy and energy efficiencies by 7% and 6.6%, respectively.Tomar et al. (2019) employed glass-to-glass PV modules in the PV-T system and found a maximum hourly-averaged electrical efficiency of ∼ 13%.In another study, Hossain et al. (2019) focused on the use of serpentine pipe PV-T collector design.The results indicated a 5.34% increase in electrical efficiency with 3.92°C reduction in the PV temperature.Employing polycrystalline and monocrystalline panels in PV-T systems, Özakin and Kaya (2019) achieved 55% and 70% increase in energy efficiency and a 70% and 30% increase in exergy efficiency compared to conventional PV-T systems, respectively.
Wavy tubes have shown great potential to improve the thermal performance compared to straight tubes.This is due to an increased heat exchange area and enhanced mixing, resulting in a more uniform temperature distribution and a greater net efficiency in the PV panel.Eisapour et al. (2020) numerically evaluated the use of wavy tubes in PV-T modules.In their best design, the energy and exergy efficiencies were increased by 6.1% and 4.3%, respectively.In another study by the same research group (Eisapour et al., 2021), the authors examined an innovative layout of wavy tubes with ascending and descending amplitudes.The authors enhanced electrical and thermal efficiencies from 10.94% and 61.04% in the case of a straight tube to 11.32% and 65.21%, respectively.The authors indicated that using SiC and microencapsulated phase change material (MPCM) (MPCM-28) as a working fluid instead of water improved the electrical efficiency by 0.4%.
The use of vortex generators to create swirling patterns and disrupt the formation of boundary layers has gained significant attention in recent years (Darbari et al., 2020;and Tian et al., 2020).Aridi et al. (2022) focused on the hydrothermal performance of a trapezoidal vortex generator in a concentric tube heat exchanger through numerical simulation.They improved the heat transfer ratios and thermal factor by 97% and 210%.Using a longitudinal vortex generator, Carpio and Valencia (2021) increased the thermal performance of the system by 52%.Arasteh et al. (2021) numerically analysed the hydrothermal performance of twisted rotational tape pitch distance inside a heat-exchanging tube.The results showed a higher performance evaluation criterion (PEC) at lower Reynolds in the case of rotating twisted tape.However, the maximum PEC was found at a Reynolds number of 1000 and a pitch distance of L/6 for a stationary twisted tape.Mashayekhi et al. (2021) corroborated the same finding of Arasteh et al. (2021) in terms of better enhancement in the thermal characteristics of rotating twisted tape at lower Reynolds numbers under the laminar regime, 32.8-39.6%increase in the Nusselt number at Reynolds number of 250.
The cooling fluid also plays a crucial role in the overall electrical performance of the PV-T system (Akbar et al., 2021;Moore & Wei, 2021).Menon et al. (2022) evaluatean unglazed PV-T system's electrical and thermal performance using water and copper oxide-based nanofluid.They reduced the panel temperature by 15°C and 23.7°C and achieved an electrical efficiency of 14.6% and 17.6% when the system was cooled with water and nanofluid, respectively.Some researchers investigated the impact of using two nanoparticles (Ghalambaz et al., 2020a;Mehryan et al., 2020), composite nanoparticles (Ghalambaz et al., 2019;Ho et al., 2020), or non-Newtonian fluids (A.Saleem et al., 2020;S. Saleem et al., 2021) on the heat transfer enhancement of PV-T systems.In this regard, Abdallah et al. (2019) reduced the temperature of the PV panel by 12°C -resulting in a maximum system efficiency of 61% -using multi-wall carbon nanotube nanofluid.Fu et al. (2021) introduced a PV-T collector that utilises two types of coolant fluids: water and MPCM slurry.The use of MPCM slurry as a coolant fluid results in a more considerable decrease in the temperature of PV-T cells, leading to 13.5% and 0.8% increase in the maximum thermal efficiency and average electrical efficiency.
Reviewing the literature highlights that many efforts have been made to develop effective thermal management systems for high-performance PV-T systems.However, despite recent advances, further research and development are required to gain sufficient insight into the combined heat transfer enhancement techniques to boost the performance of PV-T systems.While the use of twisted tapes in other heat transfer applications are emerging in recent years (e.g.Ghalambaz et al., 2020b;Mashayekhi et al., 2020), the increase of convective heat transfer coefficients of PV-T systems still remains not fully understood.Therefore, to address this research gap, the present study aims to investigate the use of twisted tape inserts in wavy tubes to evaluate the performance of PV-T systems from an energy and exergy standpoint.The present work examines the effect of twisted tape/wavy tube geometrical parameters and different types of coolant fluid on the system performance.Unlike previous studies that only considered solar radiation intensity as a constant variable, this research employs an isotropic sky module to calculate location-specific solar radiation intensity for a tilted PV-T system, allowing for a more accurate monthly performance assessment of the proposed configuration.

Model description
This study investigates a PV-T system incorporating wavy cooling tubes equipped with twisted tape inserts.Figure 1 shows the proposed PV-T system together with the PV panel, the absorber plate, cooling tubes, and the boundary conditions.In the PV-T collector, the PV module is located on the upper surface, the absorber plate is in the middle, and the cooling tubes are placed at the bottom.As illustrated in Figure 1, the boundary condition of the upper surface is considered to be the incoming solar irradiance and heat losses through convection and radiation mechanisms.
The schematic of the present module is shown in Figure 2. As seen in Figure 2(a), the PV system benefits from five wavy tubes under the panel for cooling purposes.All surfaces are considered insulated except the upper layer, which absorbs solar energy from the Sun.Due to the symmetrical condition, only 1/5 of the whole system is considered for the computational domain (Ekramian et al., 2014).The length of the absorber plate is 2000 mm, from which 1640 mm of the absorber plate is enclosed by the PV module (Eisapour et al., 2020).The absorber plate absorbs the solar thermal energy depending on the absorber's solar radiation intensity and absorptivity.
Figure 3 shows the characteristics of the proposed twisted tape inserts.In this figure, L t is the pitch of the twisted tape, W t is the width of the twisted tape, and H t is the height.Table 1 presents the characteristics of various proposed cases considered in this research.Case 1 is the system with a straight tube, which is assumed to be the base case in this study.Cases 2 and 3 are the system with wavy tubes without twisted tape inserts.Adding twisted tapes to cases 1, 2, and 3 forms cases 4, 5, and 6, respectively.

Assumptions
In order to make the PV-T numerical modelling more straightforward, the following hypotheses are applied: • The fluids are incompressible.
• The problem is assumed to be steady-state.
• The flow regimes in cases 1-3 and cases 4-6 are laminar and turbulent, respectively.The latter is due to the presence of twisted tube inserts.• The external surface of the cooling tube and the absorber plate is well insulated.• Nanofluid and MPCM slurries are in thermal equilibrium and single-phase in the calculations (Demir et al., 2011).• Convection and radiation heat losses are considered.
• No thermal resistance between the PV panel and absorber plate, leading to identical temperatures at the interface (Fontenault & Gutierrez-Miravete, 2012).

Laminar model
The continuity, momentum, and energy equations for a laminar flow are written as follows (Liaw et al., 2021):  The energy equation for the solid part is given as follows:

Turbulence model
Using twisted tape inserts inside the tube generates vortices, resulting in the transition from laminar to turbulent flow (Liaw et al., 2021).To accommodate this, the Realisable k-enhanced wall treatment model is used as below (Liaw et al., 2021): where k is the turbulent kinetic energy, is the turbulence dissipation rate, G k is the kinetic enhancement due to buoyancy, Y M is the fluctuation dilation in the compressible turbulence regime, and S k is a source term specified by the user.The values of C ε1 , C 3 , C ε2 are constants in these equations (Liaw et al., 2021).Further information on performance of different two-equation turbulence models can be found in several publications including those by the present authors (Keshmiri et al., 2008(Keshmiri et al., , 2015(Keshmiri et al., , 2016) )

Boundary conditions
The upper surface of the system absorbs solar energy from the Sun.A fraction of this energy is absorbed, transferred to the HTF, and a fraction is dissipated due to convection and radiation losses.The net heat absorbed where the plate is covered by the PV cell is calculated as (Du et al., 2016): where σ sb = 5.76 × 10 −8 (W/m 2 K) and h c = 3u a + 2.8 and u a = 0.5 m/s (Kumar & Mullick, 2010).The net heat absorbed by the sections of the absorber plate not being covered by PV cells is obtained by: The ambient conditions are considered for Shiraz, a southern city in Iran, located at 29.59 N and 52.58 E. Table 2 illustrates Shiraz's irradiance intensity and atmospheric temperature in 2020 according to the weather data file.Further information is provided in Appendix A.
The surfaces in contact are entirely adiabatic.The shared surfaces between the absorber plate and the PV module are symmetrical.Velocity inlet and pressure outlet boundary conditions are assumed at the inlet and Table 3. Absorber plate and PV panel properties (Bhattarai et al., 2012).
PV panel Absorber plate outlet of the fluid.The inlet temperature of the HTF is considered ambient for all cases (298.15K).On the tube walls, a no-slip boundary condition is employed.Table 3 displays the properties of the PV panel and absorber plate.The electrical reference efficiency of the PV cell, η ref , is evaluated when the PV temperature approaches the reference temperature.

Properties of different working fluids
The working fluids in this study include (i) pure water, (ii) Ag-water nanofluid, (iii) MPCM slurry, and (iv) MPCM nano-slurry.For Ag-water nanofluid and MPCM nanoslurry, a user-defined function (UDF) is developed to calculate the thermal conductivity of the nanofluid and the specific heat transfer coefficient of MPCM slurry as a function of temperature.

Properties Ag-water nanofluid
The properties of the nanofluids are calculated as follows (Kakaç & Pramuanjaroenkij, 2009): For φ ≤ 5%, the viscosity of nanofluid is calculated as (Mishra et al., 2014): The effective thermal conductivity can calculate as (Xuan et al., 2003): The thermophysical properties of Ag nanoparticles are listed in Table 4 (Öğüt, 2009).The properties of water are considered at 298.15 K.

MPCM slurry
The PCM material in this study is hydrocarbon noctadecane.Table 5 presents the properties of the MPCM slurry, the PCM and water composite, and the base fluid.According to the theoretical homogenous models and experimental studies, MPCM slurry properties are given as follows (Alquaity et al., 2012): The MPCM slurry viscosity is obtained from (Vand, 1945): The effective thermal conductivity is calculated from the Maxwell model (Alquaity et al., 2012): The bulk-specific heat is calculated as (Alquaity et al., 2012): It should be mentioned that melting occurs over the temperature range of 297-302 K.

MPCM nano-slurry
In this case, the Ag-water nanofluid is considered a carrier fluid combined with MPCM -containing 5% Agwater nanofluid and 15% MPCM.Eqs (12-15) and Eqs.(16-19) estimates the average properties of nanofluid and MPCM nano-slurry, respectively.The property of fluids is listed in Table 6.It is to be noted that the properties such as thermal conductivity for nanofluid and MPCM nanofluid and specific heat for MPCM slurry/nanoslurry are not known at the beginning since these parameters are the function of fluid temperature.To address this, a UDF is employed to calculate these properties.

Energy analysis
The electrical and thermal efficiencies of the PV panel are defined as (Yu et al., 2019): where p is the packing factor which is considered equal to 0.82.η pv and η I are the panel efficiency and the first law efficiency calculated as (Ji et al., 2007):

Exergy efficiency
The second thermodynamic law is adopted for the exergy analysis.The energy analysis shows the quantity of energy, while the exergy analysis shows the quality of energy.The exergy balance for the PV-T system can be calculated as (Demirel, 2013): Eq. ( 24) can be written as: The term in the bracket is useful exergy.The exergy efficiency is then given as follows (Khanjari et al.): The electrical, thermal, and solar exergy are calculated as follows (Wu et al., 2015): where The solar exergy absorbed by the PV panel is calculated as (Khanjari et al., 2016): where T sun is set to 5774 K.

Numerical method and validation
After creating the geometry using Ansys Design Modeller, Ansys Meshing software is used to generate the mesh.The SIMPLE algorithm is used to solve the pressure field, and the QUICK scheme is used to discretize the momentum and energy equations.The residual values for the convergence of continuity and momentum are 10 −6 , while for the energy equation, it is 10 −8 .The geometry is divided into blocks to achieve the structured mesh, and then the geometry is meshed using the swipe element method.Figure 4 shows the generated mesh for case 3 as an example.This figure shows that mesh size is refined in the vicinity of solid-fluid interfaces due to larger velocity/temperature gradients.The results of grid independency for case 3 are shown in Figure 5.The total number of mesh elements is varied between 830,000 and 4,560,000.The outlet and PV temperatures are chosen for the purpose of mesh independency checks.As seen in this figure, there is a marginal difference in the temperature by further refining the mesh beyond 3,418,000 grids.As such, this number of grids is chosen for all subsequent simulations.
For model validation, the results are compared with the data published in the literature for both laminar and turbulent flow.To validate the laminar flow model (cases 1-3), the numerical results of Yu et al. (2019) are used, who used straight tubes under the PV module to cool down the panel and increase the electrical/thermal efficiency of the PV cell.For this purpose, they used a microencapsulated slurry with water-based fluid.The length and width of the module considered in their research were 2000 mm and 1000 mm, of which 1640 mm was covered by the panel.The length of the cooling tube is 2000 mm, while its diameter is 8 mm.The temperature of the cooling fluid at the inlet is 298.15°C, and the solar radiation intensity is 1000 W/m 2 .Figure 6 shows the variation in PV temperature relative to the water inlet velocity for the present study and Yu et al. (2019) model.As shown, the results follow a similar trend with acceptable accuracy.
The validation of the turbulent flow (cases 4-6) is performed with the numerical results of Liaw et al. (Liaw et al., 2021).In their research, twisted tape inserts in straight and helical tubes were studied using the kmodel to solve the flow equations.The tube length was 1.258 m for the straight tube, and in the case of the helical tube, the coil pitch was 0.02 m.In addition, the width and height of the twisted tape were considered 0.4 and 10 mm, respectively.The operational conditions, including inlet temperature, wall temperature, and outlet pressure, were 298, 353 K, and zero relative pressure, respectively.Figure 7 illustrates the variation of the heat  transfer coefficient with the Reynolds number for the straight pipe, where good agreement is obtained between the results.
The accuracy of the model in terms of nanofluid flow is assessed with the numerical study of Khanjari et al. (2016).The proposed PV-T model consisted of a glass cover, a set of PV cells, a thermal sheet, and a tube collector in a header and riser configuration.To reduce the computational time, the model is simplified to include one riser tube just.In addition, the boundary condition and geometry configuration are adjusted according to Yu The MPCM slurry fluid is validated with the experimental data of Chen et al. (2008).The authors employed the MPCM slurry passing through a circular stainlesssteel tube for cooling the PV module with constant heat flux.The same physical domain is used for the model  validation.The MPCM concentration and particle diameters are 15.8% and 8.2 μm (Chen et al., 2008).Figure 9 shows the Nusselt number as a function of cylinder dimensionless length for the present work and those obtained by Chen et al. (2008) under the constant heat flux of 3 × 10 4 (W/m 2 ).The comparison of the results shows acceptable accuracy of the present model.

Results and discussions
The panel temperature and heat absorbed by the HTF for case 4 are compared against case 1 to assess the advantages of using twisted tape to enhance the cooling performance of a PV panel in a straight tube.This is followed by a detailed sensitivity analysis of geometrical parameters of twisted tape/wavy tube, types of coolant fluid, and monthly solar irradiation intensity on the overall exergy and exergy efficiencies of the PV-T systems.

Effect of twisted tape inserts in a straight tube
Figure 10 shows that the panel temperature is higher when there is no twisted tape inside the tube, reaching the maximum temperature of 326 K compared to 320 K with twisted tape.This is attributed to the twisted tape creating a rotational motion in the flow near the wall, which generates secondary flow, enhances flow mixing, and disturbs the boundary layer (Mashayekhi et al., 2022).This, in turn, leads to a reduction in the PV panel temperature and an increase in the electrical efficiency of the PV panel.
Figure 11 shows the velocity contours at different cross-sections of a straight tube with/without twisted tape at 0.1 m/s inlet velocity (cases 4 and 1).The minimum and maximum velocities are obtained near the tube wall and centre of the tube in the case without twisted tape.However, when the twisted tape is inserted inside the tube, the boundary layer is disturbed, and the thickness of the boundary layer is reduced due to enhanced mixing effects (Ghalambaz et al., 2020b).
Table 7 lists the values of the average PV temperature as well as the total heat absorption by the cooling fluid with and without twisted tape in a straight tube at an intake velocity of 0.1 m/s (cases 4 and 1).The solar radiation intensity is assumed to be 778 W/m 2 from the weather data file for Shiraz in October.The panel temperature decreases by 4.1 K, while the net heat absorbed by HTF increases by 36.3W for a case with twisted tape compared to the base case.
Figure 12(a) shows the energetic and exergetic efficiencies of the PV module for cases 4 and 1.According to Eq. ( 10), reducing the PV panel temperature is accompanied by a reduction in the temperature difference between the panel and the environment and, as a result, lower heat loss due to convection and radiation.The energy and exergy efficiencies increase when the twisted tape is inserted inside the straight tube.The thermal, electrical, and overall energy efficiencies increase by 6.8%, 4.1%, and 6.2%, while the thermal, electrical, and overall exergy efficiencies increase by 28%, 2%, and 6.3% in the case of twisted tape as compared to the base case, respectively.
The results of this study suggest that the use of a twisted tape inside a straight tube under a PV panel is an effective method for enhancing the cooling performance of the panel and increasing the overall energy and exergy efficiencies of the system.

Effect of fluid velocity
The effect of inlet velocity ranging from 0.05 m/s to 0.25 m/s is investigated in this section.Figure 13 shows the fluid temperature distribution for cases 4 and 1 under various values of inlet velocity.Inserting twisted tapes reduces the cross-sectional area of the cooling channel, increasing the flow velocity and centrifugal force -leading to a better redirection of the core flow toward the heated wall and a reduction in the fluid flow temperature.At a lower inlet velocity of 0.05 m/s, the generated centrifugal force is not strong enough to lower the PV panel temperature.However, moving to higher Reynolds numbers increases the flow velocity and momentum of the fluid -assisting the fluid in energy absorption from the heated wall and reduction in the PV panel temperature.
Table 8 lists the average PV temperature as well as the total heat absorption by the cooling fluid for different fluid velocities considering the mean solar radiation in October.It can be noticed that at higher inlet velocities, the panel temperature decreases, and the heat absorbed by the cooling fluid increases.Due to the secondary flow formation and the vortex generation, heat transfer in the twisted tape case is more significant.As a case in point, the change in the inlet velocity from 0.05 to 0.1 m/s is associated with the increase in the heat absorbed rate by 21.71% in case 4 and by 8.21% in case 1.
Figure 14(a) shows the PV-T system thermal and electrical energy efficiencies for different inlet velocities of the fluid.Increasing the inlet velocity increases both the thermal and electrical efficiencies.At 0.1 m/s, the thermal and electrical efficiencies of a straight tube with a twisted tape are 6.9% and 2% higher than that for the base case because of the higher reduction in the PV module temperature as a result of the enhancement in the heat transfer coefficient.
Figure 14(b) shows the PV-T system thermal and electrical exergy efficiencies for different velocities.Increasing the inlet velocity decreases the thermal exergy efficiency and increases the electrical exergy efficiency.The electrical exergy efficiency of a PV module with twisted tape is higher than that of a PV module without twisted tape.As velocity increases from 0.05 m/s to 25 m/s, the thermal exergy decreases as a result of higher heat transfer to the HTF.According to Eq. ( 30), the electrical exergy of the system increases due to a decrease in the PV panel temperature.Over this range of velocity, the thermal and electrical efficiencies in the case with a twisted tape drop by 44.36% and increase by 6.2%, respectively.For the same case without a twisted tube, the thermal and electrical efficiencies decrease by 52.4% and increase by 4.2%, respectively.
Figure 14(c) shows the overall energy and exergy efficiencies of the PV-T system.Based on Eqs. ( 23) and ( 27), higher energy and exergy efficiencies are obtained for the PV system equipped with a twisted tape than for the base case.The increase in the inlet velocity from 0.05 to 0.25 m/s is associated with an increase in the energy  efficiency by 20.3% and a decrease in the exergy efficiency by 5.7% for the system with a twisted tape.In contrast, for the straight tube case, the energetic and exergetic efficiencies increased by 14.1% and reduced by 8.6%, respectively.

Different types of configurations
Figure 15 shows the temperature distribution of the PV panel for different geometries at an inlet velocity of 0.1 m/s for the mean solar radiation in October.As seen in this figure, the straight tube suffers from effective cooling in the corners of the PV panel, which reduces the system performance.This can be resolved by employing wavy tubes with twisted tape.It is evident that using the wavy tube -with 200 mm wavelength and 50 mm amplitude -and twisted tape (case 6) results in a better temperature distribution than other cases due to higher heat transfer characteristics and lower stagnation regions in the fluid flow.
Figure 16(a and b) show the overall energy efficiency for different geometries.As seen in this figure, the energy efficiency for cases 2-6 increases by 5.1%, 7.8%, 6.2%, 11.5%, and 12.4% compared to the base case (case 1), respectively.The exergy efficiency in wavy-shaped tubes is higher than that of the straight tube because of the more significant heat transfer area and, consequently, more reduction in the temperature of the PV panel.Using a twisted tape in the tube also increases the exergy efficiency compared to the base case.The exergy efficiencies of cases 2-6 increased by 6.9%, 7.4%, 6.3%, 15.8%, and 16.8% compared to the base case (case 1).This is attributed to the amounts of average PV temperature and total heat absorption by the HTF for these geometries, listed in Table 9 for different fluid velocities.As seen in this table, the panel temperature decreases as the inlet velocity increases, which is more noticeable at lower velocities.Moreover, the panel temperature in cases 2-6 is 3.3, 5.2, 4.1, 7.8, and 8.3 K, cooler than the base case at the inlet velocity of 0.1 m/s.
Figure 17(a) shows the thermal and electrical energy efficiencies at different velocities.As the velocity increases, the thermal and electrical efficiencies increase.By changing the inlet velocity from 0.05 m/s to 0.25 m/s for cases 1-6, the thermal efficiency increases by 20%, 18.3%, 15.5%, 22.5%, 11.5%, and 10%, and the electrical efficiency increases by 4.1%, 5%.4.5%, 6.3%,3.4% and 3.3%, respectively.Figure 17(b) illustrates the thermal and electrical exergy efficiencies for different geometries at different velocities.As the efficiency increases, the thermal exergy efficiency decreases, and the electrical exergy increases.In cases 2, 3, 4, 5, and 6 at the inlet velocity of 0.05 m/s, the thermal exergy efficiency increases by 5.1%, 25.3%, 6.3%, 26.2%, and 30.5% compared to the base case.The electrical efficiency increases by 8.8%, 2.3%, 9.2%, 3.9%, and 4.4% for cases 2-5 compared with  the base case.In case 6, with increasing the inlet velocity, the thermal exergy efficiency decreases from 3.3% to 1.4%, and the electrical exergy increases from 9% to 9.3%.Figure 17(c) shows the overall energetic and exergetic efficiencies for different geometries at different velocities -higher velocity results in higher energy and lower exergy efficiency.Case 6 has the best performance from energy and exergy, in which a twisted tape is inserted inside a wavy tube with a wavelength of 200 mm and an amplitude of 50 mm.The overall energy efficiency for cases 2, 3, 4, 5, and 6 increases by 2.6%, 7.1%, 2. 9%, 12.8%, and 14.8% compared to the base case, while these numbers are 1.8%, 7.4%, 2%, 9% and 10.3% for exergy efficiency at the inlet velocity of 0.5 m/s.Increasing the inlet velocity from 0.05 m/s to 0.25 m/s increases the energy efficiency from 76.1% to 83.13% and decreases the exergy efficiency from 12.3% to 10.7%. Figure 17(d) shows the amount of thermal and electrical exergy.The results show that the amount of thermal exergy decreases for a higher velocity while the electrical exergy increases.In the best case (case 6), the thermal exergy decreases from 11.2 W to 4.7 W with increasing velocity, and the electrical exergy increases from 30.6 W to 31.6 W.

Different types of coolant fluid
To study the effect of fluid type on the performance of the PV-T system, eight different fluids have been used, including Ag-water nanofluid with 1%, 3%, and 5% volume fractions, n-octedane microcapsules with 5%, 10%, and 15% volume concentrations, and MPCM Nano slurry with 5% nanoparticles and 15% microcapsule.It should be noted that there is a limit for the volumetric concentration of MPCM because of the non-Newtonian behaviour of MPCM at high concentrations of MPCM;  fluid behaviour becomes non-Newtonian -which dictates the upper limit of pumping.As such, the maximum concentration of 15% is considered for the MPCM.
Figure 18(a) presents the energy efficiency of a PV-T system at an inlet velocity of 0.1 m/s for case 5, considering the mean solar radiation in October.The addition of nanoparticles increases the overall heat transfer coefficient and the performance of a PV-T system.Using microencapsulation also increases the fluid heat capacity, and, as a result, the fluid absorbs more heat.The highest thermal and electrical energy efficiencies are found in the case of water with 15% slurry due to more decrease in the

Seasonal effects
Due to the difference in the solar radiation intensity in a typical year, the performance of a PV panel varies with the geographical location and throughout the year.
Figure 19(a) shows the monthly energy efficiency of the PV panel at the inlet fluid velocity of 0.1 m/s.The highest and lowest thermal energy efficiency occurred in July and January.According to Eq. ( 10) and Table 2, absorber heat increases due to solar radiation intensity and ambient temperature enhancement.As a result, the thermal energy efficiency rises in July and decreases in January.
In contrast to the previous case, the highest and lowest electrical energy efficiency occurred in January and July due to the ambient temperature reduction presented in Table 2.The highest and lowest overall energy efficiency is in July and January.Figure 19(b) illustrates the monthly exergy efficiency of the PV panel at the inlet fluid velocity of 0.1 m/s.According to Eq. ( 28), an increase in solar radiation intensity and a reduction in ambient temperature enhance the exergy efficiency.The best thermal exergy is related to October, while the lowest is related to January.The best electrical exergy is in January, and the lowest is in October.The highest and lowest overall exergy efficiency is in October and January.The proposed system has the best and worst performance in October and January, respectively.
Table 11 shows the monthly average panel temperature and heat absorbed by the fluid.The highest panel temperature is in October (i.e.306.53 K), and the lowest is in January (i.e.302.08 K).The highest and lowest heat absorption occurs in October and January -equivalent to 255.26 and 108.25 W, respectively.

Comparison against the existing data
This section compares the results obtained from the present study with those found in the literature.Table 12 shows the energy and exergy efficiency improvement based on the proposed method in the literature and  energy and exergy performance of the system.This study indicated that changing the geometry of the tubes from a straight to a wavy shape with a twisted tape inserted under the PV panel resulted in significant improvements in energy and exergy efficiencies, corresponding to 12.4% and 16.8% increase compared to the base case, respectively.The results showed that increasing the volume concentration of nanoparticles in the coolant fluid led to an increase in energy efficiency and a decrease in exergy efficiency.Similarly, increasing the volume concentration of the MPCM coolant fluid also led to an increase in energy and exergy efficiencies.MPCM slurry with a 15% volume concentration exhibited the best energy and exergy efficiencies among the coolant fluids considered.
Compared to other papers published in the literature, the co-implementation of a wavy tube and twisted tape provides a more uniform temperature distribution along the tube wall due to the secondary flow induced in the core and near the wall regions.These dual mixing effects prohibit the formation of hot zones and boundary layers, leading to better cooling performance.This is the main reason for higher energy and exergy efficiencies obtained from this work relative to other papers published in the literature.

Conclusion
Maintaining the temperature of a photovoltaic (PV) system in the range of 35°C to 45°C requires an effective cooling system in place.The present study aimed to improve the cooling performance of a photovoltaicthermal (PV-T) system by employing twisted tape inserts in wavy tubes.Having the model validated first, a parametric study was subsequently conducted using CFD to obtain the effect of various parameters, including the inlet velocity, geometry design, different types of coolant fluid, and monthly/seasonal effects on the energy/exergy efficiencies.where nd m , is the day of the month, nd k is the number of days of the month and γ is a dimensionless parameter given by: γ = −1.498+ 1.184ζ − 27.182exp(−1.5ζ) ζ is a dimensionless parameter given by: The value of K T,min and K T,max defined as follow:  A1.
The daily diffuse radiation H d is defined by Erb's correlations as followed (Santos et al., 2003): For ω s ≤ 81.04  Finally, the solar intensity for the average day for all of the months for the city of Shiraz is presented in Figure A1.Then the moderate solar radiation intensity from 10 am to 2 pm is calculated and shown in Table 2.

Figure 1 .
Figure 1.Schematic view of the PV-T system, including the boundary conditions.

Figure 2 .
Figure 2. Schematic view of the helical wavy tube PV-T system: (a) top view, (b) front view.

Figure 3 .
Figure 3. Characteristics of twisted tape in tube: (a) isometric view, (b) front view.

Figure 4 .
Figure 4.The generated mesh for case 3 using 3.4 million cells.

Figure 5 .
Figure 5.The grid independency for case 3 for pure water with u in = 0.1 m/s in shiraz city on October 15.

Figure 6 .
Figure 6.The validation study of the PV-T system compared with Yu et al. (2019).

Figure 7 .
Figure 7.The twisted tape and turbulence model validation study compared with Liaw et al. (2021).

Figure 9 .
Figure 9. Validation of MPCM slurry thermal performance in a cylindrical cavity with Chen et al. (2008).

Figure 10 .
Figure 10.Contours of the temperature distribution in the PV-T system with (a) a straight tube and (b) a straight tube with a twisted tape.Table7.Values of the average PV temperature as well as the total heat absorption by the cooling fluid for the PV-T system with straight with and without twisted tape.

Figure 11 .
Figure 11.Cross-sectional velocity contours for straight tube (case 1) and straight tube with twisted tape insert (case 4) at u = 0.1 m/s.

Figure 12 .
Figure 12.The thermal, electrical, and overall efficiencies of the straight tube with and without twisted tape (cases 4 and 1) (a) energy efficiency and (b) exergy efficiency.

Figure 13 .
Figure 13.Temperature distribution of the PV panel at different inlet fluid velocities for the straight tube with and without twisted tape (cases 4 and 1).

Figure 14 .
Figure 14.Effect of inlet velocity on (a) thermal and electrical energy efficiencies, (b) thermal and electrical exergy efficiencies, and (c) overall energetic and exergetic efficiencies of the PV module.

Figure 15 .
Figure 15.Temperature distribution in PV-T system for different types of configurations at the inlet velocity of 0.1 m/s.

Figure 16 .
Figure 16.Effect of PV system configuration on the (a) energy efficiency and (b) exergy efficiency at the inlet velocity of 0.1 m/s.

Figure 17 .
Figure 17.Effect of inlet velocity on (a) thermal and electrical energy efficiencies, (b) thermal and electrical exergy efficiencies, (c) overall energetic and exergetic efficiencies, and (d) thermal and electrical exergy.

Figure 18 .
Figure 18.Effect of coolant fluid on the (a) energy efficiency and (b) exergy efficiency at the inlet velocity of 0.1 m/s.

Figure 19 .
Figure 19.Monthly performance of the PV system in terms of (a) energy efficiency and (b) exergy efficiency for Shiraz city at the inlet velocity of 0.1 m/s.
max = 0.6313 + 0.267 KT − 11.9( KT − 0.75) 8 (45) The solar irradiation and monthly means of the daily clear index for different months in Shiraz city of Iran are presented in Table radiation I d is obtained with the ratio of hourly diffuse to daily diffuse radiation r d as (Duffie et al.I b is calculated by subtracting I d from I.

Figure A1 .
Figure A1.Solar radiation intensity for different months of Shiraz, Iran.

Figure A2 .
Figure A2.Average low and high-temperature distribution for Shiraz city according to AccuWeather report for the year 2020.

Table 1 .
Characteristics of different cases considered in this research.

Table 2 .
The irradiance intensity and atmospheric temperature for Shiraz city in 2020.

Table 6 .
Properties of considered HTFs at 298.15 K.

Table 8 .
The average PV temperature and the total heat absorption by the cooling fluid for the PV-T system with straight with and without twisted tape for different inlet velocities.

Table 9 .
The average panel temperature and heat absorbed by HTF for the PV-T system with different cases for different velocity inlets.

Table 10 .
Average panel temperature, heat absorbed, and pressure drop for different types of coolant fluid at the inlet velocity of 0.1 m/s.

Table 11 .
Monthly average PV temperature and the total heat absorption by HTF for the PV-T system located in Shiraz at the inlet fluid velocity of 0.1 m/s.
present work relative to the base case with no design improvement.As seen in this table, the geometrical configuration of the PV-T module plays a crucial role in the

Table 12 .
comparison between the current study and other methods in terms of energy and exergy efficiency.Improvement = η with method −η without method η