E ﬀ ect of Triangular Ba ﬄ e Arrangement on Heat Transfer Enhancement of Air-Type PVT Collector

: A Photovoltaic Thermal (PVT) Collector is a device that produces electricity and simultaneously uses a heat source transmitted to back side of the Photovoltaic (PV). The PVT collector is categorized into liquid-type and air-type according to the heating medium. As an advantage, air-type PVT system is easy to manage and can be directly used for heating purposes. The performance of air-type PVT collector is determined by various factors, such as the height of air gap and air ﬂow path (by ba ﬄ es) in the collector. Ba ﬄ es are installed in the PVT collector to improve the thermal performance of the collector by generating turbulence. However, the air ﬂow that a ﬀ ects the performance of the PVT collector can vary depending on the number and placement of the ba ﬄ es. Thus, the ﬂow design using ba ﬄ es in the collector is important. In this study, the performance of an air-type PVT collector due to the arrangement of triangular ba ﬄ es and air gap height at the back of the PV module is analyzed through a simulation program. For this purpose, Computational Fluid Dynamics (CFD) analysis was performed with an NX program to compare and analyze the optimum conditions to improve the performance of the collector.


Introduction
Globally, the use of sustainable energy is increasing, and among renewable energy systems, solar energy and wind power systems are widely used. Especially, solar energy system is classified into a PV system, a solar thermal system, and a photovoltaic/thermal (PVT) system, and it operates by converting solar energy into electricity and thermal energy. However, a system that uses solar energy cannot be completely dependent on power generation and heat collection as it is not uniform. Therefore, it is applied as a method of stably supplying electricity and thermal energy to a storage device by charging it to a consumer [1,2], and a method of increasing system efficiency, such as harvesting more energy by hybridizing with a fossil fuel system such as natural gas [3].
PVT collectors are device that use heat generated at the back of a PV at the same time as the electricity produced at the front of the PV. PVT collectors are classified into air and liquid types according to the fluid used as the heat transfer medium. Air-type PVT collectors have the advantage of being easy to manage. Previous studies focused on improving the collector's own thermal and electrical efficiency through design, simulation, modeling and experimentation of air-type PVT collectors.
For instance, one review paper of previous studies considered results of air flow and single/double flow paths of various air-type PVT systems and various absorbers (i.e., fin, V-groove, round tube, etc.) [4,5]. Conclusions showed that the electrical efficiency was 10%-25% and the thermal efficiency was 40%-70%. Exergy efficiency was also in the range of 5%-25%. Riffat et al. [6] studied parameters affecting the electrical and thermal performance of various types of PVT systems. Parameters examined included the optimal flow rate inside the PVT collector, the presence or absence of a thermal absorber plate, and the thickness of the air layer. In particular, the thermal absorber plate was found to have the most influence on the thermal efficiency of the PVT system. In other studies, it was concluded that as the number of baffles increased, the thermal efficiency increased, but this increase in thermal efficiency decreased above a certain number of baffles [7,8]. So, it is necessary to set a proper number of baffles in consideration of the pressure drop.
Based on the CFD program (i.e., ANSYS Fluent), Chaube et al. [9] analyzed shape, viz. rectangular, square, chamfered, triangular and semicircle baffles after simulation in a Reynolds number range of 2900-19,500. Abuska et al. [10] examined the energy, exergy, economics and environmental performance of air-type collectors with V-groove-shaped protrusions. The thermal and exergy efficiencies ranged from 43% to 60% and 6% to 12%, respectively, and the payback period averaged 4.3-4.6 years. Experiments showed that the collector's thermal efficiency was about 6% higher than that of the flat collector. Furthermore, Fudholi et al. [11] studied the exergy and sustainability index of air-type PVT collectors with V-groove shaped protrusions. The exergy efficiency was 13.36% in theory and 12.89% according to experimental results.
In the studies by Promthaisong et al. [12], triangular ribs were applied to the thermal absorber plate and analyzed in the Reynolds number range of 3000-20,000. The influence of blockage ratio and the pitch ratio of the rib were analyzed as variables. Compared to general collectors, the coefficient of friction and Nusselt number increased by 1.01-4.93 times and 1.02-3.86 times, respectively, by promoting heat transfer. Other studies concluded that when there was a gradient baffle, the maximum enhancement in thermal and effective efficiency were 22.4% and 18.1% in the mass flow rate of 0.045 kg/s compared to general Solar Air Heater (SAH) [13]. In the studies by Bhagoria et al. [14], an air heater with a wedge (triangle) type baffle was designed, and the heat transfer and friction coefficients were analyzed through experiments. The parameters of the baffle were set as the baffle height, angle and baffle spacing; Nusselt number increased up to 2.4 times compared to that of the collector without baffles. In addition, the friction coefficient increased as the angle of the baffles increased, and the heat transfer performance was highest when the baffle angle was 10 • . Yadav et al. [15] investigated the heat transfer inside the collector by installing a triangle baffle on the absorber plate using CFD (ANSYS Fluent). Nusselt number increased with increasing Reynolds number, and was 1.4-2.7 times higher than without a baffle collector. To optimize the triangular baffle, thermal performance was compared and analyzed through CFD simulation using the baffle height, pitch, height ratio and distance ratio between the baffles as variables. The Nusselt number decreased as the distance ratio increased, and the friction coefficient increased by 3.356 times compared to the reference. For an equilateral triangle baffle with Reynolds number of 15,000, distance ratio between baffles of 7.14 and height ratio of 0.042, the Nusselt number was 3.073 times that of reference [16]. Bensaci et al. [17] performed numerical/experimental studies according to the baffle arrangement in the collector. The thermal efficiency improved with increasing number of baffles, and the coefficient of friction decreased as the Reynolds number increased. Choi et al. [18] performed CFD analysis under the same conditions by installing several resistors inside an air-type PVT collector. In terms of heat transfer performance, the intersection of triangular baffles led to an improvement of up to 1.86 times.
As described above, the application of a triangular baffle in an air-type PVT collector improves the heat transfer performance. This effect is mainly caused by the generation of a vortex in the air layer by increase of the Nusselt and Reynolds numbers. However, in order to prevent dead space and pressure drop, it is necessary to design the baffle so that the flow inside the air-type PVT collector is smooth. Consequently, baffle must be designed in consideration of the placement conditions of the baffle, so that the flow inside the collector can be smooth and the heat transfer performance can be improved.
In this study, triangular baffles were installed on a thermal absorber plate in the air gap inside an air-type PVT collector; the thermal characteristics of the baffle arrangement condition were analyzed. For this, the heat transfer, pressure drop, and thermal efficiency of the collector were scrutinized through a simulation program (i.e., NX CFD). The primary purpose was to investigate the effects according to the triangular baffle arrangement applied in the air-type PVT collector.

Model Design
The air-type PVT collector designed for this study is shown in Figure 1. The front of the collector is covered with a general PV module; the dimensions are 1011 mm × 1520 mm. The PV module has a cell covering the front surface, as in a conventional module (electrical efficiency is 17%), and consists of about 54 mono crystal cells between two layers of glass (G/G module). The air layer (height: 34 mm) in the collector has a triangular baffle (width: 82 mm, length: 75 mm, height: 16 mm, angle: 15 • ) that acts as an air flow obstruction, and has a certain spacing and arrangement inside the collector. Thus, the obstacle baffle generates turbulence and is used as an element to increase heat transfer performance.
Sustainability 2020, 12, x FOR PEER REVIEW 3 of 13 scrutinized through a simulation program (i.e., NX CFD). The primary purpose was to investigate the effects according to the triangular baffle arrangement applied in the air-type PVT collector.

Model Design
The air-type PVT collector designed for this study is shown in Figure 1. The front of the collector is covered with a general PV module; the dimensions are 1011 mm × 1520 mm. The PV module has a cell covering the front surface, as in a conventional module (electrical efficiency is 17%), and consists of about 54 mono crystal cells between two layers of glass (G/G module). The air layer (height: 34 mm) in the collector has a triangular baffle (width: 82 mm, length: 75 mm, height: 16 mm, angle: 15°) that acts as an air flow obstruction, and has a certain spacing and arrangement inside the collector. Thus, the obstacle baffle generates turbulence and is used as an element to increase heat transfer performance. The input values of lateral spacing (W) and longitudinal spacing (H1, H2) of the baffle, which affect the baffle thermal performance, are presented in Table 1. To adjust the placement of baffles in the collector, three parameters were set: W (62, 82.5, 144 mm), H1 (0, 47 mm) and H2 (83, 130, 176.3 mm). In addition, a total of 10 cases were simulated and compared, including a reference to a PVT collector with a baffle-free air layer. The input values of lateral spacing (W) and longitudinal spacing (H1, H2) of the baffle, which affect the baffle thermal performance, are presented in Table 1. To adjust the placement of baffles in the collector, three parameters were set: W (62, 82.5, 144 mm), H1 (0, 47 mm) and H2 (83, 130, 176.3 mm). In addition, a total of 10 cases were simulated and compared, including a reference to a PVT collector with a baffle-free air layer.

Energy Balance Equations
The heat transfer process along the cross section of the air-type PVT collector is depicted in Figure 2, in which it can be seen that heat is transferred in the order of the top glass, EVA, PV cell, EVA, bottom glass, air layer, baffle and insulation. However, in this simulation, the heat transfer effect of EVA was insignificant, so only the top and bottom glass layers and PV cells were considered. In Equations (1)-(9), h cv,T , h cv,B , h rd,T and h rd,B are the top and bottom surface convection and radiative heat transfer coefficients of the collector, respectively, T i is the temperature at point i, T sur is the surface temperature around the collector, T a is the ambient temperature, R i-j is the thermal resistance between points i and j, h rd,i-j is the radiative heat transfer coefficient between points i and j, S i is the amount of solar radiation absorbed at point i, P PV is the amount of power produced by the PV cells, h cv,fT and h cv,fB are the top and bottom surface heat transfer coefficients of the air layer, respectively, and q" u,T and q" u,B represent the thermal energy recovered from the top and bottom of the collector air layer. Thermal resistance between ft and fb was not calculated. The flow in the air-type PVT collector enters at the inlet, sweeps the baffles in the cavity, and passes through it. Therefore, the thermal energy recovered (q) is a factor that affects the thermal performance of the collector more than up and down thermal conduction.

Energy Balance Equations
The heat transfer process along the cross section of the air-type PVT collector is depicted in Figure 2, in which it can be seen that heat is transferred in the order of the top glass, EVA, PV cell, EVA, bottom glass, air layer, baffle and insulation. However, in this simulation, the heat transfer effect of EVA was insignificant, so only the top and bottom glass layers and PV cells were considered. In Equations (1)-(9), hcv,T, hcv,B, hrd,T and hrd,B are the top and bottom surface convection and radiative heat transfer coefficients of the collector, respectively, Ti is the temperature at point i, Tsur is the surface temperature around the collector, Ta is the ambient temperature, Ri-j is the thermal resistance between points i and j, hrd,i-j is the radiative heat transfer coefficient between points i and j, Si is the amount of solar radiation absorbed at point i, PPV is the amount of power produced by the PV cells, hcv,fT and hcv,fB are the top and bottom surface heat transfer coefficients of the air layer, respectively, and q″u,T and q″u,B represent the thermal energy recovered from the top and bottom of the collector air layer. Thermal resistance between ft and fb was not calculated. The flow in the air-type PVT collector enters at the inlet, sweeps the baffles in the cavity, and passes through it. Therefore, the thermal energy recovered (q) is a factor that affects the thermal performance of the collector more than up and down thermal conduction.  By performing energy balances in the different nodes of Figure 2, Equations (1)- (12) are obtained [19][20][21][22][23].
The heat transfer equation for the front glass top surface (point 1) is given by Equation (1).
The heat transfer equation for the front glass middle surface (point 2) is given by Equation (2).
The heat transfer equation for the PV layer (point 3) is given by Equation (3).
The heat transfer equation for the back glass middle surface (point 4) is given by Equation (4).
The heat transfer equation for the back glass bottom surface (point 5) is given by Equation (5).
The heat transfer equation for the fluid in the top portion of the air cavity (point fT) is given by Equation (6).
The heat transfer equation for the fluid in the bottom portion of the air cavity (point fB) is given by Equation (7).
The heat transfer equation for the triangle baffle top surface (point 6) is given by Equation (8).
The heat transfer equation for the insulation outer surface (point 7) is given by Equation (9).
The thermal resistances for between the top surface and middle of the front glass layer, R 1-2 , and the middle and bottom surface of the back glass layer, R 4-5 are given by Equation (10). th g , th PV , th tri-baffle and th insulation are the thickness of glazing, PV cell, triangular baffles and insulation layers, respectively. k g , k PV , k tri-baffle and k insulation are the thermal conductivity of glazing, PV cell, triangular baffles and insulation layers, respectively.
The thermal resistances for between the middle of the front glass layer and PV cell layer, R 2-3 , and the PV cell layer and middle of the bottom glass layer, R 3-4 , are given by Equation (11).
The thermal resistances for between the triangular baffles and insulation layer, R 6-7 are given by Equation (12).

Initial Modeling and Validation
The air-type PVT collector was installed on the rooftop of an educational building, located in the Cheonan campus of Kongju National University (36.85 N, 127.15 E). The collector was installed on a 2-axis tracker that can be adjusted horizontally (east-west) and vertically to introduce solar radiation into the normal plane. To evaluate the thermal performance, the collector's inlet/outlet temperature, inlet/outlet flow rate, solar radiation, outside temperature, and wind direction/wind speed were measured. The size (1011 × 1520 × 37 mm) and the triangular baffle arrangement case (W = 104 mm, H1 = 74.8 mm, H2 = 94.8 mm and B.No. = 35) of the collector used in the experiment are shown in Figure 3a,b respectively.

Initial Modeling and Validation
The air-type PVT collector was installed on the rooftop of an educational building, located in the Cheonan campus of Kongju National University (36.85 N, 127.15 E). The collector was installed on a 2-axis tracker that can be adjusted horizontally (east-west) and vertically to introduce solar radiation into the normal plane. To evaluate the thermal performance, the collector's inlet/outlet temperature, inlet/outlet flow rate, solar radiation, outside temperature, and wind direction/wind speed were Outdoor performance test was performed based on ISO 9806 (Solar energy-solar thermal collectors-test methods). Notably, the solar radiation intensity incident on the collector slope was 700 W/m 2 or higher, and the inlet and outlet flow rates were measured at 100 m 3 /h according to the test method. Experimental data were collected and analyzed for 10 min when the inlet/outlet flow rate and temperature, and the outside air temperature were in a steady state. Table 2 illustrates a comparison of experiment and numerical results of outlet temperature of the collector. Experimental and numerical data showed that the outlet temperature was 12.79 °C and 13.04 °C respectively under the same outdoor conditions, which was almost consistent with an error range of about 1.92%.
G is solar radiation (W/m2), Ta is ambient temperature (°C), Tin is inlet air temperature of PVT (°C), Tout is outlet air temperature of PVT (°C), and ΔT is Inlet and outlet temperature difference (°C) Outdoor performance test was performed based on ISO 9806 (Solar energy-solar thermal collectors-test methods). Notably, the solar radiation intensity incident on the collector slope was 700 W/m 2 or higher, and the inlet and outlet flow rates were measured at 100 m 3 /h according to the test method. Experimental data were collected and analyzed for 10 min when the inlet/outlet flow rate and temperature, and the outside air temperature were in a steady state. Table 2 illustrates a comparison of experiment and numerical results of outlet temperature of the collector. Experimental and numerical data showed that the outlet temperature was 12.79 • C and 13.04 • C respectively under the same outdoor conditions, which was almost consistent with an error range of about 1.92%.
G is solar radiation (W/m2), Ta is ambient temperature ( • C), Tin is inlet air temperature of PVT ( • C), Tout is outlet air temperature of PVT ( • C), and ∆T is Inlet and outlet temperature difference ( • C)

Modeling Conditions and Methodology
The NX program used for this study has the capability to simulate fluid flow effects, such as CFD modeling, by quickly creating flow zones for complex geometry and performing computational fluid dynamics. In addition, a combination of thermal analysis functions for conduction, convection, radiation and complex heat transfer is performed [24]. Figure 4 shows the mesh of the computational domain, and the tetrahedral method was selected for the meshing strategy. A tetrahedral mesh is set because it has high solution accuracy and has the advantage of reducing simulation time. In addition, to analyze the airflow characteristics around the triangular baffle, the mesh around the baffles was set densely unlike the mesh of the air layer in the collector.

Modeling Conditions and Methodology
The NX program used for this study has the capability to simulate fluid flow effects, such as CFD modeling, by quickly creating flow zones for complex geometry and performing computational fluid dynamics. In addition, a combination of thermal analysis functions for conduction, convection, radiation and complex heat transfer is performed [24]. Figure 4 shows the mesh of the computational domain, and the tetrahedral method was selected for the meshing strategy. A tetrahedral mesh is set because it has high solution accuracy and has the advantage of reducing simulation time. In addition, to analyze the airflow characteristics around the triangular baffle, the mesh around the baffles was set densely unlike the mesh of the air layer in the collector.  Table 3). The result values between mesh elements have a percentage difference of less than approximately 1%. Considering the calculation time, a grid of 362,677 elements was selected as an appropriate grid. To perform the CFD analysis using the NX program, 10 cases were set with different baffle arrangement as the variables (see Table 1), and the heat transfer, pressure drop, and thermal efficiency of the air-type PVT collector were compared and analyzed. To simulate the collector, the thermal energy of the collector was set to be supplied equally as 700 W/m 2 at 1.54 m 2 . Furthermore, it was input so that the outside air temperature of 7 °C with mass flow rate of 67 m 3 /h was entered into the inlet of collector. The ambient temperature was input to consider the surface heat loss of the collector, and gravity and buoyancy were applied for modeling the actual air flow. The boundary conditions are summarized in Table 4.  Table 3). The result values between mesh elements have a percentage difference of less than approximately 1%. Considering the calculation time, a grid of 362,677 elements was selected as an appropriate grid. To perform the CFD analysis using the NX program, 10 cases were set with different baffle arrangement as the variables (see Table 1), and the heat transfer, pressure drop, and thermal efficiency of the air-type PVT collector were compared and analyzed. To simulate the collector, the thermal energy of the collector was set to be supplied equally as 700 W/m 2 at 1.54 m 2 . Furthermore, it was input so that the outside air temperature of 7 • C with mass flow rate of 67 m 3 /h was entered into the inlet of collector. The ambient temperature was input to consider the surface heat loss of the collector, and gravity and buoyancy were applied for modeling the actual air flow. The boundary conditions are summarized in Table 4. The turbulence model used was the standard K-Epsilon model. It is the most common turbulence model in CFD simulation, and k is the turbulent kinetic energy and ε is the dissipation rate of turbulent energy. This model is based on the time-averaged Navier-Stokes equations, which assume that the time-varying velocities of turbulence can be divided into time-averaged velocities and velocity-dependent velocities [25]. Therefore, in the CFD evaluation for this study, air flow analysis was performed using the standard k-ε turbulence model because the analysis is not complicated and its accuracy is good in terms of computational convergence through air flow simulation data accumulated for a long time. Figure 5 shows the flow velocity distribution for cases with and without baffles inside the air-type PVT collector. Based on the results, the flow velocity of the reference collector passed to the outlet without stagnation. The collector with a baffle (W = 144 mm, H1 = 47 mm, H2 = 83 mm, B.No. = 35) tended to have a weak local flow rate at the back of the baffle, but the flow rate was significantly faster in the spaces between baffles. Therefore, compared with the reference collector, the internal flow rate of the collector with the baffle was faster and the outlet flow rate increased.  The turbulence model used was the standard K-Epsilon model. It is the most common turbulence model in CFD simulation, and k is the turbulent kinetic energy and ε is the dissipation rate of turbulent energy. This model is based on the time-averaged Navier-Stokes equations, which assume that the time-varying velocities of turbulence can be divided into time-averaged velocities and velocity-dependent velocities [25]. Therefore, in the CFD evaluation for this study, air flow analysis was performed using the standard k-ε turbulence model because the analysis is not complicated and its accuracy is good in terms of computational convergence through air flow simulation data accumulated for a long time. Figure 5 shows the flow velocity distribution for cases with and without baffles inside the airtype PVT collector. Based on the results, the flow velocity of the reference collector passed to the outlet without stagnation. The collector with a baffle (W = 144 mm, H1 = 47 mm, H2 = 83 mm, B.No. = 35) tended to have a weak local flow rate at the back of the baffle, but the flow rate was significantly faster in the spaces between baffles. Therefore, compared with the reference collector, the internal flow rate of the collector with the baffle was faster and the outlet flow rate increased.  Figure 6 shows the temperature distribution inside the air-type PVT collector for cases with and without baffles. In the collector with triangular baffles, the air temperature rose due to locally low  Figure 6 shows the temperature distribution inside the air-type PVT collector for cases with and without baffles. In the collector with triangular baffles, the air temperature rose due to locally low flow rate on the back of the baffle, but the space between the top of the baffle and the back of the PV was narrow, causing the air to sweep away quickly. Therefore, it can be deduced that air did not stagnate locally on the back of the baffle and passed quickly to the outlet. As a result, the collector with baffles had a faster outlet flow rate and a higher outlet temperature than the reference collector without baffles, so the former's heat transfer performance was advantageous.

Heat Transfer Performance
Sustainability 2020, 12, x FOR PEER REVIEW 9 of 13 flow rate on the back of the baffle, but the space between the top of the baffle and the back of the PV was narrow, causing the air to sweep away quickly. Therefore, it can be deduced that air did not stagnate locally on the back of the baffle and passed quickly to the outlet. As a result, the collector with baffles had a faster outlet flow rate and a higher outlet temperature than the reference collector without baffles, so the former's heat transfer performance was advantageous.  Figure 7 graphically illustrates the outlet velocity for each baffle arrangement condition of the air-type PVT collector with triangular baffles. The outlet velocity of the collector with the baffle was higher than that of the reference collector (0.564 m/s). It was confirmed that the smaller the left and right spacing variable (W) of the baffle was, the faster the outlet velocity and the higher the related average Reynolds number were. This is because only the arrangement was changed in a state in which the baffle size was fixed; as the length of the lateral direction (W) decreased, the flow path became narrower. However, the longer H2 was, the farther the front and back between the baffles were and the more space was created in the transverse direction. Consequently, the possibility of vortex generation by the baffle was reduced. Consequently, it was seen that the difference in the average Reynolds number range also decreased. Furthermore, the smaller the number of baffles was, the faster was the outlet velocity, but with a correspondingly lower Reynolds number. This is attributed to the flow in the air layer, which was non-uniform; the air flow rapidly passed through the flow paths in the left and right longitudinal directions in the collector with length of W.  Figure 7 graphically illustrates the outlet velocity for each baffle arrangement condition of the air-type PVT collector with triangular baffles. The outlet velocity of the collector with the baffle was higher than that of the reference collector (0.564 m/s). It was confirmed that the smaller the left and right spacing variable (W) of the baffle was, the faster the outlet velocity and the higher the related average Reynolds number were. This is because only the arrangement was changed in a state in which the baffle size was fixed; as the length of the lateral direction (W) decreased, the flow path became narrower. However, the longer H2 was, the farther the front and back between the baffles were and the more space was created in the transverse direction. Consequently, the possibility of vortex generation by the baffle was reduced. Consequently, it was seen that the difference in the average Reynolds number range also decreased. Furthermore, the smaller the number of baffles was, the faster was the outlet velocity, but with a correspondingly lower Reynolds number. This is attributed to the flow in the air layer, which was non-uniform; the air flow rapidly passed through the flow paths in the left and right longitudinal directions in the collector with length of W.  Figure 8 shows the pressure drop by baffle placement conditions of the air-type PVT collector with triangular baffles. Since the reference collector had no baffles, turbulence did not occur, resulting in a low value of pressure of 0.2302 Pa. The PVT collector with triangular baffles had a pressure drop caused by the baffles. There was a difference in pressure drop value according to the length of W; it can be seen that the longer the distance was between W, the more the vortex area decreased and the smaller the pressure drop value was. In addition, the value of number of baffles (B.No.) of 35 led to relatively lower pressure drop than did the value of 28. It can be concluded that the number of baffles was large, but stagnation areas of air flow occurred less frequently due to the relatively shorter length of H2. Therefore, it was found that the arrangement of the baffles was more favorable than the number of baffles in terms of pressure drop, and the longitudinal spacing of the baffles (H1, H2) was the main influencing factor.  Table 5 shows the thermal efficiency and heat gain according to the arrangement parameters of the air-type PVT collector with triangular baffles. The thermal efficiency of the air-type PVT collector can be calculated by Equation (13) [26], which was deduced by applying the collector's flow rate and  Figure 8 shows the pressure drop by baffle placement conditions of the air-type PVT collector with triangular baffles. Since the reference collector had no baffles, turbulence did not occur, resulting in a low value of pressure of 0.2302 Pa. The PVT collector with triangular baffles had a pressure drop caused by the baffles. There was a difference in pressure drop value according to the length of W; it can be seen that the longer the distance was between W, the more the vortex area decreased and the smaller the pressure drop value was. In addition, the value of number of baffles (B.No.) of 35 led to relatively lower pressure drop than did the value of 28. It can be concluded that the number of baffles was large, but stagnation areas of air flow occurred less frequently due to the relatively shorter length of H2. Therefore, it was found that the arrangement of the baffles was more favorable than the number of baffles in terms of pressure drop, and the longitudinal spacing of the baffles (H1, H2) was the main influencing factor.  Figure 8 shows the pressure drop by baffle placement conditions of the air-type PVT collector with triangular baffles. Since the reference collector had no baffles, turbulence did not occur, resulting in a low value of pressure of 0.2302 Pa. The PVT collector with triangular baffles had a pressure drop caused by the baffles. There was a difference in pressure drop value according to the length of W; it can be seen that the longer the distance was between W, the more the vortex area decreased and the smaller the pressure drop value was. In addition, the value of number of baffles (B.No.) of 35 led to relatively lower pressure drop than did the value of 28. It can be concluded that the number of baffles was large, but stagnation areas of air flow occurred less frequently due to the relatively shorter length of H2. Therefore, it was found that the arrangement of the baffles was more favorable than the number of baffles in terms of pressure drop, and the longitudinal spacing of the baffles (H1, H2) was the main influencing factor.  Table 5 shows the thermal efficiency and heat gain according to the arrangement parameters of the air-type PVT collector with triangular baffles. The thermal efficiency of the air-type PVT collector can be calculated by Equation (13) [26], which was deduced by applying the collector's flow rate and  Table 5 shows the thermal efficiency and heat gain according to the arrangement parameters of the air-type PVT collector with triangular baffles. The thermal efficiency of the air-type PVT collector can be calculated by Equation (13) [26], which was deduced by applying the collector's flow rate and inlet/outlet temperature. Based on Equation (10), the reference efficiency without baffles was 26.7% and the heat gain was 287Wth. The thermal efficiency of the air-type PVT collector with triangular baffles was in the range of 34.3-35%, and the heat gain was analyzed to be 369-377 Wth. The collector's thermal efficiency was up to 13.5% higher than that of the reference collector, and the heat gain was 89.4 Wth. According to the baffle arrangement parameters, the thermal efficiency and heat gain were highest when the variable conditions are W = 144 mm, H1 = 47 mm, H2 = 83 mm and B.No. = 35.

Pressure Drop and Thermal Efficiency
A PVT is surface area of the collector (m 2 ), G is solar radiation (W/m 2 ), . m is mass flow rate (kg/h), C p is specific heat of air at a constant pressure (J/kg • C), T outlet is outlet air temperature of PVT ( • C), T inlet is inlet air temperature of PVT ( • C), and η th is thermal efficiency (-).

Conclusions
To confirm the impact of the arrangement of triangular shaped baffles fitted in an air-type PVT collector, the airflow characteristics in the collector was analyzed through a validated CFD model in this study. The heat transfer characteristics, pressure drop, and thermal efficiency of the collector were examined, and the results can be summarized as follows. - Depending on the variables, the outlet temperature increased by 3.6-3.9 • C and the heat gain increased by 1.28-1.31 times compared to the collector without baffles. Therefore, thermal performance of the collector improved up to 31%. - The thermal performance of the air-type PVT collector improved when the horizontal spacing of the baffles was wider and the vertical spacing was narrower. Furthermore, it was confirmed that greater number of baffles resulted in higher thermal energy yield of the collector. - The heat transfer performance increased 1.03 times of the maximum outlet velocity and 1.05 times of the average Reynolds number; according to the triangular baffle placement, the pressure drop increased by about 6.78 Pa under these conditions. Furthermore, variable conditions in which the baffles were more evenly placed (i.e., the case of W = 144 mm, H1 = 47 mm, H2 = 83 mm and B.N = 35) resulted in relatively small pressure drop and high thermal efficiency; thus, these conditions are judged to be advantageous for improving the performance of the air-type PVT collector.
Based on the results of this study, it can be deduced that the thermal performance improvement according to the staggered arrangement spacing of triangular baffles in the air-type PVT collector was not significant, but it was advantageous in terms of heat transfer in the collector when the baffle spacing was even.