Comparative Energy Performance Analysis of Six Primary Photovoltaic Technologies in Madrid ( Spain )

There are a wealth of studies on photovoltaic cell technologies, however their performance in different climatic or geographies over an extended period is not completely established. The objective of this paper is to add to this area of study with an analysis of the principle photovoltaic technologies: monocrystalline silicon (mc-Si), polycrystalline silicon (pc-Si), tandem structure of amorphous silicon and microcrystalline silicon (a-Si/μc-Si tandem), cadmium telluride/cadmium sulfide (CdTe/CdS), copper-indium diselenide (CIS) and monocrystalline silicon with double contact back technology (mc-dc-Si), in the climatic conditions particular to the middle latitude urban environment of Madrid (Spain). To carry out this study six photovoltaic subsystems of peak power 1 kWp approximately have been installed for each selected technology on the roof of the Escuela Técnica Superior de Ingeniería y Diseño Industrial (ETSIDI) of the Universidad Politécnica de Madrid (UPM) in Spain. Each subsystem has an inverter of the same model and power for its connection to the internal electricity network of the university. This paper analyzes the energy performance of the six subsystems, calculating the reference production rates and losses, from February 2013 to December 2015. The result of the study is the extensive capture of data and detailed analysis of real time energy yields and performance ratios of key technologies resulting with patterns in line with those of other regions with comparable climatic conditions.

The electrical characteristic parameters provided by manufacturers of photovoltaic modules have been obtained under standard test conditions (STC): irradiance 1000 W/m 2 , cell temperature 25 • C, air mass (AM) 1.5 and zero incidence angle; under normal operating cell (NOC): irradiance 800 W/m 2 , ambient temperature 20 • C and wind speed 1 m/s; and in conditions of low irradiance: irradiance 200 W/m 2 , cell temperature 25 • C and air mass (AM) 1.5.This information, including efficiency, is useful for comparing different technologies, but it does not provide complete information on the energy performance of the photovoltaic module at its installation site [37,38].For this reason, IEC 61853-1:2011 [39] introduces two additional operating conditions, known as high and low temperature, high temperature condition (HTC) 1000 W/m 2 and cell temperature 75 • C and low temperature condition (LTC) 500 W/m 2 and cell temperature 15 • C.However, these last two operating conditions are not currently included in the vendor data sheets.
The reduction of the energy generated with respect to incident solar energy of a photovoltaic system can be explained by a set of factors: operating temperature of the modules [40][41][42][43][44][45][46][47][48], dirt and dust, partial shading of the modules, spatial arrangement, angular and spectral response of each technology [49][50][51][52][53][54][55], mismatch loss or connection between modules, non-compliance with the nominal power referred to STC conditions, the behavior of the inverter to work at the maximum power point of the photovoltaic generator and its loss of efficiency [56][57][58][59], the loss of power due to the degradation of the photovoltaic generator over time [60][61][62][63], ohmic drops in direct current and alternate current wiring, and by faults, breakdowns or the network connection.
Photovoltaic systems currently deployed have different energy efficiency rates depending on the cell technology, components, design and operating conditions.The conventional technologies of crystalline silicon cell (c-Si) usually have a higher temperature related power loss ~−0.45%/K with higher efficiencies in the winter than in the summer [5,41,43,44,64].In the first hours of sun exposure the c-Si suffer power degradations of 0.5-1.5% [65][66][67] and have an annual power loss of 0.5-1%/year [44,60].They are less sensitive to variation in the solar spectrum ~1-2% [49,51,55] and their angular losses can reach 3% [54].
Thin layer technologies have a more complex electrical characterization, especially for a-Si and HIT technologies [7,68].Its nominal power tolerance can reach ±12% [47] and has a lower temperature related power loss coefficient −0.21 to −0.30%/K than conventional mc-Si and pc-Si technologies.They display higher efficiencies in the summer months, as is the case with a-Si and a-Si/µc-Si tandem technology, being optimum in hot or tropical climates [3,4,41,42,44,69,70]. Thin layer technologies take advantage of diffuse irradiation on cloudy days [52], have less dependence on the angle of inclination, but are more sensitive to variations in the solar spectrum, between 2-4% [1,55,68,71,72] than c-Si.The behavior of CIS and CdTe/CdS with solar irradiation and ambient temperature is similar to the c-Si with efficiency decreasing in the summer and increasing in the winter months [7,[42][43][44]47].The energy efficiency and life cycle of CIS and CdTe/CdS modules have been studied by Raugei et al. [62] demonstrating that these technologies can be competitive with respect to conventional technologies based on polycrystalline silicon.
A study carried out in 2010 by the manufacturer SUNPOWER [73] shows that the mc-dc-Si reaches an efficiency of around to 20.4% under STC conditions, has a lower power loss coefficient −0.38%/ • C which is less than conventional c-Si, makes better use of diffuse irradiation and is less affected by the variation of the solar incidence angle and the air mass.
It is well known that thin layer technologies also undergo initial degradation in the first few hours or days of exposure to sunlight.The technology a-Si/µc-Si tandem suffers a degradation in the value of its nominal power that can reach 0.8% [69,74,75] until its stabilization.An opposite effect occurs with the CIS technology where in the first hours of operation there is a positive increase in efficiency 7-15% [68,75,76].In the case of the CdTe/CdS the first few hours of solar exposure can increase efficiency by 6-8% or suffer degradation 7-15% [69,75,77] depending on the cell design and production process.In a five year study, Rodziewicz et al. [78] found that a-Si/µc-Si tandem technology had suffered a degradation of 10% of its nominal power, this value is higher than the loss suffered by mc-Si technology in the same period of time ~7%.Cañete et al. [43] have established in Malaga (Spain) in a yearlong study, that on certain days of high ambient temperature, that the daily efficiency drops to 5.4% for the CdTe, 6.5% for a-Si/µc-Si and 7.6% for pc-Si, compared to STC.
The inverter also has a significant influence on the energy injected into the grid with maximum efficiencies of 98%.Its efficiency is related to the value of the input voltage V DC and can produce variations of ±0.005 to 0.02%/V [59,79] depending on the type of inverter [80,81].Network-connected inverters work with maximum power point tracking algorithms that try to maximize the energy produced by the photovoltaic generator [57,58].As the point of maximum power changes with irradiation, temperature and shadows, there will be times when the inverter does not work at the point of maximum power.
On the other hand, in the real operating conditions of a photovoltaic installations there are incidents, breakdowns, disconnections to the grid, etc., which can affect energy production and this is why it is important to consider operative readiness.In a study of 78 photovoltaic installations in northern and eastern Germany [82] 63% of downtime was triggered by inverters, 15% by photovoltaic modules and 22% on failures of the rest of the components of the system.
In summary, energy generation of photovoltaic systems is affected by a large number of variables which need to be observed over an extended period of time.This paper shows the operation performance data of six primary photovoltaic technologies, in the climatic conditions particular to Madrid, from February 2013 to December 2015.The result is a detailed analysis of energy yields and performance ratios of these key technologies, displaying patterns in line with those of other regions with comparable climatic conditions.
Section 2 describes the photovoltaic systems under investigation.Section 3 presents and analyzes the recorded meteorological data.In Section 4 the energy parameters of each photovoltaic technology are defined and calculated, the results are presented in tables and graphs to facilitate their comparison.In Section 5, new energy parameters are defined and calculated taking into account the availability of each subsystem.Again tables and graphs are used to display the output permitting a better comparison of the performance of each type of photovoltaic technology.Section 5 is completed with an analysis of the evolution of operational efficiency.The study finishes with a summary of the findings and is followed by the complete bibliography.

Description of the Photovoltaic and Monitoring System
The photovoltaic system under study is installed in the main building of the Escuela Técnica Superior de Ingeniería y Diseño Industrial of the Universidad Politécnica de Madrid (ETSIDI-UPM): latitude 40.4 • , longitude −3.7 • and altitude 657 m.The building is in the center of the city of Madrid, where its flat roof is well exposed to solar radiation with shading of nearby buildings reduced to positions of the sun just after sunrise and before sunset.The site has a continental climate with cold winters and hot summers.The object of the investigation consists of 6 subsystems of different cell technologies: mc-Si, pc-Si, a-Si/µc-Si tandem, CdTe/CdS, CIS and mc-dc-Si mounted on weighted fixed tilt structures.All the modules of the different technologies are coplanar with a tilt of 30 • and azimuth of 19 • east to optimize the spatial distribution according to the architectural requirements of the roof.The structure provides a separation of 20 cm between the photovoltaic modules allowing natural cooling.The photovoltaic modules have been selected using models and manufacturers representative of each cell technology, Table 1 shows their main technical characteristics.All photovoltaic modules are conventional: glass top layer, white Tedlar backsheet and aluminum frame, except CdTe/CdS which has frameless glass to glass modules and CIS technology with glass to glass modules with frame.Table 2 describes the photovoltaic subsystems and Table 3 shows the main characteristics of the installed inverter, which is identical in all subsystems to facilitate performance comparison.A photograph of the photovoltaic system is shown in Figure 1.The influence of climatic variables on the performance of the photovoltaic systems are measured in accordance with Norm IEC 61724: 1998 and the guidelines of the Joint Research Center in Ispra, Italy [83][84][85].The global solar irradiation (HI) (Wh/m 2 ) data is captured by means of a thermoelectric pyranometer (PIR) and a calibrated reference cell (CRC) [86] of polycrystalline silicon, The influence of climatic variables on the performance of the photovoltaic systems are measured in accordance with Norm IEC 61724: 1998 and the guidelines of the Joint Research Center in Ispra, Italy [83][84][85].The global solar irradiation (H I ) (Wh/m 2 ) data is captured by means of a thermoelectric pyranometer (PIR) and a calibrated reference cell (CRC) [86] of polycrystalline silicon, both coplanar with the photovoltaic modules object of the investigation.The module temperature is measured with a PT-1000 thermocouple sensor fixed to the backsheet of a central cell, of one of the central modules of the array.Ambient temperature and relative humidity are measured with a thermohygrometer while wind speed uses an anemometer (Figure 2).The technical specifications of these sensors are shown in Table 4.
Energies 2017, 10, 772 5 of 23 both coplanar with the photovoltaic modules object of the investigation.The module temperature is measured with a PT-1000 thermocouple sensor fixed to the backsheet of a central cell, of one of the central modules of the array.Ambient temperature and relative humidity are measured with a thermohygrometer while wind speed uses an anemometer (Figure 2).The technical specifications of these sensors are shown in Table 4.The electrical and meteorological variables (see Figure 3), are sampled every ten seconds with high precision Meteo Control Pro and IO web'log data loggers.The recording intervals are every five minutes, obtaining a representative average value for each interval.The measured electrical variables are taken from the inverters through the (recommended standard) RS-485 interface.Signal converters with 0.2 margin of error have been used for the meteorological variables including solar radiation, ambient temperature, relative humidity, wind speed and cells temperatures of the modules.This data can be consulted on the website CONERSA [87].
The electrical and meteorological variables (see Figure 3), are sampled every ten seconds with high precision Meteo Control Pro and IO web'log data loggers.The recording intervals are every five minutes, obtaining a representative average value for each interval.The measured electrical variables are taken from the inverters through the (recommended standard) RS-485 interface.Signal converters with 0.2 margin of error have been used for the meteorological variables including solar radiation, ambient temperature, relative humidity, wind speed and cells temperatures of the modules.This data can be consulted on the website CONERSA [87].

Solar Irradiation
The mean values are obtained every five minutes with the pyranometer and the reference solar cell with an tilt of 30° and azimuth of 19° east, are contrasted with an average meteorological year (Table 5 and Figure 4) obtained from the Photovoltaic Geographic Information System (Ispra-PVGIS) [30].The monthly values of the diffuse fraction of solar irradiation Dh/GI are taken from PVGIS and the HELIOS Solar Energy Institute (UPM) meteorological station [31] located 5 km away from the test site.Table 5 and Figure 4 show solar irradiation measurements and present small variations with respect to the average meteorological year.Clear sky conditions, low diffuse/global ratio Dh/GI ≈ 30% [7,86,88] found from June to September, produce the largest differences between the signals measured and those of HELIOS and the values of irradiation are higher in the pyranometer than the CRC.While in the winter months Dh/GI is close to ~45% the mean monthly irradiation values are similar for CRC and pyranometer (Table 5).

Solar Irradiation
The mean values are obtained every five minutes with the pyranometer and the reference solar cell with an tilt of 30 • and azimuth of 19 • east, are contrasted with an average meteorological year (Table 5 and Figure 4) obtained from the Photovoltaic Geographic Information System (Ispra-PVGIS) [30].The monthly values of the diffuse fraction of solar irradiation D h /G I are taken from PVGIS and the HELIOS Solar Energy Institute (UPM) meteorological station [31] located 5 km away from the test site.Table 5 and Figure 4 show solar irradiation measurements and present small variations with respect to the average meteorological year.Clear sky conditions, low diffuse/global ratio D h /G I ≈ 30% [7,86,88] found from June to September, produce the largest differences between the signals measured and those of HELIOS and the values of irradiation are higher in the pyranometer than the CRC.While in the winter months D h /G I is close to ~45% the mean monthly irradiation values are similar for CRC and pyranometer (Table 5).On average, the annual irradiation measured by CRC is ~2-4% lower than irradiance measured by pyranometers [86].The spectral response of the CRC is usually between 0.3-1.1 µm, while the pyranometer is 0.3-3 µm [88] because the pyranometer has a broader spectral and angular response that leads to lower PR values especially in the summer months Figure 4.The CRC measurement of irradiation is more effective on days of mixed cloud cover because it has a shorter response time ≅ 1 ms to the pyranometer 5 to 20 s, demonstrating higher irradiance values [88].The albedo effect needs to be taken into account since both devices are positioned with an inclination of 30° producing a systematic error in measured spectrum, being larger for the pyranometer.Higher performance values are obtained, especially in the summer months, when the calculations are measured by the CRC.This paper uses CRC measurements for the energy calculations.

Temperature of Operation of Modules
The energy production of each subsystem depends to a large extent on the operating cells temperature of the modules [40][41][42][43][44][45][46][47][48] which varies according to many factors: the ambient temperature (Tamb), the effective solar radiation, the wind speed [40,89], dirt that can cause hot spots and the characteristics of the installation.With the increase of the module temperature the open circuit voltage decreases and the short-circuit current rises slightly.Combining both effects results in a decrease in the peak power of the module.Figure 5 shows the mean monthly values of ambient temperature and temperature of the modules during the years 2013-2015, displaying a pattern in line with the local climatology.
Figure 5 shows that over practically the whole year the highest average module temperatures correspond to the CIS and CdTe/Cds technologies, both with glass to glass format, while the lowest correspond to the mc-cd-Si, with tedlar, although the differences are minimum (it should be noted that the mc-cd-Si modules are at the end of the rooftop and more exposed to airflows).Figure 6 shows the average monthly wind velocity values of one representative year.The anemometer is located 4 m high with respect to the position of the photovoltaic modules.The average wind speed in the summer months, with a higher solar irradiation and ambient temperature, is around 1.5 m/s so that its influence on the temperature of the modules can be considered of little relevance [89].On average, the annual irradiation measured by CRC is ~2-4% lower than irradiance measured by pyranometers [86].The spectral response of the CRC is usually between 0.3-1.1 µm, while the pyranometer is 0.3-3 µm [88] because the pyranometer has a broader spectral and angular response that leads to lower PR values especially in the summer months Figure 4.The CRC measurement of irradiation is more effective on days of mixed cloud cover because it has a shorter response time ∼ = 1 ms to the pyranometer 5 to 20 s, demonstrating higher irradiance values [88].The albedo effect needs to be taken into account since both devices are positioned with an inclination of 30 • producing a systematic error in measured spectrum, being larger for the pyranometer.Higher performance values are obtained, especially in the summer months, when the calculations are measured by the CRC.This paper uses CRC measurements for the energy calculations.

Temperature of Operation of Modules
The energy production of each subsystem depends to a large extent on the operating cells temperature of the modules [40][41][42][43][44][45][46][47][48] which varies according to many factors: the ambient temperature (T amb ), the effective solar radiation, the wind speed [40,89], dirt that can cause hot spots and the characteristics of the installation.With the increase of the module temperature the open circuit voltage decreases and the short-circuit current rises slightly.Combining both effects results in a decrease in the peak power of the module.Figure 5 shows the mean monthly values of ambient temperature and temperature of the modules during the years 2013-2015, displaying a pattern in line with the local climatology.
Figure 5 shows that over practically the whole year the highest average module temperatures correspond to the CIS and CdTe/Cds technologies, both with glass to glass format, while the lowest correspond to the mc-cd-Si, with tedlar, although the differences are minimum (it should be noted that the mc-cd-Si modules are at the end of the rooftop and more exposed to airflows).Figure 6 shows the average monthly wind velocity values of one representative year.The anemometer is located 4 m high with respect to the position of the photovoltaic modules.The average wind speed in the summer months, with a higher solar irradiation and ambient temperature, is around 1.5 m/s so that its influence on the temperature of the modules can be considered of little relevance [89].

Energy Production of Photovoltaic Installation
Table 6 shows the energy production in kWh/year and in kWh/m 2 for each subsystem in the study period.

Energy Production of Photovoltaic Installation
Table 6 shows the energy production in kWh/year and in kWh/m 2 for each subsystem in the study period.

Energy Production of Photovoltaic Installation
Table 6 shows the energy production in kWh/year and in kWh/m 2 for each subsystem in the study period.

Energy yield and Performance Ratio
In order to analyze the energy yield production indexes are calculated, according to Figure 7, of the standard IEC 61724: 98 [83] for each subsystem.Data is collated and available on databases around the world using production indexes of hundreds of photovoltaic systems located in different places and climatic conditions [7,82,90].

Energy yield and Performance Ratio
In order to analyze the energy yield production indexes are calculated, according to Figure 7, of the standard IEC 61724: 98 [83] for each subsystem.Data is collated and available on databases around the world using production indexes of hundreds of photovoltaic systems located in different places and climatic conditions [7,82,90].The available solar energy over a period of time T, is determined by the reference production index (YR) obtained dividing the global irradiance (GI) on the plane of the photovoltaic generator by the STC irradiance (GI,STC) (1 kW/m 2 ) as per Equation (1): The array yield YA and the final yield YF for the period of time T, are calculated with Equations ( 2) and (3), respectively: where: EAC = electricity generated by the photovoltaic system in kWh in period T; EDC = electricity generated at the input point of the inverter in kWh in period T; PSTC = nominal power of the photovoltaic generator STC conditions in kWp.The losses LC and LBOS can be obtained by Equations ( 4) and (5): Figures 8a,b and 9a,b show the monthly evolution of energy yields YA, YF and losses parameters LC and LBOS in kWh/kWpSTC without taking into account the operational availability for each subsystem over the studied period.
The array yield Y A and the final yield Y F for the period of time T, are calculated with Equations ( 2) and (3), respectively: where: E AC = electricity generated by the photovoltaic system in kWh in period T; E DC = electricity generated at the input point of the inverter in kWh in period T; P STC = nominal power of the photovoltaic generator STC conditions in kWp.The losses L C and L BOS can be obtained by Equations ( 4) and ( 5):  Figure 8a,b show the same evolution as YA and YF.The values of these two yields depend on certain meteorological conditions: solar radiation on the plane of the modules, ambient temperature and wind speed, which are the same for all the subsystems.These two indexes are related to the efficiency of the inverter according to Equation 9.This efficiency is very similar in all subsystems which use the same model of inverter and variations in the study period are minimal (see Section 4.3).Figure 9b follows a similar pattern, where the CdTe/CdS subsystem has declined since October 2014 and also suffered certain periods of inverter downtime.Figure 9a shows that the losses LC follow the same tendency as the temperatures of the modules, increasing in the summer months.
On the other hand, LC presents important differences between the different technologies.As the inverter is installed just behind the modules of the corresponding subsystem, losses in the DC wiring are negligible.The lowest losses are observed in CIS technology, and are even negative during the first two months of operation due to the increase at the initiation phase of the its efficiency [75], despite reaching highest module temperatures.The highest LC are observed in the CdTe/CdS technology from the outset and especially in the summer months.In the months of November and December of 2014 and January, May, June and July of 2015, there have been unexpected inverter stopages affecting drastically the value of YA.
The PR is calculated for the DC side (PRDC) and AC side (PRAC) using Equations ( 6) and (7).Normally the PR values are calculated on a monthly or annual basis [7,82,90].Calculated values for smaller intervals, such as weekly or daily, can be useful for identifying faults in the components of the installation.Figures 10 and 11 present the monthly evolution of PRDC and PRAC of all photovoltaic technologies in the years 2013, 2014 and 2015.Figure 8a,b show the same evolution as Y A and Y F .The values of these two yields depend on certain meteorological conditions: solar radiation on the plane of the modules, ambient temperature and wind speed, which are the same for all the subsystems.These two indexes are related to the efficiency of the inverter according to Equation 9.This efficiency is very similar in all subsystems which use the same model of inverter and variations in the study period are minimal (see Section 4.3).Figure 9b follows a similar pattern, where the CdTe/CdS subsystem has declined since October 2014 and also suffered certain periods of inverter downtime.Figure 9a shows that the losses L C follow the same tendency as the temperatures of the modules, increasing in the summer months.
On the other hand, L C presents important differences between the different technologies.As the inverter is installed just behind the modules of the corresponding subsystem, losses in the DC wiring are negligible.The lowest losses are observed in CIS technology, and are even negative during the first two months of operation due to the increase at the initiation phase of the its efficiency [75], despite reaching highest module temperatures.The highest L C are observed in the CdTe/CdS technology from the outset and especially in the summer months.In the months of November and December of 2014 and January, May, June and July of 2015, there have been unexpected inverter stopages affecting drastically the value of Y A .
The PR is calculated for the DC side (PR DC ) and AC side (PR AC ) using Equations ( 6) and (7).Normally the PR values are calculated on a monthly or annual basis [7,82,90].Calculated values for smaller intervals, such as weekly or daily, can be useful for identifying faults in the components of the installation.Figures 10 and 11 present the monthly evolution of PR DC and PR AC of all photovoltaic technologies in the years 2013, 2014 and 2015.Due to losses to the photovoltaic modules caused by the temperature, the PRDC values are higher in winter than in summer, it is the inverse evolution to LC except in the a-Si/µc-Si tandem technology where the PRDC value is greater during summer [41][42][43][44].Dirt in the photovoltaic modules will also affect the annual evolution of the PRDC [63].A decrease in PRDC values over the three years is observed, caused by the loss of efficiency of photovoltaic generators due to their power degradation (see Section 4.4).The lowest decrease in PRDC in the three years is in mc-Si technology with 4%, followed by mc-dc-Si 5.2%, pc-Si 5.6% and a-Si/µc-Si 5.9%.CIS technology has decreased 9% and CdTe/CdS technology stands out with a decrease of 24.1%.

Photovoltaic Generator Efficiency
The efficiency of the photovoltaic generator is the ratio of the energy generated in EDC in kWh with respect to the incident irradiation (HI) in kWh/m 2 during the same period of time T, multiplied by the capture area A in m 2 according to Equation (8).The results are shown in Table 7.
. 100 (8) The loss of efficiency in the three years studied (ΔηG) varies depending on the technology, with the extreme values for CIS 1.6% and 4.7% for CdTe/CdS.Due to losses to the photovoltaic modules caused by the temperature, the PRDC values are higher in winter than in summer, it is the inverse evolution to LC except in the a-Si/µc-Si tandem technology where the PRDC value is greater during summer [41][42][43][44].Dirt in the photovoltaic modules will also affect the annual evolution of the PRDC [63].A decrease in PRDC values over the three years is observed, caused by the loss of efficiency of photovoltaic generators due to their power degradation (see Section 4.4).The lowest decrease in PRDC in the three years is in mc-Si technology with 4%, followed by mc-dc-Si 5.2%, pc-Si 5.6% and a-Si/µc-Si 5.9%.CIS technology has decreased 9% and CdTe/CdS technology stands out with a decrease of 24.1%.

Photovoltaic Generator Efficiency
The efficiency of the photovoltaic generator is the ratio of the energy generated in EDC in kWh with respect to the incident irradiation (HI) in kWh/m 2 during the same period of time T, multiplied by the capture area A in m 2 according to Equation (8).The results are shown in Table 7.
. 100 (8) The loss of efficiency in the three years studied (ΔηG) varies depending on the technology, with the extreme values for CIS 1.6% and 4.7% for CdTe/CdS.Due to losses to the photovoltaic modules caused by the temperature, the PR DC values are higher in winter than in summer, it is the inverse evolution to L C except in the a-Si/µc-Si tandem technology where the PR DC value is greater during summer [41][42][43][44].Dirt in the photovoltaic modules will also affect the annual evolution of the PR DC [63].A decrease in PR DC values over the three years is observed, caused by the loss of efficiency of photovoltaic generators due to their power degradation (see Section 4.4).The lowest decrease in PR DC in the three years is in mc-Si technology with 4%, followed by mc-dc-Si 5.2%, pc-Si 5.6% and a-Si/µc-Si 5.9%.CIS technology has decreased 9% and CdTe/CdS technology stands out with a decrease of 24.1%.

Photovoltaic Generator Efficiency
The efficiency of the photovoltaic generator is the ratio of the energy generated in E DC in kWh with respect to the incident irradiation (H I ) in kWh/m 2 during the same period of time T, multiplied by the capture area A in m 2 according to Equation (8).The results are shown in Table 7.
The loss of efficiency in the three years studied (∆η G ) varies depending on the technology, with the extreme values for CIS 1.6% and 4.7% for CdTe/CdS.

Inverter Efficiency
Once input E DC and output E AC energy of the inverters are known, the conversion efficiency of each inverter can be calculated as per Equation ( 9): The average annual inverter efficiency obtained during operation Table 8 have remained close to 90% in the three years of study and for all subsystems, so their contribution has been homogeneous, thus facilitating the comparison of the different photovoltaic technologies.Table 8 shows that the efficiencies of the inverters are in line with the vendor data of 90.9-92.1%European efficiency rates, except for the inverters of subsystems 4 and 5 which are 0.7% and a 0.4% less, respectively.

Global Energy Efficiency
The global energy efficiency of each of the subsystems Figure 12, is calculated as the product of the efficiencies of the photovoltaic generator and the inverter according to Equation (10). Figure 12 shows the efficiency of each technology and effects of aging: Energies 2017, 10, 772 13 of 23

Inverter Efficiency
Once input EDC and output EAC energy of the inverters are known, the conversion efficiency of each inverter can be calculated as per Equation ( 9): The average annual inverter efficiency obtained during operation Table 8 have remained close to 90% in the three years of study and for all subsystems, so their contribution has been homogeneous, thus facilitating the comparison of the different photovoltaic technologies.Table 8 shows that the efficiencies of the inverters are in line with the vendor data of 90.9-92.1%European efficiency rates, except for the inverters of subsystems 4 and 5 which are 0.7% and a 0.4% less, respectively.

Global Energy Efficiency
The global energy efficiency of each of the subsystems Figure 12, is calculated as the product of the efficiencies of the photovoltaic generator and the inverter according to Equation (10). Figure 12 shows the efficiency of each technology and effects of aging: The values of the calculated global efficiencies of the photovoltaic systems depend mainly on the initial nominal efficiency ratings and their reductions are in line with the usual power degradation corresponding to each technology [43,44,60].

Operating Availability Factor, Corrected Energy Yield and Performance Ratios
There are multiple causes of lack of availability in a photovoltaic system and in some cases only a part of the system will be affected, for example, disconnection of a string of modules, while in others there can be total system shutdown caused by tripping of AC protections, absence of network, etc.In order to consider PR and availability independently, a new value of Y * R is defined which only takes into account solar irradiation when the AC power at the inverter output is different from zero.This Y * R (Equation ( 11)) allows the definition of new energy parameters of the photovoltaic subsystems L * C , PR * DC and PR * AC (Equations ( 12)-( 14)) and availability index (D) (Equation ( 15)) that eliminate the influence of the difference in the startup and stopping of the inverters and the penalties for faults outside the photovoltaics array, inverters and power grid.These new parameters are exclusively associated with the photovoltaic generator and inverter efficiencies.The availability losses affected only to L * C , PR * DC and PR * AC while that the values of Y F , Y A and L BOS remain the same as Section 4.
where G * I and H * I are the global solar irradiance (W/m 2 ) and the global solar radiation (Wh/m 2 ) respectively on the plane of the photovoltaic generator over a period of time T for alternating power above zero.
The new values of the production and loss indexes are shown in the Tables 9 and 10.How the values of PR * DC and PR * AC have increased with respect to those initially calculated can be seen.The D factor of all the subsystems has remained high in the years 2013 and 2014.In the year 2015 it has fallen ~3% in all the subsystems due to power cuts of the general network, outside the photovoltaic system, except in the subsystem 4 that has more downtime.CdTe/CdS technology continues to have the greatest decrease in the value of PR    The corrected annual operational efficiency ( * ) and aggregated loss of efficiency in the three years studied (Δ * ) of the photovoltaic array of each technology has been calculated (Equation ( 16)) again from this new scenario (Table 11).
 * (%) = * .× 100 (16)       The corrected annual operational efficiency ( * ) and aggregated loss of efficiency in the three years studied (Δ * ) of the photovoltaic array of each technology has been calculated (Equation ( 16)) again from this new scenario (Table 11).
 * (%) = * .× 100 (16)   The corrected annual operational efficiency (η * G ) and aggregated loss of efficiency in the three years studied (∆η * G ) of the photovoltaic array of each technology has been calculated (Equation ( 16)) again from this new scenario (Table 11).The annual loss of efficiency in all technologies are very similar to the values obtained in studies carried out with the same technologies in other locations [41][42][43][44].Figure 15 shows the evolution of the monthly-corrected operational efficiency of all technologies.In the first month of solar exposure, February 2013, there is a significant loss of efficiency in all technologies with respect to the values provided by the vendor under STC conditions, except in CIS technology, being 0.8% mc-Si, 1.50% pc-Si, 1.4% a-Si/µc-Si, 1.1% CdTe/CdS, −0.1% CIS and 1.3% mc-dc-Si.Figure 14 shows the annual operational efficiency losses.The technologies pc-Si and CdTe/CdS reduced 1.9% and 1.8% respectively in the first year with CIS technology losing only 0.4%.This annual efficiency loss has been lower in the years 2014 and 2015 as can be seen in Figure 16.
The annual loss of efficiency in all technologies are very similar to the values obtained in studies carried out with the same technologies in other locations [41][42][43][44].Figure 15 shows the evolution of the monthly-corrected operational efficiency of all technologies.In the first month of solar exposure, February 2013, there is a significant loss of efficiency in all technologies with respect to the values provided by the vendor under STC conditions, except in CIS technology, being 0.8% mc-Si, 1.50% pc-Si, 1.4% a-Si/µc-Si, 1.1% CdTe/CdS, −0.1% CIS and 1.3% mc-dc-Si.Figure 14 shows the annual operational efficiency losses.The technologies pc-Si and CdTe/CdS reduced 1.9% and 1.8% respectively in the first year with CIS technology losing only 0.4%.This annual efficiency loss has been lower in the years 2014 and 2015 as can be seen in Figure 16.The annual loss of efficiency in all technologies are very similar to the values obtained in studies carried out with the same technologies in other locations [41][42][43][44].Figure 15 shows the evolution of the monthly-corrected operational efficiency of all technologies.In the first month of solar exposure, February 2013, there is a significant loss of efficiency in all technologies with respect to the values provided by the vendor under STC conditions, except in CIS technology, being 0.8% mc-Si, 1.50% pc-Si, 1.4% a-Si/µc-Si, 1.1% CdTe/CdS, −0.1% CIS and 1.3% mc-dc-Si.Figure 14 shows the annual operational efficiency losses.The technologies pc-Si and CdTe/CdS reduced 1.9% and 1.8% respectively in the first year with CIS technology losing only 0.4%.This annual efficiency loss has been lower in the years 2014 and 2015 as can be seen in Figure 16.Table 12 provides a summary of the study.The percentage variations of the loss of efficiency in the first month, the first year and the study period are shown with respect nominal efficiency STC of each photovoltaic technology.Also included are the percentage changes in global efficiency and energy production indexes PR * AC and PR * DC .The efficiency losses in the first month are similar in all technologies except mc-Si that presents minor losses and CIS which has a slightly positive balance positive due to the increase in efficiency during the first hours of sun exposure.The trend continues throughout the first year except for the a-Si/µc-Si technology which has a small yearly increase resulting in a reduction of efficiency loss compared to the first month.The variations in photovoltaic efficiency over the study period is similar for the crystalline silicon technologies, with a-Si/µc-Si showing intermediate values while CdTe/CdS and CIS are at the highest and lowest end of the spectrum respectively.Moreover, CdTe/CdS displays the largest decrease of the global efficiency.Regarding the variations of the PR * DC and PR * AC , CdTe/CdS and CIS are highest while the rest of the technologies are very similar.

Conclusions
A photovoltaic installation on the rooftop of the university campus has permitted a comparative study of the energy production rates of six selected photovoltaic technologies connected to the internal electricity network of the university using the same model of inverter under the same physical and climatic conditions, over a period of three years with the following conclusions.The solar irradiation measurements over the study period present small variations with respect to the average meteorological year.The ambient temperature has followed the usual pattern of local climate.The influence of the wind speed can be considered of little relevance.The use of the availability index allows the energy comparative analysis of the technologies for the photovoltaic generator and inverter efficiencies.
The mc-Si, pc-Si, CIS and mc-cd-Si technologies reach an average value of PR * AC above 80%, and a-Si/µc-Si and CdTe/CdS remain at 74.5% and 64.3%, respectively.The loss of efficiency in all technologies during the first month is evident, except for the CIS technology because it initially achieves a gain in efficiency.
The conventional technologies mc-Si and pc-Si displayed very similar thermal and energy behavior.The decrease in the PR * DC and PR * AC values during the three years were lower than for the other technologies, which indicates a more stable behavior, being the values of the pc-Si technology, which had a corrected capture loss of 9%, the highest, while for mc-Si technology it was 12.4%, despite having suffered a major power degradation during the first month.
With respect to thin-film technologies, CIS technology (subsystem 5), is the one that reaches a higher temperature in the summer months.The losses in CIS technology, PR * AC value and global efficiency performance in the three years were 6.3% and 0.79%, respectively, and the generator efficiency loss was 1.3%.Its initial degradation during the first year (0.4%) was the lowest.The mean values of PR * DC and PR * AC during the three years are the highest, with 94.5% and 84.6%, respectively.The corrected capture losses have been the lowest of all technologies at around 6%.
The behavior of the a-Si/µc-Si tandem technology obtains better PR * AC value results during the summer months than the other technologies, having a coefficient of loss of power with the low temperature, confirming it is a more appropriate technology for warm climates.Moreover, the loss of efficiency of the generator in the first year and over the entire number of years has been the lowest, with the exception of CIS technology.The average value of PR * AC is 74.5% and its corrected capture losses reach 20.6%.
CdTe/CdS technology is the one that has had the worst performance.Its loss of efficiency in the first year was similar to that of the pc-Si technology, but during the following two years it suffered a degradation of 2%.Mean values of PR * DC and PR * AC during the three years are the lowest of all technologies, with 72% and 64.3%, respectively.Its corrected capture losses reach a value of 28.6%.This behavior and degradation has been confirmed in previous studies.
The high-efficiency mc-dc-Si technology has an initial loss of efficiency similar to those of the mc-Si and pc-Si technologies.The module operating temperature during the summer months is the lowest out of all technologies.It is the most efficient technology in STC conditions, which translates into the greater value of overall efficiency, but its decrease during the investigation period was 0.8% higher than mc-Si and pc-Si, with 0.53% and 0.54%, respectively.The average PR * AC and corrected capture losses were 80.6% and 10.5%, respectively.
This study expands the performance database of the principal photovoltaic technologies for middle latitude continental climates over an extended period of time, thus enriching the data available to calculate energy produced in the long term, key for determining the Levelized Cost of Energy (LCOE) [91] and the Energy Payback Time (EPBT) [92] of photovoltaic systems.
The study concentrates on side-by-side comparisons and analysis between different commercial PV technologies in the same urban location in Madrid (Spain).The emphasis is placed on the operational availability of each of the platforms as seen by the use of the corrected performance ratio values, permitting an in depth exhaustive study of the concerned technologies.The results obtained add to the body of photovoltaic performance data available worldwide.

Figure 1 .
Figure 1.PV systems on the rooftop of the ETSIDI-UPM in Madrid.

Figure 2 .
Figure 2. Measuring sensors placed on the rooftop of the ETSIDI-UPM.(a) Thermohygrometer and anemometer.(b) PIR and CRC installed with a tilt of 30° and azimuth of 19° east.

Figure 2 .
Figure 2. Measuring sensors placed on the rooftop of the ETSIDI-UPM.(a) Thermohygrometer and anemometer.(b) PIR and CRC installed with a tilt of 30 • and azimuth of 19 • east.

2 )Figure 3 .
Figure 3. Monitored electrical and meteorological variables of each of the photovoltaic subsystems.

Figure 3 .
Figure 3. Monitored electrical and meteorological variables of each of the photovoltaic subsystems.

igure 4 .
J S N J M M J S N J M M J S N Evolution of global solar radiation H I (kWh/m 2 ), Years 2013-2015.

Figure 5 .
Figure 5. Evolution of average ambient and modules temperatures.

Figure 7 .
Figure 7. Energy yields and losses parameters.

Figure 7 .
Figure 7. Energy yields and losses parameters.

Figures 8a,b and
Figures 8a,b and 9a,b show the monthly evolution of energy yields Y A , Y F and losses parameters L C and L BOS in kWh/kWp STC without taking into account the operational availability for each subsystem over the studied period.
show the mean values corresponding to three years of * , YA, YF, * , LBOS, * , * and factor D.
show the mean values corresponding to three years of * , YA, YF, * , LBOS, * , * and factor D.

Figure 15 .
Figure 15.Evolution of the operating PV generators efficiency.

Figure 15 .
Figure 15.Evolution of the operating PV generators efficiency.

Table 1 .
Technical parameters of photovoltaic modules.

Table 1 .
Technical parameters of photovoltaic modules.
Figure 1.PV systems on the rooftop of the ETSIDI-UPM in Madrid.

Table 4 .
Measuring sensors: manufacturer and specifications.

Table 4 .
Measuring sensors: manufacturer and specifications.

Table 5 .
Average monthly values of global solar irradiation HI (kWh/m 2 ) and diffuse irradiation fraction (Dh/GI) measured with a tilt of 30° and azimuth of 19° east.

Table 5 .
Average monthly values of global solar irradiation H I (kWh/m 2 ) and diffuse irradiation fraction (D h /G I ) measured with a tilt of 30 • and azimuth of 19 • east.

Table 6 .
Annual energy production of each subsystem.

Table 6 .
Annual energy production of each subsystem.

Table 6 .
Annual energy production of each subsystem.

Table 7 .
Nominal and operational efficiency of the photovoltaic generators.

Table 7 .
Nominal and operational efficiency of the photovoltaic generators.

Table 10 .
Comparison between values of PR DC , PR * DC (%) and PR AC , PR * AC (%) and operative availability factor D (%) in three years.
AC of 15.5%, ending with 56.2% in 2015.The CIS technology continues to have the highest PR * AC value, with 81.9% in 2015.Figures 13 and 14 show the mean values corresponding to three years of Y * R , Y A , Y F , L * C , L BOS , PR * DC , PR * AC and factor D. have increased with respect to those initially calculated can be seen.The D factor of all the subsystems has remained high in the years 2013 and 2014.In the year 2015 it has fallen ~3% in all the subsystems due to power cuts of the general network, outside the photovoltaic system, except in the subsystem 4 that has more downtime.CdTe/CdS technology continues to have the greatest decrease in the value of * of 15.5%, ending with 56.2% in 2015.

Table 11 .
Corrected operational efficiency and cumulative losses of the PV array.
have increased with respect to those initially calculated can be seen.The D factor of all the subsystems has remained high in the years 2013 and 2014.In the year 2015 it has fallen ~3% in all the subsystems due to power cuts of the general network, outside the photovoltaic system, except in the subsystem 4 that has more downtime.CdTe/CdS technology continues to have the greatest decrease in the value of * of 15.5%, ending with 56.2% in 2015.

Table 11 .
Corrected operational efficiency and cumulative losses of the PV array.

Table 11 .
Corrected operational efficiency and cumulative losses of the PV array.

Table 12 .
Percentage changes produced in generators efficiency and performance ratios.